40†years of the homeobox: mechanisms of Hox spatial-temporal collinearity in vertebrates.

Animal body plans are established during embryonic development by the Hox genes. This patterning process relies on the differential expression of Hox genes along the head-to-tail axis. Hox spatial collinearity refers to the relationship between the organization of Hox genes in clusters and the differential Hox expression, whereby the relative order of the Hox genes within a cluster mirrors the spatial sequence of expression in the developing embryo. In vertebrates, the cluster organization is also associated with the timing of Hox activation, which harmonizes Hox expression with the progressive emergence of axial tissues. Thereby, in vertebrates, Hox temporal collinearity is intimately linked to Hox spatial collinearity. Understanding the mechanisms contributing to Hox temporal and spatial collinearity is thus key to the comprehension of vertebrate patterning. Here, we provide an overview of the main discoveries pertaining to the mechanisms of Hox spatial-temporal collinearity.

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes.

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Mining genic resources regulating nitrogen-use efficiency based on integrative biological analyses and their breeding applications in maize and other crops.

Nitrogen (N) is an essential factor for limiting crop yields, and cultivation of crops with low nitrogen-use efficiency (NUE) exhibits increasing environmental and ecological risks. Hence, it is crucial to mine valuable NUE improvement genes, which is very important to develop and breed new crop varieties with high NUE in sustainable agriculture system. Quantitative trait locus (QTL) and genome-wide association study (GWAS) analysis are the most common methods for dissecting genetic variations underlying complex traits. In addition, with the advancement of biotechnology, multi-omics technologies can be used to accelerate the process of exploring genetic variations. In this study, we integrate the substantial data of QTLs, quantitative trait nucleotides (QTNs) from GWAS, and multi-omics data including transcriptome, proteome, and metabolome and further analyze their interactions to predict some NUE-related candidate genes. We also provide the genic resources for NUE improvement among maize, rice, wheat, and sorghum by homologous alignment and collinearity analysis. Furthermore, we propose to utilize the knowledge gained from classical cases to provide the frameworks for improving NUE and breeding N-efficient varieties through integrated genomics, systems biology, and modern breeding technologies.

Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock.

Collinear regulation of Hox genes in space and time has been an outstanding question ever since the initial work of Ed Lewis in 1978. Here we discuss recent advances in our understanding of this phenomenon in relation to novel concepts associated with large-scale regulation and chromatin structure during the development of both axial and limb patterns. We further discuss how this sequential transcriptional activation marks embryonic stem cell-like axial progenitors in mammals and, consequently, how a temporal genetic system is further translated into spatial coordinates via the fate of these progenitors. In this context, we argue the benefit and necessity of implementing this unique mechanism as well as the difficulty in evolving an alternative strategy to deliver this critical positional information.

Transcriptome analysis reveals that PbMYB61 and PbMYB308 are involved in the regulation of lignin biosynthesis in pear fruit stone cells.

Pear fruit stone cells have thick walls and are formed by the secondary deposition of lignin in the primary cell wall of thin-walled cells. Their content and size seriously affect fruit characteristics related to edibility. To reveal the regulatory mechanism underlying stone cell formation during pear fruit development and to identify hub genes, we examined the stone cell and lignin contents of 30 'Shannongsu' pear flesh samples and analyzed the transcriptomes of 15 pear flesh samples collected at five developmental stages. On the basis of the RNA-seq data, 35 874 differentially expressed genes were detected. Additionally, two stone cell-related modules were identified according to a WGCNA. A total of 42 lignin-related structural genes were subsequently obtained. Furthermore, nine hub structural genes were identified in the lignin regulatory network. We also identified PbMYB61 and PbMYB308 as candidate transcriptional regulators of stone cell formation after analyzing co-expression networks and phylogenetic relationships. Finally, we experimentally validated and characterized the candidate transcription factors and revealed that PbMYB61 regulates stone cell lignin formation by binding to the AC element in the PbLAC1 promoter to upregulate expression. However, PbMYB308 negatively regulates stone cell lignin synthesis by binding to PbMYB61 to form a dimer that cannot activate PbLAC1 expression. In this study, we explored the lignin synthesis-related functions of MYB family members. The results presented herein are useful for elucidating the complex mechanisms underlying lignin biosynthesis during pear fruit stone cell development.

Genome-wide identification and molecular evolution of Dof gene family in Camellia oleifera.

DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes' promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera.

Genome-wide identification and molecular evolution of Dof transcription factors in Cyperus esculentus.

Dof transcription factor family in Cyperus esculentus genome was identified and analyzed using bioinformatics. The analysis results revealed that C.esculentus genome contains 29 Dof genes (CesDof), all of which are located in the nucleus according to subcellular localization prediction. CesDof proteinrs have a range of 124 to 512 amino acids, with most being basic proteins. Their secondary structure was mainly irregular curl. The promoter sequence of CesDof genes contains cis-acting elements that respond to light, drought, hormones, low temperature, and circadian rhythm. Codon preference analysis showed that CesDof genes' codon preference ends in T/A. Collinearity analysis revealed that C.esculentus had three pairs of collinear CesDof genes. Additionally, there were 15 pairs of collinear genes between C.esculentus and Arabidopsis thaliana. The genetic relationship between C.esculentus and Rhynchospora pubera was found to be the closest. Phylogenetic tree analysis revealed that 29 CesDof genes of C.esculentus can be classified into 4 subgroups. Additionally, 144 miRNAs were predicted to target these CesDof genes. Furthermore, protein interaction analysis indicated that 15 Dof proteins in C.esculentus had interactions. The qRT-PCR verification results of drought stress and salt stress treatment experiments showed that most CesDof genes were involved in drought stress and salt stress responses, and the gene expression trends under drought stress and salt stress conditions were consistent. These results lay a theoretical foundation for further studying the molecular functions of Dof gene family in C.esculentus and its molecular mechanisms in regulating the life activities of C.esculentus.

The genomes of Australian wild limes.

Australian wild limes occur in highly diverse range of environments and are a unique genetic resource within the genus Citrus. Here we compare the haplotype-resolved genome assemblies of six Australian native limes, including four new assemblies generated using PacBio HiFi and Hi-C sequencing data. The size of the genomes was between 315 and 391 Mb with contig N50s from 29.5 to 35 Mb. Gene completeness of the assemblies was estimated to be from 98.4 to 99.3% and the annotations from 97.7 to 98.9% based upon BUSCO, confirming the high contiguity and completeness of the assembled genomes. High collinearity was observed among the genomes and the two haplotype assemblies for each species. Gene duplication and evolutionary analysis demonstrated that the Australian citrus have undergone only one ancient whole-genome triplication event during evolution. The highest number of species-specific and expanded gene families were found in C. glauca and they were primarily enriched in purine, thiamine metabolism, amino acids and aromatic amino acids metabolism which might help C. glauca to mitigate drought, salinity, and pathogen attacks in the drier environments in which this species is found. Unique genes related to terpene biosynthesis, glutathione metabolism, and toll-like receptors in C. australasica, and starch and sucrose metabolism genes in both C. australis and C. australasica might be important candidate genes for HLB tolerance in these species. Expanded gene families were not lineage specific, however, a greater number of genes related to plant-pathogen interactions, predominantly disease resistant protein, was found in C. australasica and C. australis.

Genome-Wide Identification and Expression Analysis of the COL Gene Family in Hemerocallis citrina Baroni.

Hemerocallis citrina Baroni (H. citrina) is an important specialty vegetable that is not only edible and medicinal but also has ornamental value. However, much remains unknown about the regulatory mechanisms associated with the growth, development, and flowering rhythm of this plant. CO, as a core regulatory factor in the photoperiod pathway, coordinates light and circadian clock inputs to transmit flowering signals. We identified 18 COL genes (HcCOL1-HcCOL18) in the H. citrina cultivar 'Mengzihua' and studied their chromosomal distribution, phylogenetic relationships, gene and protein structures, collinearity, and expression levels in the floral organs at four developmental stages. The results indicate that these genes can be classified into three groups based on phylogenetic analysis. The major expansion of the HcCOL gene family occurred via segmental duplication, and the Ka/Ks ratio indicated that the COL genes of Arabidopsis thaliana, Oryza sativa, Phalaenopsis equestris, and H. citrina were under purifying selection. Many cis-elements, including light response elements, abiotic stress elements, and plant hormone-inducible elements, were distributed in the promoter sequences of the HcCOL genes. Expression analysis of HcCOL genes at four floral developmental stages revealed that most of the HcCOL genes were expressed in floral organs and might be involved in the growth, development, and senescence of the floral organs of H. citrina. This study lays a foundation for the further elucidation of the function of the HcCOL gene in H. citrina and provides a theoretical basis for the molecular design breeding of H. citrina.

