Gene duplication and subsequent diversification strongly affect phenotypic evolvability and robustness

Identification of Apple Flower Development-Related Gene Families and Analysis of Transcriptional Regulation

Apple (Malus domestica Borkh.) stands out as a globally significant fruit tree with considerable economic importance. Nonetheless, the orchard production of 'Fuji' apples faces significant challenges, including delayed flowering in young trees and inconsistent annual yields in mature trees, ultimately resulting in suboptimal fruit yield due to insufficient flower bud formation. Flower development represents a pivotal process influencing plant adaptation to environmental conditions and is a crucial determinant of successful plant reproduction. The three gene or transcription factor (TF) families, C2H2, DELLA, and FKF1, have emerged as key regulators in plant flowering regulation; however, understanding their roles during apple flowering remains limited. Consequently, this study identified 24 MdC2H2, 6 MdDELLA, and 6 MdFKF1 genes in the apple genome with high confidence. Through phylogenetic analyses, the genes within each family were categorized into three distinct subgroups, with all facets of protein physicochemical properties and conserved motifs contingent upon subgroup classification. Repetitive events between these three gene families within the apple genome were elucidated via collinearity analysis. qRT-PCR analysis was conducted and revealed significant expression differences among MdC2H2-18, MdDELLA1, and MdFKF1-4 during apple bud development. Furthermore, yeast two-hybrid analysis unveiled an interaction between MdC2H2-18 and MdDELLA1. The genome-wide identification of the C2H2, DELLA, and FKF1 gene families in apples has shed light on the molecular mechanisms underlying apple flower bud development.

Identification and Analysis of the Expression of the PIP5K Gene Family in Tomatoes

To explore the function of phosphatidylinositol 4-phosphate 5-kinase (PIP5K) in tomatoes, members of the tomato PIP5K family were identified and characterized using bioinformatic methods, and their expression patterns were also analyzed under salt stress and in different tissues. Twenty-one PIP5K members-namely, SlPIP5K1-SlPIP5K21-were identified from ten chromosomes, and these were divided into three groups according to a phylogenetic analysis. Further bioinformatic analysis showed four pairs of collinear relationships and fragment replication events among the SlPIP5K family members. To understand the possible roles of the SlPIP5Ks, a cis-acting element analysis was conducted, which indicated that tomato PIP5Ks could be associated with plant growth, hormones, and stress responses. We further validated the results of the in silico analysis by integrating RNA-seq and qRT-PCR techniques for salt- and hormone-treated tomato plants. Our results showed that SlPIP5K genes exhibited tissue- and treatment-specific patterns, and some of the SlPIP5Ks exhibited significantly altered expressions after our treatments, suggesting that they might be involved in these stresses. We selected one of the SlPIP5Ks that responded to our treatments, SlPIP5K2, to further understand its subcellular localization. Our results showed that SlPIP5K2 was located on the membrane. This study lays a foundation for the analysis of the biological functions of the tomato SlPIP5K genes and can also provide a theoretical basis for the selection and breeding of new tomato varieties and germplasm innovation, especially under salt stress.

The non-deterministic genotype-phenotype map of RNA secondary structure

Selection and variation are both key aspects in the evolutionary process. Previous research on the mapping between molecular sequence (genotype) and molecular fold (phenotype) has shown the presence of several structural properties in different biological contexts, implying that these might be universal in evolutionary spaces. The deterministic genotype-phenotype (GP) map that links short RNA sequences to minimum free energy secondary structures has been studied extensively because of its computational tractability and biologically realistic nature. However, this mapping ignores the phenotypic plasticity of RNA. We define a GP map that incorporates non-deterministic (ND) phenotypes, and take RNA as a case study; we use the Boltzmann probability distribution of folded structures and examine the structural properties of ND GP maps for RNA sequences of length 12 and coarse-grained RNA structures of length 30 (RNAshapes30). A framework is presented to study robustness, evolvability and neutral spaces in the ND map. This framework is validated by demonstrating close correspondence between the ND quantities and sample averages of their deterministic counterparts. When using the ND framework we observe the same structural properties as in the deterministic GP map, such as bias, negative correlation between genotypic robustness and evolvability, and positive correlation between phenotypic robustness and evolvability.

