1
|
Ferguson S, Jones A, Murray K, Andrew RL, Schwessinger B, Bothwell H, Borevitz J. Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus. Gigascience 2024; 13:giae029. [PMID: 38869149 PMCID: PMC11170218 DOI: 10.1093/gigascience/giae029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Revised: 03/11/2024] [Accepted: 05/14/2024] [Indexed: 06/14/2024] Open
Abstract
Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.
Collapse
Affiliation(s)
- Scott Ferguson
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory, 2600 Australia
| | - Ashley Jones
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory, 2600 Australia
| | - Kevin Murray
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory, 2600 Australia
- Department of Molecular Biology, Max Planck Institute for Biology Tübingen, Tübingen, 72076 Germany
| | - Rose L Andrew
- Botany & N.C.W. Beadle Herbarium, School of Environmental and Rural Science, University of New England, Armidale, NSW 2351, Australia
| | - Benjamin Schwessinger
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory, 2600 Australia
| | - Helen Bothwell
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory, 2600 Australia
- Warnell School of Forestry & Natural Resources, University of Georgia, Athens 30602 GA, United States
| | - Justin Borevitz
- Research School of Biology, Australian National University, Canberra, Australian Capital Territory, 2600 Australia
| |
Collapse
|
2
|
Hong-Le T, Crouse WL, Keele GR, Holl K, Seshie O, Tschannen M, Craddock A, Das SK, Szalanczy AM, McDonald B, Grzybowski M, Klotz J, Sharma NK, Geurts AM, Key CCC, Hawkins G, Valdar W, Mott R, Solberg Woods LC. Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats. Diabetes 2023; 72:135-148. [PMID: 36219827 PMCID: PMC9797320 DOI: 10.2337/db22-0252] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 10/04/2022] [Indexed: 01/21/2023]
Abstract
Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.
Collapse
Affiliation(s)
- Thu Hong-Le
- Genetics Institute, University College London, London, U.K
| | - Wesley L. Crouse
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | | | - Katie Holl
- Medical College of Wisconsin, Milwaukee, WI
| | - Osborne Seshie
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | | | - Ann Craddock
- Department of Biochemistry, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Swapan K. Das
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Alexandria M. Szalanczy
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Bailey McDonald
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | | | | | - Neeraj K. Sharma
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | | | - Chia-Chi Chuang Key
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| | - Gregory Hawkins
- Department of Biochemistry, Wake Forest University School of Medicine, Winston-Salem, NC
| | - William Valdar
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Richard Mott
- Genetics Institute, University College London, London, U.K
| | - Leah C. Solberg Woods
- Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC
| |
Collapse
|
3
|
Fruciano C, Colangelo P, Castiglia R, Franchini P. Does divergence from normal patterns of integration increase as chromosomal fusions increase in number? A test on a house mouse hybrid zone. Curr Zool 2020; 66:527-538. [PMID: 33293931 PMCID: PMC7705516 DOI: 10.1093/cz/zoaa035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2020] [Accepted: 07/03/2020] [Indexed: 11/14/2022] Open
Abstract
Chromosomal evolution is widely considered an important driver of speciation because it can promote the establishment of reproductive barriers. Karyotypic reorganization is also expected to affect the mean phenotype, as well as its development and patterns of phenotypic integration, through processes such as variation in genetic linkage between quantitative trait loci or between regulatory regions and their targets. Here we explore the relationship between chromosomal evolution and phenotypic integration by analyzing a well-known house mouse parapatric contact zone between a highly derived Robertsonian (Rb) race (2n = 22) and populations with standard karyotype (2n = 40). Populations with hybrid karyotypes are scattered throughout the hybrid zone connecting the two parental races. Using mandible shape data and geometric morphometrics, we test the hypothesis that patterns of integration progressively diverge from the “normal” integration pattern observed in the standard race as they accumulate Rb fusions. We find that the main pattern of integration observed between the posterior and anterior part of the mandible can be largely attributed to allometry. We find no support for a gradual increase in divergence from normal patterns of integration as fusions accumulate. Surprisingly, however, we find that the derived Rb race (2n = 22) has a distinct allometric trajectory compared with the standard race. Our results suggest that either individual fusions disproportionately affect patterns of integration or that there are mechanisms which “purge” extreme variants in hybrids (e.g. reduced fitness of hybrid shape).
Collapse
Affiliation(s)
- Carmelo Fruciano
- Institut de Biologie de l'Ecole Normale Supérieure (IBENS), Ecole Normale Supérieure, CNRS, PSL Université Paris, Paris, 75005, France.,School of Biological Sciences, University of Portsmouth, Portsmouth, PO1 2DY, UK
| | - Paolo Colangelo
- National Research Council, Research Institute on Terrestrial Ecosystems, Montelibretti (RM), 00010, Italy
| | - Riccardo Castiglia
- Department of Biology and Biotechnology "Charles Darwin", "La Sapienza" University of Rome, Rome, 00161, Italy
| | - Paolo Franchini
- Department of Biology, Lehrstuhl für Zoologie und Evolutionsbiologie, University of Konstanz, Konstanz, 78457, Germany
| |
Collapse
|
4
|
Scott MF, Ladejobi O, Amer S, Bentley AR, Biernaskie J, Boden SA, Clark M, Dell'Acqua M, Dixon LE, Filippi CV, Fradgley N, Gardner KA, Mackay IJ, O'Sullivan D, Percival-Alwyn L, Roorkiwal M, Singh RK, Thudi M, Varshney RK, Venturini L, Whan A, Cockram J, Mott R. Multi-parent populations in crops: a toolbox integrating genomics and genetic mapping with breeding. Heredity (Edinb) 2020; 125:396-416. [PMID: 32616877 PMCID: PMC7784848 DOI: 10.1038/s41437-020-0336-6] [Citation(s) in RCA: 68] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Revised: 06/16/2020] [Accepted: 06/16/2020] [Indexed: 11/21/2022] Open
Abstract
Crop populations derived from experimental crosses enable the genetic dissection of complex traits and support modern plant breeding. Among these, multi-parent populations now play a central role. By mixing and recombining the genomes of multiple founders, multi-parent populations combine many commonly sought beneficial properties of genetic mapping populations. For example, they have high power and resolution for mapping quantitative trait loci, high genetic diversity and minimal population structure. Many multi-parent populations have been constructed in crop species, and their inbred germplasm and associated phenotypic and genotypic data serve as enduring resources. Their utility has grown from being a tool for mapping quantitative trait loci to a means of providing germplasm for breeding programmes. Genomics approaches, including de novo genome assemblies and gene annotations for the population founders, have allowed the imputation of rich sequence information into the descendent population, expanding the breadth of research and breeding applications of multi-parent populations. Here, we report recent successes from crop multi-parent populations in crops. We also propose an ideal genotypic, phenotypic and germplasm 'package' that multi-parent populations should feature to optimise their use as powerful community resources for crop research, development and breeding.
Collapse
Affiliation(s)
| | | | - Samer Amer
- University of Reading, Reading, RG6 6AH, UK
- Faculty of Agriculture, Alexandria University, Alexandria, 23714, Egypt
| | - Alison R Bentley
- The John Bingham Laboratory, NIAB, 93 Lawrence Weaver Road, Cambridge, CB3 0LE, UK
| | - Jay Biernaskie
- Department of Plant Sciences, University of Oxford, South Parks Road, Oxford, OX1 3RB, UK
| | - Scott A Boden
- School of Agriculture, Food and Wine, University of Adelaide, Glen Osmond, SA, 5064, Australia
| | | | | | - Laura E Dixon
- Faculty of Biological Sciences, University of Leeds, Leeds, LS2 9JT, UK
| | - Carla V Filippi
- Instituto de Agrobiotecnología y Biología Molecular (IABIMO), INTA-CONICET, Nicolas Repetto y Los Reseros s/n, 1686, Hurlingham, Buenos Aires, Argentina
| | - Nick Fradgley
- The John Bingham Laboratory, NIAB, 93 Lawrence Weaver Road, Cambridge, CB3 0LE, UK
| | - Keith A Gardner
- The John Bingham Laboratory, NIAB, 93 Lawrence Weaver Road, Cambridge, CB3 0LE, UK
| | - Ian J Mackay
- SRUC, West Mains Road, Kings Buildings, Edinburgh, EH9 3JG, UK
| | | | | | - Manish Roorkiwal
- Center of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad, India
| | - Rakesh Kumar Singh
- International Center for Biosaline Agriculture, Academic City, Dubai, United Arab Emirates
| | - Mahendar Thudi
- Center of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad, India
| | - Rajeev Kumar Varshney
- Center of Excellence in Genomics and Systems Biology, International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Hyderabad, India
| | | | - Alex Whan
- CSIRO, GPO Box 1700, Canberra, ACT, 2601, Australia
| | - James Cockram
- The John Bingham Laboratory, NIAB, 93 Lawrence Weaver Road, Cambridge, CB3 0LE, UK
| | - Richard Mott
- UCL Genetics Institute, Gower Street, London, WC1E 6BT, UK
| |
Collapse
|
5
|
Luan MW, Zhang XM, Zhu ZB, Chen Y, Xie SQ. Evaluating Structural Variation Detection Tools for Long-Read Sequencing Datasets in Saccharomyces cerevisiae. Front Genet 2020; 11:159. [PMID: 32211024 PMCID: PMC7075250 DOI: 10.3389/fgene.2020.00159] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 02/11/2020] [Indexed: 01/01/2023] Open
Abstract
Structural variation (SV) represents a major form of genetic variations that contribute to polymorphic variations, human diseases, and phenotypes in many organisms. Long-read sequencing has been successfully used to identify novel and complex SVs. However, comparison of SV detection tools for long-read sequencing datasets has not been reported. Therefore, we developed an analysis workflow that combined two alignment tools (NGMLR and minimap2) and five callers (Sniffles, Picky, smartie-sv, PBHoney, and NanoSV) to evaluate the SV detection in six datasets of Saccharomyces cerevisiae. The accuracy of SV regions was validated by re-aligning raw reads in diverse alignment tools, SV callers, experimental conditions, and sequencing platforms. The results showed that SV detection between NGMLR and minimap2 was not significant when using the same caller. The PBHoney was with the highest average accuracy (89.04%) and Picky has the lowest average accuracy (35.85%). The accuracy of NanoSV, Sniffles, and smartie-sv was 68.67%, 60.47%, and 57.67%, respectively. In addition, smartie-sv and NanoSV detected the most and least number of SVs, and SV detection from the PacBio sequencing platform was significantly more than that from ONT (p = 0.000173).
Collapse
Affiliation(s)
- Mei-Wei Luan
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants (Ministry of Education), Hainan Key Laboratory for Biology of Tropical Ornamental Plant Germplasm, College of Forestry, Hainan University, Haikou, China
| | - Xiao-Ming Zhang
- College of Grassland, Resources and Environment, Inner Mongolia Agricultural University, Huhhot, China
| | - Zi-Bin Zhu
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants (Ministry of Education), Hainan Key Laboratory for Biology of Tropical Ornamental Plant Germplasm, College of Forestry, Hainan University, Haikou, China
| | - Ying Chen
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China
| | - Shang-Qian Xie
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants (Ministry of Education), Hainan Key Laboratory for Biology of Tropical Ornamental Plant Germplasm, College of Forestry, Hainan University, Haikou, China
| |
Collapse
|
6
|
Goel M, Sun H, Jiao WB, Schneeberger K. SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies. Genome Biol 2019; 20:277. [PMID: 31842948 PMCID: PMC6913012 DOI: 10.1186/s13059-019-1911-0] [Citation(s) in RCA: 251] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2019] [Accepted: 12/02/2019] [Indexed: 01/27/2023] Open
Abstract
Genomic differences range from single nucleotide differences to complex structural variations. Current methods typically annotate sequence differences ranging from SNPs to large indels accurately but do not unravel the full complexity of structural rearrangements, including inversions, translocations, and duplications, where highly similar sequence changes in location, orientation, or copy number. Here, we present SyRI, a pairwise whole-genome comparison tool for chromosome-level assemblies. SyRI starts by finding rearranged regions and then searches for differences in the sequences, which are distinguished for residing in syntenic or rearranged regions. This distinction is important as rearranged regions are inherited differently compared to syntenic regions.
Collapse
Affiliation(s)
- Manish Goel
- Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany
| | - Hequan Sun
- Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany
| | - Wen-Biao Jiao
- Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany
| | - Korbinian Schneeberger
- Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany
- Faculty of Biology, LMU Munich, 82152 Planegg-Martinsried, Germany
| |
Collapse
|
7
|
Components of the ribosome biogenesis pathway underlie establishment of telomere length set point in Arabidopsis. Nat Commun 2019; 10:5479. [PMID: 31792215 PMCID: PMC6889149 DOI: 10.1038/s41467-019-13448-z] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 11/08/2019] [Indexed: 12/20/2022] Open
Abstract
Telomeres cap the physical ends of eukaryotic chromosomes to ensure complete DNA replication and genome stability. Heritable natural variation in telomere length exists in yeast, mice, plants and humans at birth; however, major effect loci underlying such polymorphism remain elusive. Here, we employ quantitative trait locus (QTL) mapping and transgenic manipulations to identify genes controlling telomere length set point in a multi-parent Arabidopsis thaliana mapping population. We detect several QTL explaining 63.7% of the total telomere length variation in the Arabidopsis MAGIC population. Loss-of-function mutants of the NOP2A candidate gene located inside the largest effect QTL and of two other ribosomal genes RPL5A and RPL5B establish a shorter telomere length set point than wild type. These findings indicate that evolutionarily conserved components of ribosome biogenesis and cell proliferation pathways promote telomere elongation. Major effect loci controlling natural, heritable variation in telomere length are not known. Here, the authors use QTL mapping and transgenic manipulations in Arabidopsis to implicate the rRNA-processing genes NOP2A and RPL5 in telomere length set point regulation in this model species.
Collapse
|
8
|
Lyra DH, Galli G, Alves FC, Granato ÍSC, Vidotti MS, Bandeira E Sousa M, Morosini JS, Crossa J, Fritsche-Neto R. Modeling copy number variation in the genomic prediction of maize hybrids. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2019; 132:273-288. [PMID: 30382311 DOI: 10.1007/s00122-018-3215-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2018] [Accepted: 10/20/2018] [Indexed: 06/08/2023]
Abstract
Our study indicates that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids. Moreover, predicting hybrid phenotypes by combining additive-dominance effects with copy variants has the potential to be a viable predictive model. Non-additive effects resulting from the actions of multiple loci may influence trait variation in single-cross hybrids. In addition, complementation of allelic variation could be a valuable contributor to hybrid genetic variation, especially when crossing inbred lines with higher contents of copy gains. With this in mind, we aimed (1) to study the association between copy number variation (CNV) and hybrid phenotype, and (2) to compare the predictive ability (PA) of additive and additive-dominance genomic best linear unbiased prediction model when combined with the effects of CNV in two datasets of maize hybrids (USP and HELIX). In the USP dataset, we observed a significant negative phenotypic correlation of low magnitude between copy number loss and plant height, revealing a tendency that more copy losses lead to lower plants. In the same set, when CNV was combined with the additive plus dominance effects, the PA significantly increased only for plant height under low nitrogen. In this case, CNV effects explicitly capture relatedness between individuals and add extra information to the model. In the HELIX dataset, we observed a pronounced difference in PA between additive (0.50) and additive-dominance (0.71) models for predicting grain yield, suggesting a significant contribution of dominance. We conclude that copy variants may play an essential role in the phenotypic variation of complex traits in maize hybrids, although the inclusion of CNVs into datasets does not return significant gains concerning PA.
Collapse
Affiliation(s)
- Danilo Hottis Lyra
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil.
- Department of Computational and Analytical Sciences, Rothamsted Research, West Common, Harpenden, AL52JQ, UK.
| | - Giovanni Galli
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| | - Filipe Couto Alves
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| | - Ítalo Stefanine Correia Granato
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| | - Miriam Suzane Vidotti
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| | - Massaine Bandeira E Sousa
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| | - Júlia Silva Morosini
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| | - José Crossa
- Biometrics and Statistics Unit, International Maize and Wheat Improvement Center (CIMMYT), 06600, Texcoco, D.F, Mexico
| | - Roberto Fritsche-Neto
- Department of Genetics, Luiz de Queiroz College of Agriculture, University of São Paulo (ESALQ/USP), Piracicaba, São Paulo, Brazil
| |
Collapse
|
9
|
Menard GN, Bryant FM, Kelly AA, Craddock CP, Lavagi I, Hassani-Pak K, Kurup S, Eastmond PJ. Natural variation in acyl editing is a determinant of seed storage oil composition. Sci Rep 2018; 8:17346. [PMID: 30478395 PMCID: PMC6255774 DOI: 10.1038/s41598-018-35136-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2018] [Accepted: 10/26/2018] [Indexed: 01/09/2023] Open
Abstract
Seeds exhibit wide variation in the fatty acid composition of their storage oil. However, the genetic basis of this variation is only partially understood. Here we have used a multi-parent advanced generation inter-cross (MAGIC) population to study the genetic control of fatty acid chain length in Arabidopsis thaliana seed oil. We mapped four quantitative trait loci (QTL) for the quantity of the major very long chain fatty acid species 11-eicosenoic acid (20:1), using multiple QTL modelling. Surprisingly, the main-effect QTL does not coincide with FATTY ACID ELONGASE 1 and a parallel genome wide association study suggested that LYSOPHOSPHATIDYLCHOLINE ACYLTRANSFERASE 2 (LPCAT2) is a candidate for this QTL. Regression analysis also suggested that LPCAT2 expression and 20:1 content in seeds of the 19 MAGIC founder accessions are related. LPCAT is a key component of the Lands cycle; an acyl editing pathway that enables acyl-exchange between the acyl-Coenzyme A and phosphatidylcholine precursor pools used for microsomal fatty acid elongation and desaturation, respectively. We Mendelianised the main-effect QTL using biparental chromosome segment substitution lines and carried out complementation tests to show that a single cis-acting polymorphism in the LPCAT2 promoter causes the variation in seed 20:1 content, by altering the LPCAT2 expression level and total LPCAT activity in developing siliques. Our work establishes that oilseed species exhibit natural variation in the enzymic capacity for acyl editing and this contributes to the genetic control of storage oil composition.
Collapse
Affiliation(s)
- Guillaume N Menard
- Department of Plant Science, Rothamsted Research, Harpenden, Hertfordshire, AL5 2JQ, UK
| | - Fiona M Bryant
- Department of Plant Science, Rothamsted Research, Harpenden, Hertfordshire, AL5 2JQ, UK
| | - Amélie A Kelly
- Georg-August-University, Albrecht-von-Haller-Institute for Plant Sciences, Justus-von-Liebig Weg 11, 37077, Göttingen, Germany
| | - Christian P Craddock
- Mt. San Jacinto College, Menifee Valley Campus, 28237 La Piedra Road, Menifee, CA, 92584, USA
| | - Irene Lavagi
- Department of Microbiology and Plant Pathology, University of California Riverside, Riverside, CA, 92521, USA
| | - Keywan Hassani-Pak
- Department of Plant Science, Rothamsted Research, Harpenden, Hertfordshire, AL5 2JQ, UK
| | - Smita Kurup
- Department of Plant Science, Rothamsted Research, Harpenden, Hertfordshire, AL5 2JQ, UK
| | - Peter J Eastmond
- Department of Plant Science, Rothamsted Research, Harpenden, Hertfordshire, AL5 2JQ, UK.
| |
Collapse
|
10
|
Perlaza-Jiménez L, Walther D. A genome-wide scan for correlated mutations detects macromolecular and chromatin interactions in Arabidopsis thaliana. Nucleic Acids Res 2018; 46:8114-8132. [PMID: 29986106 PMCID: PMC6144803 DOI: 10.1093/nar/gky576] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2018] [Accepted: 06/14/2018] [Indexed: 01/05/2023] Open
Abstract
The concept of exploiting correlated mutations has been introduced and applied successfully to identify interactions within and between biological macromolecules. Its rationale lies in the preservation of physical interactions via compensatory mutations. With the massive increase of available sequence information, approaches based on correlated mutations have regained considerable attention. We analyzed a set of 10 707 430 single nucleotide polymorphisms detected in 1135 accessions of the plant Arabidopsis thaliana. To measure their covariance and to reveal the global genome-wide sequence correlation structure of the Arabidopsis genome, the adjusted mutual information has been estimated for each possible pair of polymorphic sites. We developed a series of filtering steps to account for genetic linkage and lineage relations between Arabidopsis accessions, as well as transitive covariance as possible confounding factors. We show that upon appropriate filtering, correlated mutations prove indeed informative with regard to molecular interactions, and furthermore, appear to reflect on chromosomal interactions. Our study demonstrates that the concept of correlated mutations can also be applied successfully to within-species sequence variation and establishes a promising approach to help unravel the complex molecular interactions in A. thaliana and other species with broad sequence information.
Collapse
Affiliation(s)
- Laura Perlaza-Jiménez
- Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, 14476 Potsdam-Golm, Germany
| | - Dirk Walther
- Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, 14476 Potsdam-Golm, Germany
| |
Collapse
|
11
|
Sedlazeck FJ, Rescheneder P, Smolka M, Fang H, Nattestad M, von Haeseler A, Schatz MC. Accurate detection of complex structural variations using single-molecule sequencing. Nat Methods 2018; 15:461-468. [PMID: 29713083 PMCID: PMC5990442 DOI: 10.1038/s41592-018-0001-7] [Citation(s) in RCA: 836] [Impact Index Per Article: 139.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2017] [Accepted: 03/16/2018] [Indexed: 02/08/2023]
Abstract
Structural variations are the greatest source of genetic variation, but they remain poorly understood because of technological limitations. Single-molecule long-read sequencing has the potential to dramatically advance the field, although high error rates are a challenge with existing methods. Addressing this need, we introduce open-source methods for long-read alignment (NGMLR; https://github.com/philres/ngmlr ) and structural variant identification (Sniffles; https://github.com/fritzsedlazeck/Sniffles ) that provide unprecedented sensitivity and precision for variant detection, even in repeat-rich regions and for complex nested events that can have substantial effects on human health. In several long-read datasets, including healthy and cancerous human genomes, we discovered thousands of novel variants and categorized systematic errors in short-read approaches. NGMLR and Sniffles can automatically filter false events and operate on low-coverage data, thereby reducing the high costs that have hindered the application of long reads in clinical and research settings.
Collapse
Affiliation(s)
- Fritz J Sedlazeck
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA.
| | - Philipp Rescheneder
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna, Medical University of Vienna, Vienna, Austria
| | - Moritz Smolka
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna, Medical University of Vienna, Vienna, Austria
| | - Han Fang
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Maria Nattestad
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Arndt von Haeseler
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna, Medical University of Vienna, Vienna, Austria
- Bioinformatics and Computational Biology, Faculty of Computer Science, University of Vienna, Vienna, Austria
| | - Michael C Schatz
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA.
- Departments of Computer Science and Biology, Johns Hopkins University, Baltimore, MD, USA.
| |
Collapse
|
12
|
Rabanal FA, Mandáková T, Soto-Jiménez LM, Greenhalgh R, Parrott DL, Lutzmayer S, Steffen JG, Nizhynska V, Mott R, Lysak MA, Clark RM, Nordborg M. Epistatic and allelic interactions control expression of ribosomal RNA gene clusters in Arabidopsis thaliana. Genome Biol 2017; 18:75. [PMID: 28464948 PMCID: PMC5414317 DOI: 10.1186/s13059-017-1209-z] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2016] [Accepted: 04/06/2017] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND Ribosomal RNA (rRNA) accounts for the majority of the RNA in eukaryotic cells, and is encoded by hundreds to thousands of nearly identical gene copies, only a subset of which are active at any given time. In Arabidopsis thaliana, 45S rRNA genes are found in two large ribosomal DNA (rDNA) clusters and little is known about the contribution of each to the overall transcription pattern in the species. RESULTS By taking advantage of genome sequencing data from the 1001 Genomes Consortium, we characterize rRNA gene sequence variation within and among accessions. Notably, variation is not restricted to the pre-rRNA sequences removed during processing, but it is also present within the highly conserved ribosomal subunits. Through linkage mapping we assign these variants to a particular rDNA cluster unambiguously and use them as reporters of rDNA cluster-specific expression. We demonstrate that rDNA cluster-usage varies greatly among accessions and that rDNA cluster-specific expression and silencing is controlled via genetic interactions between entire rDNA cluster haplotypes (alleles). CONCLUSIONS We show that rRNA gene cluster expression is controlled via complex epistatic and allelic interactions between rDNA haplotypes that apparently regulate the entire rRNA gene cluster. Furthermore, the sequence polymorphism we discovered implies that the pool of rRNA in a cell may be heterogeneous, which could have functional consequences.
Collapse
Affiliation(s)
- Fernando A Rabanal
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria.
| | - Terezie Mandáková
- Central European Institute of Technology (CEITEC), Masaryk University, Brno, Czech Republic
| | - Luz M Soto-Jiménez
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria
| | | | - David L Parrott
- Department of Biology, University of Utah, Salt Lake City, UT, USA
| | - Stefan Lutzmayer
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria
| | - Joshua G Steffen
- Department of Natural Sciences, Colby-Sawyer College, New London, NH, USA
| | - Viktoria Nizhynska
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria
| | - Richard Mott
- Genetics Institute, University College London (UCL), Gower Street, London, WC1E 6BT, UK
| | - Martin A Lysak
- Central European Institute of Technology (CEITEC), Masaryk University, Brno, Czech Republic
| | - Richard M Clark
- Department of Biology, University of Utah, Salt Lake City, UT, USA
- Center for Cell and Genome Science, University of Utah, Salt Lake City, UT, USA
| | - Magnus Nordborg
- Gregor Mendel Institute (GMI), Austrian Academy of Sciences, Vienna Biocenter (VBC), Dr. Bohr-Gasse 3, 1030, Vienna, Austria.
| |
Collapse
|
13
|
Rabanal FA, Nizhynska V, Mandáková T, Novikova PY, Lysak MA, Mott R, Nordborg M. Unstable Inheritance of 45S rRNA Genes in Arabidopsis thaliana. G3 (BETHESDA, MD.) 2017; 7:1201-1209. [PMID: 28188182 PMCID: PMC5386868 DOI: 10.1534/g3.117.040204] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/16/2016] [Accepted: 02/07/2017] [Indexed: 12/30/2022]
Abstract
The considerable genome size variation in Arabidopsis thaliana has been shown largely to be due to copy number variation (CNV) in 45S ribosomal RNA (rRNA) genes. Surprisingly, attempts to map this variation by means of genome-wide association studies (GWAS) failed to identify either of the two likely sources, namely the nucleolus organizer regions (NORs). Instead, GWAS implicated a trans-acting locus, as if rRNA gene CNV was a phenotype rather than a genotype. To explain these results, we investigated the inheritance and stability of rRNA gene copy number using the variety of genetic resources available in A. thaliana - F2 crosses, recombinant inbred lines, the multiparent advanced-generation inter-cross population, and mutation accumulation lines. Our results clearly show that rRNA gene CNV can be mapped to the NORs themselves, with both loci contributing equally to the variation. However, NOR size is unstably inherited, and dramatic copy number changes are visible already within tens of generations, which explains why it is not possible to map the NORs using GWAS. We did not find any evidence of trans-acting loci in crosses, which is also expected since changes due to such loci would take very many generations to manifest themselves. rRNA gene copy number is thus an interesting example of "missing heritability"-a trait that is heritable in pedigrees, but not in the general population.
Collapse
Affiliation(s)
- Fernando A Rabanal
- Gregor Mendel Institute, Austrian Academy of Sciences, Vienna Biocenter (VBC), 1030, Austria
| | - Viktoria Nizhynska
- Gregor Mendel Institute, Austrian Academy of Sciences, Vienna Biocenter (VBC), 1030, Austria
| | - Terezie Mandáková
- Central European Institute of Technology, Masaryk University, 625 00 Brno, Czech Republic
| | - Polina Yu Novikova
- Gregor Mendel Institute, Austrian Academy of Sciences, Vienna Biocenter (VBC), 1030, Austria
| | - Martin A Lysak
- Central European Institute of Technology, Masaryk University, 625 00 Brno, Czech Republic
| | - Richard Mott
- University College London Genetics Institute, WC1E 6BT, UK
| | - Magnus Nordborg
- Gregor Mendel Institute, Austrian Academy of Sciences, Vienna Biocenter (VBC), 1030, Austria
| |
Collapse
|
14
|
Menard GN, Moreno JM, Bryant FM, Munoz-Azcarate O, Kelly AA, Hassani-Pak K, Kurup S, Eastmond PJ. Genome Wide Analysis of Fatty Acid Desaturation and Its Response to Temperature. PLANT PHYSIOLOGY 2017; 173:1594-1605. [PMID: 28108698 PMCID: PMC5338679 DOI: 10.1104/pp.16.01907] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Accepted: 01/20/2017] [Indexed: 05/22/2023]
Abstract
Plants modify the polyunsaturated fatty acid content of their membrane and storage lipids in order to adapt to changes in temperature. In developing seeds, this response is largely controlled by the activities of the microsomal ω-6 and ω-3 fatty acid desaturases, FAD2 and FAD3. Although temperature regulation of desaturation has been studied at the molecular and biochemical levels, the genetic control of this trait is poorly understood. Here, we have characterized the response of Arabidopsis (Arabidopsis thaliana) seed lipids to variation in ambient temperature and found that heat inhibits both ω-6 and ω-3 desaturation in phosphatidylcholine, leading to a proportional change in triacylglycerol composition. Analysis of the 19 parental accessions of the multiparent advanced generation intercross (MAGIC) population showed that significant natural variation exists in the temperature responsiveness of ω-6 desaturation. A combination of quantitative trait locus (QTL) analysis and genome-wide association studies (GWAS) using the MAGIC population suggests that ω-6 desaturation is largely controlled by cis-acting sequence variants in the FAD2 5' untranslated region intron that determine the expression level of the gene. However, the temperature responsiveness of ω-6 desaturation is controlled by a separate QTL on chromosome 2. The identity of this locus is unknown, but genome-wide association studies identified potentially causal sequence variants within ∼40 genes in an ∼450-kb region of the QTL.
Collapse
Affiliation(s)
- Guillaume N Menard
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Jose Martin Moreno
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Fiona M Bryant
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Olaya Munoz-Azcarate
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Amélie A Kelly
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Keywan Hassani-Pak
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Smita Kurup
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| | - Peter J Eastmond
- Department of Plant Biology and Crop Science, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, United Kingdom (G.N.M., J.M.M., F.M.B., O.M.-A., K.H.-P., S.K., P.J.E.); and
- School of Life Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom (A.A.K.)
| |
Collapse
|