The pattern of genetic variability in a core collection of 2,021 cowpea accessions.

Cowpea is a highly drought-adapted leguminous crop with great promise for improving agricultural sustainability and food security. Here, we report analyses derived from array-based genotyping of 2,021 accessions constituting a core subset of the world's largest cowpea collection, held at the International Institute of Tropical Agriculture (IITA) in Ibadan, Nigeria. We used this dataset to examine genetic variation and population structure in worldwide cowpea. We confirm that the primary pattern of population structure is two geographically defined subpopulations originating in West and East Africa, respectively, and that population structure is associated with shifts in phenotypic distribution. Furthermore, we establish the cowpea core collection as a resource for genome-wide association studies by mapping the genetic basis of several phenotypes, with a focus on seed coat pigmentation patterning and color. We anticipate that the genotyped IITA Cowpea Core Collection will serve as a powerful tool for mapping complex traits, facilitating the acceleration of breeding programs to enhance the resilience of this crop in the face of rapid global climate change.

Effects of internal cooling on physical performance, physiological and perceptional parameters when exercising in the heat: A systematic review with meta-analyses.

Background: An elevated core temperature (Tcore) increases the risk of performance impairments and heat-related illness. Internal cooling (IC) has the potential to lower Tcore when exercising in the heat. The aim of the review was to systematically analyze the effects of IC on performance, physiological, and perceptional parameters. Methods: A systematic literature search was performed in the PubMed database on 17 December 2021. Intervention studies were included assessing the effects of IC on performance, physiological, or perceptional outcomes. Data extraction and quality assessment were conducted for the included literature. The standardized mean differences (SMD) and 95% Confidence Intervals (CI) were calculated using the inverse-variance method and a random-effects model. Results: 47 intervention studies involving 486 active subjects (13.7% female; mean age 20-42 years) were included in the meta-analysis. IC resulted in significant positive effects on time to exhaustion [SMD (95% CI) 0.40 (0.13; 0.67), p < 0.01]. IC significantly reduced Tcore [-0.19 (22120.34; -0.05), p < 0.05], sweat rate [-0.20 (-0.34; -0.06), p < 0.01], thermal sensation [-0.17 (-0.33; -0.01), p < 0.05], whereas no effects were found on skin temperature, blood lactate, and thermal comfort (p > 0.05). IC resulted in a borderline significant reduction in time trial performance [0.31 (-0.60; -0.02), p = 0.06], heart rate [-0.13 (-0.27; 0.01), p = 0.06], rate of perceived exertion [-0.16 (-0.31; -0.00), p = 0.05] and borderline increased mean power output [0.22 (0.00; 0.44), p = 0.05]. Discussion: IC has the potential to affect endurance performance and selected physiological and perceptional parameters positively. However, its effectiveness depends on the method used and the time point of administration. Future research should confirm the laboratory-based results in the field setting and involve non-endurance activities and female athletes. Systematic review registration: https://www.crd.york.ac.uk/PROSPERO/, identifier: CRD42022336623.

On the hidden temporal dynamics of plant adaptation.

Adaptation to a wide range of environments is a major driver of plant diversity. It is now possible to catalog millions of potential adaptive genomic differences segregating between environments within a plant species in a single experiment. Understanding which of these changes contributes to adaptive phenotypic divergence between plant populations is a major goal of evolutionary biologists and crop breeders. In this review, we briefly highlight the approaches frequently used to understand the genetic basis of adaptive phenotypes in plants, and we discuss some of the limitations of these methods. We propose that direct observation of the process of adaptation using multigenerational studies and whole genome sequencing is a crucial missing component of recent studies of plant adaptation because it complements several shortcomings of sampling-based techniques.

A KARRIKIN INSENSITIVE2 paralog in lettuce mediates highly sensitive germination responses to karrikinolide.

Karrikins (KARs) are chemicals in smoke that can enhance germination of many plants. Lettuce (Lactuca sativa) cv. Grand Rapids germinates in response to nanomolar karrikinolide (KAR1). Lettuce is much less responsive to KAR2 or a mixture of synthetic strigolactone analogs, rac-GR24. We investigated the molecular basis of selective and sensitive KAR1 perception in lettuce. The lettuce genome contains two copies of KARRIKIN INSENSITIVE2 (KAI2), which in Arabidopsis (Arabidopsis thaliana) encodes a receptor that is required for KAR responses. LsKAI2b is more highly expressed than LsKAI2a in dry achenes and during early stages of imbibition. Through cross-species complementation assays in Arabidopsis, we found that an LsKAI2b transgene confers robust responses to KAR1, but LsKAI2a does not. Therefore, LsKAI2b likely mediates KAR1 responses in lettuce. We compared homology models of KAI2 proteins from lettuce and a fire-follower, whispering bells (Emmenanthe penduliflora). This identified pocket residues 96, 124, 139, and 161 as candidates that influence the ligand specificity of KAI2. Further support for the importance of these residues was found through a broader comparison of pocket residues among 281 KAI2 proteins from 184 asterid species. Almost all KAI2 proteins had either Tyr or Phe identity at position 124. Genes encoding Y124-type KAI2 are more broadly distributed in asterids than in F124-type KAI2. Substitutions at residues 96, 124, 139, and 161 in Arabidopsis KAI2 produced a broad array of responses to KAR1, KAR2, and rac-GR24. This suggests that the diverse ligand preferences observed among KAI2 proteins in plants could have evolved through relatively few mutations.

