Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding

Modeling 0.6 million genes for the rational design of functional cis-regulatory variants and de novo design of cis-regulatory sequences

Rational design of plant cis-regulatory DNA sequences without expert intervention or prior domain knowledge is still a daunting task. Here, we developed PhytoExpr, a deep learning framework capable of predicting both mRNA abundance and plant species using the proximal regulatory sequence as the sole input. PhytoExpr was trained over 17 species representative of major clades of the plant kingdom to enhance its generalizability. Via input perturbation, quantitative functional annotation of the input sequence was achieved at single-nucleotide resolution, revealing an abundance of predicted high-impact nucleotides in conserved noncoding sequences and transcription factor binding sites. Evaluation of maize HapMap3 single-nucleotide polymorphisms (SNPs) by PhytoExpr demonstrates an enrichment of predicted high-impact SNPs in cis-eQTL. Additionally, we provided two algorithms that harnessed the power of PhytoExpr in designing functional cis-regulatory variants, and de novo creation of species-specific cis-regulatory sequences through in silico evolution of random DNA sequences. Our model represents a general and robust approach for functional variant discovery in population genetics and rational design of regulatory sequences for genome editing and synthetic biology.

Development of near homozygous lines for diploid hybrid TPS breeding in potatoes

Diploid inbred-based F1 hybrid True Potato Seed (DHTPS) breeding is a novel technique to transform potato breeding and cultivation across the globe. Significant efforts are being made to identify elite diploids, dihaploids and develop diploid inbred lines for heterosis exploitation in potatoes. Self-incompatibility is the first obstacle for developing inbred lines in diploid potatoes, which necessitates the introgression of a dominant S locus inhibitor gene (Sli) for switching self-incompatibility to self-compatibility. We evaluated a set of 357 diploid clones in different selfing generations for self-compatibility and degree of homozygosity using Kompetitive Allele Specific PCR (KASP) Single Nucleotide Polymorphism (SNP) markers. A subset of 10 KASP markers of the Sli candidate region on chromosome 12 showed an association with the phenotype for self-compatibility. The results revealed that the selected 10 KASP markers for the Sli gene genotype could be deployed for high throughput rapid screening of self-compatibility in diploid populations and to identify new sources of self-compatibility. The homozygosity assessed through 99 KASP markers distributed across all the chromosomes of the potato genome was 20-78 % in founder diploid clones, while different selfing generations, i.e., S0, S1, S2 and S3 observed 36.1-80.4, 56.9-82.8, 59.5-85.4 and 73.7-87.8 % average homozygosity, respectively. The diploid plants with ∼80 % homozygosity were also observed in the first selfing generation, which inferred that homozygosity assessment in the early generations itself could identify the best plants with high homozygosity to speed up the generation of diploid inbred lines.

Cotton transposon-related variome reveals roles of transposon-related variations in modern cotton cultivation

INTRODUCTION

Transposon plays a vital role in cotton genome evolution, contributing to the expansion and divergence of genomes within the Gossypium genus. However, knowledge of transposon activity in modern cotton cultivation is limited.

OBJECTIVES

In this study, we aimed to construct transposon-related variome within Gossypium genus and reveal role of transposon-related variations during cotton cultivation. In addition, we try to identify valuable transposon-related variations for cotton breeding.

METHODS

We utilized graphical genome construction to build up the graphical transposon-related variome. Based on the graphical variome, we integrated t-test, eQTL analysis and Mendelian Randomization (MR) to identify valuable transposon activities and elite genes. In addition, a convolutional neural network (CNN) model was constructed to evaluate epigenomic effects of transposon-related variations.

RESULTS

We identified 35,980 transposon activities among 10 cotton genomes, and the diversity of genomic and epigenomic features was observed among 21 transposon categories. The graphical cotton transposon-related variome was constructed, and 9,614 transposon-related variations with plasticity in the modern cotton cohort were used for eQTL, phenotypic t-test and Mendelian Randomization. 128 genes were identified as gene resources improving fiber length and strength simultaneously. 4 genes were selected from 128 genes to construct the elite gene panel whose utility has been validated in a natural cotton cohort and 2 accessions with phenotypic divergence. Based on the eQTL analysis results, we identified transposon-related variations involved in cotton's environmental adaption and human domestication, providing evidence of their role in cotton's adaption-domestication cooperation.

CONCLUSIONS

The cotton transposon-related variome revealed the role of transposon-related variations in modern cotton cultivation, providing genomic resources for cotton molecular breeding.

A survey of the Sli gene in wild and cultivated potato

