Recursive Test of Hardy-Weinberg Equilibrium in Tetraploids.

Testing for deviations from Hardy-Weinberg equilibrium (HWE) can provide fundamental information about genetic variation and evolutionary processes in natural populations. In contrast to diploids, where genotype frequencies remain constant after a single episode of random mating, polyploids, characterized by polysomic inheritance, approach HWE gradually. Here, we mathematically show the asymptotic trajectory of tetraploid equilibrium from any initial genotype frequencies. We formulate a statistical framework to test and estimate the degree of deviation from HWE at individual loci in allotetraploids and autotetraploids. Knowledge about HWE test fills an important gap in population genetic studies of tetraploids related to their evolution and ecology.

Epidemiological characteristics of respiratory syncytial virus infection in pediatric patients before, during the COVID-19 pandemic and after easing of COVID-19 restrictive measures in China.

We aimed to assess the epidemiological characteristics of respiratory syncytial virus (RSV) infection in Chinese children at different phases of the coronavirus disease 2019 (COVID-19) pandemic, that is, before, during the pandemic and after easing of restrictive measures. We included 123 623 patients aged 0-18 years with respiratory infection symptoms who were suspected with RSV infection from January 1, 2019 to June 30, 2023 in Hangzhou Children's Hospital. Clinical information and RSV test result were extracted from the laboratory information system. We calculated the positive rate of RSV detection by age groups, gender, seasons, types of patients and phases of COVID-19 pandemic. Nonlinear associations between age and risk of RSV infection in three phases of pandemic were assessed by restricted cubic spline regression models. Among 123 623 patients, 3875 (3.13%) were tested as positive. The highest positive rate was observed in children aged 0-28 days (i.e., 12.28%). RSV infection was most prevalent in winter (6.04%), and followed by autumn (2.52%). Although there is no statistical significance regarding the positive rate at three phases of the pandemic, we observed that the rate was lowest during the pandemic and increased after easing the measures in certain age groups (p < 0.05), which was consisted with results from the nonlinear regression analyses. In addition, regression analyses suggested that the age range of children susceptible to RSV got wider, that is, 0-3.5 years, after easing all restrictive measures compared with that before (i.e., 0-3 years) and during the pandemic (i.e., 0-1 year). Based on our findings, we called for attention from health professionals and caregivers on the new epidemiological characteristics of RSV infection in the post-pandemic era after easing the restrictive measures.

Epidemiological characteristics of norovirus infection in pediatric patients during the COVID-19 pandemic.

To assess the epidemiological characteristics of norovirus infection. We included 5564 patients under the age of 18 years who visited the hospital in which the study took place from December 2020 to November 2022 with a primary diagnosis of acute diarrhea. Clinical information was extracted from the electronic health record system. We calculated the prevalence of norovirus infection by age, gender, season, year, and type of patients. A nonlinear association between age and prevalence rates was assessed using a restricted cubic spline regression model. A total of 5564 patients completed the test for human norovirus, among whom 1442 (25.9%) tested positive. The prevalence of norovirus infection was significantly lower in 2022 than in 2021 (35.9% vs. 53.7%, p < 0.001), and the highest prevalence was observed in winter (35.1%) and then followed by autumn (27.5%). Regarding the age pattern, the highest rate was seen in children aged 1-3 years (37.5%). Children at age 1.5 years may have the highest risk of having norovirus infection (Pnonlinear < 0.001). The prevalence of norovirus infection of norovirus during the COVID-19 pandemic was similar to that before the pandemic shown in literatures. A relatively high rate was observed in cool seasons and in younger children (i.e., 1-3 years).

Design, synthesis, and biological evaluation of deuterated indolepropionic acid derivatives as novel long-acting pan PPARα/γ/δ agonists.

