Search

BackgroundSalidroside (SAL), the main component of Rhodiola rosea extract, is a flavonoid with biological activities, such as antioxidative stress, anti-inflammatory, and hypolipidemic. In this study, the potential therapeutic targets and mechanisms of SAL against oxidative stress in retinal ganglion cells (RGCs) were investigated on the basis of in-vitro experiments, network pharmacology, and molecular docking techniques.MethodsRGC oxidative stress models were constructed, and cell activity, reactive oxygen species (ROS), and apoptosis levels were examined for differences. The genes corresponding to rhodopsin, RGCs, and oxidative stress were screened from GeneCards, TCMSP database, and an analysis platform. The intersection of the three was taken, and a Venn diagram was drawn. Protein interactions, GO functional enrichment, and KEGG pathway enrichment data were analyzed by STRING database, Cytohubba plugin, and Metascape database. The key factors in the screening pathway were validated using qRT-PCR. Finally, molecular docking prediction was performed using MOE 2019 software, molecular dynamic simulations was performed using Gromacs 2018 software.ResultsIn the RGC oxidative stress model in vitro, the cell activity was enhanced, ROS was reduced, and apoptosis was decreased after SAL treatment. A total of 16 potential targets of oxidative stress in SAL RGCs were obtained, and the top 10 core targets were screened by network topology analysis. GO analysis showed that SAL retinal oxidative stress treatment mainly involved cellular response to stress, transcriptional regulatory complexes, and DNA-binding transcription factor binding. KEGG analysis showed that most genes were mainly enriched in multiple cancer pathways and signaling pathways in diabetic complications, nonalcoholic fatty liver, and lipid and atherosclerosis. Validation by PCR, molecular docking and molecular dynamic simulations revealed that SAL may attenuate oxidative stress and reduce apoptosis in RGCs by regulating SIRT1, NRF2, and NOS3.ConclusionThis study initially revealed the antioxidant therapeutic effects and molecular mechanisms of SAL on RGCs, providing a theoretical basis for subsequent studies.

BackgroundNumerous observational studies have investigated on the correlation of whole, semi-skimmed, and skimmed milk with coronary artery disease (CAD) and myocardial infarction (MI) risk; However, no consensus has been reached and evidence on any causal links between these exposures and outcomes remains unclear. This study aimed to conduct univariate and multivariate Mendelian randomization (MR) analyses, using publicly released genome-wide association study summary statistics (GWAS) from the IEU GWAS database, to ascertain the causal association of milk with various fat content with CAD and MI risk.MethodsFor the exposure data, 29, 15, and 30 single-nucleotide polymorphisms for whole milk, semi-skimmed milk, and skimmed milk, respectively, obtained from 360,806 Europeans, were used as instrumental variables. CAD and MI comprised 141,217 and 395,795 samples, respectively. We used inverse variance weighted (IVW), weighted median, MR-Egger regression, and MR Pleiotropy Residual Sum and Outlier analyses to determine whether pleiotropy and heterogeneity could skew the MR results. Sensitivity tests were conducted to verify the robustness of the results.ResultsAfter adjusting for false discovery rates (FDR), we discovered proof that skimmed milk intake is a genetically predicted risk factor for CAD (odds ratio [OR] = 5.302; 95% confidence interval [CI] 2.261-12.432; P < 0.001; FDR-corrected P < 0.001) and MI (OR = 2.287; 95% CI 1.218-4.300; P = 0.010; FDR-corrected P = 0.009). Most sensitivity assessments yielded valid results. Multivariable MR for CAD and MI produced results consistent with those obtained using the IVW method. There was no causal relationship between whole or semi-skimmed milk, and CAD or MI.ConclusionOur findings indicate that the consumption of skimmed milk may increase the risk of CAD and MI. This evidence may help inform dietary recommendations for preventing cardiovascular disease. Further studies are required to elucidate the underlying mechanisms.