Genome-wide analysis of NAC transcription factors and exploration of candidate genes regulating selenium metabolism in Broussonetia papyrifera.

MAIN CONCLUSION: Genome-wide identification revealed 79 BpNAC genes belonging to 16 subfamilies, and their gene structures and evolutionary relationships were characterized. Expression analysis highlighted their importance in plant selenium stress responses. Paper mulberry (Broussonetia papyrifera), a deciduous arboreal plant of the Moraceae family, is distinguished by its leaves, which are abundant in proteins, polysaccharides, and flavonoids, positioning it as a novel feedstock. NAC transcription factors, exclusive to plant species, are crucial in regulating growth, development, and response to biotic and abiotic stress. However, extensive characterization of the NAC family within paper mulberry is lacking. In this study, 79 BpNAC genes were identified from the paper mulberry genome, with an uneven distribution across 13 chromosomes. A comprehensive, genome-wide analysis of BpNACs was performed, including investigating gene structures, promoter regions, and chromosomal locations. Phylogenetic tree analysis, alongside comparisons with Arabidopsis thaliana NACs, allowed for categorizing these genes into 16 subfamilies in alignment with gene structure and motif conservation. Collinearity analysis suggested a significant homologous relationship between the NAC genes of paper mulberry and those in Morus notabilis, Ficus hispida, Antiaris toxicaria, and Cannabis sativa. Integrating transcriptome data and Se content revealed that 12 BpNAC genes were associated with selenium biosynthesis. Subsequent RT-qPCR analysis corroborated the correlation between BpNAC59, BpNAC62 with sodium selenate, and BpNAC55 with sodium selenite. Subcellular localization experiments revealed the nuclear functions of BpNAC59 and BpNAC62. This study highlights the potential BpNAC transcription factors involved in selenium metabolism, providing a foundation for strategically breeding selenium-fortified paper mulberry.

Genome-wide identification, phylogenetic, structural and functional evolution of the core components of ABA signaling in plant species: a focus on rice.

MAIN CONCLUSION: A genome-wide analysis had identified 642 ABA core component genes from 20 plant species, which were further categorized into three distinct subfamilies. The gene structures and evolutionary relationships of these genes had been characterized. PP2C_1, PP2C_2, and SnRK2_1 had emerged as key players in mediating the ABA signaling transduction pathway, specifically in rice, in response to abiotic stresses. The plant hormone abscisic acid (ABA) is essential for growth, development, and stress response, relying on its core components, pyrabactin resistance, pyrabactin resistance-like, and the regulatory component of ABA receptor (PYR/PYL/RCAR), 2C protein phosphatase (PP2C), sucrose non-fermenting-1-related protein kinase 2 (SnRK2). However, there's a lack of research on their structural evolution and functional differentiation across plants. Our study analyzed the phylogenetic, gene structure, homology, and duplication evolution of this complex in 20 plant species. We found conserved patterns in copy number and homology across subfamilies. Segmental and tandem duplications drove the evolution of these genes, while whole-genome duplication (WGD) expanded PYR/PYL/RCAR and PP2C subfamilies, enhancing environmental adaptation. In rice and Arabidopsis, the PYR/PYL/RCAR, PP2C, and SnRK2 genes showed distinct tissue-specific expression and responded to various stresses. Notably, PP2C_1 and PP2C_2 interacted with SnRK2_1 and were crucial for ABA signaling in rice. These findings offered new insights into ABA signaling evolution, interactions, and integration in green plants, benefiting future research in agriculture, evolutionary biology, ecology, and environmental science.

A systematic genome-wide analysis and screening of the NAC family genes related to wood formation in Cinnamomum camphora.