Metabolic syndrome is a complex disease with diverse symptoms, but current pharmacological interventions have limited efficacy. Indeglitazar, a pan-agonist targeting the three-peroxisome proliferator activated receptors (PPAR), exhibits significant therapeutic effects on both diabetic and fatty liver animal models. However, its short half-life limits the in vivo efficacy, which might be attributed to the ß-oxidation of indolepropionic acid at Indeglitazar. To overcome this metabolic instability, two deuterium atoms were introduced to the α-position of indolepropionic acid to block the ß-oxidation. In this study, several deuterated derivatives were found to sustain PPARs activity and extend the half-life of liver microsomes. In oral glucose tolerance tests, I-1 exhibited the strongest glucose-lowering effect on ob/ob mice in this series. In db/db mice, I-1 reduced lipid levels, liver steatosis and promoted UCP1 expression in white adipose tissue. Mechanistic studies further revealed that I-1 exerts stronger effects than Indeglitazar on the regulation of genes related to lipid metabolism, mitochondrial function, and oxidative stress. Furthermore, I-1 significantly reduced liver steatosis, hepatocellular ballooning, inflammation, and fibrosis in NASH model induced by HFD + CCl4, and even exerted better therapeutic effect than that of Indeglitazar. With the above attractive efficacy, deuterated derivative I-1 is considered as a promising treatment for metabolic syndrome.

Genome-Wide Identification of the MAPK and MAPKK Gene Families in Response to Cold Stress in Prunus mume.

Protein kinases of the MAPK cascade family (MAPKKK-MAPKK-MAPK) play an essential role in plant stress response and hormone signal transduction. However, their role in the cold hardiness of Prunus mume (Mei), a class of ornamental woody plant, remains unclear. In this study, we use bioinformatic approaches to assess and analyze two related protein kinase families, namely, MAP kinases (MPKs) and MAPK kinases (MKKs), in wild P. mume and its variety P. mume var. tortuosa. We identify 11 PmMPK and 7 PmMKK genes in the former species and 12 PmvMPK and 7 PmvMKK genes in the latter species, and we investigate whether and how these gene families contribute to cold stress responses. Members of the MPK and MKK gene families located on seven and four chromosomes of both species are free of tandem duplication. Four, three, and one segment duplication events are exhibited in PmMPK, PmvMPK, and PmMKK, respectively, suggesting that segment duplications play an essential role in the expansion and evolution of P. mume and its gene variety. Moreover, synteny analysis suggests that most MPK and MKK genes have similar origins and involved similar evolutionary processes in P. mume and its variety. A cis-acting regulatory element analysis shows that MPK and MKK genes may function in P. mume and its variety's development, modulating processes such as light response, anaerobic induction, and abscisic acid response as well as responses to a variety of stresses, such as low temperature and drought. Most PmMPKs and PmMKKs exhibited tissue-specifific expression patterns, as well as time-specific expression patterns that protect them through cold. In a low-temperature treatment experiment with the cold-tolerant cultivar P. mume 'Songchun' and the cold-sensitive cultivar 'Lve', we find that almost all PmMPK and PmMKK genes, especially PmMPK3/5/6/20 and PmMKK2/3/6, dramatically respond to cold stress as treatment duration increases. This study introduces the possibility that these family members contribute to P. mume's cold stress response. Further investigation is warranted to understand the mechanistic functions of MAPK and MAPKK proteins in P. mume development and response to cold stress.

The chromosome-level genome provides insight into the molecular mechanism underlying the tortuous-branch phenotype of Prunus mume.

Plant with naturally twisted branches is referred to as a tortuous-branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome-level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi-C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi-omics in deciphering the properties of P. mume plant architecture.

HDACs Gene Family Analysis of Eight Rosaceae Genomes Reveals the Genomic Marker of Cold Stress in Prunus mume.