It is well documented that statins protect atherosclerotic patients from inflammatory changes and plaque instability in coronary arteries. However, the underlying mechanisms are not fully understood. Using a previously established mouse model for vulnerable atherosclerotic plaque, we investigated the effect of atorvastatin (10 mg/kg/day) on plaque morphology. Atorvastatin did not lower plasma total cholesterol levels or affect plaque progression at this dosage; however, vulnerable plaque numbers were significantly reduced in the atorvastatin-treated group compared to control. Detailed examinations revealed that atorvastatin significantly decreased macrophage infiltration and subendothelial lipid deposition, reduced intimal collagen content, and elevated collagenase activity and expression of matrix metalloproteinases (MMPs). Because vascular inflammation is largely driven by changes in monocyte/macrophage numbers in the vessel wall, we speculated that the anti-inflammatory effect of atorvastatin may partially result from decreased monocyte recruitment to the endothelium. Further experiments showed that atorvastatin downregulated expression of the chemokines monocyte chemoattractant protein (MCP)-1, chemokine (C-X3-C motif) ligand 1 (CX3CL1) and their receptors CCR2 and, CX3CR1, which are mainly responsible for monocyte recruitment. In addition, levels of the plasma inflammatory markers C-reactive protein (CRP) and tumor necrosis factor (TNF)-α were also significantly decrease in atorvastatin-treated mice. Collectively, our results demonstrate that atorvastatin can improve plaque stability in mice independent of plasma cholesterol levels. Given the profound inhibition of macrophage infiltration into atherosclerotic plaques, we propose that statins may partly exert protective effects by modulating levels of chemokines and their receptors. These findings elucidate yet another atheroprotective mechanism of statins.

Efficient membrane fusion has been successfully mimicked in vitro using artificial membranes and a number of cellular proteins that are currently known to participate in membrane fusion. However, these proteins are not sufficient to promote efficient fusion between biological membranes, indicating that critical fusogenic factors remain unidentified. We have recently identified a TIP30 protein complex containing TIP30, acyl-CoA synthetase long-chain family member 4 (ACSL4) and Endophilin B1 (Endo B1) that promotes the fusion of endocytic vesicles with Rab5a vesicles, which transport endosomal acidification enzymes vacuolar (H⁺)-ATPases (V-ATPases) to the early endosomes in vivo. Here, we demonstrate that the TIP30 protein complex facilitates the fusion of endocytic vesicles with Rab5a vesicles in vitro. Fusion of the two vesicles also depends on arachidonic acid, coenzyme A and the synthesis of arachidonyl-CoA by ACSL4. Moreover, the TIP30 complex is able to transfer arachidonyl groups onto phosphatidic acid (PA), producing a new lipid species that is capable of inducing close contact between membranes. Together, our data suggest that the TIP30 complex facilitates biological membrane fusion through modification of PA on membranes.

Activity-dependent transcription of brain-derived neurotrophic factor (BDNF) has been studied as an important model to elucidate the mechanisms underlying numerous aspects of neuroplasticity. It has been extensively emphasized that Ca(2+) influx through different routes may have significantly different effects on BDNF transcription. Here, we examined the regulatory property of the major calcium responsive elements (CaRE) in BDNF promoter IV in cultured rat cortical neurons. BDNF promoter IV, as well as CaRE1 and CaRE3, was significantly activated by Ca(2+) influx through L-type voltage-gated calcium channel (L-VGCC) or NMDA receptor (NMDAR). However, the L-VGCC- and NMDAR-mediated activation of CaRE was differentially regulated by different Ca(2+)-stimulated protein kinases. Specifically, PKA, CaMKI, and CaMKIV activity were required for L-VGCC-, but not NMDAR-mediated CaRE1 activation. CaMKI activity was required for NMDAR- but not L-VGCC-mediated CaRE3 activation. Surprisingly, the activation of CaRF, a previously identified transcription factor for CaRE1, was stimulated via L-VGCC but not NMDAR, and required MEK, PI3K and CaMKII activity. These results suggest a new working model that activity-dependent BDNF IV up-regulation may be coordinately mediated by CaRE1 and CaRE3 activity, which show different responses to Ca(2+)-stimulated kinases. Our data also explain how the individual cis-element in BDNF promoter is distinctively coupled to different Ca(2+) routes.