Inbred-hybrid breeding of diploid potatoes necessitates breeding lines that are self-compatible. One way of incorporating self-compatibility into incompatible cultivated potato (Solanum tuberosum) germplasm is to introduce the S-locus inhibitor gene (Sli), which functions as a dominant inhibitor of gametophytic self-incompatibility. To learn more about Sli diversity and function in wild species relatives of cultivated potato, we obtained Sli gene sequences that extended from the 5'UTR to the 3'UTR from 133 individuals from 22 wild species relatives of potato and eight diverse cultivated potato clones. DNA sequence alignment and phylogenetic trees based on genomic and protein sequences show that there are two highly conserved groups of Sli sequences. DNA sequences in one group contain the 533 bp insertion upstream of the start codon identified previously in self-compatible potato. The second group lacks the insertion. Three diploid and four polyploid individuals of wild species collected from geographically disjointed localities contained Sli with the 533 bp insertion. For most of the wild species clones examined, however, Sli did not have the insertion. Phylogenetic analysis indicated that Sli sequences with the insertion, in wild species and in cultivated clones, trace back to a single origin. Some diploid wild potatoes that have Sli with the insertion were self-incompatible and some wild potatoes that lack the insertion were self-compatible. Although there is evidence of positive selection for some codon positions in Sli, there is no evidence of diversifying selection at the gene level. In silico analysis of Sli protein structure did not support the hypothesis that amino acid changes from wild-type (no insertion) to insertion-type account for changes in protein function. Our study demonstrated that genetic factors besides the Sli gene must be important for conditioning a switch in the mating system from self-incompatible to self-compatible in wild potatoes.

Fingerprint Finder: Identifying Genomic Fingerprint Sites in Cotton Cohorts for Genetic Analysis and Breeding Advancement

Genomic data in Gossypium provide numerous data resources for the cotton genomics community. However, to fill the gap between genomic analysis and breeding field work, detecting the featured genomic items of a subset cohort is essential for geneticists. We developed FPFinder v1.0 software to identify a subset of the cohort's fingerprint genomic sites. The FPFinder was developed based on the term frequency-inverse document frequency algorithm. With the short-read sequencing of an elite cotton pedigree, we identified 453 pedigree fingerprint genomic sites and found that these pedigree-featured sites had a role in cotton development. In addition, we applied FPFinder to evaluate the geographical bias of fiber-length-related genomic sites from a modern cotton cohort consisting of 410 accessions. Enriching elite sites in cultivars from the Yangtze River region resulted in the longer fiber length of Yangze River-sourced accessions. Apart from characterizing functional sites, we also identified 12,536 region-specific genomic sites. Combining the transcriptome data of multiple tissues and samples under various abiotic stresses, we found that several region-specific sites contributed to environmental adaptation. In this research, FPFinder revealed the role of the cotton pedigree fingerprint and region-specific sites in cotton development and environmental adaptation, respectively. The FPFinder can be applied broadly in other crops and contribute to genetic breeding in the future.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these two fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we derive a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., Genome-Wide Association Studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur using analytical theory and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate this by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study of this, we re-examine an analysis testing for co-evolution of expression levels between genes across a fungal phylogeny, and show that including covariance matrix eigenvectors as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Enhancing cold resistance in Banana (Musa spp.) through EMS-induced mutagenesis, L-Hyp pressure selection: phenotypic alterations, biomass composition, and transcriptomic insights

BACKGROUND

The cultivation of bananas encounters substantial obstacles, particularly due to the detrimental effects of cold stress on their growth and productivity. A potential remedy that has gained attention is the utilization of ethyl mesylate (EMS)-induced mutagenesis technology, which enables the creation of a genetically varied group of banana mutants. This complex procedure entails subjecting the mutants to further stress screening utilizing L-Hyp in order to identify those exhibiting improved resistance to cold. This study conducted a comprehensive optimization of the screening conditions for EMS mutagenesis and L-Hyp, resulting in the identification of the mutant cm784, which exhibited remarkable cold resistance. Subsequent investigations further elucidated the physiological and transcriptomic responses of cm784 to low-temperature stress.

RESULTS

EMS mutagenesis had a substantial effect on banana seedlings, resulting in modifications in shoot and root traits, wherein a majority of seedlings exhibited delayed differentiation and limited elongation. Notably, mutant leaves displayed altered biomass composition, with starch content exhibiting the most pronounced variation. The application of L-Hyp pressure selection aided in the identification of cold-resistant mutants among seedling-lethal phenotypes. The mutant cm784 demonstrated enhanced cold resistance, as evidenced by improved survival rates and reduced symptoms of chilling injury. Physiological analyses demonstrated heightened activities of antioxidant enzymes and increased proline production in cm784 when subjected to cold stress. Transcriptome analysis unveiled 946 genes that were differentially expressed in cm784, with a notable enrichment in categories related to 'Carbohydrate transport and metabolism' and 'Secondary metabolites biosynthesis, transport, and catabolism'.

CONCLUSION

The present findings provide insights into the molecular mechanisms that contribute to the heightened cold resistance observed in banana mutants. These mechanisms encompass enhanced carbohydrate metabolism and secondary metabolite biosynthesis, thereby emphasizing the adaptive strategies employed to mitigate the detrimental effects induced by cold stress.