Genotype-phenotype mapping with polyominos made from DNA origami tiles

The correlation between genetic information and characteristics of a living cell-its genotype and its phenotype-constitutes the basis of genetics. Here, we experimentally realize a primitive form of genotype-phenotype mapping with DNA origami. The DNA origami can polymerize into two-dimensional lattices (phenotype) via blunt-end stacking facilitated by edge staples at the seam of the planar DNA origami. There are 80 binding positions for edge staples, which allow us to translate an 80-bit long binary code (genotype) onto the DNA origami. The presence of an edge staple thus corresponds to a "1" and its absence to a "0." The interactions of our DNA-based system can be reproduced by a polyomino model. Polyomino growth simulations qualitatively reproduce our experimental results. We show that not only the absolute number of base stacks but also their sequence position determine the cluster size and correlation length of the orientation of single DNA origami within the cluster. Importantly, the mutation of a few bits can result in major morphology changes of the DNA origami cluster, while more often, major sequence changes have no impact. Our experimental realization of a correlation between binary information ("genotype") and cluster morphology ("phenotype") thus reproduces key properties of genotype-phenotype maps known from living systems.

Characterization of the WRKY gene family reveals its contribution to the adaptability of almond (Prunus dulcis)

Background

WRKY (WRKY DNA-binding domain) transcription factors an important gene family that widely regulates plant resistance to biological and abiotic stresses, such as drought, salt and ion stresses. However, research on the WRKY family in almond has not yet been reported. Almond is an economically important fruit tree in Xinjiang that have strong resistance to various stresses.

Results

A total of 62 PdWRKY genes were identified (including six pairs of homologous genes), and the phylogenetic tree was divided into three groups according to the WRKY domain and zinc finger motifs. The members of each group had a significant number of conserved motifs and exons/introns distributed unevenly across eight chromosomes, as well as 24 pairs of fragment duplicates and nine pairs of tandem duplicates. Moreover, the synteny and Ka/Ks analyses of the WRKY genes among almond and distinct species provided more detailed evidence for PdWRKY genes evolution. The examination of different tissue expression patterns showed that PdWRKY genes have tissue-specific expression characteristics. The qRT-PCR results showed that PdWRKY genes participate in the resistance of almond to the effects of low-temperature, drought and salt stress and that the expression levels of these genes change over time, exhibiting spatiotemporal expression characteristics. It is worth noting that many genes play a significant role in low-temperature stress resistance. In addition, based on the conserved WRKY motif, 321 candidate target genes were identified as having functions in multiple pathways.

Conclusions

We conducted systematic bioinformatics analysis and abiotic stress research on the WRKY gene family in almond, laying the foundation for future PdWRKY genes research and improvements to almond production and breeding.

Fast free-energy-based neutral set size estimates for the RNA genotype-phenotype map

The genotype–phenotype (GP) map of RNA secondary structure links each RNA sequence to its corresponding secondary structure. Previous research has shown that the large-scale structural properties of GP maps, such as the size of neutral sets in genotype space, can influence evolutionary outcomes. In order to use neutral set sizes, efficient and accurate computational methods are needed to compute them. Here, we propose a new method, which is based on free energy estimates and is much faster than existing sample-based methods. Moreover, this approach can give insight into the reasons behind neutral set size variations, for example, why structures with fewer stacks tend to have larger neutral set sizes. In addition, we generalize neutral set size calculations from the previously studied many-to-one framework, where each sequence folds into a single energetically preferred structure, to a fuller many-to-many framework, where several low-energy structures are included. We find that structures with high neutral sets in one framework also tend to have large neutral sets in the other framework for a range of parameters and thus the choice of GP map does not fundamentally affect which structures have the largest neutral set sizes.