A sensitive assay to identify biomarkers using non-invasively collected clinical specimens is ideal for breast cancer detection. While there are other studies showing disease biomarkers in saliva for breast cancer, our study tests the hypothesis that there are breast cancer discriminatory biomarkers in saliva using de novo discovery and validation approaches. This is the first study of this kind and no other study has engaged a de novo biomarker discovery approach in saliva for breast cancer detection. In this study, a case-control discovery and independent preclinical validations were conducted to evaluate the performance and translational utilities of salivary transcriptomic and proteomic biomarkers for breast cancer detection.Salivary transcriptomes and proteomes of 10 breast cancer patients and 10 matched controls were profiled using Affymetrix HG-U133-Plus-2.0 Array and two-dimensional difference gel electrophoresis (2D-DIGE), respectively. Preclinical validations were performed to evaluate the discovered biomarkers in an independent sample cohort of 30 breast cancer patients and 63 controls using RT-qPCR (transcriptomic biomarkers) and quantitative protein immunoblot (proteomic biomarkers). Transcriptomic and proteomic profiling revealed significant variations in salivary molecular biomarkers between breast cancer patients and matched controls. Eight mRNA biomarkers and one protein biomarker, which were not affected by the confounding factors, were pre-validated, yielding an accuracy of 92% (83% sensitive, 97% specific) on the preclinical validation sample set.Our findings support that transcriptomic and proteomic signatures in saliva can serve as biomarkers for the non-invasive detection of breast cancer. The salivary biomarkers possess discriminatory power for the detection of breast cancer, with high specificity and sensitivity, which paves the way for prediction model validation study followed by pivotal clinical validation.

Figs and fig pollinators have co-evolved species-specific systems of mutualism. So far, it was unknown how visual opsin genes of pollinators have evolved in the light conditions inside their host figs. We cloned intact full-length mRNA sequences of four opsin genes from a species of fig pollinator, Ceratosolen solmsi, and tested for selective pressure and expressional plasticity of these genes. Molecular evolutionary analysis indicated that the four opsin genes evolved under different selective constraints. Subsets of codons in the two long wavelength sensitive opsin (LW1, LW2) genes were positively selected in ancestral fig pollinators. The ultraviolet sensitive opsin (UV) gene was under strong purifying selection, whereas a relaxation of selective constrains occurred on several amino acids in the blue opsin. RT-qPCR analysis suggested that female and male fig pollinators had different expression patterns possibly due to their distinct lifestyles and different responses to light within the syconia. Co-evolutionary history with figs might have influenced the evolution and expression plasticity of opsin genes in fig pollinators. [ABSTRACT FROM AUTHOR]

Oxidative stress has recently been considered as a pivotal player in the pathogenesis of diabetic gastrointestinal dysfunction. We therefore investigated the role of 2, 3, 5, 4'-tetrahydroxystilbene-2-O-beta-D-glucoside (THSG) that has a strong anti-oxidant property, in diabetic gastrointestinal dysmotility as well as the underlying protective mechanisms. THSG restored the delayed gastric emptying and the increased intestinal transit in streptozotocin (STZ)-induced diabetic mice. Loss of neuronal nitric oxide synthase (nNOS) expression and impaired nonadrenergic, noncholinergic (NANC) relaxations in diabetic mice were relieved by long-term preventive treatment with THSG. Meanwhile, THSG (10-7,10-4 mol/L) enhanced concentration-dependently NANC relaxations of isolated colons in diabetic mice. Diabetic mice displayed a significant increase in Malondialdehyde (MDA) level and decrease in the activity of glutathione peroxidase (GSH-Px), which were ameliorated by THSG. Inhibition of caspase-3 and activation of ERK phosphorylation related MAPK pathway were involved in prevention of enhanced apoptosis in diabetes afforded by THSG. Moreover, THSG prevented the significant decrease in PPAR-γ and SIRT1 expression in diabetic ileum. Our study indicates that THSG improves diabetic gastrointestinal disorders through activation of MAPK pathway and upregulation of PPAR-γ and SIRT1. [ABSTRACT FROM AUTHOR]

Figs and fig wasps form a peculiar closed community in which the Ficus tree provides a compact syconium (inflorescence) habitat for the lives of a complex assemblage of Chalcidoid insects. These diverse fig wasp species have intimate ecological relationships within the closed world of the fig syconia. Previous surveys of Wolbachia, maternally inherited endosymbiotic bacteria that infect vast numbers of arthropod hosts, showed that fig wasps have some of the highest known incidences of Wolbachia amongst all insects. We ask whether the evolutionary patterns of Wolbachia sequences in this closed syconium community are different from those in the outside world. In the present study, we sampled all 17 fig wasp species living on Ficus benjamina, covering 4 families, 6 subfamilies, and 8 genera of wasps. We made a thorough survey of Wolbachia infection patterns and studied evolutionary patterns in wsp (Wolbachia Surface Protein) sequences. We find evidence for high infection incidences, frequent recombination between Wolbachia strains, and considerable horizontal transfer, suggesting rapid evolution of Wolbachia sequences within the syconium community. Though the fig wasps have relatively limited contact with outside world, Wolbachia may be introduced to the syconium community via horizontal transmission by fig wasps species that have winged males and visit the syconia earlier. [ABSTRACT FROM AUTHOR]