Histone deacetylases (HDACs) play important roles in plant growth, development, and stress response. However, the pattern of how they are expressed in response to cold stress in the ornamental woody plant Prunus mume is poorly understood. Here, we identify 121 RoHDACs from eight Rosaceae plants of which 13 PmHDACs genes are from P. mume. A phylogenetic analysis suggests that the RoHDACs family is classified into three subfamilies, HDA1/RPD3, HD2, and SIR2. We identify 11 segmental duplication gene pairs of RoHDACs and find, via a sequence alignment, that the HDACs gene family, especially the plant-specific HD2 family, has experienced gene expansion and contraction at a recent genome evolution history. Each of the three HDACs subfamilies has its own conserved domains. The expression of PmHDACs in mei is found to be tissue-specific or tissue-wide. RNA-seq data and qRT-PCR experiments in cold treatments suggest that almost all PmHDACs genes-especially PmHDA1/6/14, PmHDT1, and PmSRT1/2-significantly respond to cold stress. Our analysis provides a fundamental insight into the phylogenetic relationship of the HDACs family in Rosaceae plants. Expression profiles of PmHDACs in response to cold stress could provide an important clue to improve the cold hardiness of mei.

The human brain acetylome reveals that decreased acetylation of mitochondrial proteins associates with Alzheimer's disease.

Metabolic changes that correlate to cognitive changes are well-known in Alzheimer's disease (AD). Metabolism is often linked to functional changes in proteins by post-translational modifications. The importance of the regulation of transcription by acetylation is well documented. Advanced mass spectrometry reveals hundreds of acetylated proteins in multiple tissues, but the acetylome of human brain, its functional significance, and the changes with disease are unknown. Filling this gap is critical for understanding the pathophysiology and development of therapies. To fill this gap, we assessed the human brain acetylome in human brain and its changes with AD. More than 5% of the 4,442 proteins from the human brain global proteome were acetylated. Acetylated proteins were primarily found in the cytosol (148), mitochondria (100), nucleus (91), and plasma membrane (58). The comparison of the brain acetylome in controls to that of patients with AD revealed striking and selective differences in terms of its abundances of acetylated peptides/sites. Acetylation of 18 mitochondrial proteins decreased, while acetylation of two cytosolic proteins, tau and GFAP, increased. Our experiments demonstrate that acetylation at some specific lysine sites alters enzyme function. The results indicate that general activation of de-acetylases (i.e., sirtuins) is not an appropriate therapeutic approach for AD.

Stapled, Long-Acting Xenopus GLP-1-Based Dual GLP-1/Glucagon Receptor Agonists with Potent Therapeutic Efficacy for Metabolic Disease.

Novel peptidic glucagon receptor (GCGR) and glucagon-like peptide 1 receptor (GLP-1R) dual agonists are reported to have increased efficacy over GLP-1R monoagonists for the treatment of diabetes and obesity. We identified a novel Xenopus GLP-1-based dual GLP-1R/GCGR agonist (xGLP/GCG-13) designed with a proper activity ratio favoring the GLP-1R versus the GCGR. However, the clinical utility of xGLP/GCG-13 is limited by its short in vivo half-life. Starting from xGLP/GCG-13, dual Cys mutation was performed, followed by covalent side-chain stapling and serum albumin binder incorporation, resulting in a stabilized secondary structure, enhanced agonist potency at GLP-1R and GCGR, and improved stability. The lead peptide 2c (stapled xGLP/GCG-13 analogue with a palmitic acid albumin binder) exhibits balanced GLP-1R and GCGR activations and potent, long-lasting effects on in vivo glucose control. 2c was further explored pharmacologically in diet-induced obesity and db/db rodent models. Chronic administration of 2c potently induced body weight loss and hypoglycemic effects, improved glucose tolerance, increased energy expenditure, and normalized lipid metabolism and adiposity in relevant animal models. These results indicated that 2c has potential for development as a novel antidiabetic and/or antiobesity drug. Furthermore, we propose that the incorporation of a proper serum protein-binding motif into a di-Cys staple is an effective method for improving the stabilities and bioactivities of peptides. This approach is likely applicable to other therapeutic peptides, such as glucose-dependent insulin-tropic peptide receptor (GIPR) and GLP-1R dual agonists or GLP-1R/GCGR/GIPR triagonists.