The Fate of Deleterious Variants in a Barley Genomic Prediction Population.

Targeted identification and purging of deleterious genetic variants has been proposed as a novel approach to animal and plant breeding. This strategy is motivated, in part, by the observation that demographic events and strong selection associated with cultivated species pose a "cost of domestication." This includes an increase in the proportion of genetic variants that are likely to reduce fitness. Recent advances in DNA resequencing and sequence constraint-based approaches to predict the functional impact of a mutation permit the identification of putatively deleterious SNPs (dSNPs) on a genome-wide scale. Using exome capture resequencing of 21 barley lines, we identified 3855 dSNPs among 497,754 total SNPs. We generated whole-genome resequencing data of Hordeum murinum ssp. glaucum as a phylogenetic outgroup to polarize SNPs as ancestral vs. derived. We also observed a higher proportion of dSNPs per synonymous SNPs (sSNPs) in low-recombination regions of the genome. Using 5215 progeny from a genomic prediction experiment, we examined the fate of dSNPs over three breeding cycles. Adjusting for initial frequency, derived alleles at dSNPs reduced in frequency or were lost more often than other classes of SNPs. The highest-yielding lines in the experiment, as chosen by standard genomic prediction approaches, carried fewer homozygous dSNPs than randomly sampled lines from the same progeny cycle. In the final cycle of the experiment, progeny selected by genomic prediction had a mean of 5.6% fewer homozygous dSNPs relative to randomly chosen progeny from the same cycle.

Comprehensive lifestyle intervention vs soy protein-based meal regimen in non-alcoholic steatohepatitis.

BACKGROUND: Non-alcoholic steatohepatitis (NASH) has become one of the leading causes of liver disease in the western world. In obese patients weight reduction is recommended. Up to now there are no specific guidelines for weight loss in order to reduce hepatic fat content. AIM: To investigate the effects of a 24-wk guided lifestyle intervention program compared to a meal replacement regimen based on soy protein. METHODS: Twenty-six subjects with NASH participated in a randomized single-center study. They were randomly assigned to either meal replacement group (MR-G) with soy-yogurt-honey preparation or to guided lifestyle change group (LC-G) with endurance activity and nutrition counselling. Serum alanine transaminase (ALT), aspartate transaminase (AST), lipid parameters, and adipokines were measured. Liver fat content and lipid composition were determined by magnetic resonance imaging and magnetic resonance spectroscopy. Body fat mass and lean body mass were assessed using Bod Pod® device. Pre- and post-intervention monitoring of parameters was performed. Statistical analyses were conducted with SPSS software, results were expressed as median (interquartile range). RESULTS: Twenty-two subjects (MR-G, n = 11 and LC-G, n = 11) completed the study (9 women, 13 men; age 52.1 (15.0) years, body mass index (BMI) 32.3 (3.3) kg/m²). In both groups a significant weight loss was achieved (MR-G: -6.4 (3.6) kg, P < 0.01; LC-G: -9.1 (10.4) kg, P < 0.01). BMI dropped in both groups (MR-G: -2.3 (1.5) kg/m2, P = 0.003; LC-G: -3.0 (3.4) kg/m2, P = 0.006). Internal fat and hepatic lipid content were markedly reduced in both groups in comparable amount. There was a strong correlation between reduction in liver fat and decrease in ALT. Likewise, both groups showed an improvement in glycemic control and lipid profile. Changes in adipokines, particularly in adiponectin and leptin were closely related to intrahepatic lipid changes. CONCLUSION: Comprehensive lifestyle intervention and meal replacement regimen have comparable effects on body and liver fat, as well as decrease in markers of hepatic inflammation among NASH patients.

Long-term balancing selection drives evolution of immunity genes in Capsella.

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses.

The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

Genetic architecture of nonadditive inheritance in Arabidopsis thaliana hybrids.