Neurofilament, light polypeptide (NEFL) was demonstrated to be ectopically expressed in breast cancer tissues and decreased in lymph node metastases compared to the paired primary breast cancers in our previous study. Moreover, in several studies, NEFL was regarded as a tumor suppressor gene, and its loss of heterozygosity (LOH) was related to carcinogenesis and metastasis in several types of cancer. To explore the role of NEFL in the progression of breast cancer and to evaluate its clinical significance, we detected the NEFL mRNA level in normal breast tissues, primary breast cancer samples and lymph node metastases, and then analyzed the association between the NEFL expression level and several clinicopathological parameters and disease-free survival (DFS). NEFL mRNA was found to be expressed in 92.3% of breast malignancies and down-regulated in lymph node metastases compared to the paired primary tumors. NEFL mRNA level was lower in primary breast cancers with positive lymph nodes than in cancers with negative lymph nodes. Moreover, a low expression level of NEFL mRNA indicated a poor five-year DFS for early-stage breast cancer patients. Thus, NEFL mRNA is ectopically expressed in breast malignancies and could be a potential prognostic factor for early-stage breast cancer patients. [ABSTRACT FROM AUTHOR]

Among the Chalcidoids, hymenopteran parasitic wasps that have diversified lifestyles, a partial mitochondrial genome has been reported only from Nasonia. This genome had many unusual features, especially a dramatic reorganization and a high rate of evolution. Comparisons based on more mitochondrial genomic data from the same superfamily were required to reveal weather these unusual features are peculiar to Nasonia or not. In the present study, we sequenced the nearly complete mitochondrial genomes from the species Philotrypesis. pilosa and Philotrypesis sp., both of which were associated with Ficus hispida. The acquired data included all of the protein-coding genes, rRNAs, and most of the tRNAs, and in P. pilosa the control region. High levels of nucleotide divergence separated the two species. A comparison of all available hymenopteran mitochondrial genomes (including a submitted partial genome from Ceratosolen solmsi) revealed that the Chalcidoids had dramatic mitochondrial gene rearrangments, involved not only the tRNAs, but also several protein-coding genes. The AT-rich control region was translocated and inverted in Philotrypesis. The mitochondrial genomes also exhibited rapid rates of evolution involving elevated nonsynonymous mutations. [ABSTRACT FROM AUTHOR]

SMLS (Sitobion miscanthi L type symbiont) is a newly reported aphid secondary symbiont. Phylogenetic evidence from molecular markers indicates that SMLS belongs to the Rickettsiaceae and has a sibling relationship with Orientia tsutsugamushi. A comparative analysis of coxA nucleotide sequences further supports recognition of SMLS as a new genus in the Rickettsiaceae. In situ hybridization reveals that SMLS is housed in both sheath cells and secondary bacteriocytes and it is also detected in aphid hemolymph. The population dynamics of SMLS differ from those of Buchnera aphidicola and titer levels of SMLS increase in older aphids. A survey of 13 other aphids reveals that SMLS only occurs in wheat-associated species. [ABSTRACT FROM AUTHOR]

Cryptic and polymorphic species can complicate traditional taxonomic research and both of these concerns are common in fig wasp communities. Species identification is very difficult, despite great effort and the ecological importance of fig wasps. Herein, we try to identify all chalcidoid wasp species hosted by one species of fig, using both morphological and molecular methods. We compare the efficiency of four different DNA regions and find that ITS2 is highly effective for species identification, while mitochondrial COI and Cytb regions appear less reliable, possibly due to the interference signals from either nuclear copies of mtDNA, i.e. NUMTs, or the effects of Wolbachia infections. The analyses suggest that combining multiple markers is the best choice for inferring species identifications as any one marker may be unsuitable in a given case. [ABSTRACT FROM AUTHOR]