Many processes, such as growth, aging, and adaptation to abiotic stress, are regulated in plants by NAC transcription factors. In woody plants, NAC transcription factors acts as a primary switch that regulates secondary xylem development by activating various downstream transcription factors and modulating expression levels of genes involved in the synthesis of the secondary cell wall. Our team had previously sequenced the whole genome of the camphor tree (Cinnamomum camphora). Here, we performed a detailed analysis of the NAC gene family of C. camphora and examined its evolutionary history. The genomic sequences of 121 NAC genes of C. camphora were identified and classified into 20 subfamilies in 2 major classes based on the phylogenetic analysis and structural features. Expansion of the CcNAC gene family occurred mainly by fragment replication and was influenced by the purifying selection. By analyzing predicted interactions of the homologous AtNAC proteins, we identified five CcNACs that potentially regulate xylem development in C. camphora. RNA sequencing revealed distinct expression profiles of CcNACs in seven different plant tissues. Subcellular localization prediction revealed that 120, 3, and 2 CcNACs have biological functions in the nucleus, cytoplasm, and chloroplast, respectively. Furthermore, we verified expression patterns of five CcNACs (CcNAC012, CcNAC028, CcNAC055, CcNAC080, and CcNAC119) in various tissue types using qRT-PCR. Our results will facilitate further in-depth studies of the molecular mechanisms by which CcNAC transcription factors regulate wood formation and other processes in C. camphora.

Understanding the Consequences of Collinearity for Multilevel Models: The Importance of Disaggregation Across Levels.

In multilevel models, disaggregating predictors into level-specific parts (typically accomplished via centering) benefits parameter estimates and their interpretations. However, the importance of level-specificity has been sparsely addressed in multilevel literature concerning collinearity. In this study, we develop novel insights into the interactivity of centering and collinearity in multilevel models. After integrating the broad literatures on centering and collinearity, we review level-specific and conflated correlations in multilevel data. Next, by deriving formal relationships between predictor collinearity and multilevel model estimates, we demonstrate how the consequences of collinearity change across different centering specifications and identify data characteristics that may exacerbate or mitigate those consequences. We show that when all or some level-1 predictors are uncentered, slope estimates can be greatly biased by collinearity. Disaggregation of all predictors eliminates the possibility that fixed effect estimates will be biased due to collinearity alone; however, under some data conditions, collinearity is associated with biased standard errors and random effect (co)variance estimates. Finally, we illustrate the importance of disaggregation for diagnosing collinearity in multilevel data and provide recommendations for the use of level-specific collinearity diagnostics. Overall, the necessity of disaggregation for identifying and managing collinearity's consequences in multilevel models is clarified in novel ways.

Molecular Characteristics and Expression Patterns of Carotenoid Cleavage Oxygenase Family Genes in Rice (Oryza sativa L.).

Carotenoid cleavage oxygenases (CCOs) cleave carotenoid molecules to produce bioactive products that influence the synthesis of hormones such as abscisic acid (ABA) and strigolactones (SL), which regulate plant growth, development, and stress adaptation. Here, to explore the molecular characteristics of all members of the OsCCO family in rice, fourteen OsCCO family genes were identified in the genome-wide study. The results revealed that the OsCCO family included one OsNCED and four OsCCD subfamilies. The OsCCO family was phylogenetically close to members of the maize ZmCCO family and the Sorghum SbCCO family. A collinearity relationship was observed between OsNCED3 and OsNCED5 in rice, as well as OsCCD7 and OsNCED5 between rice and Arabidopsis, Sorghum, and maize. OsCCD4a and OsCCD7 were the key members in the protein interaction network of the OsCCO family, which was involved in the catabolic processes of carotenoids and terpenoid compounds. miRNAs targeting OsCCO family members were mostly involved in the abiotic stress response, and RNA-seq data further confirmed the molecular properties of OsCCO family genes in response to abiotic stress and hormone induction. qRT-PCR analysis showed the differential expression patterns of OsCCO members across various rice organs. Notably, OsCCD1 showed relatively high expression levels in all organs except for ripening seeds and endosperm. OsNCED2a, OsNCED3, OsCCD1, OsCCD4a, OsCCD7, OsCCD8a, and OsCCD8e were potentially involved in plant growth and differentiation. Meanwhile, OsNCED2a, OsNCED2b, OsNCED5, OsCCD8b, and OsCCD8d were associated with reproductive organ development, flowering, and seed formation. OsNCED3, OsCCD4b, OsCCD4c, OsCCD8b, and OsCCD8c were related to assimilate transport and seed maturation. These findings provide a theoretical basis for further functional analysis of the OsCCO family.

Phylogeny of the plant receptor-like kinase (RLK) gene family and expression analysis of wheat RLK genes in response to biotic and abiotic stresses.