The ubiquity of nonparental hybrid phenotypes, such as hybrid vigor and hybrid inferiority, has interested biologists for over a century and is of considerable agricultural importance. Although examples of both phenomena have been subject to intense investigation, no general model for the molecular basis of nonadditive genetic variance has emerged, and prediction of hybrid phenotypes from parental information continues to be a challenge. Here we explore the genetics of hybrid phenotype in 435 Arabidopsis thaliana individuals derived from intercrosses of 30 parents in a half diallel mating scheme. We find that nonadditive genetic effects are a major component of genetic variation in this population and that the genetic basis of hybrid phenotype can be mapped using genome-wide association (GWA) techniques. Significant loci together can explain as much as 20% of phenotypic variation in the surveyed population and include examples that have both classical dominant and overdominant effects. One candidate region inherited dominantly in the half diallel contains the gene for the MADS-box transcription factor AGAMOUS-LIKE 50 (AGL50), which we show directly to alter flowering time in the predicted manner. Our study not only illustrates the promise of GWA approaches to dissect the genetic architecture underpinning hybrid performance but also demonstrates the contribution of classical dominance to genetic variance.

A sister of PIN1 gene in tomato (Solanum lycopersicum) defines leaf and flower organ initiation patterns by maintaining epidermal auxin flux.

The spatiotemporal localization of the plant hormone auxin acts as a positional cue during early leaf and flower organogenesis. One of the main contributors to auxin localization is the auxin efflux carrier PIN-FORMED1 (PIN1). Phylogenetic analysis has revealed that PIN1 genes are split into two sister clades; PIN1 and the relatively uncharacterized Sister-Of-PIN1 (SoPIN1). In this paper we identify entire-2 as a loss-of-function SlSoPIN1a (Solyc10g078370) mutant in Solanum lycopersicum. The entire-2 plants are unable to specify proper leaf initiation leading to a frequent switch from the wild type spiral phyllotactic pattern to distichous and decussate patterns. Leaves in entire-2 are large and less complex and the leaflets display spatial deformities in lamina expansion, vascular development, and margin specification. During sympodial growth in entire-2 the specification of organ position and identity is greatly affected resulting in variable branching patterns on the main sympodial and inflorescence axes. To understand how SlSoPIN1a functions in establishing proper auxin maxima we used the auxin signaling reporter DR5: Venus to visualize differences in auxin localization between entire-2 and wild type. DR5: Venus visualization shows a widening of auxin localization which spreads to subepidermal tissue layers during early leaf and flower organogenesis, showing that SoPIN1 functions to focus auxin signaling to the epidermal layer. The striking spatial deformities observed in entire-2 help provide a mechanistic framework for explaining the function of the SoPIN1 clade in S.lycopersicum.

Improving the Annotation of Arabidopsis lyrata Using RNA-Seq Data.

Gene model annotations are important community resources that ensure comparability and reproducibility of analyses and are typically the first step for functional annotation of genomic regions. Without up-to-date genome annotations, genome sequences cannot be used to maximum advantage. It is therefore essential to regularly update gene annotations by integrating the latest information to guarantee that reference annotations can remain a common basis for various types of analyses. Here, we report an improvement of the Arabidopsis lyrata gene annotation using extensive RNA-seq data. This new annotation consists of 31,132 protein coding gene models in addition to 2,089 genes with high similarity to transposable elements. Overall, ~87% of the gene models are corroborated by evidence of expression and 2,235 of these models feature multiple transcripts. Our updated gene annotation corrects hundreds of incorrectly split or merged gene models in the original annotation, and as a result the identification of alternative splicing events and differential isoform usage are vastly improved.

Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella.

The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.

Beyond the thale: comparative genomics and genetics of Arabidopsis relatives.

For decades a small number of model species have rightly occupied a privileged position in laboratory experiments, but it is becoming increasingly clear that our knowledge of biology is greatly improved when informed by a broader diversity of species and evolutionary context. Arabidopsis thaliana has been the primary model organism for plants, benefiting from a high-quality reference genome sequence and resources for reverse genetics. However, recent studies have made a group of species also in the Brassicaceae family and closely related to A. thaliana a focal point for comparative molecular, genomic, phenotypic and evolutionary studies. In this Review, we emphasize how such studies complement continued study of the model plant itself, provide an evolutionary perspective and summarize our current understanding of genetic and phenotypic diversity in plants.

Evolution of DNA methylation patterns in the Brassicaceae is driven by differences in genome organization.

DNA methylation is an ancient molecular modification found in most eukaryotes. In plants, DNA methylation is not only critical for transcriptionally silencing transposons, but can also affect phenotype by altering expression of protein coding genes. The extent of its contribution to phenotypic diversity over evolutionary time is, however, unclear, because of limited stability of epialleles that are not linked to DNA mutations. To dissect the relative contribution of DNA methylation to transposon surveillance and host gene regulation, we leveraged information from three species in the Brassicaceae that vary in genome architecture, Capsella rubella, Arabidopsis lyrata, and Arabidopsis thaliana. We found that the lineage-specific expansion and contraction of transposon and repeat sequences is the main driver of interspecific differences in DNA methylation. The most heavily methylated portions of the genome are thus not conserved at the sequence level. Outside of repeat-associated methylation, there is a surprising degree of conservation in methylation at single nucleotides located in gene bodies. Finally, dynamic DNA methylation is affected more by tissue type than by environmental differences in all species, but these responses are not conserved. The majority of DNA methylation variation between species resides in hypervariable genomic regions, and thus, in the context of macroevolution, is of limited phenotypic consequence.