Hedgehog (Hh) genes play major roles in animal development and studies of their evolution, expression and function point to major differences among chordates. Here we focused on Hh genes in lampreys in order to characterize the evolution of Hh signalling at the emergence of vertebrates. Screening of a cosmid library of the river lamprey Lampetra fluviatilis and searching the preliminary genome assembly of the sea lamprey Petromyzon marinus indicate that lampreys have two Hh genes, named Hha and Hhb. Phylogenetic analyses suggest that Hha and Hhb are lamprey-specific paralogs closely related to Sonic/Indian Hh genes. Expression analysis indicates that Hha and Hhb are expressed in a Sonic Hh-like pattern. The two transcripts are expressed in largely overlapping but not identical domains in the lamprey embryonic brain, including a newly-described expression domain in the nasohypophyseal placode. Global alignments of genomic sequences and local alignment with known gnathostome regulatory motifs show that lamprey Hhs share conserved non-coding elements (CNE) with gnathostome Hhs albeit with sequences that have significantly diverged and dispersed. Functional assays using zebrafish embryos demonstrate gnathostome-like midline enhancer activity for CNEs contained in intron2. We conclude that lamprey Hh genes are gnathostome Shh-like in terms of expression and regulation. In addition, they show some lamprey-specific features, including duplication and structural (but not functional) changes in the intronic/regulatory sequences. [ABSTRACT FROM AUTHOR]

Lymphoproliferative diseases (LPDs) of the ocular adnexa encompass the majority of orbital diseases and include reactive follicular hyperplasia (RFH), atypical lymphoid hyperplasia (ALH), and mucosa-associated lymphoid tissue lymphoma (MALToma). Lymphoid follicles (LFs) are usually observed during the histological examination of LPDs. Currently, because there is a lack of specific clinical signs and diagnostic immunohistochemical biomarkers, it is difficult for pathologists to distinguish MALToma from ocular RFH and ALH, which makes the clinical management of these conditions difficult. Here, we analyzed the clinical features of patients with ocular adnexal LPDs (n = 125) and investigated the structure of LFs in paraffin-embedded tissue samples using anti-CD23 and anti-IgD immunochemistry. We found that some clinical features including age, sex, and laterality were different among RFH, LFH, and MALToma. Additionally, immunohistochemistry revealed that the expression of IgD and CD23 was higher in RFH patients and decreased in patients with ALH and MALToma. Moreover, LFs in RFH were intact, whereas the structures of most LFs were disrupted in ALH. In MALToma specimens, few intact LFs were observed. In a further investigation, we combined the results for CD23/IgD immunohistochemistry and the structure of LFs to establish a scoring method for the differential diagnosis of LPDs. According to the BIOMED-2 protocol, we further detected IgH gene monoclonal rearrangement in 73 cases (35 RFH, 17 ALH, and 21 MALToma cases). The sensitivity of our scoring method, based on a comparison with the results of IgH gene monoclonal rearrangement detection, was 85.7% (18/21) for MALToma and 35.3% (6/17) for ALH. Our study provides a method that may be useful for the differential diagnosis of RFH, ALH, and MALToma.

Basic/helix-loop-helix (bHLH) proteins comprise one of the largest transcription factor families and play important roles in diverse cellular and molecular processes. Comprehensive analyses of the composition and evolution of the bHLH family in cotton are essential to elucidate their functions and the molecular basis of cotton development. By searching bHLH homologous genes in sequenced diploid cotton genomes (Gossypium raimondii and G. arboreum), a set of cotton bHLH reference genes containing 289 paralogs were identified and named as GobHLH001-289. Based on their phylogenetic relationships, these cotton bHLH proteins were clustered into 27 subfamilies. Compared to those in Arabidopsis and cacao, cotton bHLH proteins generally increased in number, but unevenly in different subfamilies. To further uncover evolutionary changes of bHLH genes during tetraploidization of cotton, all genes of S5a and S5b subfamilies in upland cotton and its diploid progenitors were cloned and compared, and their transcript profiles were determined in upland cotton. A total of 10 genes of S5a and S5b subfamilies (doubled from A- and D-genome progenitors) maintained in tetraploid cottons. The major sequence changes in upland cotton included a 15-bp in-frame deletion in GhbHLH130D and a long terminal repeat retrotransposon inserted in GhbHLH062A, which eliminated GhbHLH062A expression in various tissues. The S5a and S5b bHLH genes of A and D genomes (except GobHLH062) showed similar transcription patterns in various tissues including roots, stems, leaves, petals, ovules, and fibers, while the A- and D-genome genes of GobHLH110 and GobHLH130 displayed clearly different transcript profiles during fiber development. In total, this study represented a genome-wide analysis of cotton bHLH family, and revealed significant changes in sequence and expression of these genes in tetraploid cottons, which paved the way for further functional analyses of bHLH genes in the cotton genus.

