Students enrolled in AIOU’s Code 419 educational program in the spring semester face a number of challenges that require careful consideration and application of educational theory. These assignments are designed to assess students’ understanding of various educational concepts, teaching methods, and their ability to critically review relevant literature.
Addressing these challenges not only contributes to overall academic performance but also increases students’ practical knowledge, preparing them for the challenges they face in the field of education. Topics covered in the spring semester may include curriculum design, educational psychology, teaching strategies, and assessment methods.
AIOU Course Code 419 Solved Assignments: The AIOU Course Code 419, titled Education (ایجوکیشن), is an extensive 6-credit hour course comprising four assignments. This course delves into various educational theories, practices, and policies. It aims to provide students with a thorough understanding of the education system, teaching methodologies, curriculum development, and the role of education in society. By studying this course, students will gain valuable insights into how educational processes can be improved and tailored to meet diverse learning needs.
At Asif Brain Academy, we offer a wide range of resources to assist students in mastering Course Code 419. Our comprehensive support includes solved assignments that simplify complex educational concepts, past papers that help students understand the exam patterns, and guess papers that offer focused preparation for exams. Additionally, we provide access to essential educational materials such as date sheets, notes, and textbooks.
ASIFBRAINACADEMY provide latest educational information.
Save my name, email, and website in this browser for the next time I comment.
Welcome to aiou master acadmey, aiou(علم التعلیم ) 419 solved assignment semester spring 2024 pdf.
aiou(علم التعلیم ) 419 solved assignment semester spring 2024 pdf
Assignment no# 1
Click to download
Assignment no# 2
Assignment no# 3
Assignment no# 4
To download your assignment, simply click the download button above, and the PDF will be saved to your computer or mobile device within seconds.
Remember, submitting assignments is vital for academic success as it showcases your comprehension of the course material and enables instructors to evaluate your progress accurately. Failure to submit assignments can significantly impact your overall grade and may even lead to exam failure. To avoid this, it's essential to plan and prioritize effectively, setting realistic deadlines and allocating sufficient time for each task. Additionally, don't hesitate to seek assistance from your instructor or classmates if you encounter difficulties with an assignment, as timely support can often resolve challenges and enhance learning outcomes. Furthermore, maintaining open communication with your instructor is crucial; if you face any obstacles or unforeseen circumstances that may affect your ability to submit assignments on time, be proactive in notifying your instructor and seeking guidance on possible accommodations or extensions.
Post a comment, popular posts, 436 (: سیرت طیبہ )solved assignment spring 2024 pdf july 07, 2024, 437( اسلامیات اختیاری ) solved assignment spring 2024 pdf july 07, 2024, 426(پاکستانی ادب ) solved assignment spring 2024 pdf july 07, 2024, 484 solved assignment spring 2024 pdf. july 07, 2024.
Download and view past papers for Education code number 419 (AIOU) for the semester of spring and autumn. It is a 6 credit hours book.Old papers will provide a clear understanding of important areas or parts of the subjects. Also previous examinations papers will provide a guess for the upcoming examination.Below table provide past papers. Click on the link to view:
Serial No | Year/Semester | View |
---|---|---|
1 | ||
2 | ||
3 | ||
4 | ||
5 |
Serial No | Year/Semester | View |
---|---|---|
1 | ||
2 |
Serial No | Name | View |
---|---|---|
1 |
Home » AIOU » Assignments » Bachelor » 419 Education Bachelor Free Solved Assignments
Download 419 education bachelor free solved assignments semester spring 2024.
Firstly welcome to download page 419 Education free for the Bachelor class. That is to say, that you can download 419 totally free for the semester Spring 2024.
Firstly we will guide you through the download process. So, you could download 419 Education Assignments for Bachelor class in the semester Spring 2024. Allama Iqbal open university has more than four hundred codes for all classes, but all codes have specific subjects.
That is to say, that It is quite easy to download 419 Education free solved assignments from aiou studio 9. Certainly, you need to scroll all the way down and click the download button it will open a page where you could see all codes select one, and click download pdf. We also have prepared a video for you have a great look for any confusion do comment us.
https://www.youtube.com/watch?v=Re-7KnV-zOQ&t=21s
Make sure, you keep yourself up to date, with 419 Education assignments submission dates, so that you submit them on a time. So, you make sure to submit these assignments using the date. Consequently exceeding these dates could cause failure in assignments, in other words, failure in papers. So dear, don’t wait for the last dates and submit our assignments at your earliest.
We know students do take a lot of tension for their papers that is why aiou studio 9 is presenting 419 free guess papers. You can download Education 419 guess papers for the semester Spring 2024.
Download Guess Papers
Please click the below button and it will redirect you to the download page after that you can download Education 419 Free solved Assignments.
Matric Free Solved Assignments Download
New Date sheet – Updated
BA Free solved Assignments Download – New
FA I Com Free Solved Assignments Download
Guess Papers Download Free
AiouStudio9.com
Allama Iqbal Open University has a system of Assignments. One must have to send Assignments the Hand-Written assignments to the AIOU Tutor assigned to him/her.
Now you can download all updated AIOU Solved Assignments for BA, BCom in PDF for semester Spring and Autumn.
Read Also: Check AIOU Tutor Address Online by Roll Number
Allama Iqbal Open University now providing the facility of PDF books in Soft Form for its students. If one has not received its books then he/she can download the AIOU books in PDF .
As we above mentioned that it is necessary to submit assignments because of AIOU count Assignments Marks in Final Result . Submission of Assignments is necessary and students do not have proper guides to solve the questions. That’s why they look for the pre-solved assignments of AIOU . Or, they contact in contact to get pre-hand-written assignments.
Let’s dive in the AIOU Solved Assignments B.A PDF.
Follow us on our Facebook Page or Join our Facebook group for AIOU updates
Follow us on Twitter
402 | Economics | ||||
404 | Urdu | ||||
405 | Iqbaliat | ||||
406 | Economics of Pakistan | ||||
407 | Modern Muslim World | ||||
408 | Arabic | ||||
409 | Commercial Geography | ||||
411 | Sociology – I | ||||
412 | Social & Cultural Anthropology | ||||
414 | Demography | ||||
419 | Education | ||||
422 | Organizing Library Resources | ||||
423 | Library Services | ||||
426 | Pakistani Adab – II | – | |||
429 | Mass Education / Literacy | ||||
430 | Principles of Journalism | ||||
436 | Seerat – e – Tayyaba | ||||
437 | Islamiat (E) | ||||
445 | History of Urdu Adab | – | – | – | |
458 | Community Development | ||||
464 | Islamic Fiqh | ||||
465 | Population and Development | – | – | – | – |
466 | History of Libraries with Reference to Pakistan | ||||
467 | Classification and Cataloguing | ||||
472 | Quran – E – Hakim | ||||
473 | Hadith |
And, here is the complete details about deadline of the submission of these assignments.
Read Also: Download FREE AIOU Solved Assignments in PDF
Assignment No 1 | |
Assignment No 2 | |
Assignment No 3 | 15-03-2021 |
Assignment No 4 | 15-04-2021 |
413 | Sociology – II | ||
416 | Islamiyat | ||
417 | Pakistan Studies | ||
427 | Pakistani – Adab – I | ||
431 | Reporting | ||
434 | Iqbal’s Urdu Prose | – | |
438 | Principles of Accounting | ||
439 | Advertising and Sales Promotion | – | |
444 | Advance Accounting | ||
447 | Marketing Management | – | |
449 | Magazine Journalism | ||
451 | Public Relations | ||
452 | Mass Communication | ||
453 | Radio Broadcasting | ||
454 | TV Broadcasting | ||
455 | Book Editing | ||
456 | Business Taxation | ||
460 | Mercantile Law | ||
461 | Advertising | – | |
462 | Cost Accounting | ||
463 | Fundamentals of Business | ||
470 | Principles of Marketing | ||
481 | Auditing | ||
482 | Food Microbiology | ||
484 | Food and Nutrition | ||
485 | Health and Nutrition | ||
487 | Child Development | ||
1414 | Fundamental of Money and Banking | – | |
1416 | Business Communication | – | |
1421 | Introduction to Environment | ||
1422 | Environmental Pollution | ||
1423 | English – I | ||
1424 | English – II | ||
1425 | Basics of Technical Writing | ||
1426 | English Literature | ||
1428 | Commercial Geography | – | |
1429 | Business Mathematics | ||
1430 | Business Statistics | ||
1431 | Basics of ICT (English) | ||
1431 | Basics of ICT (Urdu) |
Assignment No 1 | |
Assignment No 2 | 15-04-2021 |
These all books are present in PDF, you can easily avail these in your any device which can open the PDF file.
What is the last date of aiou ba assignments submission.
Allama Iqbal Open University has extended the date of Assignments. You can submit your first three assignments by March 15, 2021 . And, the last date for all assignments submission is April 15, 2021 .
Well, it is pretty much easy to check AIOU tutor online. You just need your Roll Number to check your assigned tutor. Open this link and here enter your Roll Number to check your Tutor Information online.
Check your AIOU BA result online by entering your Roll Number. Open this link to check your AIOU BA result online. AIOU has successfully announced BA results. Moreover, you can check the AIOU BA result of past years from 2002 to 2020 also by entering your Roll Number.
You can easily ask your questions in the comments below. Our team will try their best to promptly help you and will answer your questions as soon as possible
Aiou previous result of all years | 2022 to 2021 history of result, aiou extends the last date of admissions (spring 2022), aiou admission 2022 apply online: semester spring & autumn, [2022] download aiou books pdf (all programs), download aiou roll no slips 2022 • all programs, 158 comments.
Plz send me the 1416 assessment semester 3
AOA… Respected Sir… Kindly Sir/Medam upload B.com 1415&1426 assignment as soon as possible…
Save my name, email, and website in this browser for the next time I comment.
arXiv's Accessibility Forum starts next month!
Help | Advanced Search
Title: multi-agent target assignment and path finding for intelligent warehouse: a cooperative multi-agent deep reinforcement learning perspective.
Abstract: Multi-agent target assignment and path planning (TAPF) are two key problems in intelligent warehouse. However, most literature only addresses one of these two problems separately. In this study, we propose a method to simultaneously solve target assignment and path planning from a perspective of cooperative multi-agent deep reinforcement learning (RL). To the best of our knowledge, this is the first work to model the TAPF problem for intelligent warehouse to cooperative multi-agent deep RL, and the first to simultaneously address TAPF based on multi-agent deep RL. Furthermore, previous literature rarely considers the physical dynamics of agents. In this study, the physical dynamics of the agents is considered. Experimental results show that our method performs well in various task settings, which means that the target assignment is solved reasonably well and the planned path is almost shortest. Moreover, our method is more time-efficient than baselines.
Subjects: | Artificial Intelligence (cs.AI); Multiagent Systems (cs.MA) |
Cite as: | [cs.AI] |
(or [cs.AI] for this version) | |
Focus to learn more arXiv-issued DOI via DataCite |
Access paper:.
Code, data and media associated with this article, recommenders and search tools.
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs .
This site displays a prototype of a “Web 2.0” version of the daily Federal Register. It is not an official legal edition of the Federal Register, and does not replace the official print version or the official electronic version on GPO’s govinfo.gov.
The documents posted on this site are XML renditions of published Federal Register documents. Each document posted on the site includes a link to the corresponding official PDF file on govinfo.gov. This prototype edition of the daily Federal Register on FederalRegister.gov will remain an unofficial informational resource until the Administrative Committee of the Federal Register (ACFR) issues a regulation granting it official legal status. For complete information about, and access to, our official publications and services, go to About the Federal Register on NARA's archives.gov.
The OFR/GPO partnership is committed to presenting accurate and reliable regulatory information on FederalRegister.gov with the objective of establishing the XML-based Federal Register as an ACFR-sanctioned publication in the future. While every effort has been made to ensure that the material on FederalRegister.gov is accurately displayed, consistent with the official SGML-based PDF version on govinfo.gov, those relying on it for legal research should verify their results against an official edition of the Federal Register. Until the ACFR grants it official status, the XML rendition of the daily Federal Register on FederalRegister.gov does not provide legal notice to the public or judicial notice to the courts.
Design Updates: As part of our ongoing effort to make FederalRegister.gov more accessible and easier to use we've enlarged the space available to the document content and moved all document related data into the utility bar on the left of the document. Read more in our feature announcement .
A Rule by the Defense Department , the General Services Administration , and the National Aeronautics and Space Administration on 08/29/2024
This document has been published in the Federal Register . Use the PDF linked in the document sidebar for the official electronic format.
This table of contents is a navigational tool, processed from the headings within the legal text of Federal Register documents. This repetition of headings to form internal navigation links has no substantive legal effect.
Supplementary information:, list of subjects in 48 cfr parts 4 and 19, part 4—administrative and information matters, part 19—small business programs.
Comments are no longer being accepted. See DATES for details.
FederalRegister.gov retrieves relevant information about this document from Regulations.gov to provide users with additional context. This information is not part of the official Federal Register document.
This document is also available in the following formats:
More information and documentation can be found in our developer tools pages .
This PDF is the current document as it appeared on Public Inspection on 08/28/2024 at 8:45 am.
It was viewed 1 times while on Public Inspection.
If you are using public inspection listings for legal research, you should verify the contents of the documents against a final, official edition of the Federal Register. Only official editions of the Federal Register provide legal notice of publication to the public and judicial notice to the courts under 44 U.S.C. 1503 & 1507 . Learn more here .
Document headings vary by document type but may contain the following:
See the Document Drafting Handbook for more details.
General services administration, national aeronautics and space administration.
Department of Defense (DoD), General Services Administration (GSA), and National Aeronautics and Space Administration (NASA).
Final rule.
This document makes administrative changes to 48 CFR parts 4 and 19 . The date change is to provide additional time to implement the policy addressing the assignment of North American Industry Classification System codes to orders placed under multiple-award contracts, as covered by changes made by FAR Case 2014-002, Set-Asides Under Multiple-Award Contracts, 85 FR 11746 .
Effective: September 30, 2024.
Ms. Lois Mandell, Regulatory Secretariat Division (MVCB), at 202-501-4755 or [email protected] . Please cite FAC 2024-07, Technical Amendments.
This document makes editorial changes to 48 CFR parts 4 and 19 .
William F. Clark,
Director, Office of Government-Wide Acquisition Policy, Office of Acquisition Policy, Office of Government-Wide Policy.
Therefore, DoD, GSA, and NASA amend 48 CFR parts 4 and 19 as set forth below:
1. The authority citation for 48 CFR parts 4 and 19 continues to read as follows:
Authority: 40 U.S.C. 121(c) ; 10 U.S.C. chapter 4 and 10 U.S.C. chapter 137 legacy provisions (see 10 U.S.C. 3016 ); and 51 U.S.C. 20113 .
4.1202[Amended]
2. Amend section 4.1202 by removing from paragraph (a) introductory text the date “October 1, 2025” and adding the date “October 1, 2028” in its place.
3. Amend section 19.102 by removing from paragraphs (b)(2)(i) and (b)(2)(ii) introductory text the date “October 1, 2025” and adding the date “October 1, 2028” in their places, respectively.
4. Amend section 19.309 by removing from paragraphs (a)(3) and (c)(2) the date “October 1, 2025” and adding the date “October 1, 2028” in their places, respectively.
[ FR Doc. 2024-19174 Filed 8-28-24; 8:45 am]
BILLING CODE 6820-6EP-P
Information.
Aiou solved assignments 1 code 419 autumn 2022 | learning ki dunya.
Hi! Students are you looking for aiou 419 assignments 2021. Here all 419 assignment 1 in pdf are available. You can download aiou assignments.
Education code 419 assignment number 1.
Education code 419 assignment number 3, education code 419 assignment number 4.
First and foremost Assignments are an essential part of the students. Because this plays a very important role. Moreove r , you can’t clear your semester without assignments. Most Importantly send your assignments in time. However, you can not pass exams without assignments. In other words, your attendance in exams is wasteless.
You will be considered to fail without submitting your assignment. In other words, even if you get 95% marks in the examination you will be considered as fail. If you clear your assignments, you have the chance to appear in the examination. Therefore you can have 3 attempts to appear.
Following are the passing marks for each assignment
In contrast, subjects having 2 assignments follow the above same percentage.
Above All AIOU Master Academy strives to help AIOU students. Moreover, we provide exam preparation notes like Key books , Guess papers , and 5 years papers.
First of all our mission is to educate our young generation. Secondly Academy tries its best to build civilized and sophisticated mankind in our country. besides, we are not indulged in any kind of miscellaneous activity that detracts students to take any shortcut in their education. subsequently, we always promote legal education as a priority, as well as help students self-learning and build skills to perform best in the future.
Our basic earnings are from selling hardcopy notes like Key books , Guess papers , and 5 years’ Papers to students at a very reasonable price. The Academy doesn’t earn by providing any kind of consultancy services. We assist students on Call or WhatsApp and we have premium support like others to earn money from needy students.
We also provide free Guess Papers , 5 Years Papers , and Soft Copy of Assignments to students on our website. Assignments are only available on the website for guidelines and awareness about the conten t and pattern of related courses.
Despite this, our earnings are from the internet, YouTube , and Google . Similarly, They pay us because of the ads they publish on our YouTube channel or on our website.
Book editing code 455 assignments of autumn 2023, tv broadcasting code 454 assignments of autumn 2023, radio broadcasting code 453 assignments of autumn 2023, advance of accounting code 444 assignments of autumn 2023, leave a reply cancel reply.
Your email address will not be published. Required fields are marked *
Save my name, email, and website in this browser for the next time I comment.
Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Nature Ecology & Evolution ( 2024 ) Cite this article
796 Accesses
4 Altmetric
Metrics details
Closely related species often use the same genes to adapt to similar environments. However, we know little about why such genes possess increased adaptive potential and whether this is conserved across deeper evolutionary lineages. Adaptation to climate presents a natural laboratory to test these ideas, as even distantly related species must contend with similar stresses. Here, we re-analyse genomic data from thousands of individuals from 25 plant species as diverged as lodgepole pine and Arabidopsis (~300 Myr). We test for genetic repeatability based on within-species associations between allele frequencies in genes and variation in 21 climate variables. Our results demonstrate significant statistical evidence for genetic repeatability across deep time that is not expected under randomness, identifying a suite of 108 gene families (orthogroups) and gene functions that repeatedly drive local adaptation to climate. This set includes many orthogroups with well-known functions in abiotic stress response. Using gene co-expression networks to quantify pleiotropy, we find that orthogroups with stronger evidence for repeatability exhibit greater network centrality and broader expression across tissues (higher pleiotropy), contrary to the ‘cost of complexity’ theory. These gene families may be important in helping wild and crop species cope with future climate change, representing important candidates for future study.
Is evolution repeatable? This question, captured by Stephen Jay Gould’s ‘replaying the tape of life’ metaphor 1 , has been the subject of decades of empirical research (reviewed in ref. 2 ). The answer appears to be context-dependent, with variation in repeatability among taxa, populations within species and genes in the genome, which leads to new questions about which evolutionary forces govern the continuum between identical parallel change and divergent responses 3 , 4 . The relative importance of deterministic versus stochastic and contingent explanations has been the subject of great interest and has been explored extensively at short evolutionary timescales among populations (for example, in stickleback 5 ) or closely related species (for example, Brassicaceae 6 , Arabidopsis spp. 7 and Heliconius spp. 8 ). These systems have highlighted predominant roles of the availability of common adaptive variation through shared inheritance 9 , 10 , 11 , gene flow 8 or recurrent mutation 12 . Experimental evolution has further emphasized the inverse associations between genetic repeatability and the complexity of adaptive trait architecture 13 . We know much less, however, about the processes and contingencies shaping genetic repeatability at much longer evolutionary distances or whether such distances impose hard limits on observing any repeatability at all. The Mc1r gene provides a textbook example of repeatability across deep time, driving adaptive colour polymorphism in distantly related vertebrates, from fish to mammoths 14 . Some have speculated that the widespread re-use of this gene is due to minimal interactions with other genes, facilitating its modification while incurring minimal pleiotropic disruption 15 . We have little understanding of whether the kinds of factors affecting repeatability for candidate genes such as Mc1r also generalize to drive repeatability at genome-wide scales or to adaptive phenotypes beyond simple traits such as colour.
Climatic variation across the species ranges of plants is a ubiquitous selection pressure among distantly related species, presenting a natural laboratory to study repeated adaptation. Such variation exerts strong selection pressure for local adaptation and genotypic responses 16 , demonstrable through common garden experiments or provenance trials and reciprocal transplants 17 . Numerous candidate genes for drought and thermal tolerance have been identified in, for example, Arabidopsis thaliana 18 , 19 , Panicum hallii 20 and conifers 21 . Plants may use a wide array of strategies to adapt to climate, including phenological shifts, such as flowering later to avoid frost or modification of structures, such as smaller leaves to mitigate the effect of air temperature in hotter climates 22 . While such studies have advanced our understanding of climate adaptation within individual species and have provided examples of a few individual candidate genes with evidence of adaptiveness in several species, there has been little combined analysis of whole-genome patterns across many species and large phylogenetic distances. There is reason to believe that phylogenetic constraints might exist for repeatability among plants, if, for example, mutational target size differs with genome size 23 . Alternatively, diverging lifestyles may limit common selection pressures or functional solutions 24 . Repeatability in adaptive genetic responses to climate across independently evolving species is expected if biological adaptation is limited to conserved functions, whereas a lack of repeatability might indicate highly polygenic, alternative adaptive strategies to the problem of climate adaptation 25 . The survival of plant populations challenged by anthropogenic climate change may depend on adaptation, which itself may depend on the maintenance of adaptive genetic variation within the metapopulation of a species through local adaptation 26 , 27 . Thus, studying repeatability of climate adaptation will give insights into the flexibility of such genetic responses, which will both deepen our understanding of evolution and help predict how species may respond to changes in the future.
Here, we analyse sequencing data from thousands of individuals from 25 plant species (Fig. 1 ) across a distribution spanning diverse biomes across four continents (Fig. 1a ) to investigate repeatability in the genetic basis of adaptation to climate variation. By processing all raw data through a common pipeline, we compiled by far the largest population genomics dataset to date in terms of phylogenetic breadth and sequencing effort to examine the phenomenon of repeated local adaptation across the genome. Our aims are to: (1) identify gene families that exhibit repeated evolutionary associations with climate across several species; (2) test whether genes driving repeated adaptation tend to fall within particular functional groups, suggestive of repeatability at the level of gene function; and (3) test whether properties of these gene families, specifically pleiotropy and duplication history, may facilitate their repeated evolution and significance for climate adaptation across the plant kingdom.
a , Sampling locations from all 29 datasets (25 species) and global annual mean temperature. b , An example of six bioclim variables across a single dataset with ten sampling locations from ref. 65 . c , An example for the same dataset showing how climate change variables were calculated, taking monthly climate data from the 1960s and 2010s and calculating the effect size of the difference between decades. d , The species tree derived from the OrthoFinder2 analysis for the 17 reference genomes used in this study; species where a reference genome assembly was used to map related species but that species was not included in analyses, are marked with asterisks. The reference genome assembly at each tip came from the species not listed in brackets and species listed in brackets below had their data mapped onto that reference genome. The age of high-confidence nodes is shown with values pulled from TimeTree (Extended Data Fig. 1 ). Species analysed here are placed according to the reference genome used and the reference’s position in the tree. e , f , The per genome distribution of the number of paralogues per orthogroup ( e ) and number of representative genomes per orthogroup ( f ) are depicted as stacked bars, with bar fill scaling from light-to-dark according to most- to least-statistically tractable value. Bars are ordered vertically in line with the tree in d . Ma, million years ago.
Assembling datasets and defining orthology across species.
Twenty-nine datasets were downloaded from the NCBI sequencing read archive (SRA) or the EBI European nucleotide archive (ENA) repositories, covering 25 unique species (Supplementary Table 1 and Extended Data Fig. 1 ). These data were individual whole-genome sequencing (WGS), capture-based sequencing (CAPTURE) or pool-sequencing (POOL) data, each with a minimum of five sampling locations ( Methods gives full selection criteria). Raw reads were processed using a common single-nucleotide polymorphism (SNP)-calling pipeline (Supplementary Fig. 1 ) to minimize the influence of bioinformatic technical artefacts. The number of SNPs and individuals sampled per dataset ranged from 173,119 to 23,406,976 and 46 to ~1,300, respectively (Supplementary Table 1 ). To enable comparisons of several species (Fig. 1d ), we reconstructed orthology relationships using OrthoFinder2 (ref. 28 ) and classified genes into orthologous sets across species (orthogroups), which could include several copies of a gene in a given species due to gene duplication. In total, 44,861 orthogroups were classified, 14,328 were species-specific and 92.8% of all genes from all genomes were assigned to an orthogroup. Importantly, orthogroups were predominantly characterized by low rates of paralogy (Fig. 1e ), with rates of single-copy genes ranging from 39.8% to 62.4% (Supplementary Table 2 ) and high levels of presence across species (Fig. 1f ). From a statistical point of view, low paralogy is beneficial in aiding interpretation and comparison of genes among species and reducing confounding effects, whereas high presence across species increases statistical power to detect repeatability.
To identify genes with signatures of local adaptation within species, we performed genotype–environment association (GEA) scans by testing the association among population allele frequencies and population climatic variation (bioclim variables 1–19; Fig. 1b ) taken from worldclim (v.2.1, 2.5 arcmin resolution) 29 using non-parametric Kendall’s τ correlations ( Methods ). We also defined and quantified two variables that capture change in local climate (maximum temperature and precipitation) over the last 50 years. Briefly, these represent the effect size per sampling site of the difference between climate, quantified across a decade’s worth of data from either the 1960s or 2010s (example, in Fig. 1c ; Methods ). We then combined per SNP P values to calculate per gene P values using the weighted- Z analysis (WZA) 30 correcting for associations between SNP count and WZA variance. This WZA GEA method exhibits increased power and reduced error for identifying adaptive genes compared with other commonly used methods 30 .
To identify orthogroups with repeated signatures of association across several species, we applied PicMin 31 ( Methods ) across 8,470 orthogroups with at least 20 species represented for each of the 21 climate variables. These focal orthogroups were slightly more likely to include genes with stronger associations with climate relative to untested orthogroups (Supplementary Results 1 ). For orthogroups with several paralogues within a given species, we include the paralogue with the strongest evidence of association to test for repeatability after correction for multiple testing. This approach, therefore, tests for repeated adaptation driven by any member of a gene family.
Across all climate variables, we identified 141 repeatedly associated orthogroups (hereafter, RAOs) at a lenient false discovery rate (FDR) of 50%. These were made up of 108 unique orthogroups, with some showing significant hits for several climate variables, caused in part by covariance among climate variables. This set of RAOs represents an approximately threefold enrichment relative to expectations by chance, so although this list probably includes many individual RAOs that are false positives given the lenient threshold, it is highly unlikely to observe this many RAOs at FDR < 0.5 under a null model of no repeatability (1,000 random permutations, median expected number of RAOs was 36 across all variables, P < 0.001; Fig. 2a,b ). Full details of randomizations and justification of thresholds are in Supplementary Methods 1 . In line with our first aim, these results demonstrate robust evidence for genetic repeatability across ~300 Myr of plant evolution and yield a suite of candidate orthogroups for further analysis. These 108 candidate RAOs will include false positives but these should only add noise to analyses of gene properties. A summary of statistical tests for all 141 RAOs is provided in Supplementary Table 3 .
a , b , Enrichment of the total ( a ) and unique ( b ) number of RAOs observed at FDR < 0.5 (red lines) relative to random permutations (distribution and mean dashed line). c , The number of orthogroups exhibiting significant repeatability at different FDR thresholds per climate variable. Details about the interpretation of FDR values are provided in the main text.
RAOs were observed across most climate variables, although the number of RAOs per climate variable and strength of statistical support varied. Generally, temperature variables exhibited more RAOs with stronger evidence of repeatability, in particular maximum temperature in the warmest month (Fig. 2c ). This indicates that the adaptive molecular response to temperature variation across plants may be more repeatable at the level of individual genes, compared with precipitation. This might reflect adaptive constraints underlying temperature adaptation or the added complexity of how precipitation interacts with soil to modulate drought effects. Variability in the number of RAOs identified across climate variables for a given species was not linked to inferred GEA power or relative niche breadth among species (as a proxy for selection variability) (Supplementary Results 2 and Supplementary Figs. 2 – 4 ). For downstream analyses, we focus on two sets of orthogroups: those with FDR < 0.5 to explore general trends in the distribution of RAOs among species and climate variables and to discuss specific candidate genes; and those with a PicMin P < 0.005 (based on the top decile cutoff of strongest evidence of repeatability across all climate variables per orthogroup) within each climate variable to explore general properties of RAOs.
For RAOs (FDR < 0.5), we identified species driving the signature of repeatability on the basis of their per orthogroup P value ( Methods ). Visualizing the magnitude of these contributions by species and climate variables (Fig. 3a ) and among pairs of species (Fig. 3b ) demonstrates that no specific species, or cluster of species, contributes excessively to the repeatability signatures that we observe. In agreement, we failed to observe any evidence of phylogenetic signal of GEA results in our RAOs with respect to random orthogroups (Supplementary Results 3 and Extended Data Fig. 2 ). We reproduced Fig. 3a,b using RAOs detected at a more stringent FDR < 0.3 threshold and observed a comparable lack of phylogenetic signal (Extended Data Fig. 3 ). This shows that our identification of significant orthogroups is not driven by groups of closely related species but rather a signature observed broadly across the phylogenetic tree. This stands in contrast to expectations that repeatability declines as a function of relatedness among species 24 , 32 . Our result may be linked to the necessity of looking at conserved gene families, which may be less likely to have functionally diverged among species. In line with this, when repeating PicMin analyses only within Brassicaceae, we observed that statistical signals of repeatability tended to be stronger in our tested set of orthogroups (those with at least 20 species overall) relative to untested orthogroups with lower conservation (orthogroups with fewer than 20 species overall but including all Brassicaceae; Supplementary Results 4 and Supplementary Fig. 5 ).
a , b , Heatmaps show the contribution of individual species to orthogroup repeatability at FDR < 0.5 for different climate variables ( a ) and among pairs of species ( b ), with species ordered phylogenetically. In each case, the fill of each cell represents the proportion of orthogroups where a given species contributes towards the signature of repeatability based on its minimum GEA P value. c , Heatmap showing the number of times a species contributes towards repeatability for a given RAO (FDR < 0.5). In c , cell fill denotes the number of climate variables where a species contributes a low P value to a given RAO, with grey being 0. The vertical green line separates gymnosperms and angiosperms. Row-wise summations are shown as bars, where the blue bar shows the number of species associated with a given orthogroup (species n ) and the red bar shows the total number of species and climate variables (total n ). RAOs associated with a single climate variable only have a blue bar and those with stacked red bars are RAOs associated with multiple climate variables. The strongest statistical support for repeatability is also shown per orthogroup as −log 10 -transformed FDR. Only 73 RAOs with at least five contributing species are shown.
As an additional exploration into sources of variation within our repeatability analyses, we conducted leave-one-out cross-validations by removing each species and re-running PicMin (Supplementary Results 5 ). These analyses highlighted variation among species in terms of either increasing or decreasing the number of RAOs at FDR < 0.5 and <0.3. At FDR < 0.5, removing either Eucalyptus albens or E. sideroxylon reduced the number of RAOs the most, whereas removing Amaranthus tuberculatus caused the largest increase (Extended Data Fig. 4 ). Comparing these against features of the original SNP datasets indicated that reductions in RAOs were associated with removing datasets in which the sampling covered a larger amount of the global range of the species, geographically and in terms of climate breadth (Extended Data Fig. 4 ). Other factors, such as environmental similarity 5 , may also explain heterogeneity in repeatability among pairs of species (Fig. 3b ). The variability of dataset quality could confound signals among close relatives if, for example, they were sampled in different environments or to different extents of their ranges in our dataset.
Sixteen orthogroups were associated with repeatability across several climate variables. The orthogroup including the A. thaliana genes AT5G60100 ( PRR3 ) and AT5G02810 ( PRR7 ), with instrumental functional roles in circadian rhythm 33 and flowering time 34 , was repeatedly adaptive across 10 of 21 variables (most strongly with mean temperature in the dry quarter, FDR = 0.102; Extended Data Fig. 5 ) and across several variables within the same species (difference between red and blue bars in Fig. 3c ). The role of this gene family in circadian rhythm may contribute to its repeated association with several climate variables if these also vary with latitude. Another notable RAO includes a family of four genes encoding ubiquitin-related proteins: AT1G31340 ( RUB1 ), AT1G11980 ( RUB3 ), AT2G35635 ( UBQ7 ) , AT1G11970 , which was associated with six climate variables and most strongly with mean temperature (FDR = 0.051; Fig. 3c and Supplementary Table 3 ). The related-to-ubiquitin (RUB)-conjugation pathway has been implicated in the auxin response, embryo development and growth 35 . A summary of genes in RAOs that are associated with phenotypes known to be adaptive under climatic stress 36 is shown in Table 1 (full details of RAO gene contents in Supplementary Tables 4 – 6 ).
The functions of these genes reflect the expected responses to climatic variation, such as phenological avoidance of drought or frost through changes to flowering time or seed dormancy 22 or modification to root hair number and structure in response to temperature changes 37 . Changes to growth or stomatal function through hormone signalling 38 meanwhile may facilitate tolerance to drought by reducing water loss and genes associated with salt stress may be involved in surviving salt accumulation in soils due to aridity.
We found only a single RAO associated with our two climate change variables at FDR < 0.5; harbouring the A. thaliana genes AT5G53480 ( ATKPNB1 ), AT3G08943 and AT3G08947 . ATKPNB1 is sensitive to abscisic acid and is involved in drought tolerance through stomatal closure 39 . The few RAOs here may be due to the relatively short amount of time that our climate change variables are calculated over (~50 years) and the limited time to respond to selection subsequently, particularly in longer-lived species. Three species contributed to repeatability in this RAO (Supplementary Table 3 ), two of which are short-lived. However, given associations between our climate change variables and general bioclim variables (Extended Data Fig. 6 ), we cannot rule out contributions by longer-lived species.
To examine repeatability beyond the gene level we characterized the functions of RAOs, to assess cases where several genes from within the same molecular pathway are used for adaptation across several species. Examples of ‘functional repeatability’ have been documented in adaptation to whole-genome duplication in Arabidopsis 40 and highland adaptation in maize 41 . To explore this phenomenon in accordance with our second aim, we used the STRING database 42 (v.11.5) to provide a network-based representation of protein–protein interactions (direct and indirect associations compiled from genomic context, experimental evidence such as co-expression and text-mining of literature) and tested whether RAOs formed networks with more interactions than expected ( Methods ). We grouped RAOs all together and into temperature- and precipitation-related groups and tested each group to see whether RAOs as a group contained genes that were more likely to interact with one another than random orthogroups.
Each group of RAOs tended to include more protein–protein interactions than random orthogroups and this was statistically significant for RAOs identified through precipitation variables, where ~2.4× more connections were observed than expected by chance (10,000 random permutations, average observed interactions was 19.06, expected number of interactions was 8.04, P = 0.015) (Fig. 4a ). Gene ontology (GO) enrichment over the same groups of RAOs reflected this, as precipitation-related RAOs were highly enriched (hypergeometric test, one-tailed FDR < 0.1) for biological processes with clear adaptive roles (Fig. 4b,c and Supplementary Table 7 ). Orthogroups associated with ‘root hair tip’ were enriched across all RAOs (hypergeometric test, observed = 0.028, expected = 0.001, FDR = 0.046) and temperature RAOs were enriched for several processes, notably ‘regulation of photoperiodism, flowering’ (hypergeometric test, observed = 0.056, expected = 0.004, FDR = 0.016) and ‘brassinosteroid biosynthetic process’ (hypergeometric test, observed = 0.028, expected = 0.002, FDR = 0.099), which are known to be involved in thermotolerance 22 , 43 . These results suggest that independently identified RAOs contain genes involved in similar functional processes, which implies that repeated adaptation is also occurring beyond the level of the gene.
a , Mean number of STRING interactions within single-gene per orthogroup networks based on grouped RAOs (triangles) compared to results from permutation tests. Circles and lines represent the mean and 5/95% quantiles STRING interactions among 10,000 random gene sets sampled for single genes in the same way. The colour of circles shows the number of orthogroups in each group of RAOs. b , The network derived from all genes in precipitation-related RAOs, with orange highlighting which genes are members of four enriched GO terms within the same network. c , Enriched GO terms (FDR < 0.1) showing the number of contributing orthogroups from each species, with −log 10 -transformed FDR-adjusted significance shown as adjacent bars. GO terms are ordered on the y axis based on semantic similarity clustering, to group GO terms associated with similar genes together. Cell fill shows the number of orthogroups per species and GO term, with grey being 0.
The above patterns of enrichment may occur because a given adaptive response requires the coordinated modification of several functionally related genes in all species involved. Alternatively, this signal could also be driven by different subsets of functionally related genes contributing to adaptation in different subsets of species (for example, if genes A, B, C and D are functionally related, species 1, 2 and 3 adapt via genes A and D while species 4, 5 and 6 adapt via genes B and C). In either case, our results show that certain pathways or functional groups of genes are particularly important for adaptation to these climatic stressors, particularly with regards to adaptation to precipitation variation.
We next ask whether the gene families identified that contribute to repeated adaptation share particular characteristics with respect to their degree of pleiotropy. Pleiotropy is a fundamental attribute of a gene describing the number of traits it affects. On the basis of Fisher’s model of universal pleiotropy 44 , the ‘cost of complexity’ hypothesis 45 posits a reduced adaptive potential for genes with greater pleiotropy, as constraint increases with organismal ‘complexity’. In keeping with this, greater fitness consequences are predicted by the degree of pleiotropy in yeast 46 . However, empirical evidence from mice, nematodes and yeast suggests that this cost may be counteracted by a greater mutational effect size per trait observed for genes with greater pleiotropy 47 , 48 . In line with this, others 49 found genes repeatedly involved in stickleback adaptation exhibited elevated levels of pleiotropy.
To test the importance of pleiotropy by several definitions in our dataset, we used public databases of gene expression for A. thaliana and Medicago truncatula genes extracted from Expression Atlas 50 and ATTED-II 51 . We explored pleiotropy by two definitions: tissue specificity 52 and condition-independent co-expression with other genes 53 . Tissue specificity of gene expression is inversely associated with pleiotropy, has previously been linked to increased rates of evolution 54 and was estimated here according to the τ metric 55 (Fig. 5a ). The τ metric describes tissue specificity, the inverse of pleiotropy, so to avoid confusion we will describe changes to the breadth of expression, which we define as lower τ . Contrary to the cost of complexity prediction, we found that RAOs with the strongest evidence of repeatability were strongly associated with increased expression breadth (Stouffer’s Z , P = 5.44 × 10 −4 ). Expression breadth also tended to decrease in subsets of orthogroups with increasingly weaker evidence of repeatability, such that orthogroups with the weakest evidence of repeatability were enriched for genes with high specificity (Stouffer’s Z , P = 4.74 × 10 −6 ; Fig. 5b ).
a , An example of breadth of tissue expression. Each cell in the heatmap shows tissue-specific gene expression scaled by the maximum across all tissues for five genes, with their τ statistic as a bar. b , Stouffer’s Z calculated for orthogroups grouped into deciles on the basis of the strongest evidence of repeatability across all climate variables (green, P < 0.05 more pleiotropic than expected; red, P < 0.05 less pleiotropic than expected). c , An example network of 20 genes with four centrality measures calculated across nodes. Node size and colour (small-blue, low; large-red, high) denotes centrality by each metric. d , Stouffer’s Z for node degree and tissue specificity metrics within climate variables, focussing on orthogroups with at least one PicMin P < 0.005. e , Non-parametric correlations of per orthogroup pleiotropy and duplication metrics, along with a randomized control variable. At, A. thaliana ; Mt, M. truncatula .
Alternatively, pleiotropic constraint can be considered as various node centrality statistics within a gene co-expression network (Fig. 5c ), where the ‘distance’ between two gene nodes is lower when co-expression of those genes increases (for example, across experimental treatments or tissue types). Co-expression centrality measures have been inversely linked to rate of evolution in several eukaryotic protein–protein networks and changes to genes with high centrality have a higher chance of being lethal 56 , 57 , indicative of evolutionary constraint. However, in contrast to these negative associations with evolvability, we observed clear positive associations between evidence of repeatability and co-expression centrality for node closeness, degree and strength across both co-expression networks (Fig. 5b ). Similar to results for expression breadth, centrality was significantly greater (Stouffer’s Z , P < 0.05) in orthogroups with the strongest evidence of repeatability and significantly lower (Stouffer’s Z , P < 0.05) in orthogroups with the weakest evidence. These clear trends across both specificity of tissue expression and co-expression networks highlight a robust association between increased pleiotropy and evidence of adaptive repeatability across all tested orthogroups. Associated measures of pleiotropy based on A. thaliana estimates for all orthogroups are included in Supplementary Table 8 .
The association we find between repeatability of local adaptation and increased pleiotropy stands in apparent contrast to previous findings of increased contributions to adaptation by genes with reduced pleiotropy 54 , 56 , 57 . This difference may arise because of the scale at which adaptation is occurring; here, we have focused on local adaptation, which involves a tension between migration and spatially divergent selection that tends to favour the contribution of alleles of large effect, as they can overcome migration swamping 25 , 58 . As the phenotypic effect size of mutations tends to increase with pleiotropy 47 , 48 , increased pleiotropy may therefore be favoured by natural selection during local adaptation. By contrast, when a species adapts to a temporal change in environment across its whole range, there is no tension between migration and selection and no further advantage for alleles of larger effect 25 . This may explain the reduced pleiotropy previously observed in rapidly evolving genes 54 , 56 , 57 . Our results suggest a robust impact of pleiotropy on local adaptation across several plant species, consistent with similar observations in stickleback 49 , ragweed 59 and A. thaliana 60 . It is unknown whether the association between pleiotropy and repeatability is monotonic or if intermediate pleiotropy promotes repeatability and extreme pleiotropy remains constraining, as suggested by refs. 49 , 60 . In the A. thaliana tissue specificity and node degree data, RAO sets were diminished for the least pleiotropic genes but also enriched for the most pleiotropic (Extended Data Fig. 7 ). It is important to make clear that we cannot rule out an association between repeatability and pleiotropy because of an increased likelihood of detecting large-effect alleles with GEA methods. A comparable analysis centred on global adaptation through selective sweeps could distinguish these, as selective sweep methods are similarly biased towards large-effect alleles but there is no assumption of a biological role for increased pleiotropy facilitating global adaptation at the species level 44 , 45 .
We were also interested in whether pleiotropy enrichment was limited to specific climate variables. We therefore repeated enrichment analyses for orthogroups exhibiting the strongest evidence of repeatability within climate variables (PicMin P < 0.005) for tissue expression specificity and A. thaliana node degree (Fig. 5d ). Most of these sets of orthogroups exhibited elevated pleiotropy. Notably, orthogroups with the strongest evidence of repeatability associated with our climate change variables were highly enriched for pleiotropic genes by both measures (precipitation change—expression breadth P = 0.0201, A. thaliana node degree P = 0.0003; maximum temperature change—expression breadth P = 0.0129, A. thaliana node degree P = 0.0112). Given that these variables only capture environmental change over ~50 years, the genes with the strongest evidence of repeatability associated with these variables may be highly pleiotropic as a result of genes with greater effect sizes facilitating rapid adaptation to shifting fitness optima 61 .
Our orthogroup-level estimates of pleiotropy were also negatively associated with the number of gene duplication events within an orthogroup (Fig. 5e ). This was not an artefact of how per gene pleiotropy estimates were condensed to orthogroup-level estimates (Supplementary Results 6 ). Variability in gene duplication rate may also influence the likelihood of repeatability among gene families 62 . Indeed, we tended to see stronger evidence of repeatability in orthogroups with fewer duplication events, even when accounting for the effect of increased multiple-corrections (Extended Data Fig. 8 ). This negative association between duplications and repeatability is the opposite trend to that observed between conifer species by ref. 62 . This difference may stem from interactions between gene duplication and functional divergence over long periods of time. For example, between closely related species gene duplication may promote repeatability by alleviating functional constraints through neofunctionalization. However, across deep-time, gene families with higher rates of duplication may be expected to functionally diverge, inhibiting repeatability through a degradation in the common mapping of genotype to phenotype. Owing to inverse associations between duplication and pleiotropy estimates, however, it is difficult to disentangle these two properties of gene families (Supplementary Results 6 ).
Across plant species separated by >300 Myr of evolution, there is repeatability in the genetic basis of local adaptation to climate variation. The gene families identified in our analyses will be of substantial adaptive relevance as plants find themselves under increased climatic selection due to anthropogenic global change. We show that the adaptive responses to such changing selection involves conserved, core functional responses. We have also demonstrated that, contrary to expectations that pleiotropic constraint impedes adaptation, the genes identified as repeatedly adaptive bear signatures of increased pleiotropy. These results support the model of pleiotropy driving increasing phenotypic effect sizes and such large effects being a principal driver of local adaptation in the genome. Whether these results extend beyond the dominant biomes of boreal and temperate forest found in the Global North to species in the tropics and Global South is a question of great importance.
We selected 29 datasets covering 25 species from 21 studies 6 , 62 , 63 , 64 , 65 , 66 , 67 , 68 , 69 , 70 , 71 , 72 , 73 , 74 , 75 , 76 , 77 , 78 , 79 , 80 , 81 for which sequence data were available with WGS, POOL or CAPTURE (Supplementary Table 1 ). Such datasets provide broad genomic representation and maximize the resolution of the number of genes for each species. We limited our study to datasets with at least five populations distributed along latitudinal, longitudinal and/or altitudinal gradients. Here, the minimum number of populations was set to allow for a minimum baseline of statistical power for Kendall’s τ correlations between environment and allele frequency (see below). The list of included datasets is not exhaustive but rather was selected to achieve a desired number of approximately 20–30 species for which power to detect repeated adaptation has been demonstrated for the methods used here 82 . Datasets also had to provide locality information, either for individual samples or groups of individuals. Species were selected to include both gymnosperms and angiosperms, setting the most recent common ancestor of all species sampled at approximately 300 million years ago (Extended Data Fig. 1 ) 83 .
We used two SNP-calling pipelines depending on whether sequencing data came from individuals (WGS and CAPTURE) or from pools of individuals (POOL). These pipelines were necessarily different but were based on similar approaches to reduce bioinformatic discrepancies between the data types. For all data types, raw fastq data were retrieved from either SRA or ENA. Accession codes for all data are provided in Supplementary Table 1 . The reference genomes used for each species, or closely related species if genomes were not available, are also in Supplementary Table 1 . SNP pipelines were designed to balance computational time and low false-positive rates 84 .
The pipeline for individual-based data was as follows, note that selfing and outcrossing species were processed using the same pipeline. Raw fastq files were cleaned and trimmed for adaptor sequences using fastp (v.0.20.1) 85 before being aligned against the reference genome using bwa-mem (v.0.7.17-r1188) 86 . BAM files were generated, sorted and indexed using samtools (v.1.16.1) 87 , skipping alignments with MAPQ < 10 (-q 10). We then collected quality metrics with Picard Tools (v.2.26.3) 88 based on alignment summary (CollectAlignmentSummaryMetrics), insert size metrics (CollectInsertSizeMetrics) and coverage (CollectWgsMetricsWithNonZeroCoverage). We then marked and removed duplicates using Picard’s MarkDuplicates and used AddOrReplaceReadGroups to amend read groups. In some cases, datasets split sequencing data from individual samples across several technical replicates, so we then merged BAM files within samples with Picard’s MergeSamFiles. We ran a realignment of the cleaned, merged BAM files by running the RealignerTargetCreator and IndelRealigner from the Genome Analysis Tool Kit (GATK v.3.8) 89 and repeated the aforementioned quality metrics on final BAM files. To identify SNPs, we generated genotype likelihoods using BCFtools 87 mpileup by specifying a minimum mapping quality of five for an alignment to be used and retained further annotation information such as allelic depth. From there, individual pileups were converted into SNP VCFs by using the BCFtools call programme. We set sample ploidy (-S) information to match the known ploidy of a species and called for genotype quality to be reported while excluding any group samples (-G -) information. Finally, we filtered raw VCF files with VCFtools 90 by removing sites with quality value below 30 (--minQ 30), Genotype quality below 20 (--minGQ 20), minimum read depth of 5 (--minDP 5), before finally retaining only biallelic (--max-alleles 2) genotypes present in >70% of individuals (--max-missing 0.7). For downstream analyses, we performed more filtering on the basis of minor allele frequency (maf) and minor allele count (mac), retaining only sites with maf > 0.05 and mac > 5, whichever was most stringent.
Pooled data were processed using a similarly structured workflow. Raw fastq files were cleaned and trimmed with fastp and aligned to references with bwa-mem using the additional flag to mark shorter split hits as secondary (-M). BAM files were generated, sorted and indexed using samtools, skipping alignments with MAPQ < 20 (-q 20) and bedfiles were generated from indexed BAM files using BEDtools (v.2.27.1) 91 . Duplicates were then marked and removed with MarkDuplicates before indel realignment was performed with Picard’s RealignerTargetCreator and IndelRealigner. SNP-calling was then performed using mpileup followed by the VarScan (v.2.4.2) mpileup2cns programme. Variants were called on the basis of minimum read depth of 8 (--in-coverage 8), a P value threshold of 0.05 (--p-value 0.05), a minimum frequency for calling homozygotes of 80% (--min-freq-for-hom 0.8), ignoring variants with >90% support on one strand (--strand-filter 1), minimum base quality at position of 20 (--min-avg-qual 20) and a minimum variant allele frequency of 0 (--min-var-freq = 0). To extract allele frequencies, pooled SNPs were converted to SNP tables using GATK (v.4.1.0.0) 92 VariantsToTable, splitting multi-allelic sites across multiple rows (--split-multi-allelic) and extracting the AF field. Final SNPTables for POOL data were filtered for indels, retaining only biallelic sites and a minor allele frequency cutoff of 0.05.
We grouped genes from each reference genome into orthogroups (sets of genes across species that are descended from a single gene in the last common ancestor 28 ) to facilitate comparisons among species. To construct orthogroups, we built proteomes for each reference genomes by extracting the protein sequence of the longest isoform for each gene using the AGAT (v.1.0.0) scripts agat_sp_keep_longest_isoform.pl and agat_sp_extract_sequences.pl. Proteomes were parsed as input to OrthoFinder2 (v.2.5.2) 28 using default settings. The species tree (Fig. 1d ) based on the reference genomes was inferred using all orthogroups where all genomes were represented ( n = 5,003; Fig. 1f ) 93 , 94 .
Resulting orthogroups were then filtered for those that would yield the most statistical power, that is low rates of paralogy and high rates of species representation. Low rates of paralogy were necessary owing to multiple-comparison corrections performed among paralogues within species (see below), thus we removed a genome from an orthogroup if it had more than ten paralogues and retained only orthogroups with at least 20 representative species (this was dropped down to 19 species when analysing isothermality, as no GEA data were calculable for Eucalyptus magnificata for this variable). Of the 44,861 classified orthogroups, 8,470 were retained as high-quality for downstream analyses. The focal 8,470 high-quality orthogroups (20–25 species represented, maximum paralogue n per species of 10; Supplementary Fig. 6 ) covered 37.1%–61.6% of genes within each genome (Supplementary Table 2 ). The species tree from OrthoFinder was in good agreement with the species tree from TimeTree 83 (Supplementary Methods 2 and Supplementary Fig. 7 ). Only 355 orthogroups exhibited one-to-one orthology over 20 or more species, representing a substantial limitation of scope if only these were to be analysed.
For each unique latitude–longitude co-ordinate pair, we extracted bioclimatic variable (bioclim) data from the worldclim database (v.2.1) 29 at 2.5 arcmin resolution for all 19 bioclim variables. To describe recent climate change, we added two variables that quantified change in maximum temperature (tmax) and precipitation between all worldclim data collected monthly for the 1960s compared with the 2010s. As a measure of change, we estimated a non-parametric effect size between the distribution for each decade at each sampling location using the rstatix::wilcox_effsize() function 95 . Thus, a larger effect size is indicative of a greater change in recent climate (example shown in Fig. 1c ). We then calculated associations between individual allele frequencies and each of the 21 bioclim variables individually, without accounting for expected covariance among them (Supplementary Methods 3 ).
We first converted individual-based data into allele frequencies by combining individuals labelled as sampled from the same location and calculating per sampling site allele frequencies. These should not be considered as populations in the traditional population genetics sense, in that we do not assign populations on the basis of shared ancestry among individuals. In doing so, we trade off the risk of allele frequency sampling error when only a handful of individuals, or even a single individual, are sampled at a given location, against increased sampling power in downstream GEAs.
For each dataset, we combined per sampling location allele frequencies with our 21 climate variables and performed GEAs on the basis of non-parametric Kendall’s τ correlations (that is, no correction for population structure; justification of methodology is in Supplementary Methods 4 ). Each Kendalls’ τ correlation yields a P value for the correlation between allele frequency and environment. Because we are interested in repeatability at the gene level, we retained all SNPs within annotated gene boundaries (plus 500 base pair flanking regions) and converted P values to empirical P values (e P values: the rank of the P value divided by the number of P values) within each environmental variable. This yielded 21 uniform distributions of e P values across all SNPs for our 21 climate variables. To get gene-level GEA results, we combined evidence of all SNPs within each gene using the WZA 30 . This approach exploits the signal of elevated linkage among nearby SNPs within genes that is representative of local adaptation and has been shown to have similar or improved power and error relative to other GEA approaches including BayPass, RDA and LFMM 30 . The approach combines e P values from all SNPs within genes and flanks while weighting each SNP relative to its expected heterozygosity, resulting in a theoretically normal distribution of per gene Z -scores ( Z WZA ). We performed additional corrections to the Z WZA to account for heteroscedasiticity linked to variable SNP count per gene (Supplementary Methods 5 ).
If we had several datasets for an individual species ( A. thaliana , n = 3; A. lyrata , n = 2; A. halleri , n = 2), we combined Z WZA P values within genes using Fisher’s approach, assuming each dataset represents an independent test of association between gene and environment, given datasets from the same species do not overlap spatially.
Because genes are not directly comparable among genomes, we combined per gene Z WZA e P values (e P value WZA ) with orthogroup assignments. Owing to low numbers of simple one-to-one orthogroups, we included species in orthogroups with a maximum of ten paralogues. For each species, we then retained the strongest signal of association (minimum e P value WZA ) within an orthogroup, correcting for the number of paralogues with a Dunn–Šidák correction, akin to Tippett’s method 96 . This approach of taking the strongest GEA among paralogues within the same orthogroup is based on the assumption that modification of any paralogue, rather than all paralogues, may be sufficient for adaptation if they are functionally similar 97 . This was preferable in contrast to simply combining evidence of association among all paralogues.
For each climate variable, we tested orthogroups with e P values WZA from at least 20 species for statistical clustering of low e P values WZA using PicMin 82 . PicMin was run in separate configurations depending on whether an orthogroup contained any of 20–25 species, with results combined at the end. Our null model for the among-species Tippett correction was that adaptation associated with an orthogroup is observed in one or no species. This test is performed by removing the minimum e P value WZA and testing each of the remaining n − 1 e P values WZA against the expected probability density functions of n − 1 random uniform draws. This generates a parametric P value for each observed e P value WZA based on how much lower the observed is relative to the beta expectation, which are combined according to Tippett’s method, that is taking the minimum P value corrected for the number of tests with a Dunn–Šidák correction and the expected correlation structure. This combined P value provides an estimate of the strength of evidence for the orthogroup contributing to repeated adaptation, with the rank of the minimum P value providing an estimation of the number of species driving the repeatability signature. A feature of PicMin under this null model is that as the number of species being tested increases, the multiple-comparison corrections applied increases to the point at which the distribution of PicMin P values across randomly distributed uniform P values is upwardly biased towards larger P values (Supplementary Fig. 8 ). This leads to an expected marginal loss of power. To address this, we simulated randomly distributed uniform P values from 20–25 species for 1,000,000 orthogroups and ran these through PicMin to generate an empirical distribution of PicMin P values under our random null model. We used this empirical null distribution to correct our observed PicMin P values using the qvalue::empPvals() function. For each climate variable analysed, the resulting corrected per orthogroup PicMin P values were FDR corrected and orthogroups with q < 0.5 were considered as showing evidence of repeatability (RAOs). Justification of this threshold is in Supplementary Methods 1 .
We tested for repeatability at a functional level by exploring functional interactions among RAOs. To do this, we took the A. thaliana genes from RAOs associated with each climate variable (FDR < 0.5) and asked whether these sets of genes were enriched for protein–protein interactions among orthogroups. We explored interactions within three groups: (1) all unique RAOs across all 21 climate variables; (2) unique RAOs related to temperature variables (bioclim 1–11 and tmax_clim_change); and (3) unique RAOs related to precipitation variables (bioclim 12–19 and prec_clim_change). To quantify interactions among genes, we used the STRING database 42 . We proposed that gene sets that were enriched for particular functions would include different orthogroups that included different, but functionally similar, genes that would interact with one another.
To explicitly measure interactions among orthogroups, as opposed to within, we first repeatedly sampled a single random gene from RAOs and quantified STRING interactions among these random single-gene per orthogroup gene sets. This was done to not count interactions among genes from the same orthogroup, which would be more likely to be functionally similar and potentially interact than random genes. Although this approach is conservative, single-copy genes are the most common in our orthogroups (Fig. 1e ). In the A. thaliana data used here, 58.5% of tested orthogroups were single copy, with two paralogues the next most common orthogroup size, so in most cases the loss of power incurred through randomly sampling a single gene should be minimal. With a single random gene per orthogroup, we then constructed the STRING network across these genes and counted the number of interactions across the network with medium or greater support (>0.4). We took the mean number of interactions across 1,000 networks (each time drawing a random single gene per orthogroup) for each of the three groups of RAOs as the ‘observed’ number of interactions (Fig. 4a ). We finally compared these observed values against 10,000 random draws, in which for each iteration we drew a random set of non-RAOs equivalent in size to each of the three groups. From this, we identified which sets of RAOs were associated with significantly more interactions among orthogroups than expected for random gene sets of equivalent size. An excess of interactions suggests that different orthogroups identified as repeatedly adaptive across different sets of species probably contain genes that are performing functionally similar roles but are not orthologous across species. Consequently, these networks may be particularly helpful at identifying more general biological processes that are associated with adaptation to climate if such processes are enriched within highly interactive networks.
We therefore asked what functions were enriched within all RAOs, temperature RAOs and precipitation RAOs. To do this, we collapsed GO assignments for all A. thaliana genes within a given orthogroup and removed duplicated GO terms. For example, if an orthogroup included two paralogues and paralogue 1 has GO terms GO1, GO2, GO3 whereas paralogue 2 has GO terms GO2, GO3 and GO4, the GO terms we retained for that orthogroup would be GO1, GO2, GO3 and GO4 (each once). This is important as our analysis is non-specific with regard to paralogues within orthogroups, so we cannot know which GO terms among paralogues may be relevant. We did this for all 8,470 tested orthogroups to produce a custom GO background and then assessed enrichment of GO terms within each set of RAOs. Enrichment of GO terms was determined on the basis of the hypergeometric expectation and P values were FDR corrected.
We estimated pleiotropy using two separate approaches: specificity of tissue expression and connectivity within co-expression networks. To calculate tissue specificity, we downloaded tissue expression data for A. thaliana from Expression Atlas 50 (accession no. E-MTAB-7978; ref. 98 ). This dataset comprises tissue expression (transcripts per million, TPM) across developmental stages, tissue types and subtissue type. Because we are interested specifically in specificity across different tissue types, we took the mean TPM across all developmental stages and subtissue types within the tissue type field. This resulted in mean TPM within each of the 23 tissue types (example in Fig. 5a ). The tissue specificity metric τ was calculated following ref. 55 as:
where, for a given gene, x i corresponds to the mean TPM for a given tissue type normalized by the maximum mean TPM across n tissue types. This yielded an estimate of τ for 30,074 A. thaliana genes. Of these, 25,831 could be matched to orthogroups and of these 15,472 were in the orthogroups tested for repeatability. To condense these to single orthogroup estimates, we converted τ estimates to rank-based e P values (least specific gives lowest e P value), took the lowest e P value within each orthogroup and corrected for the number of paralogues with a Dunn–Šidák correction. Finally, we transformed per orthogroup e P values to Z -scores with a mean of 0 and s.d. of 1 across all orthogroups. Alternative approaches were explored and are discussed in Supplementary Methods 6 .
As a complement to tissue specificity, we explored connectivity of genes in co-expression networks. We built two co-expression networks using co-expression data from ATTED-II 51 for A. thaliana and M. truncatula . Co-expression gene tables were downloaded for each species: A. thaliana , Ath-u.c3-0; M. truncatula , Mtr-u.c3-0. We discarded all edges with −5 < Z < 2.33 following the recommendations for significant negative or positive co-expression 51 . The A. thaliana network included 18,570 genes or 13,424 after retaining only genes in orthogroups tested for repeatability. Similarly, the M. truncatula network included 17,786 genes and 12,558 genes for the same groups. Networks were produced using the igraph package in R. Node betweenness and closeness were calculated using the estimate_betweenness() and closeness() functions, respectively. Node degree and strength were calculated as the number and absolute sum of edges respectively. The same approach was repeated for the co-expression network derived from M. truncatula . Orthogroup-level Z -scores were calculated for co-expression metrics as for expression specificity.
To assess how tissue expression specificity and co-expression centrality are associated with RAOs, we grouped orthogroups into deciles based on the strongest evidence for repeatability (minimum PicMin P value) observed for each orthogroup across all 21 climate variables. This was to rank orthogroups from ‘most repeatable’ to ‘least repeatable’ taking account of all tests performed. We then combined Z -scores within each decile based on Stouffer’s approach under a null hypothesis that if there is no association between repeatability and pleiotropy estimates, each decile should draw Z -scores randomly from the total distribution and yield Stouffer’s combined Z -scores of approximately 0. Stouffer’s Z -scores were finally converted to two-sided P values using the qnorm() function. We also performed the same analysis within each climate variable (Fig. 5d ).
Information on duplications within orthogroups was obtained from the OrthoFinder2 outputs. We retained duplication events that had support >0.7 and occurred at nodes which included species that contributed GEA P values. This total number of duplications is therefore associated to the number of paralogues within orthogroups but counting the duplication events as opposed to the number of paralogues avoids counting the same duplication event several times. As well as counting the number of duplication events within each orthogroup gene tree, we also counted specifically the number of duplications that occur within species, that is all gene-tree tips downstream of the duplication node that include only the genome of a single species. We refer to these as species-specific duplications and were interested in these because of the potential for sub- and neo-functionalization to occur within species. We also quantified the number of single-copy genes per orthogroup. To examine associations of each of these duplication statistics with repeatability, we again split orthogroups into deciles on the basis of their strongest evidence of adaptive repeatability across the 21 climate variables and assessed how each per decile mean duplication metric varied from strongest to weakest evidence of repeatability.
To rule out the possibility that orthogroup structure and the number of paralogues drove associations between duplications and repeatability we used a randomization procedure (Supplementary Methods 7 ).
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
The SRA codes for all raw sequencing data are in Supplementary Table 1 . All VCFs and sample location data are available via Dryad at https://doi.org/10.5061/dryad.15dv41p57 (ref. 99 ).
All scripts for SNP-calling and analyses are available and documented at https://github.com/JimWhiting91/RepAdapt and available via Zenodo at https://doi.org/10.5281/zenodo.12680122 (ref. 100 ).
Gould, S. J. Wonderful Life: The Burgess Shale and the Nature of History (WW Norton & Company, 1990).
Blount, Z. D., Lenski, R. E. & Losos, J. B. Contingency and determinism in evolution: replaying life’s tape. Science 362 , eaam5979 (2018).
Article PubMed Google Scholar
Bolnick, D. I., Barrett, R. D. H., Oke, K. B., Rennison, D. J. & Stuart, Y. E. (Non)parallel evolution. Annu. Rev. Ecol. Evol. Syst. 49 , 303–330 (2018).
Article Google Scholar
Arendt, J. & Reznick, D. Convergence and parallelism reconsidered: what have we learned about the genetics of adaptation? Trends Ecol. Evol. 23 , 26–32 (2008).
Magalhaes, I. S. et al. Intercontinental genomic parallelism in multiple three-spined stickleback adaptive radiations. Nat. Ecol. Evol. 5 , 251–261 (2021).
Rellstab, C. et al. Genomic signatures of convergent adaptation to Alpine environments in three Brassicaceae species. Mol. Ecol. 29 , 4350–4365 (2020).
Article CAS PubMed PubMed Central Google Scholar
Bohutínská, M. et al. Genomic basis of parallel adaptation varies with divergence in Arabidopsis and its relatives. Proc. Natl Acad. Sci. USA 118 , e2022713118 (2021).
Article PubMed PubMed Central Google Scholar
Montejo-Kovacevich, G. et al. Repeated genetic adaptation to altitude in two tropical butterflies. Nat. Commun. 13 , 4676 (2022).
Waters, J. M. & McCulloch, G. A. Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution. Mol. Ecol. 30 , 4162–4172 (2021).
Konečná, V. et al. Parallel adaptation in autopolyploid Arabidopsis arenosa is dominated by repeated recruitment of shared alleles. Nat. Commun. 12 , 4979 (2021).
Louis, M. et al. Selection on ancestral genetic variation fuels repeated ecotype formation in bottlenose dolphins. Sci. Adv. 7 , eabg1245 (2021).
Zou, D. et al. Vulture genomes reveal molecular adaptations underlying obligate scavenging and low levels of genetic diversity. Mol. Biol. Evol. 38 , 3649–3663 (2021).
Long, A., Liti, G., Luptak, A. & Tenaillon, O. Elucidating the molecular architecture of adaptation via evolve and resequence experiments. Nat. Rev. Genet. 16 , 567–582 (2015).
Manceau, M., Domingues, V. S., Linnen, C. R., Rosenblum, E. B. & Hoekstra, H. E. Convergence in pigmentation at multiple levels: mutations, genes and function. Philos. Trans. R. Soc. Lond. B 365 , 2439–2450 (2010).
Article CAS Google Scholar
Mundy, N. I. A window on the genetics of evolution: MC1R and plumage colouration in birds. Proc. Biol. Sci. 272 , 1633–1640 (2005).
CAS PubMed PubMed Central Google Scholar
Exposito-Alonso, M. et al. Natural selection on the Arabidopsis thaliana genome in present and future climates. Nature 573 , 126–129 (2019).
Article CAS PubMed Google Scholar
Lortie, C. J. & Hierro, J. L. A synthesis of local adaptation to climate through reciprocal common gardens. J. Ecol. 110 , 1015–1021 (2022).
Monroe, J. G. et al. Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function. eLife 7 , e41038 (2018).
Exposito-Alonso, M. et al. Genomic basis and evolutionary potential for extreme drought adaptation in Arabidopsis thaliana . Nat. Ecol. Evol. 2 , 352–358 (2017).
Lovell, J. T. et al. The genomic landscape of molecular responses to natural drought stress in Panicum hallii . Nat. Commun. 9 , 5213 (2018).
MacLachlan, I. R. et al. Genome-wide shifts in climate-related variation underpin responses to selective breeding in a widespread conifer. Proc. Natl Acad. Sci. USA 118 , e2016900118 (2021).
James, M. E., Brodribb, T., Wright, I. J., Rieseberg, L. H. & Ortiz-Barrientos, D. Replicated evolution in plants. Annu. Rev. Plant Biol. 74 , 697–725 (2023).
Mei, W., Stetter, M. G., Gates, D. J., Stitzer, M. C. & Ross-Ibarra, J. Adaptation in plant genomes: bigger is different. Am. J. Bot. 105 , 16–19 (2018).
Bohutínská, M. & Peichel, C. L. Divergence time shapes gene reuse during repeated adaptation. Trends Ecol. Evol. 39 , 396–407 (2024).
Yeaman, S. Evolution of polygenic traits under global vs local adaptation. Genetics 220 , iyab134 (2022).
Anderson, J. & Song, B.-H. Plant adaptation to climate change—where are we? J. Syst. Evol. 58 , 533–545 (2020).
Christmas, M. J., Breed, M. F. & Lowe, A. J. Constraints to and conservation implications for climate change adaptation in plants. Conserv. Genet. 17 , 305–320 (2016).
Emms, D. M. & Kelly, S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20 , 238 (2019).
Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37 , 4302–4315 (2017).
Booker, T. R., Yeaman, S., Whiting, J. R. & Whitlock, M. C. The WZA: a window-based method for characterizing genotype–environment association. Mol. Ecol. Resour. 24 , e13768 (2024).
Booker, T. R., Yeaman, S. & Whitlock, M. C. Using genome scans to identify genes used repeatedly for adaptation. Evolution 77 , 801–811 (2023).
Conte, G. L., Arnegard, M. E., Peichel, C. L. & Schluter, D. The probability of genetic parallelism and convergence in natural populations. Proc. Biol. Sci. 279 , 5039–5047 (2012).
PubMed PubMed Central Google Scholar
Liu, T., Carlsson, J., Takeuchi, T., Newton, L. & Farré, E. M. Direct regulation of abiotic responses by the Arabidopsis circadian clock component PRR7. Plant J. 76 , 101–114 (2013).
Guan, J. et al. Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding. Nat. Commun. 13 , 5707 (2022).
Bostick, M., Lochhead, S. R., Honda, A., Palmer, S. & Callis, J. Related to ubiquitin 1 and 2 are redundant and essential and regulate vegetative growth, auxin signaling, and ethylene production in Arabidopsis . Plant Cell 16 , 2418–2432 (2004).
Gray, S. B. & Brady, S. M. Plant developmental responses to climate change. Dev. Biol. 419 , 64–77 (2016).
Calleja-Cabrera, J., Boter, M., Oñate-Sánchez, L. & Pernas, M. Root growth adaptation to climate change in crops. Front. Plant Sci. 11 , 544 (2020).
Acharya, B. R. & Assmann, S. M. Hormone interactions in stomatal function. Plant Mol. Biol. 69 , 451–462 (2009).
Luo, Y. et al. An Arabidopsis homolog of importin β1 is required for ABA response and drought tolerance. Plant J. 75 , 377–389 (2013).
Bohutínská, M. et al. Novelty and convergence in adaptation to whole genome duplication. Mol. Biol. Evol. 38 , 3910–3924 (2021).
Wang, L. et al. Molecular parallelism underlies convergent highland adaptation of maize landraces. Mol. Biol. Evol. 38 , 3567–3580 (2021).
Szklarczyk, D. et al. The STRING database in 2021: customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49 , D605–D612 (2021).
Li, N. et al. Plant hormone-mediated regulation of heat tolerance in response to global climate change. Front. Plant Sci. 11 , 627969 (2020).
Fisher, R. A. The Genetical Theory of Natural Selection Vol. 272 (Clarendon, 1930).
Orr, H. A. Adaptation and the cost of complexity. Evolution 54 , 13–20 (2000).
Cooper, T. F., Ostrowski, E. A. & Travisano, M. A negative relationship between mutation pleiotropy and fitness effect in yeast. Evolution 61 , 1495–1499 (2007).
Wang, Z., Liao, B.-Y. & Zhang, J. Genomic patterns of pleiotropy and the evolution of complexity. Proc. Natl Acad. Sci. USA 107 , 18034–18039 (2010).
Wagner, G. P. et al. Pleiotropic scaling of gene effects and the ‘cost of complexity’. Nature 452 , 470–472 (2008).
Rennison, D. J. & Peichel, C. L. Pleiotropy facilitates parallel adaptation in sticklebacks. Mol. Ecol. 31 , 1476–1486 (2022).
Papatheodorou, I. et al. Expression Atlas: gene and protein expression across multiple studies and organisms. Nucleic Acids Res. 46 , D246–D251 (2018).
Obayashi, T., Hibara, H., Kagaya, Y., Aoki, Y. & Kinoshita, K. ATTED-II v11: a plant gene coexpression database using a sample balancing technique by subagging of principal components. Plant Cell Physiol. 63 , 869–881 (2022).
Watanabe, K. et al. A global overview of pleiotropy and genetic architecture in complex traits. Nat. Genet. 51 , 1339–1348 (2019).
Proulx, S. R., Promislow, D. E. L. & Phillips, P. C. Network thinking in ecology and evolution. Trends Ecol. Evol. 20 , 345–353 (2005).
Mack, K. L., Phifer-Rixey, M., Harr, B. & Nachman, M. W. Gene expression networks across multiple tissues are associated with rates of molecular evolution in wild house mice. Genes 10 , 225 (2019).
Yanai, I. et al. Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification. Bioinformatics 21 , 650–659 (2005).
Hahn, M. W. & Kern, A. D. Comparative genomics of centrality and essentiality in three eukaryotic protein-interaction networks. Mol. Biol. Evol. 22 , 803–806 (2005).
Mähler, N. et al. Gene co-expression network connectivity is an important determinant of selective constraint. PLoS Genet. 13 , e1006402 (2017).
Lenormand, T. Gene flow and the limits to natural selection. Trends Ecol. Evol. 17 , 183–189 (2002).
Hämälä, T., Gorton, A. J., Moeller, D. A. & Tiffin, P. Pleiotropy facilitates local adaptation to distant optima in common ragweed ( Ambrosia artemisiifolia ). PLoS Genet. 16 , e1008707 (2020).
Frachon, L. et al. Intermediate degrees of synergistic pleiotropy drive adaptive evolution in ecological time. Nat. Ecol. Evol. 1 , 1551–1561 (2017).
Jain, K. & Stephan, W. Rapid adaptation of a polygenic trait after a sudden environmental shift. Genetics 206 , 389–406 (2017).
Yeaman, S. et al. Convergent local adaptation to climate in distantly related conifers. Science 353 , 1431–1433 (2016).
Kreiner, J. M. et al. Multiple modes of convergent adaptation in the spread of glyphosate-resistant Amaranthus tuberculatus . Proc. Natl Acad. Sci. USA 116 , 21076–21084 (2019).
Fischer, M. C. et al. Estimating genomic diversity and population differentiation—an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri . BMC Genomics 18 , 69 (2017).
Kubota, S. et al. A genome scan for genes underlying microgeographic-scale local adaptation in a wild Arabidopsis species. PLoS Genet. 11 , e1005361 (2015).
Willi, Y., Fracassetti, M., Zoller, S. & Van Buskirk, J. Accumulation of mutational load at the edges of a species range. Mol. Biol. Evol. 35 , 781–791 (2018).
Mattila, T. M., Tyrmi, J., Pyhäjärvi, T. & Savolainen, O. Genome-wide analysis of colonization history and concomitant selection in Arabidopsis lyrata . Mol. Biol. Evol. 34 , 2665–2677 (2017).
Frachon, L. et al. A genomic map of climate adaptation in Arabidopsis thaliana at a micro-geographic scale. Front. Plant Sci. 9 , 967 (2018).
1001 Genomes Consortium. 1,135 genomes reveal the global pattern of polymorphism in Arabidopsis thaliana . Cell 166 , 481–491 (2016).
Wang, B. et al. Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta . Genome Biol. 20 , 126 (2019).
Koenig, D. et al. Long-term balancing selection drives evolution of immunity genes in Capsella . eLife 8 , e43606 (2019).
Murray, K. D. et al. Landscape drivers of genomic diversity and divergence in woodland Eucalyptus . Mol. Ecol. 28 , 5232–5247 (2019).
Todesco, M. et al. Massive haplotypes underlie ecotypic differentiation in sunflowers. Nature 584 , 602–607 (2020).
Yoder, J. B. et al. Genomic signature of adaptation to climate in Medicago truncatula . Genetics 196 , 1263–1275 (2014).
Gould, B. A. et al. Population genomics and climate adaptation of a C4 perennial grass, Panicum hallii (Poaceae). BMC Genomics 19 , 792 (2018).
Karunarathne, P., Zhou, Q., Lascoux, M. & Milesi, P. Hybridization mediated range expansion and climate change resilience in two conifers. Glob. Change Biol. https://doi.org/10.1111/gcb.17262 (2024).
Tyrmi, J. S. et al. Genomics of clinal local adaptation in Pinus sylvestris under continuous environmental and spatial genetic setting. G3 10 , 2683–2696 (2020).
Fahrenkrog, A. M. et al. Population genomics of the eastern cottonwood ( Populus deltoides ). Ecol. Evol. 7 , 9426–9440 (2017).
Wang, J. et al. A major locus controls local adaptation and adaptive life history variation in a perennial plant. Genome Biol. 19 , 72 (2018).
Holliday, J. A., Zhou, L., Bawa, R., Zhang, M. & Oubida, R. W. Evidence for extensive parallelism but divergent genomic architecture of adaptation along altitudinal and latitudinal gradients in Populus trichocarpa . New Phytol. 209 , 1240–1251 (2016).
Leroy, T. et al. Adaptive introgression as a driver of local adaptation to climate in European white oaks. New Phytol. 226 , 1171–1182 (2020).
Booker, T. R., Yeaman, S. & Whitlock, M. C. Using genome scans to identify genes used repeatedly for adaptation. Evolution 77 , 801–811 (2022).
Kumar, S., Stecher, G., Suleski, M. & Hedges, S. B. TimeTree: a resource for timelines, timetrees, and divergence times. Mol. Biol. Evol. 34 , 1812–1819 (2017).
Jasper, R. J. et al. Evaluating the accuracy of variant calling methods using the frequency of parent–offspring genotype mismatch. Mol. Ecol. Resour. 22 , 2524–2533 (2022).
Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34 , i884–i890 (2018).
Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at https://doi.org/10.48550/arXiv.1303.3997 (2013).
Danecek, P. et al. Twelve years of SAMtools and BCFtools. Gigascience 10 , giab008 (2021).
Picard (Broad Institute, 2019); http://broadinstitute.github.io/picard/
DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43 , 491–498 (2011).
Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27 , 2156–2158 (2011).
Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26 , 841–842 (2010).
Van der Auwera, G. A. & O’Connor, B. D. Genomics in the Cloud: Using Docker, GATK, and WDL in Terra (O’Reilly Media, 2020).
Emms, D. M. & Kelly, S. STAG: species tree inference from all genes. Preprint at bioRxiv https://doi.org/10.1101/267914 (2018).
Emms, D. M. & Kelly, S. STRIDE: Species tree root inference from gene duplication events. Mol. Biol. Evol. 34 , 3267–3278 (2017).
Kassambara, A. rstatix: Pipe-friendly framework for basic statistical tests. R package version 0.6.0 (2023).
Tippett, L. H. C. The Methods of Statistics: An Introduction Mainly for Workers in the Biological Sciences (Williams & Norgate, 1931).
Gabaldón, T. & Koonin, E. V. Functional and evolutionary implications of gene orthology. Nat. Rev. Genet. 14 , 360–366 (2013).
Mergner, J. et al. Mass-spectrometry-based draft of the Arabidopsis proteome. Nature 579 , 409–414 (2020).
Yeaman, S. & Whiting, J. R. The genetic architecture of repeated local adaptation to climate in distantly related plants. Dryad https://doi.org/10.5061/dryad.15dv41p57 (2024).
Whiting, J. R. et al. JimWhiting91/RepAdapt: Accepted Manuscript Release (v1.0.0). Zenodo https://doi.org/10.5281/zenodo.12680122 (2024).
Download references
We would like to thank O. Savolainen, T. Hämälä, T. Mitchell-Olds, D. Lowry, M. Kirst, P. Tiffin, S. Kubota and A. Widmer who were involved in the collection and preservation of data used in this study. Funding was provided by NSERC Discovery (RGPIN/03310-2023; S.Y.) and Alberta Innovates (212201729; S.Y.), with computational resources and support provided by the Digital Research Alliance of Canada (S.Y.).
Authors and affiliations.
Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada
James R. Whiting, Clément Rougeux, Brandon M. Lind, Pooja Singh, Mengmeng Lu & Sam Yeaman
Department of Zoology, Faculty of Science, University of British Columbia, Vancouver, British Colombia, Canada
Tom R. Booker & Michael C. Whitlock
Department of Forest and Conservation Sciences, Faculty of Forestry, University of British Columbia, Vancouver, British Columbia, Canada
Tom R. Booker, Brandon M. Lind & Sally N. Aitken
Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
Pooja Singh
EAWAG, Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland
Department of Biological Sciences, University of Notre Dame, Notre Dame, IN, USA
Mengmeng Lu
Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
Kaichi Huang, Julia M. Kreiner, Loren H. Rieseberg & Marco Todesco
School of Environmental and Rural Science, University of New England, Armidale, New South Wales, Australia
Rose L. Andrew, Jeremy J. Bruhl, Momena Khandaker & Ian R. H. Telford
Research School of Biology, Australian National University, Canberra, Australian Capital Territory, Australia
Justin O. Borevitz & Kevin D. Murray
Department of Planning and Environment, Queanbeyan, New South Wales, Australia
Timothy L. Collins
Department of Climate Change, Energy, the Environment and Water, Queanbeyan, New South Wales, Australia
ETH Zurich: Institute of Integrative Biology (IBZ), ETH Zurich, Zurich, Switzerland
Martin C. Fischer
School of Biological Sciences, Monash University, Melbourne, Victoria, Australia
Kathryn A. Hodgins
Department of Forest Resources and Environmental Conservation, Virginia Tech, Blacksburg, VA, USA
Jason A. Holliday & Lecong Zhou
Department of Plant Biology, Swedish University of Agricultural Sciences, Uppsala, Sweden
Pär K. Ingvarsson
Biology Department, Vancouver Island University, Nanaimo, British Columbia, Canada
Jasmine K. Janes
Department of Ecosystem Science and Management, University of Northern British Columbia, Prince George, British Columbia, Canada
Species Survival Commission, Orchid Specialist Group, IUCN North America, Washington, DC, USA
Department of Botany and Plant Sciences, University of California, Riverside, CA, USA
Daniel Koenig
Institute for Integrative Genome Biology, University of California, Riverside, CA, USA
Department of Ecology & Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada
Julia M. Kreiner, John R. Stinchcombe & Stephen I. Wright
UMR BIOGECO, INRAE, Université de Bordeaux; 69 Route d’Arcachon, Cestas, France
Antoine Kremer
Program in Plant Ecology and Evolution, Department of Ecology and Genetics, Evolutionary Biology Centre and Science for Life Laboratory, Uppsala University, Uppsala, Sweden
Martin Lascoux & Pascal Milesi
GenPhySE, Université de Toulouse, INRAE, ENVT, Castanet Tolosan, France
Thibault Leroy
Department of Molecular Biology, Max Planck Institute for Biology Tübingen, Tübingen, Germany
Kevin D. Murray & Detlef Weigel
Department of Forest Sciences, University of Helsinki, Helsinki, Finland
Tanja Pyhäjärvi
Viikki Plant Science Centre, University of Helsinki, Helsinki, Finland
Swiss Federal Research Institute WSL, Birmensdorf, Switzerland
Christian Rellstab
Laboratoire des Interactions Plantes-Microbes-Environnement, Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement, CNRS, Université de Toulouse, Castanet-Tolosan, France
Fabrice Roux
Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia, Canada
Marco Todesco
Department of Biology, University of British Columbia, Kelowna, British Columbia, Canada
Department of Ecology and Genetics, University of Oulu, Oulu, Finland
Jaakko S. Tyrmi
South China National Botanical Garden, Guangzhou, China
Baosheng Wang
Department of Environmental Sciences, University of Basel, Basel, Switzerland
Yvonne Willi
You can also search for this author in PubMed Google Scholar
J.R.W. carried out the analysis and wrote the manuscript. S.Y. conceived the study, helped design the methods and analysis, and contributed to writing the manuscript. C. Rougeux organized the initial data collection and solicited contributions from collaborators. T.R.B. and M.C.W. helped design the methods. B.M.L., P.S., M. Lu and K. Huang contributed to bioinformatic methods and data processing. All other co-authors contributed data and edited the manuscript.
Correspondence to James R. Whiting or Sam Yeaman .
Competing interests.
The authors declare no competing interests.
Peer review information.
Nature Ecology & Evolution thanks Magdalena Bohutinska and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data fig. 1 timetree phylogeny for the focal species studied here..
TimeTree 83 phylogeny for the species studied here. The phylogeny at the top provides estimates of node ages that are assessed across multiple studies that have dated the splits among clades. Asterisks denote cases where a substitute species was selected by TimeTree to be used as data on the target species was unavailable. The substituted species are not the same as those analysed in our study. Substituted species are: Amaranthus hybridus ( A. tuberculatus ); Eucalyptus erythrocorys ( E. albens ); Quercus rubra ( Quercus petraea ). A number of species are missing that were analysed, but congeneric representatives are shown (missing species: Helianthus petiolaris , E. magnificata , E. sideroxylon , Pinus contorta , Picea glauca x engelmannii ).
The species tree that was used for phylogenetic tests is shown in panel A . Panel B shows the distribution of mean ‘K’ values across 1,000 random draws, each of 141 orthogroups. The mean ‘K’ observed in true RAOs is shown as a red line, and the mean of the random distribution is shown as a dashed black line.
Heatmaps show the contribution of individual species to orthogroup repeatability at FDR < 0.3 for different climate variables ( A ) and among pairs of species ( B ), with species ordered phylogenetically. In each case, the fill of each cell represents the proportion of orthogroups where a given species contributes towards the signature of repeatability based on its minimum GEA p-value.
Panels A (FDR <0.5) and B (FDR <0.3) show stacked bars for RAOs identified when removing one species and testing the remaining 24, alongside the full 25 species dataset (vertical dashed line). The heatmap in C shows how the change in RAO number varied by species (LOO CV Change), alongside other features of species datasets including the breadth of sampling (geographically and climatically) relative to the total species range, and technical features related to genome sequencing and sample size (see Supplementary Results 5 ). The association between dataset features and the cross-validation results are shown as correlation coefficients in panel D . Negative correlation coefficients imply that removing datasets with lower dataset feature values increases the number of RAOs, and vice versa.
Heatmap in ( A ) shows the per species -log 10 -transformed GEA p-value for each species and climate variable. Note that isothermality is absent for Eucalyptus magnificata as there was no climate variation here. Alongside the heatmap, the -log 10 -transformed PicMin FDR values are shown, and rows are ordered according to the most significant to least significant. Individual GEA p-value vectors are plotted as histograms in ( B ), with the most significant in terms of PicMin FDR shown in top-left through to least significant in bottom-right.
Each boxplot shows the non-parametric correlation coefficients calculated across all individual datasets (N = 29 biological replicates) between either precipitation climate change or maximum temperature climate change. Each box shows the median, quartiles, standard range (1.5 x IQR) and points show outliers beyond the standard range. Deviations from the central x = 0 are indicative of persistent association between a climate change variable and a given bioclim variable.
Each pair of bars shows the proportion of Arabidopsis thaliana genes belonging to the relevant decile based on either specificity of tissue expression ( A ) or co-expression node degree ( B ), relative to the random expectation (red bars). Deciles are ordered 1-10 from least to most pleiotropic.
Panel A shows a simplified gene tree with 4 species (S1-S4), 3 duplications (D1-3), and 2 species-specific duplications (D2-3). Panel B shows depletion of duplications in orthogroups grouped according to their strongest evidence of repeatability. The mean duplication value is plotted for each repeatability decile. The black points and lines show a comparable analysis over 100 randomisations of all per gene ep-values WZA . The black point shows the mean duplication metric per decile and the lines show the maximum and minimum values across the 100 randomisations. Panel C shows duplication results exclusively looking at contributing species within RAOs, compared with species with low ep-values WZA (<0.1) within 10,000 draws of random orthogroup-climate sets. Lines show the observed mean relative to the distribution of random means, where green lines show permuted p-value < 0.05. The one-sided p-values for permutation tests = 0.153 (Mean Duplication N) and 0.035 (Mean Proportion of Duplicated Genes).
Supplementary information.
Supplementary Methods 1–7, Results 1–7 and Figs. 1–12.
Peer review file, supplementary tables.
Supplementary Tables 1–8. Table legends are included with each sheet in the workbook.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ .
Reprints and permissions
Cite this article.
Whiting, J.R., Booker, T.R., Rougeux, C. et al. The genetic architecture of repeated local adaptation to climate in distantly related plants. Nat Ecol Evol (2024). https://doi.org/10.1038/s41559-024-02514-5
Download citation
Received : 03 November 2023
Accepted : 22 July 2024
Published : 26 August 2024
DOI : https://doi.org/10.1038/s41559-024-02514-5
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.
IMAGES
VIDEO
COMMENTS
Spring Semester Topics: Students enrolled in AIOU's Code 419 educational program in the spring semester face a number of challenges that require careful consideration and application of educational theory. These assignments are designed to assess students' understanding of various educational concepts, teaching methods, and their ability to ...
AIOU Course Code 419 Solved Assignments: The AIOU Course Code 419, titled Education (ایجوکیشن), is an extensive 6-credit hour course comprising four assignments. This course delves into various educational theories, practices, and policies. It aims to provide students with a thorough understanding of the education system, teaching ...
AIOU Solved Assignments Spring 2023: bit.ly/3JtCoiWSkilling Foundation | All About EducationWebsite: Skilling.pkYouTube: http://bit.ly/2DNz53ZFacebook: https...
Education Code 419 Assignment Number 4. Important Note: 1. Allama Iqbal Open University ASSIGNMENTS :-. First and foremost Assignments are an essential part of the students. Because this plays a very important role, moreover, you can't clear your semester without assignments. Most Importantly send your assignments on time.
Click to download. Assignment no# 4. Click to download. To download your assignment, simply click the download button above, and the PDF will be saved to your computer or mobile device within seconds. Remember, submitting assignments is vital for academic success as it showcases your comprehension of the course material and enables instructors ...
⏩ Download Link:⏩ AIOU Code 419 Solved Assignment No.1 & 2 Autumn 2022 https://youtu.be/SiXmc--bM9U⏩ AIOU Code 419 Solved Assignment No.3 & 4 Spring 2022 htt...
aiou ba semester spring 2022 code 419 solved assignment number 1 | easy method of educationAIOU Course Code 419Solved Assignment No.1The Main points that we...
The download process of the AIOU Bachelor BA, B.COM, and associate degree AD assignment is very simple. We have listed all assignments in the shape of a list against each code just scroll down, and click 1st, 2nd, 3rd, or 4th assignment against any code to start downloading. Download Matric FA class assignments for semester Spring 2024 from below:
Download and view past papers for Education code number 419 (AIOU) for the semester of spring and autumn. It is a 6 credit hours book.Old papers will provide a clear understanding of important areas or parts of the subjects. Also previous examinations papers will provide a guess for the upcoming examination.Below table provide past papers.
Here you can find AIOU Assignments Code 419. Assignments for B.A to get a better insight into Exams. Five Years Past Papers for AIOU Matric, F.A, B.A, BS, B.ed, M.ed of all semesters. Assignments for a wide range of undergraduate, graduate, Master.
AIOU علم التعلیم Code 419 Assignment Number 4. 1. Allama Iqbal Open University ASSIGNMENTS :-. First and foremost Assignments are an essential part of the students. Because this plays a very important role. Moreover, you can't clear your semester without assignments. Most Importantly send your assignments in time.
How to Download Education Solved Assignments Spring 2024 for free. That is to say, that It is quite easy to download 419 Education free solved assignments from aiou studio 9. Certainly, you need to scroll all the way down and click the download button it will open a page where you could see all codes select one, and click download pdf.
Code 419 Solved Assignment No.1 Autumn 2021 https://youtu.be/XevKK4xEoTcCode 419 Solved Assignment No.2 Autumn 2021 https://youtu.be/EzDrZIw7FNwCode 419 Solv...
Allama Iqbal Open University. Allama Iqbal Open University has a system of Assignments. One must have to send Assignments the Hand-Written assignments to the AIOU Tutor assigned to him/her. Now you can download all updated AIOU Solved Assignments for BA, BCom in PDF for semester Spring and Autumn. AIOU Solved Assignments of FA, ICom 2021.
Aiou 419 solved past paper 2020 in pdf. Hi dear students, you can download aiou 416 paper in pdf and get help in solving your paper. But remember that it is not recommended that to copy all these available solutions. You only need to get help from here. This site is not responsible if you get fail an exam due to copyrighted material.
Education Code 419 Assignment Number 4. Important Note: 1. Allama Iqbal Open University ASSIGNMENTS :-. First and foremost Assignments are an essential part of the students. Because this plays a very important role, moreover, you can't clear your semester without assignments. Most Importantly send your assignments on time.
419 - Data Science Class IX & X - 2023-2024 Page 2 of 2 CBSE | DEPARTMENT OF SKILL EDUCATION DATA SCIENCE (SUBJECT CODE - 419) CLASS - X (SESSION 2023-2024) Total Marks: 100 (Theory-50 + Practical-50) UNITS MAXIMUM MARKS FOR THEORY AND PRACTICAL EMPLOYABILITY SKILLS Unit 1: Communication Skills-II 2 Unit 2: Self-Management Skills-II 2
Here you can find AIOU Assignments B.A Code 419. Assignments for B.A to get a better insight of B.A exams. Five Years Past Papers for AIOU B.A of all semesters. Assignments for a wide range of undergraduate, graduate, Master's.
CODE OF CONDUCT ... ($419). o The iPad Air or 11" iPad Pro are acceptable, but not required. o The 12.9" iPad Pro is not recommended as its size will be difficult to use in the clinical setting. ... Loss of credit for work involved in violation of code (assignment, examination, quiz) ...
⏩ Download Link:⏩ AIOU Code 419 Solved Assignment No.1 Autumn 2023 https://youtu.be/CtpYgDOWMjg⏩ AIOU Code 419 Solved Assignment No.2 Autumn 2023 https://you...
Multi-agent target assignment and path planning (TAPF) are two key problems in intelligent warehouse. However, most literature only addresses one of these two problems separately. In this study, we propose a method to simultaneously solve target assignment and path planning from a perspective of cooperative multi-agent deep reinforcement learning (RL). To the best of our knowledge, this is the ...
AIOU علم التعلیم Code 419 Assignments of BA B.COM Class Autumn 2021. 20/12/2021. 0 18,130. There are 4 assignments for this course. You will find all the assignments below.
The date change is to provide additional time to implement the policy addressing the assignment of North American Industry Classification System codes to orders placed under multiple-award contracts, as covered by changes made by FAR Case 2014-002, Set-Asides Under Multiple-Award Contracts, 85 FR 11746. DATES: Effective: September 30, 2024.
⏩ Download Link:⏩ AIOU Code 419 Solved Assignment No.1 Spring 2022 https://youtu.be/BRPfrJihax0⏩ AIOU Code 419 Solved Assignment No.2 Spring 2022 https://you...
Here you can find AIOU Assignments B.A Code 419. Assignments for B.A to get better insight into B.A exams. Five Years Past Papers for AIOU B.A of all semesters. Assignments for a wide range of undergraduate, graduate, Masters. M.Phil and Doctoral programs are also available at learningkidunya. Course:Education (419) Assignment#1.
Education Code 419 Assignment Number 4. Important Note: 1. Allama Iqbal Open University ASSIGNMENTS :-. First and foremost Assignments are an essential part of the students. Because this plays a very important role. Moreover, you can't clear your semester without assignments. Most Importantly send your assignments in time.
Accession codes for all data are provided in Supplementary Table 1. ... Orthology assignment. ... Dev. Biol. 419, 64-77 (2016).