Mating system shifts and transposable element evolution in the plant genus Capsella.

BACKGROUND: Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora. RESULTS: We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species. CONCLUSION: Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

The genome of the stress-tolerant wild tomato species Solanum pennellii.

Solanum pennellii is a wild tomato species endemic to Andean regions in South America, where it has evolved to thrive in arid habitats. Because of its extreme stress tolerance and unusual morphology, it is an important donor of germplasm for the cultivated tomato Solanum lycopersicum. Introgression lines (ILs) in which large genomic regions of S. lycopersicum are replaced with the corresponding segments from S. pennellii can show remarkably superior agronomic performance. Here we describe a high-quality genome assembly of the parents of the IL population. By anchoring the S. pennellii genome to the genetic map, we define candidate genes for stress tolerance and provide evidence that transposable elements had a role in the evolution of these traits. Our work paves a path toward further tomato improvement and for deciphering the mechanisms underlying the myriad other agronomic traits that can be improved with S. pennellii germplasm.

Truncation of LEAFY COTYLEDON1 protein is required for asexual reproduction in Kalanchoë daigremontiana.

Kalanchoë daigremontiana reproduces asexually by generating numerous plantlets on its leaf margins. The formation of plantlets requires the somatic initiation of organogenic and embryogenic developmental programs in the leaves. However, unlike normal embryogenesis in seeds, leaf somatic embryogenesis bypasses seed dormancy to form viable plantlets. In Arabidopsis (Arabidopsis thaliana), seed dormancy and embryogenesis are initiated by the transcription factor LEAFY COTYLEDON1 (LEC1). The K. daigremontiana ortholog of LEC1 is expressed during leaf somatic embryo development. However, KdLEC1 encodes for a LEC1-type protein that has a unique B domain, with 11 unique amino acids and a premature stop codon. Moreover, the truncated KdLEC1 protein is not functional in Arabidopsis. Here, we show that K. daigremontiana transgenic plants expressing a functional, chimeric KdLEC1 gene under the control of Arabidopsis LEC1 promoter caused several developmental defects to leaf somatic embryos, including seed dormancy characteristics. The dormant plantlets also behaved as typical dormant seeds. Transgenic plantlets accumulated oil bodies and responded to the abscisic acid biosynthesis inhibitor fluridone, which broke somatic-embryo dormancy and promoted their normal development. Our results indicate that having a mutated form of LEC1 gene in K. daigremontiana is essential to bypass dormancy in the leaf embryos and generate viable plantlets, suggesting that the loss of a functional LEC1 promotes viviparous leaf somatic embryos and thus enhances vegetative propagation in K. daigremontiana. Mutations resulting in truncated LEC1 proteins may have been of a selective advantage in creating somatic propagules, because such mutations occurred independently in several Kalanchoë species, which form plantlets constitutively.

The Capsella rubella genome and the genomic consequences of rapid mating system evolution.

The shift from outcrossing to selfing is common in flowering plants, but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible <200,000 years ago. We report a C. rubella reference genome sequence and compare RNA expression and polymorphism patterns between C. rubella and its outcrossing progenitor Capsella grandiflora. We found a clear shift in the expression of genes associated with flowering phenotypes, similar to that seen in Arabidopsis, in which self fertilization evolved about 1 million years ago. Comparisons of the two Capsella species showed evidence of rapid genome-wide relaxation of purifying selection in C. rubella without a concomitant change in transposable element abundance. Overall we document that the transition to selfing may be typified by parallel shifts in gene expression, along with a measurable reduction of purifying selection.

Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato.

Although applied over extremely short timescales, artificial selection has dramatically altered the form, physiology, and life history of cultivated plants. We have used RNAseq to define both gene sequence and expression divergence between cultivated tomato and five related wild species. Based on sequence differences, we detect footprints of positive selection in over 50 genes. We also document thousands of shifts in gene-expression level, many of which resulted from changes in selection pressure. These rapidly evolving genes are commonly associated with environmental response and stress tolerance. The importance of environmental inputs during evolution of gene expression is further highlighted by large-scale alteration of the light response coexpression network between wild and cultivated accessions. Human manipulation of the genome has heavily impacted the tomato transcriptome through directed admixture and by indirectly favoring nonsynonymous over synonymous substitutions. Taken together, our results shed light on the pervasive effects artificial and natural selection have had on the transcriptomes of tomato and its wild relatives.