A unified DNA sequence and non-DNA sequence mapping model of complex traits.

Increasing evidence shows that quantitative inheritance is based on both DNA sequence and non-DNA sequence variants. However, how to simultaneously detect these variants from a mapping study has been unexplored, hampering our effort to illustrate the detailed genetic architecture of complex traits. We address this issue by developing a unified model of quantitative trait locus (QTL) mapping based on an open-pollinated design composed of randomly sampling maternal plants from a natural population and their half-sib seeds. This design forms a two-level hierarchical platform for a joint linkage-linkage disequilibrium analysis of population structure. The EM algorithm was implemented to estimate and test DNA sequence-based effects and non-DNA sequence-based effects of QTLs. We applied this model to analyze genetic mapping data from the OP design of a gymnosperm coniferous species, Torreya grandis, identifying 25 significant DNA sequence and non-DNA sequence QTLs for seedling height and diameter growth in different years. Results from computer simulation show that the unified model has good statistical properties and is powerful for QTL detection. Our model enables the tests of how a complex trait is affected differently by DNA-based effects and non-DNA sequence-based transgenerational effects, thus allowing a more comprehensive picture of genetic architecture to be charted and quantified.

A dissection model for mapping complex traits.

Many quantitative traits are composites of other traits that contribute differentially to genetic variation. Quantitative trait locus (QTL) mapping of these composite traits can benefit by incorporating the mechanistic process of how their formation is mediated by the underlying components. We propose a dissection model by which to map these interconnected components traits under a joint likelihood setting. The model can test how a composite trait is determined by pleiotropic QTLs for its component traits or jointly by different sets of QTLs each responsible for a different component. The model can visualize the pattern of time-varying genetic effects for individual components and their impacts on composite traits. The dissection model was used to map two composite traits, stemwood volume growth decomposed into its stem height, stem diameter and stem form components for Euramerican poplar adult trees, and total lateral root length constituted by its average lateral root length and lateral root number components for Euphrates poplar seedlings. We found the pattern of how QTLs for different components contribute to phenotypic variation in composite traits. The detailed understanding of the genetic machineries of composite traits will not only help in the design of molecular breeding in plants and animals, but also shed light on the evolutionary processes of quantitative traits under natural selection.

Functional mapping of N deficiency-induced response in wheat yield-component traits by implementing high-throughput phenotyping.

As overfertilization leads to environmental concerns and the cost of N fertilizer increases, the issue of how to select crop cultivars that can produce high yields on N-deficient soils has become crucially important. However, little information is known about the genetic mechanisms by which crops respond to environmental changes induced by N signaling. Here, we dissected the genetic architecture of N-induced phenotypic plasticity in bread wheat (Triticum aestivum L.) by integrating functional mapping and semiautomatic high-throughput phenotyping data of yield-related canopy architecture. We identified a set of quantitative trait loci (QTLs) that determined the pattern and magnitude of how wheat cultivars responded to low N stress from normal N supply throughout the wheat life cycle. This analysis highlighted the phenological landscape of genetic effects exerted by individual QTLs, as well as their interactions with N-induced signals and with canopy measurement angles. This information may shed light on our mechanistic understanding of plant adaptation and provide valuable information for the breeding of N-deficiency tolerant wheat varieties.

How trees allocate carbon for optimal growth: insight from a game-theoretic model.

How trees allocate photosynthetic products to primary height growth and secondary radial growth reflects their capacity to best use environmental resources. Despite substantial efforts to explore tree height-diameter relationship empirically and through theoretical modeling, our understanding of the biological mechanisms that govern this phenomenon is still limited. By thinking of stem woody biomass production as an ecological system of apical and lateral growth components, we implement game theory to model and discern how these two components cooperate symbiotically with each other or compete for resources to determine the size of a tree stem. This resulting allometry game theory is further embedded within a genetic mapping and association paradigm, allowing the genetic loci mediating the carbon allocation of stemwood growth to be characterized and mapped throughout the genome. Allometry game theory was validated by analyzing a mapping data of stem height and diameter growth over perennial seasons in a poplar tree. Several key quantitative trait loci were found to interpret the process and pattern of stemwood growth through regulating the ecological interactions of stem apical and lateral growth. The application of allometry game theory enables the prediction of the situations in which the cooperation, competition or altruism is an optimal decision of a tree to fully use the environmental resources it owns.