Telomere-to-telomere genome of the allotetraploid legume Sesbania cannabina reveals transposon-driven subgenome divergence and mechanisms of alkaline stress tolerance

Alkaline soils pose an increasing problem for agriculture worldwide, but using stress-tolerant plants as green manure can improve marginal land. Here, we show that the legume Sesbania cannabina is very tolerant to alkaline conditions and, when used as a green manure, substantially improves alkaline soil. To understand genome evolution and the mechanisms of stress tolerance in this allotetraploid legume, we generated the first telomere-to-telomere genome assembly of S. cannabina spanning ∼2,087 Mb. The assembly included all centromeric regions, which contain centromeric satellite repeats, and complete chromosome ends with telomeric characteristics. Further genome analysis distinguished A and B subgenomes, which diverged approximately 7.9 million years ago. Comparative genomic analysis revealed that the chromosome homoeologs underwent large-scale inversion events (>10 Mb) and a significant, transposon-driven size expansion of the chromosome 5A homoeolog. We further identified four specific alkali-induced phosphate transporter genes in S. cannabina; these may function in alkali tolerance by relieving the deficiency in available phosphorus in alkaline soil. Our work highlights the significance of S. cannabina as a green tool to improve marginal lands and sheds light on subgenome evolution and adaptation to alkaline soils.

Unlikely heroes on the long and winding road to potato inbreeding

Conversion of potato from a tetraploid, heterozygous, vegetatively propagated crop to a diploid F1 hybrid, propagated via botanical seed, would constitute a considerable advance for global agriculture, but faces multiple challenges. One such challenge is the difficulty in inbreeding potato, which involves purging deleterious alleles from its genome. This commentary discusses possible reasons for this difficulty and highlights a recent sequence-based effort to classify SNP variation, in potato germplasm, according to its deleterious potential. Tools and strategies connected to this database may facilitate development of F1 hybrids.

Dissected Leaf 1 encodes an MYB transcription factor that controls leaf morphology in potato

KEY MESSAGE

The transcription factor StDL1 regulates dissected leaf formation in potato and the genotype frequency of recessive Stdl1/Stdl1, which results in non-dissected leaves, has increased in cultivated potatoes. Leaf morphology is a key trait of plants, influencing plant architecture, photosynthetic efficiency and yield. Potato (Solanum tuberosum L.), the third most important food crop worldwide, has a diverse leaf morphology. However, despite the recent identification of several genes regulating leaf formation in other plants, few genes involved in potato leaf development have been reported. In this study, we identified an R2R3 MYB transcription factor, Dissected Leaf 1 (StDL1), regulating dissected leaf formation in potato. A naturally occurring allele of this gene, Stdl1, confers non-dissected leaves in young seedlings. Knockout of StDL1 in a diploid potato changes the leaf morphology from dissected to non-dissected. Experiments in N. benthamiana and yeast show that StDL1 is a transcriptional activator. Notably, by calculating the genotype frequency of the Stdl1/Stdl1 in 373-potato accessions, we found that it increases significantly in cultivated potatoes. This work reveals the genetic basis of dissected leaf formation in potato and provides insights into plant leaf morphology.

Unraveling Prevalence and Effects of Deleterious Mutations in Maize Elite Lines across Decades of Modern Breeding

Future breeding is likely to involve the detection and removal of deleterious alleles, which are mutations that negatively affect crop fitness. However, little is known about the prevalence of such mutations and their effects on phenotypic traits in the context of modern crop breeding. To address this, we examined the number and frequency of deleterious mutations in 350 elite maize inbred lines developed over the past few decades in China and the United States. Our findings reveal an accumulation of weakly deleterious mutations and a decrease in strongly deleterious mutations, indicating the dominant effects of genetic drift and purifying selection for the two types of mutations, respectively. We also discovered that slightly deleterious mutations, when at lower frequencies, were more likely to be heterozygous in the developed hybrids. This is consistent with complementation as a potential explanation for heterosis. Subsequently, we found that deleterious mutations accounted for more of the variation in phenotypic traits than nondeleterious mutations with matched minor allele frequencies, especially for traits related to leaf angle and flowering time. Moreover, we detected fewer deleterious mutations in the promoter and gene body regions of differentially expressed genes across breeding eras than in nondifferentially expressed genes. Overall, our results provide a comprehensive assessment of the prevalence and impact of deleterious mutations in modern maize breeding and establish a useful baseline for future maize improvement efforts.

Unloading potatoes: Potato breeding moves forward with only half the genome

Potato is one of the most important food crops in the world and, in contrast to other staples, has not seen large improvements in yield. Agha, Shannon, and Morrell preview an article recently published in Cell, "Phylogenomic discovery of deleterious mutations facilitates hybrid potato breeding," which advances potato breeding strategies via a genetic approach.