Brown cotton fiber is the major raw material for colored cotton industry. Previous studies have showed that the brown pigments in cotton fiber belong to proanthocyanidins (PAs). To clarify the details of PA biosynthesis pathway in brown cotton fiber, gene expression profiles in developing brown and white fibers were compared via digital gene expression profiling and qRT-PCR. Compared to white cotton fiber, all steps from phenylalanine to PA monomers (flavan-3-ols) were significantly up-regulated in brown fiber. Liquid chromatography mass spectrometry analyses showed that most of free flavan-3-ols in brown fiber were in 2, 3-trans form (gallocatechin and catechin), and the main units of polymeric PAs were trihydroxylated on B ring. Consistent with monomeric composition, the transcript levels of flavonoid 3', 5'-hydroxylase and leucoanthocyanidin reductase in cotton fiber were much higher than their competing enzymes acting on the same substrates (dihydroflavonol 4-reductase and anthocyanidin synthase, respectively). Taken together, our data revealed a detailed PA biosynthesis pathway wholly activated in brown cotton fiber, and demonstrated that flavonoid 3', 5'-hydroxylase and leucoanthocyanidin reductase represented the primary flow of PA biosynthesis in cotton fiber.

Bioactive gibberellins (GAs) comprise an important class of natural plant growth regulators and play essential roles in cotton fiber development. To date, the molecular base of GAs' functions in fiber development is largely unclear. To address this question, the endogenous bioactive GA levels in cotton developing fibers were elevated by specifically up-regulating GA 20-oxidase and suppressing GA 2-oxidase via transgenic methods. Higher GA levels in transgenic cotton fibers significantly increased micronaire values, 1000-fiber weight, cell wall thickness and cellulose contents of mature fibers. Quantitative RT-PCR and biochemical analysis revealed that the transcription of sucrose synthase gene GhSusA1 and sucrose synthase activities were significantly enhanced in GA overproducing transgenic fibers, compared to the wild-type cotton. In addition, exogenous application of bioactive GA could promote GhSusA1 expression in cultured fibers, as well as in cotton hypocotyls. Our results suggested that bioactive GAs promoted secondary cell wall deposition in cotton fibers by enhancing sucrose synthase expression.

Crop agronomic parameters (leaf area index (LAI), nitrogen (N) uptake, total chlorophyll (Chl) content ) are very important for the prediction of crop growth. The objective of this experiment was to investigate whether the wheat LAI, N uptake, and total Chl content could be accurately predicted using spectral indices collected at different stages of wheat growth. Firstly, the product of the optimized soil-adjusted vegetation index and wheat biomass dry weight (OSAVI×BDW) were used to estimate LAI, N uptake, and total Chl content; secondly, BDW was replaced by spectral indices to establish new spectral indices (OSAVI×OSAVI, OSAVI×SIPI, OSAVI×CIred edge, OSAVI×CIgreen mode and OSAVI×EVI2); finally, we used the new spectral indices for estimating LAI, N uptake, and total Chl content. The results showed that the new spectral indices could be used to accurately estimate LAI, N uptake, and total Chl content. The highest R(2) and the lowest RMSEs were 0.711 and 0.78 (OSAVI×EVI2), 0.785 and 3.98 g/m(2) (OSAVI×CIred edge) and 0.846 and 0.65 g/m(2) (OSAVI×CIred edge) for LAI, nitrogen uptake and total Chl content, respectively. The new spectral indices performed better than the OSAVI alone, and the problems of a lack of sensitivity at earlier growth stages and saturation at later growth stages, which are typically associated with the OSAVI, were improved. The overall results indicated that this new spectral indices provided the best approximation for the estimation of agronomic indices for all growth stages of wheat.