HpQTL: a geometric morphometric platform to compute the genetic architecture of heterophylly.

Heterophylly, i.e. morphological changes in leaves along the axis of an individual plant, is regarded as a strategy used by plants to cope with environmental change. However, little is known of the extent to which heterophylly is controlled by genes and how each underlying gene exerts its effect on heterophyllous variation. We described a geometric morphometric model that can quantify heterophylly in plants and further constructed an R-based computing platform by integrating this model into a genetic mapping and association setting. The platform, named HpQTL, allows specific quantitative trait loci mediating heterophyllous variation to be mapped throughout the genome. The statistical properties of HpQTL were examined and validated via computer simulation. Its biological relevance was demonstrated by results from a real data analysis of heterophylly in a wood plant, mei (Prunus mume). HpQTL provides a powerful tool to analyze heterophylly and its underlying genetic architecture in a quantitative manner. It also contributes a new approach for genome-wide association studies aimed to dissect the programmed regulation of plant development and evolution.

The genetic architecture of heterochrony as a quantitative trait: lessons from a computational model.

Heterochrony is known as a developmental change in the timing or rate of ontogenetic events across phylogenetic lineages. It is a key concept synthesizing development into ecology and evolution to explore the mechanisms of how developmental processes impact on phenotypic novelties. A number of molecular experiments using contrasting organisms in developmental timing have identified specific genes involved in heterochronic variation. Beyond these classic approaches that can only identify single genes or pathways, quantitative models derived from current next-generation sequencing data serve as a more powerful tool to precisely capture heterochronic variation and systematically map a complete set of genes that contribute to heterochronic processes. In this opinion note, we discuss a computational framework of genetic mapping that can characterize heterochronic quantitative trait loci that determine the pattern and process of development. We propose a unifying model that charts the genetic architecture of heterochrony that perceives and responds to environmental perturbations and evolves over geologic time. The new model may potentially enhance our understanding of the adaptive value of heterochrony and its evolutionary origins, providing a useful context for designing new organisms that can best use future resources.

Novel glucagon- and OXM-based peptides acting through glucagon and GLP-1 receptors with body weight reduction and anti-diabetic properties.

Oxyntomodulin (OXM) is an endogenous gastrointestinal hormone, which activates both the Glucagon-like peptide-1 receptor (GLP-1R) and the glucagon receptor (GCGR). However, OXM has shortcomings including poor GLP-1R agonism to control glycemia, short half-life and others. Inspired from the sequence relationship between OXM and glucagon, in this study, we introduced different C-terminus residues of GLP-1, exenatide and OXM to glucagon to get a series of hybrid peptides with enhanced GLP-1R activation. The formed glucagon-exenatide hybrid peptide shows higher GLP-1R activation properties than OXM. Then the peptides based on the glucagon-exenatide hybrid peptide were coupled with fatty acid side chains to prolong their half-lives. As a result, the most potent compound 16a could stimulate insulin secretion and maintain blood glucose in normal level for ~42.6 h in diabetic mice. 16a exhibited reduced HbA1c level in diabetic mice, lowered body weight significantly in obesity mice on chronic treatment assay. 16a, combined efficient GCGR/GLP-1R activity, is potential as novel treatment for obesity and diabetes. This finding provides new insights into balancing GLP-1/GCGR potency of glucagon-exenatide hybrid peptide and is helpful for discovery of novel anti-diabetic and bodyweight-reducing drugs.

Integrating Evolutionary Game Theory into Mechanistic Genotype-Phenotype Mapping.

