Multiple Sources of Introduction of North American Arabidopsis thaliana from across Eurasia

Biologically meaningful genome interpretation models to address data underdetermination for the leaf and seed ionome prediction in Arabidopsis thaliana

Genome interpretation (GI) encompasses the computational attempts to model the relationship between genotype and phenotype with the goal of understanding how the first leads to the second. While traditional approaches have focused on sub-problems such as predicting the effect of single nucleotide variants or finding genetic associations, recent advances in neural networks (NNs) have made it possible to develop end-to-end GI models that take genomic data as input and predict phenotypes as output. However, technical and modeling issues still need to be fixed for these models to be effective, including the widespread underdetermination of genomic datasets, making them unsuitable for training large, overfitting-prone, NNs. Here we propose novel GI models to address this issue, exploring the use of two types of transfer learning approaches and proposing a novel Biologically Meaningful Sparse NN layer specifically designed for end-to-end GI. Our models predict the leaf and seed ionome in A.thaliana, obtaining comparable results to our previous over-parameterized model while reducing the number of parameters by 8.8 folds. We also investigate how the effect of population stratification influences the evaluation of the performances, highlighting how it leads to (1) an instance of the Simpson's Paradox, and (2) model generalization limitations.

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Genetic adaptations in the population history of Arabidopsis thaliana

A population encounters a variety of environmental stresses, so the full source of its resilience can only be captured by collecting all the signatures of adaptation to the selection of the local environment in its population history. Based on the multiomic data of Arabidopsis thaliana, we constructed a database of phenotypic adaptations (p-adaptations) and gene expression (e-adaptations) adaptations in the population. Through the enrichment analysis of the identified adaptations, we inferred a likely scenario of adaptation that is consistent with the biological evidence from experimental work. We analyzed the dynamics of the allele frequencies at the 23,880 QTLs of 174 traits and 8,618 eQTLs of 1,829 genes with respect to the total SNPs in the genomes and identified 650 p-adaptations and 3,925 e-adaptations [false discovery rate (FDR) = 0.05]. The population underwent large-scale p-adaptations and e-adaptations along 4 lineages. Extremely cold winters and short summers prolonged seed dormancy and expanded the root system architecture. Low temperatures prolonged the growing season, and low light intensity required the increased chloroplast activity. The subtropical and humid environment enhanced phytohormone signaling pathways in response to the biotic and abiotic stresses. Exposure to heavy metals selected alleles for lower heavy metal uptake from soil, lower growth rate, lower resistance to bacteria, and higher expression of photosynthetic genes were selected. The p-adaptations are directly interpretable, while the coadapted gene expressions reflect the physiological requirements for the adaptation. The integration of this information characterizes when and where the population has experienced environmental stress and how the population responded at the molecular level.

Microbe-associated molecular pattern recognition receptors have little effect on endophytic Arabidopsis thaliana microbiome assembly in the field

Plant microbiome structure affects plant health and productivity. A limited subset of environmental microbes successfully establishes within plant tissues, but the forces underlying this selectivity remain poorly characterized. Transmembrane pattern recognition receptors (PRRs), used by plants to detect microbe-associated molecular patterns (MAMPs), are strong candidates for achieving this selectivity because PRRs can potentially interact with many members of the microbiome. Indeed, MAMPs found in many microbial taxa, including beneficials and commensals, can instigate a robust immune response that affects microbial growth. Surprisingly, we found that MAMP-detecting PRRs have little effect on endophytic bacterial and fungal microbiome structure in the field. We compared the microbiomes of four PRR knockout lines of Arabidopsis thaliana to wild-type plants in multiple tissue types over several developmental stages and detected only subtle shifts in fungal, but not bacterial, β-diversity in one of the four PRR mutants. In one developmental stage, lore mutants had slightly altered fungal β-diversity, indicating that LORE may be involved in plant-fungal interactions in addition to its known role in detecting certain bacterial lipids. No other effects of PRRs on α-diversity, microbiome variability, within-individual homogeneity, or microbial load were found. The general lack of effect suggests that individual MAMP-detecting PRRs are not critical in shaping the endophytic plant microbiome. Rather, we suggest that MAMP-detecting PRRs must either act in concert and/or are individually maintained through pleiotropic effects or interactions with coevolved mutualists or pathogens. Although unexpected, these results offer insights into the role of MAMP-detecting PRRs in plant-microbe interactions and help direct future efforts to uncover host genetic elements that control plant microbiome assembly.

Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

Genetic diversity pattern reveals the primary determinant of burcucumber (Sicyos angulatus L.) invasion in Korea

Biological invasion is a complex process associated with propagule pressure, dispersal ability, environmental constraints, and human interventions, which leave genetic signatures. The population genetics of an invasive species thus provides invaluable insights into the patterns of invasion. Burcucumber, one of the most detrimental weeds for soybean production in US, has recently colonized Korea and rapidly spread posing a great threat to the natural ecosystem. We aim to infer the determinants of the rapid burcucumber invasion by examining the genetic diversity, demography, and spread pattern with advanced genomic tools. We employed 2,696 genome-wide single-nucleotide polymorphisms to assess the level of diversity and the spatial pattern associated with the landscape factors and to infer the demographic changes of 24 populations (364 genotypes) across four major river basins with the east coastal streams in South Korea. Through the approximate Bayesian computation, we inferred the likely invasion scenario of burcucumber in Korea. The landscape genetics approach adopting the circuit theory and MaxEnt model was applied to determine the landscape contributors. Our data suggested that most populations have experienced population bottlenecks, which led to lowered within-population genetic diversity and inflated population divergences. Burcucumber colonization in Korea has strongly been affected by demographic bottlenecks and multiple introductions, whereas environmental factors were not the primary determinant of the invasion. Our work highlighted the significance of preventing secondary introductions, particularly for aggressive weedy plants such as the burcucumber.

Complex population history affects admixture analyses in nine-spined sticklebacks

Introgressive hybridization is an important process in evolution but challenging to identify, undermining the efforts to understand its role and significance. On the contrary, many analytical methods assume direct descent from a single common ancestor, and admixture among populations can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome sequences of nine-spined sticklebacks (Pungitius pungitius) revealed a complex pattern of population ancestry involving multiple waves of gene flow and introgression across northern Europe. The two recognized lineages were found to have drastically different histories, and their secondary contact zone was wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious deviations from the expected pattern. Interestingly, the freshwater isolates provided peeks into the past and helped to understand the intermediate states of evolutionary processes. Our analyses and findings paint a detailed picture of the complex colonization history of northern Europe and provide backdrop against which introgression and its role in evolution can be investigated. However, they also expose the challenges in analyses of admixed populations and demonstrate how hidden admixture and colonization history misleads the estimation of admixture proportions and population split times.