BACKGROUND: The receptor-like kinase (RLK) gene families in plants contains a large number of members. They are membrane proteins with an extracellular receptor domain and participate in biotic and abiotic stress responses. RESULTS: In this study, we identified RLKs in 15 representative plant genomes, including wheat, and classified them into 64 subfamilies by using four types of phylogenetic trees and HMM models. Conserved exon‒intron structures with conserved exon phases in the kinase domain were found in many RLK subfamilies from Physcomitrella patens to Triticum aestivum. Domain distributions of RLKs were also diagrammed. Collinearity events and tandem gene clusters suggested that polyploidization and tandem duplication events contributed to the member expansions of T. aestivum RLKs. Global expression pattern analysis was performed by using public transcriptome data. These analyses were involved in T. aestivum, Aegilops tauschii and Brachypodium distachyon RLKs under biotic and abiotic stresses. We also selected 9 RLKs to validate the transcriptome prediction by using qRT‒PCR under drought treatment and with Fusarium graminearum infection. The expression trends of these 9 wheat RLKs from public transcriptome data were consistent with the results of qRT‒PCR, indicating that they might be stress response genes under drought or F. graminearum treatments. CONCLUSION: In this study, we identified, classified, evolved, and expressed RLKs in wheat and related plants. Thus, our results will provide insights into the evolutionary history and molecular mechanisms of wheat RLKs.

Genome wide characterization of R2R3 MYB transcription factor from Apocynum venetum revealed potential stress tolerance and flavonoid biosynthesis genes.

MYB transcription factors are crucial in regulating stress tolerance and expression of major genes involved in flavonoid biosynthesis. The functions of MYBs is well explored in a number of plants, yet no study is reported in Apocynum venetum. We identified a total of 163 MYB candidates, that comprised of 101 (61.96%) R2R3, 6 3R, 1 4R and 55 1R. Syntenic analysis of A. venetum R2R3 (AvMYBs) showed highest orthologous pairs with Vitis vinifera MYBs followed by Arabidopsis thaliana among the four species evaluated. Thirty segmental duplications and 6 tandem duplications were obtained among AvMYB gene pairs signifying their role in the MYB gene family expansion. Nucleotide substitution analysis (Ka/Ks) showed the AvMYBs to be under the influence of strong purifying selection. Expression analysis of selected AvMYBs under low temperature and cadmium stresses resulted in the identification of AvMYB48, AvMYB97, AvMYB8, AvMYB4 as potential stress responsive genes and AvMYB10 and AvMYB11 in addition, proanthocyanidin biosynthesis regulatory genes which is consistent with their annotated homologues in Arabidopsis. Tissue specific expression profile analysis of the AvMYBs further supported the qPCR analysis result. MYBs with higher transcript levels in root, stem and leaf like AvMYB4 for example, was downregulated under the stresses and such with low transcript level such as AvMYB48 which had low transcript in the leaf was upregulated under both stresses. Transcriptome and phylogenetic analyses suggested AvMYB42 as a potential regulator of anthocyanin biosynthesis. Thus, this study provided valuable information on AvR2R3-MYB gene family with respect to stress tolerance and flavonoid biosynthesis.

The BIDIAP index: a clinical, analytical and ultrasonographic score for the diagnosis of acute appendicitis in children.

BACKGROUND: Pediatric acute appendicitis (PAA) continues to be a diagnostic challenge today. The diagnostic performance of classical indices is only moderate, especially in pediatric population. This study aimed to define a clinical, radiological and analytical index for the diagnosis of PAA. MATERIALS AND METHODS: This prospective study included 151 patients divided into two groups: (1) 53 patients with non-surgical abdominal pain (NSAP) and (2) 98 patients with a confirmed PAA. Sociodemographic and clinical characteristics were compared between groups using the Mann-Whitney U test and the Fisher exact test. To identify the predictors of PAA, we performed a multivariable logistic regression using a forward stepwise analysis and we assigned multiples of integer values to the selected variables. The diagnostic performance of the index was assessed by calculating the area under the receiver operating characteristic curve. Intra-cohort calibration was assessed with the Hosmer-Lemeshow test. RESULTS: We developed the BIDIAP index (BIomarkers for the DIagnosis of Appendicitis in Pediatrics), which included three variables that independently predicted higher odds of PAA: appendiceal caliber (≥ 6.9 mm), systemic immune-inflammation index (≥ 890) and peritoneal irritation, which scored 4, 3 and 2 points, respectively. Mean (SD) score of the participants was 2.38 (2.06) in group 1 and 7.89 (1.50) in group 2. The area under the ROC was 0.97 (95% CI 0.95-0.99). The cut-off point was established at 4 points, resulting in a sensitivity of 98.98% and a specificity of 77.78%. CONCLUSIONS: The BIDIAP index has an exceptional diagnostic performance in PAA. The importance of these results lies in its novelty and in the simplicity of the index. Although external validation will be necessary, initial results look promising.