Natural selection has shaped the evolution of organisms toward optimizing their structural and functional design. However, how this universal principle can enhance genotype-phenotype mapping of quantitative traits has remained unexplored. Here we show that the integration of this principle and functional mapping through evolutionary game theory gains new insight into the genetic architecture of complex traits. By viewing phenotype formation as an evolutionary system, we formulate mathematical equations to model the ecological mechanisms that drive the interaction and coordination of its constituent components toward population dynamics and stability. Functional mapping provides a procedure for estimating the genetic parameters that specify the dynamic relationship of competition and cooperation and predicting how genes mediate the evolution of this relationship during trait formation.

Revealing conformational dynamics of 2'-O-methyl-RNA guanine modified G-quadruplex by replica exchange molecular dynamics.

Thrombin-binding DNA aptamer (TBA) can fold into an antiparallel unimolecular G-quadruplex (G4) structure. Different types of modifications lead to various effects on the structure and stability of the G4 structure. Previous study has shown that a modified TBA (mTBA) that 2'-deoxy guanine (dG) at positions 10 and 11 in the TBA sequence were replaced by 2'-O-methyl-RNA guanine (2'OMe-G) can't fold into a well-defined G4 structure. In order to investigate the detailed structural information and probe the instability factors, we successfully employed the replica exchange molecular dynamics (REMD) to characterize the large conformational variations of the mTBA and systemically describe the influences of the 2'OMe-G on the mTBA in terms of conformation variations and the probability distributions of Hoogsteen hydrogen bonds, dihedral, sugar pucker and glycosyl torsion angle. Replacing position 10 with the 2'OMe-G (2'OMe-G10) induced a strong destabilization of the aptamer, while the 2'OMe-G at position 11(2'OMe-G11) was less destabilizing. More importantly, the glycosyl torsion angle and sugar pucker of 2'OMe-G10 were the most critical destabilization factors. These results were in good agreement with the theoretical and experimental results. Moreover, the structure information can be used as guidelines for the further design of modifications on G4 structure.

Structural and mechanistic insights into modified G-quadruplex thrombin-binding DNA aptamers.

Thrombin-binding aptamer (TBA) can fold into a G-quadruplex structure necessary for interacting with thrombin. When one thymidine residue of the TGT loop at position 7 is replaced with unlocked uracil (UNA), d-isothymidine (D-isoT) or l-isothymidine (L-isoT), these modified sequences display different activities. To date, the mechanisms of how D/L-isoT and UNA influence the biological properties of TBA have not been illustrated in the literature. In this paper, we fill this gap by probing the structure variations and binding modes of these modified TBAs via molecular dynamics (MD) simulation and free energy calculation. Comparative structural analyses demonstrated that both D-IsoT and UNA changed the local conformation of TGT loop and formed stronger interactions with the target protein. Particularly, D-IsoT and UNA adopted similar conformation which can well explain their similar biological activities. In addition, the flexibility of the two TT loops were described clearly. In contrast, L-IsoT at position 7 led to an obvious tendency to unfold. Free energy calculation and the analysis of key residues energy contributions eventually provide a clear picture of interactions for further understanding of the structure-activity relationships. Collectively, our findings open the way for a rational design of modified aptamers.

A computing platform to map ecological metabolism by integrating functional mapping and the metabolic theory of ecology.

Whole-organism metabolic rate co-varies allometrically with body mass, and is also affected by temperature through different biochemical mechanisms. Here we implement a computational platform to map specific quantitative trait loci (QTLs) that govern the dependence of metabolic rate on size and temperature. The model was formulated within settings of genetic mapping or genome-wide association studies through a mapping population genotyped by a set of molecular markers throughout the genome and phenotyped for metabolic parameters over a range of temperature. The model, estimated by a maximum-likelihood approach, allows a genome-wide search for the underlying metabolic QTLs and the estimation of genotype-specific parameters that specify the metabolism of an organism. Our model provides a tool to detect pleiotropy and epistasis that cause the size- and temperature-dependent change of metabolic rate.