Collinearity Issues in Autoregressive Models with Time-Varying Serially Dependent Covariates.

First-order autoregressive models are popular to assess the temporal dynamics of a univariate process. Researchers often extend these models to include time-varying covariates, such as contextual factors, to investigate how they moderate processes' dynamics. We demonstrate that doing so has implications for how well one can estimate the autoregressive and covariate effects, as serial dependence in the variables can imply predictor collinearity. This is a noteworthy contribution, since in current practice serial dependence in a time-varying covariate is rarely considered important. We first recapitulate the role of predictor collinearity for estimation precision in an ordinary least squares context, by discussing how it affects estimator variances, covariances and correlations. We then derive a general formula detailing how predictor collinearity in first-order autoregressive models is impacted by serial dependence in the covariate. We provide a simulation study to illustrate the implications of the formula for different types of covariates. The simulation results highlight when the collinearity issue becomes severe enough to hamper interpretation of the effects. We also show that the effect estimates can be biased in small samples (i.e., 50 time points). Implications for study design, the use of time as a predictor, and related model variants are discussed.

Genome-Wide Identification of MADS-Box Genes in Taraxacum kok-saghyz and Taraxacum mongolicum: Evolutionary Mechanisms, Conserved Functions and New Functions Related to Natural Rubber Yield Formation.

MADS-box transcription regulators play important roles in plant growth and development. However, very few MADS-box genes have been isolated in the genus Taraxacum, which consists of more than 3000 species. To explore their functions in the promising natural rubber (NR)-producing plant Taraxacum kok-saghyz (TKS), MADS-box genes were identified in the genome of TKS and the related species Taraxacum mongolicum (TM; non-NR-producing) via genome-wide screening. In total, 66 TkMADSs and 59 TmMADSs were identified in the TKS and TM genomes, respectively. From diploid TKS to triploid TM, the total number of MADS-box genes did not increase, but expansion occurred in specific subfamilies. Between the two genomes, a total of 11 duplications, which promoted the expansion of MADS-box genes, were identified in the two species. TkMADS and TmMADS were highly conserved, and showed good collinearity. Furthermore, most TkMADS genes exhibiting tissue-specific expression patterns, especially genes associated with the ABCDE model, were preferentially expressed in the flowers, suggesting their conserved and dominant functions in flower development in TKS. Moreover, by comparing the transcriptomes of different TKS lines, we identified 25 TkMADSs related to biomass formation and 4 TkMADSs related to NR content, which represented new targets for improving the NR yield of TKS.

Genome-Wide Characterization and Analysis of bHLH Transcription Factors Related to Anthocyanin Biosynthesis in Cinnamomum camphora ('Gantong 1').

Cinnamomum camphora is one of the most commonly used tree species in landscaping. Improving its ornamental traits, particularly bark and leaf colors, is one of the key breeding goals. The basic helix-loop-helix (bHLH) transcription factors (TFs) are crucial in controlling anthocyanin biosynthesis in many plants. However, their role in C. camphora remains largely unknown. In this study, we identified 150 bHLH TFs (CcbHLHs) using natural mutant C. camphora 'Gantong 1', which has unusual bark and leaf colors. Phylogenetic analysis revealed that 150 CcbHLHs were divided into 26 subfamilies which shared similar gene structures and conserved motifs. According to the protein homology analysis, we identified four candidate CcbHLHs that were highly conserved compared to the TT8 protein in A. thaliana. These TFs are potentially involved in anthocyanin biosynthesis in C. camphora. RNA-seq analysis revealed specific expression patterns of CcbHLHs in different tissue types. Furthermore, we verified expression patterns of seven CcbHLHs (CcbHLH001, CcbHLH015, CcbHLH017, CcbHLH022, CcbHLH101, CcbHLH118, and CcbHLH134) in various tissue types at different growth stages using qRT-PCR. This study opens a new avenue for subsequent research on anthocyanin biosynthesis regulated by CcbHLH TFs in C. camphora.