Shea tree (Vitellaria paradoxa) is an important fruit tree crop because of its oil used for cooking and industrial manufacture of cosmetics. Despite its many benefits, quantitative trait loci linked to the economic traits have not yet been studied. In this study, we performed association mapping on a panel of 374 shea tree accessions using 7,530 single-nucleotide polymorphisms (SNPs) markers for oil yield and seed related traits. Twenty three markers that were significantly (–log10 (p) = 4.87) associated to kernel oil content, kernel length; width and weight were identified. The kernel dry matter oil content and kernel width had the most significant Marker Trait Association (MTA) on chromosomes 1 and 8 respectively. Sixteen candidate genes that condition early induction of flower buds and somatic embryos, seed growth and development, substrate binding, transport, lipid biosynthesis, metabolic processes during seed germination and disease resistance and abiotic stress adaptation were identified. The presence of these genes suggest their role in promoting shea bioactive functions that condition high oil synthesis. This study provides insights into the important marker-linked seed traits with genes controlling them, useful for molecular breeding for improving oil yield in the species.

Machine learning, including deep learning, reinforcement learning, and generative artificial intelligence are revolutionising every area of our lives when data is made available. With the help of these methods, we can decipher information from larger datasets while addressing the complex nature of biological systems in a more efficient way. Although machine learning methods have been introduced to human genetic epidemiological research as early as 2004, those were never used to their full capacity. In this review, we outline some of the main applications of machine learning to assigning human genetic loci to health outcomes. We summarise widely used methods and discuss their advantages and challenges. We also identify several tools, such as Combi, GenNet and GMSTool, specifically designed to integrate these methods for hypothesis-free analysis of genetic variation data. We elaborate on the additional value and limitations of these tools from a geneticist’s perspective. Finally, we discuss the fast-moving field of foundation models and large multi-modal omics biobank initiatives.

Genome-wide association studies (GWAS) have successfully revealed many disease-associated genetic variants. For a case-control study, the adequate power of an association test can be achieved with a large sample size, although genotyping large samples is expensive. A cost‐effective strategy to boost power is to integrate external control samples with publicly available genotyped data. However, the naïve integration of external controls may inflate the type I error rates if ignoring the systematic differences (batch effect) between studies, such as the differences in sequencing platforms, genotype calling procedures, population stratification, and so forth. To account for the batch effect, we propose an approach by integrating External Controls into the Association Test by Regression Calibration (iECAT-RC) in case-control association studies. Extensive simulation studies show that iECAT-RC not only can control type I error rates but also can boost statistical power in all models. We also apply iECAT-RC to the UK Biobank data for M72 Fibroblastic disorders by considering genotype calling as the batch effect. Four SNPs associated with Fibroblastic disorders have been detected by iECAT-RC and the other two comparison methods. However, our method has a higher probability of identifying these significant SNPs in the scenario of an unbalanced case-control association study.

Sorghum (Sorghum bicolor L.) ranks fifth as the most crucial cereal crop globally, yet its seed morphology remains relatively unexplored. This study investigated seed morphology in sorghum based on 115 mini core and 130 Senegalese germplasms. Eight seed morphology traits encompassing size, shape, and color parameters were assessed. Statistical analyses explored potential associations between these traits and resistance to three major sorghum diseases: anthracnose, head smut, and downy mildew. Furthermore, genome-wide association studies (GWAS) were conducted using phenotypic data from over 24,000 seeds and over 290,000 publicly available single nucleotide polymorphisms (SNPs) through the Genome Association and Prediction Integrated Tool (GAPIT) R package. Significant SNPs associated with various seed morphology traits were identified and mapped onto the reference sorghum genome to identify novel candidate defense genes.

Pasmo is one of the most widespread diseases threatening flax production. To identify genetic regions associated with pasmo resistance (PR), a genome-wide association study was performed on 370 accessions from the flax core collection. Evaluation of pasmo severity was performed in the field from 2012 to 2016 in Morden, MB, Canada. Genotyping-by-sequencing has identified 258,873 single nucleotide polymorphisms (SNPs) distributed on all 15 flax chromosomes. Marker-trait associations were identified using ten different statistical models. A total of 692 unique quantitative trait nucleotides (QTNs) associated with 500 putative quantitative trait loci (QTL) were detected from six phenotypic PR datasets (five individual years and average across years). Different QTNs were identified with various statistical models and from individual PR datasets, indicative of the complementation between analytical methods and/or genotype × environment interactions of the QTL effects. The single-locus models tended to identify large-effect QTNs while the multi-loci models were able to detect QTNs with smaller effects. Among the putative QTL, 67 had large effects (3-23%), were stable across all datasets and explained 32-64% of the total variation for PR in the various datasets. Forty-five of these QTL spanned 85 resistance gene analogs including a large toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (TNL) type gene cluster on chromosome 8. The number of positive effect QTL (NPQTL) in accessions was significantly correlated with PR (R²=0.55), suggesting additive effects. NPQTL was also significantly associated with morphotype (R²=0.52) and major positive effect QTL were present in the fiber type accessions. The 67 large effect QTL are suited for marker-assisted selection and the 500 QTL for effective genomic prediction in PR molecular breeding.

To evaluate genetic variability and seedling plant response to a dominated Warrior (-) race of yellow rust in Northern and Central European germplasm, we used a population of 229 winter wheat cultivars and breeding lines for genome-wide association study (GWAS). A wide variation in yellow rust disease severity (based on infection types 1-9) was observed in this panel. Four breeding lines TS049 (from Austria), TS111, TS185, and TS229 (from Germany) and one cultivar TS158 (KWS Talent) from Germany were found to be resistant to Warrior (-) FS 53/20 and Warrior (-) G 23/19. GWAS identified five significant SNPs associated with yellow rust on chromosomes 1B, 2A, 5B, and 7A for Warrior (-) FS 53/20, while one SNP on chromosome 5B was associated with disease for Warrior (-) G 23/19. For Warrior (-) FS 53/20, we discovered a new QTL for yellow rust resistance associated with the marker Kukri_c5357_323 on chromosome 1B. The resistant allele G at the marker locus Kukri_c5357_323 on chromosome 1B and the susceptible allele T at the marker locus Excalibur_c17489_804 on chromosome 5B showed the largest effects (1.21 and 0.81, respectively) on the severity of yellow rust detected in Warrior (-) FS 53/20 and Warrior (-) G 23/19. Among 144 putative genes within the flanking sequence of the significant SNPs detected by GWAS, the function of the best candidate genes was determined as protein kinase activity and oxidoreductase activity. Our results provide the basis for knowledge-based resistance breeding in the face of the enormous impact of the Warrior (-) race on wheat production in Europe.

Oligopeptide transporters (OPTs) are a group of membrane localized proteins that have a broad range of substrate transport capabilities and contribute to numerous biological processes. However, limited information has been reported on OPTs in higher plants. In this study, a comprehensive analysis of the OPT gene family in Medicago truncatula was performed. A total of 26 OPT genes (MtOPT01-MtOPT26) have been identified in the Medicago truncatula genome. Phylogenetic analyses indicated that MtOPTs consisted of two distinct subgroups, 12 MtOPTs belonged to the peptide transport subgroup (PT-OPT) based on their predicted amino acid sequences containing the two highly conserved motifs (NPG and KIPPR) and 14 MtOPTs belonged to yellow stripe subgroup (YS-OPT). The MtOPTs distributed on each of 8 chromosomes in Medicago truncatula. Sequence analysis verified that MtOPTs significant similar to those in other plants. The copy number of MtOPTs was low and the multiply of MtOPTs was simple relatively. Gene structure analysis showed that most of the MtOPTs have various numbers of introns. The multiple of MtOPTs might play different biological roles which were supported by the fact that MtOPTs have a distinct tissue-specific expression pattern. The data obtained in this study will help to better understand the complexity of the MtOPTs gene family and provide new evidence for the function and evolution of the OPT gene family in higher plants.

Owing to the rising demand for vegetable soybean products, there is an increasing need for high-yield soybean varieties to raise the output. However, complex correlation patterns among quantitative traits with genetic interactions bring a challenge for vegetable soybean breeding. Here, a genome-wide association study (GWAS) was applied for six yield-related traits with 188 vegetable soybean accessions. A total of 81 SNPs were identified for pod length, plant height, pods number, pod thickness, pod width, and seed weight per plant by Blink model. Among these significant SNPs, 79 were novel SNPs, while 2 overlapped with previously reported SNPs. A total of 220 genes were found in 100 kb upstream and downstream regions of significant SNPs, of which 17 genes were functional proteins. Among these function proteins, four candidate genes, Glyma.13G109100, Glyma.03G183200, Glyma.09G102200, and Glyma.09G102300 were analyzed as significant haplotype variations, which encode MYB-related transcription factor, auxin-responsive family protein, F-box protein and CYP450, respectively. In addition, the relative expression with four candidate genes in the V030 and V071 (Among them, the plant height, pod number and fresh pod weight of V030 were lower than V071 strains) vegetable soybean were significantly different and these genes could be involved in plant growth and development via various pathways. Therefore, we proposed four genes from vegetable soybean germplasm that could be the key candidate genes for pod yield and plant height and these studies will accelerate the process of high-yield vegetable soybean breeding.

(1) Background: The goal of our study was to identify SNPs, metabolic pathways (KEGG), and gene ontology (GO) terms significantly associated with calving and workability; (2) Methods: Based on the EuroGenomics reference data set, we analyzed direct (DCE) and maternal (MCE) calving ease, direct (DSB), and maternal (MSB) stillbirth, milking speed (MSP), and temperament (TEM). We estimated SNP effects using a multi-SNP mixed-model. Further, SNP positions were mapped to genes, and GO terms/KEGG pathways of the corresponding genes were assigned. The estimation of GO term/KEGG pathway effects was based on a mixed-model using SNP effects as dependent variables; (3) Results: The number of significant SNPs comprised 59 for DCE, 25 for DSB and MSP, 17 for MCE and MSB, and 7 for TEM. Significant KEGG pathways were found for MSB (2), TEM (2), and MSP (1), while 11 GO terms were significant for MSP, 10 for DCE, 8 for DSB and TEM, 5 for MCE, and 3 for MSB.; (4) Conclusions: From the perspective of a better understanding of the genomic background of the phenotypes, traits with low heritabilities suggest that the focus should be moved from single genes to metabolic pathways or gene ontologies significant for the phenotype.

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-21, 2021-22, and 2022-23), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K SNP Infinium iSelect SNP assay to obtain genome-wide SNP marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 MTAs located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p≤0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-21 and 2021-22 environments. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interactions related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins, were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Tiller angle (TA) is one of the most important agronomic traits which directly affected plant architecture, photosynthetic efficiency, and planting density of rice. Although many quantitative trait loci (QTL) and associated genes for tiller angle have been identified, the detection of QTL for TA still required emphasis. In this study, 238 micro-core germplasm populations were used for genome-wide association analysis (GWAS) of the tiller angle. The frequency distribution of TA was shown to be continuous in 2018 and 2019, five QTL designated as qTA1, qTA2, qTA6, qTA9, qTA12, and two QTL designated as qTA1, qTA1-1 were detected in 2018 and 2019, respectively. These QTL explained phenotypic variation ranged from 5.94% to 12.38%. Among them, qTA2 and qTA6 are new QTL for TA, and only qTA1 was detected repeatedly in two years. 253 germplasm from the 3K Rice Genome Project were used for candidate gene and haplotypes analysis, the candidate gene and haplotypes analyzed showed LOC_Os01g67770, LOC_Os02g35180 and LOC_Os09g38130 were the candidate gene of qTA1, qTA2, and qTA9, respectively. The results of this study provide insight into the genetic improvement of plant architecture in rice.

The development of modern genomic tools has sped up the breeding program of food crops. More recently, genomic resources are being developed for tropical forages and will be key resources for developing climate resilient high yielding forage varieties. In this study, we present a genome-wide association study for biomass yield and feed quality traits in buffel grass (Cenchrus ciliaris L. aka Pennisetum ciliare L.,). Genome-wide markers generated using the DArTSeq platform and mapped onto the reference genome of Setaria italica were used for the genome-wide association study. The results revealed several markers associated with biomass yield and feed quality traits. A total of 78 marker trait associations were identified with R2 values ranging from 0.138 to 0.236. The marker trait associations were distributed across different chromosomes. Of these associations, the most marker trait associations (23) were observed on Chr9, followed by Chr5 with 12. The fewest number of marker trait associations were observed on Chr4 with 2. In terms of traits, 17 markers were associated with biomass yield, 24 with crude protein, 26 with TDN, 14 with ADF, 10 with NDF and six with DMI. Twenty of the identified markers were associated with at least two traits. The identified marker trait associations provide a useful genomic resource for future improvement and breeding of buffel grass.

In order to integrate genomics in breeding and development of drought tolerant groundnut genotypes, identification of genomic regions/genetic markers for drought surrogate traits is essential. We used SNP markers for a genetic analysis of the ICRISAT groundnut minicore collection for genome wide marker-trait association for some physiological traits and to determine the magnitude of linkage disequilibrium (LD) present in the genetic resources. The LD analysis showed that about 36% of loci pairs were in significant LD (P < 0.05 and r2 > 0.2) and 3.14% of the pairs were in complete LD. There was rapid decline in LD with distance and the LD was <0.2 at a distance of 41635 bp. The marker trait association (MTAs) studies revealed 20 significant MTAs (p <0.001) with 11 markers for leaf area index (4), canopy temperature (13), chlorophyll content (1) and NDVI (2). The markers explained 2 to 21% of the phenotypic variation observed. Most of the MTAs identified on the A subgenome were also identified on the respective homeologous chromosome on the B subgenome. The duplications of effect observed could be due to common ancestor of the A and B genome which explains the linkage detected between markers lying on different chromosomes seen in the current study. The present study identified a total of 20 highly significant marker trait associations with 11 markers for four physiological traits of importance in groundnut; LAI, CT, SCMR and NDVI. The markers identified in this study can serve as useful genomic resources to initiate marker-assisted selection and trait introgression of groundnut for drought tolerance. The identified markers in this study may be useful for marker assisted selection after further validation.

Ningxiang pigs exhibit a diverse array of fatty acids, making them an intriguing model for exploring the genetic underpinnings of Fatty acid metabolism. In this study, we conducted a genome-wide association study using a dataset comprising 50,697 SNPs and samples from over 600 Ningxiang pigs. Our investigation yielded novel candidate genes linked to five saturated fatty acids (SFAs), four monounsaturated fatty acids (MUFAs), and five polyunsaturated fatty acids (PUFAs). Re-markably, 37, 21, and 16 SNPs showed significant associations with SFAs, MUFAs, and PUFAs, respectively. Notably, ALGA0047587, H3GA0046208, DRGA0016063, and ALGA0031262 demonstrated substantial phenotypic variance (≥ 30%) in Arachidic acid (C20:0), Elaidic Acid (C18:1n-9(t)), and Arachidonic acid (C20:4n-6), respectively. Several significant SNPs were posi-tioned proximally to previously reported genes. In total, we identified 11, 6, and 5 candidate genes associated with SFAs, MUFAs, and PUFAs, respectively. These findings hold great promise for advancing breeding strategies aimed at optimizing meat quality and enhancing lipid metabolism within the intramuscular fat (IMF) of Ningxiang pigs.

At least 50% of factors predisposing to alcohol dependence (AD) are genetic and women affected with this disorder present with more psychiatric comorbidities, probably indicating different genetic factors involved. We aimed to run a genome-wide association study (GWAS) followed by a bioinformatic functional annotation of associated genomic regions in male and female patients with AD and eight related clinical measures. A genome-wide significant association of rs220677 with AD (p-value = 1.33×10^-8 calculated with the Yates-corrected Chi-square test under the assumption of dominant inheritance) was discovered in female patients. Associations of AD and related clinical measures with seven other single nucleotide polymorphisms listed in previous GWAS of psychiatric and addiction traits were differently replicated in male and female patients. The bioinformatic analysis showed that regulatory elements in the eight associated linkage disequilibrium blocks define the expression of 80 protein-coding genes. Nearly 68% of these and of 120 previously published coding genes associated with alcohol phenotypes directly interact in a single network. This study indicates that a number of genes behind the pathogenesis of AD are different in male and female patients, but implicated molecular mechanisms are functionally connected. The results also suggest the genetic basis of sex-specific psychiatric comorbidities of AD.

Control of Tropical Theileriosis, a tick-borne disease with a strong impact on cattle breeding, can be facilitated by using marker-assisted selection in breeding programs. Genome-wide association studies (GWAS) using high-density arrays are extremely important for the ongoing process of identifying genomic variants associated with resistance to Theileria annulata infection. In this work, Single Nucleotide Polymorphisms (SNPs) were analyzed in the Portuguese autochthonous cattle breeds Alentejana and Mertolenga. In total, 14 SNPs suggestive of significance (p ≤ 10-4) were identified for Alentejana cattle and 20 SNPs for Mertolenga cattle. The genomic regions around these SNPs were further investigated for annotated genes and Quantitative trait loci (QTLs) previously described by other authors. Regarding the Alentejana breed, the MAP3K1, CMTM7, SSFA2, and ATG13 genes are located near suggestive SNPs and appear as candidate genes for resistance to Tropical Theileriosis, considering its action the immune response and resistance to other diseases. On the other hand, in the Mertolenga breed, the UOX gene is also a candidate gene due to its apparent link to the pathogenesis of the disease. These results may represent a first step towards the possibility of including genetic markers for resistance to Tropical Theileriosis in current breed selection programs.

Meloidogyne javanica causing root-knot nematode in soybean is an important problem in soybean areas, leading to several yield losses. Some accessions have been identified carrying resistance loci to this nematode specie. In this study, a set of 317 soybean accessions were characterized for resistance to M. javanica. Genome-wide association study (GWAS) was performed using SNPs from genotyping-by-sequencing (GBS), and a region of 29.2 Kbp on chromosome 13 was identified. The haplotype analysis showed that SNPs were able to discriminate susceptible and resistant accessions, leading to 25 accessions sharing the resistance locus. Furthermore, 5 accessions may be new M. javanica resistance sources. The screening of the SNPs in the USDA soybean germplasm showed that several accessions previous reported as resistance to other nematodes also showed the resistance haplotype on chromosome 13. High levels of concordance among the phenotypes of Brazilian cultivars and the SNPs in chromosome 13 were observed. A in silico analysis of the mapped region on soybean genome revealed a presence of 5 genes with structural similarity with major resistance genes. The expression levels of the candidate genes in the interval demonstrated a potential pseudogene, and other two model genes up-regulated in the resistance source after pathogen infection. The SNPs associated to the region conferring resistance is a important tool for introgression of the resistance by marker-assisted selection in soybean breeding programs.

Several herbivores feed on maize in field and storage setups making the development of multiple-insect resistance a critical breeding target. In this study, an association mapping panel of 341 tropical maize lines was evaluated in three field environments for resistance to FAW (fall armyworm) whilst bulked grains were subjected to MW (maize weevil) bioassay, genotyped with Diversity Array Technologies single nucleotide polymorphisms (SNPs) markers. A multi-locus genome-wide association study (GWAS) revealed 62 quantitative trait nucleotides (QTNs) associated with FAW and MW resistance traits on all 10 maize chromosomes, of which, 47 and 31 were discovered at stringent Bonferroni genome-wide significance level of 0.05 and 0.01, respectively, and located within or close to multiple-insect resistance genomic regions (MIRGRs) concerning FAW, SB, and MW. Sixteen QTNs influenced multiple-traits of which six were associated with resistance to both FAW and MW suggesting a pleiotropic genetic control. Functional prioritization of candidate genes (CGs) located within 10-30kb of the QTNs revealed 64 putative GWAS-based CGs (GbCGs) showing evidence of involvement in plant defense mechanisms. Only one GbCG was associated with each of five of the six combined-resistance QTNs, thus, reinforcing the pleiotropy hypothesis. In addition, through In-silico co-functional network inferences, an additional 107 Network-based CGs (NbCGs), biologically connected to the 64 GbCGs, differentially expressed under biotic or abiotic stress were revealed within MIRGRs. The provided multiple-insect resistance physical map should contribute to the development of combined-insect resistance in maize.

A genome-wide association study (GWAS) was performed on a set of 260 lines which belong to three different bi-parental flax mapping populations. These lines were sequenced to an averaged genome coverage of 19× using the Illumina Hi-Seq platform. Phenotypic data for 11 seed yield and oil quality traits were collected in eight year/location environments. A total of 17,288 single nucleotide polymorphisms were identified, which explained more than 80% of the phenotypic variation for days to maturity (DTM), iodine value (IOD), palmitic (PAL), stearic, linoleic (LIO) and linolenic (LIN) acid contents. Twenty-three unique genomic regions associated with 33 QTL for the studied traits were detected, thereby validating four genomic regions previously identified. The 33 QTL explained 48-73% of the phenotypic variation for oil content, IOD, PAL, LIO and LIN but only 8-14% for plant height, DTM and seed yield. A genome-wide selective sweep scan for selection signatures detected 114 genomic regions that accounted for 7.82% of the flax pseudomolecule and overlapped with the 11 GWAS-detected genomic regions associated with 18 QTL for 11 traits. The results demonstrate the utility of GWAS combined with selection signatures for dissection of the genetic structure of traits and for pinpointing genomic regions for breeding improvement.

Fat has a high energy density and excessive fatness has been recognized as a problem for egg production and the welfare of chickens. The identification of a genetic polymorphism controlling fat deposition would be helpful to select against excessive fatness in the laying hen. This study aimed to estimate genomic heritability and identify the genetic architecture of abdominal fat deposition in a population of chickens from a Dongxiang blue-shelled local breed crossbred with the White Leghorn. A genome-wide association study wase conducted on abdominal fat percentage, egg production and body weights using a sample of 1534 hens genotyped with a 600K Chicken Genotyping Array. The analysis yielded a heritability estimate of 0.19 ± 0.04 for abdominal fat percentage, 0.56 ± 0.04 for body weight at 72 weeks, 0.11 ± 0.03 for egg production and 0.24 ± 0.04 for body weight gain. The genetic correlation of abdominal fat percentage with egg production between 60 and 72 weeks of age was 0.35 ± 0.18. This implies a potential trade-off between these two traits related to allocation of resource. Strong positive genetic correlations were found between fat deposition and weight traits. A promising locus close to COL12A1 on chromosome 3, associated with abdominal fat percent, was found in the present study. The other region located around HTR2A on chromosome 1, where allele substitution was predicted to be associated with body weight gain, accounted for 2.9% of phenotypic variance. Another region located on chromosome 1, but close to SOX5, was associated with egg production. These results may be used to influence the balanced genetic selection for laying hens.

Japanese apricot (Prunus mume) is a culturally and economically important fruit tree in East Asia. ‘Nanko’, the primary cultivar of the Japanese apricot well known in Japan, usually suffers from scab, a disease caused by Venturia carpophila. Scab results in the development of black spots on the fruit surface, which considerably reduces its commercial value. Despite being an economic concern for Japanese apricot, there have been few reports on the phenotypic variation in scab resistance/susceptibility, the underlying genetic factors, and the development of DNA markers. Therefore, in this study, we performed a comprehensive trait analysis for scab resistance and susceptibility for four years and conducted a genome-wide association study (GWAS) to identify the associated loci. The results showed that significant GWAS peaks were detected in 2017 and 2018. Promising candidate gene encoding domains associated with disease resistance were found at the nine single-nucleotide polymorphisms (SNPs) identified in 2017 and 2018. These SNPs were thought to be associated with scab susceptibility of ‘Nanko’ lineages. These findings shed light on the mechanisms of scab resistance in P. mume and will assist future breeding programs to improve scab resistance in ‘Nanko’ lineages.

There are rich proper noun derivative semantic relationships in geographic entity vector data. The identification of traditional proper noun derivative semantic relationships usually requires toponym researchers to examine the source of derived toponyms by consulting a large number of toponym documents. This method cannot adapt to the mining of proper noun derivative semantic relationships for large-scale geographical entity. In this regard, this paper proposes a method of identifying the proper noun derivative semantic relationships for geographical entity by combining prompt learning and XGBoost model. This method transforms the extraction problem of proper noun derivative semantic relationships into a binary classification problem by analyzing the characteristics of the derivation relationship of toponyms in the form of toponyms and the spatial and temporal distribution of geography. Firstly, prompt learning and gaussian mixture model are used to extract and cluster the derivation words of toponyms, so as to realize the recognition of derived toponyms collections. Then, enumerate the possible proper noun derivations in the toponym set, and then the XGBoost model is used to identify the semantic relations derived from the proper noun of geographical entities. In the experiment, the recognition method of proper noun derivative semantic relationships based on prompt learning and XGBoost proposed in this paper has achieved good relationship recognition effect in the experiment. This method has certain practical value in the fields of toponym management and toponym translation.

Lodging poses a significant challenge to rice yield, prompting the need to identify elite alleles for lodging resistance traits to improve cultivated rice varieties. In this study, a natural population of 518 rice accessions was examined to identify elite alleles associated with plant height (PH), stem diameter (SD), stem anti-thrust (AT/S), and various internode lengths: (first (FirINL), second (SecINL), third (ThirINL), fourth (ForINL), and fifth (FifINL) internode lengths). 262 SSR markers linked to these traits were uncovered through association mapping in two environmental conditions. Phenotypic evaluations revealed striking differences among cultivars, and genetic diversity assessments showed polymorphisms across the accessions. Favorable alleles were identified for PH, SD, AT/S, and 1 to 5 internode lengths, with specific alleles displaying considerable effects. Noteworthy alleles include RM6811-160 bp on chromosome 6 (which reduces PH) and RM161-145 bp on chromosome 5 (which increases SD). The study identified a total of 42 novel QTLs. Specifically, seven QTLs were identified for PH, four for SD, five for AT/S, five for FirINL, six for SecINL, five for ThirINL, six for ForINL, and four for FifINL. QTLs qAT/S-2, qPH2.1, qForINL2.1, and qFifINL exhibited the most significant phenotypic variance (PVE) of 3.99% for the stem lodging trait. AT/S, PH, ForINL and FifINL had additive effects of 5.31 kPa, 5.42 cm, 4.27cm and 4.27cm, respectively, offering insights into eight distinct cross-combinations for enhancing each trait. This research suggests the potential for cross-breeding superior parents based on stacked alleles, promising improved rice cultivars with enhanced lodging resistance to meet market demands.

The ZIPs (ZRT/IRT-like proteins) play critical roles in the absorption, transport, and intracellular balance of metal ions essential for various physiological processes in plants. However, little is known about the pan-genomic characteristics and properties of ZIP genes in cucumber (Cucumis sativus L.). In this study, we identified 10 CsZIP genes from the pan-genome of 13 C. sativus accessions. Among them, only CsZIP10 showed no variation in protein sequence length. We analyzed the gene structure, conserved domains, promoter cis-elements, and phylogenetic relationships of these 10 CsZIP genes derived from "9930". Based on phylogenetic analysis, the CsZIP genes were classified into three branches. Amino acid sequence comparison revealed the presence of conserved histidine residues in the ZIP proteins. Analysis of promoter cis-elements showed that most promoters contained elements responsive to plant hormones. Expression profiling in different tissues showed that most CsZIP genes were expressed at relatively low levels in C. sativus leave, stem, and tendril, and CsZIP7 and CsZIP10 were specifically expressed in root, indicating their potential involvement in the absorption and transport of metal ions. Transcriptomic data indicated that these 10 ZIP genes displayed responses to both downy mildew and powdery mildew, and CsZIP1 were significantly downregulated after both salt and heat treatments. In conclusion, this study deepens our understanding of the ZIP gene family and enhances our knowledge of the biological functions of CsZIP genes in C. sativus.

For the different harvest targets, the requirement for the prolificacy trait of maize was also different, so prolificacy is of great significance for modern production. Although some QTLs and genes associated with prolificacy in teosinte have been reported, the genetic mechanism of prolificacy in maize has not been fully elucidated. In this study, two RIL populations and GWAS population were used to genetic research of prolificacy trait in maize, with multi-environment. Combine linkage analysis and Genome-wide association study has identified a total of 13 QTLs and 8 significant SNPs. There were two genes related to tissue differentiation in the stable QTL qP9-2, and two significant SNPs corresponding to three genes were in QTL qP5-1 and QTL qP7-1, respectively. Four candidate genes GRMZM2G317262, GRMZM2G317584, GRMZM5G882364 and GRMZM2G141679 were finally screened out by qRT-PCR analysis. Based on the function of candidate genes, ethylene signaling pathway plays an important role in the formation of prolificacy in maize. It has deepened our understanding of the formation mechanism of prolificacy and laid a foundation for breeding new varieties with various prolificacy in maize.

New flaxseed cultivars differing in seed mucilage content (MC) with low hull content (HC) represent an attractive option to simultaneously target the food and feed markets. Here, a genome-wide association study (GWAS) was conducted for MC and HC in 200 diverse flaxseed accessions genotyped with 1.7 million SNP markers. The data obtained for MC and HC indicated a broad phenotypic variation and high (~70%) and a moderate (~49%) narrow sense heritability, respectively. MC and HC did not differ statistically between fiber and oil morphotypes, but yellow-seeded accessions had 2.7% less HC than brown-seeded ones. The genome wide linkage disequilibrium (LD) decayed to r2 = 0.1 at a physical distance of ~100 Kb. Seven and four QTL were identified for MC and HC, respectively. Promising candidate genes included Linum usitatissimum orthologs of the Arabidopsis thaliana genes TRANSPARENT TESTA 8, SUBTILISIN-LIKE SERINE PROTEASE, GALACTUROSYL TRANSFERASE-LIKE 5, MUCILAGE-MODIFIED 4, AGAMOUS-LIKE MADS-BOX PROTEIN AGL62, GLYCOSYL HYDROLASE FAMILY 17 and UDP-GLUCOSE FLAVONOL 3-O-GLUCOSYLTRANSFERASE that have been shown to play a role in mucilage synthesis and release, seed coat development and anthocyanin biosynthesis in A. thaliana were identified. The favorable alleles will be useful in flaxseed breeding towards the goal of achieving the ideal flaxseed cultivars for food and feed by genomic-based breeding.

Diacylglycerol kinase (DGK) is a lipid kinase that phosphorylates diacylglycerol (DAG) to generate phosphatidic acid (PA). Based on converting one important signaling molecule (DAG) to another (PA), DGK plays an important role in plant responses to abiotic stress including waterlogging stress. However, no studies have been reported on the characterization of the DGK gene family in the waterlogging-tolerant kiwifruit germplasm, Actinidia valvata Dunn. Here, we identified 18 AvDGK genes in the A. valvata genome. The phylogenetic analysis showed that AvDGKs can be classified into three clusters, and members within the same cluster have similar domain distribution, exon-intron structures and conserved motif compositions. Chromosome localization analysis revealed that all the AvDGK genes are located across 18 different chromosomes. There were 29 duplicated gene pairs in kiwifruit and all had undergone purifying selection during evolution. Promoter cis-element analysis revealed that the cis-elements within AvDGK genes are associated with multiple functions, including phytohormone signal transduction, stress response, and plant growth and development. The expression pattern analyses indicated that AvDGK play an important role in the fruit development and plant response to waterlogging stress. AvDGKs gene family in the tetraploid A. valvata genome might promote PA synthesis and subsequent signal transduction both under short-term and long-term waterlogging stress. These results provide information regarding the structural characteristics and potential function of DGK genes within kiwifruit and lay a fundamental basis for further research into breeding for enhancing the kiwifruit’s tolerance to waterlogging stress.

Mercury (Hg) pollution not only poses a threat to the environment but also adversely affects the growth and development of plants, with potential repercussions for animals and humans through bioaccumulation in the food chain. Maize, a crucial source of food, industrial materials, and livestock feed, requires special attention in understanding the genetic factors influencing mercury accumulation. Developing maize varieties with low mercury accumulation is vital for both maize production and human health. In this study, a comprehensive genome-wide association study (GWAS) was conducted using an enlarged SNP panel comprising 1.25 million single nucleotide polymorphisms (SNPs) in 230 maize inbred lines across three environments. The analysis identified 111 significant SNPs within 78 quantitative trait loci (QTL), involving 169 candidate genes under the Q model. Compared to the previous study, the increased marker density and optimized statistical model led to the discovery of 74 additional QTL, demonstrating improved statistical power. Gene Ontology (GO) enrichment analysis revealed that most genes participate in arsenate reduction and stress responses. Notably, GRMZM2G440968, associated with the significant SNP chr6.S_155668107 of axis tissue, encodes a cysteine proteinase inhibitor, suggesting its potential role in alleviating Hg toxicity by inhibiting cysteine. Haplotype analyses provided further insights, indicating that lines carrying hap3 exhibited the lowest mercury content compared to other haplotypes. In summary, our study significantly enhances the statistical power of GWAS, identifying additional genes related to mercury accumulation and metabolism. These findings offer valuable insights into unraveling the genetic basis of mercury content in maize and contribute to the development of maize varieties with low mercury accumulation.

Eating disorders (ED) are characterized by alterations in eating behavior. The genetic factors shared between ED diagnoses have been underexplored. The present study aimed to perform a genome-wide association study on individuals with disordered eating behaviors in the Mexican population, blood methylation quantitative trait loci (blood-meQTL) analysis, and in silico function prediction by different algorithms. The analysis included a total of 1803 individuals. Genome-wide association study and blood-meQTL analysis were performed by logistic and linear regression. In silico functional variant prediction, phenome-wide, and transcriptome-wide association studies by different algorithms were analyzed. In the genome-wide association study, we identified 44 single-nucleotide polymorphisms (SNP) associated at a nominal value and 7 blood-meQTL at a genome-wide umbral. The SNPs were enriched in genome-wide associations of the metabolic and immunologic domains. In the in silico analysis, the SNP rs10419198 located on an enhancer mark could change the expression of PRR12 on blood, adipocytes, and brain areas that regulate food intake. The present study supports the previous associations of genetic variation in the metabolic domain with ED.

Gibberella ear rot (GER) caused by Fusarium graminearum (teleomorph Gibberella zeae) is one of the most destructive diseases in maize that severely reduces grain yield and contaminates several potential mycotoxins. However, few efforts had been devoted to dissect the genetic basis of maize GER resistance. In the present study, a genome-wide association study (GWAS) was conducted in a maize association panel consisting of 303 diverse inbred lines. The phenotypes of GER severity were evaluated using kernel bioassay across multiple time points in the laboratory. Then, three models including fixed and random model circulating probability unification model (FarmCPU), general linear model (GLM) and mixed linear model (MLM), were conducted simultaneously in GWAS to identify single-nucleotide polymorphisms (SNPs) significantly associated with GER resistance. A total of four individual significant association SNPs with the phenotypic variation explained (PVE) ranging from 3.51 to 6.42% were obtained. Interestingly, the peak SNP (PUT-163a-71443302-3341) with the greatest PVE value, was co-localized in all models. Subsequently, 12 putative genes were captured from the peak SNP that several of these genes were directly or indirectly involved in disease resistance. Overall, these findings contribute to understand the complex plant-pathogen interactions in maize GER resistance. The regions and genes identified herein provide a list of candidate targets for further investigation, in addition to the kernel bioassay that can be used for evaluating and selecting elite germplasm resources with GER resistance in maize.

Wheat dwarf virus (WDV) causes an important vector transmitted virus disease, which leads to significant yield losses in barley production. Due to the fact, that at the moment no plant protection products are approved to combat the vector Psammotettix alienus and this disease cannot be controlled by chemical means, therefore the use of WDV resistant or tolerant genotypes is the most efficient method to control and reduce negative effects of WDV on barley growth and production. In this study, a set of 480 barley genotypes were screened to identify genotypic differences in response to WDV and five traits were assessed under infected and non-infected conditions. In total, 32 genotypes showed resistance or tolerance against WDV. Subsequently, the phenotypic data of 191 out of 480 genotypes combined with 34,408 single nucleotide polymorphisms (SNPs) in a genome wide association study to identify quantitative trait loci (QTLs) and markers linked to resistance/tolerance to WDV, using Tassel, GAPIT, and FARM CPU. In total, 1, 3, 2, 2 and 1 significantly associated markers on chromosomes 3H, 4H, 5H and 7H identified by all three methods for ELISA-60, relative performance of total grain weight, plant height, number of ears per plant and thousand grain weight, respectively.

In this work we looked for a stable and accurate procedure to perform feature selection in datasets with a much higher number of predictors than individuals, as in Genome-Wide Association Studies. Due to the instability in feature selection when many potential predictors are measured, a variable selection procedure is proposed that combines several replications of shrinkage regression models. A weighted formulation is used to define the final predictors. The procedure is applied to investigate Single Nucleotide Polymorphisms (SNPs) predictors associated to Alzheimer’s disease in the Alzheimer’s disease Neuroimaging Initiative (ADNI) dataset. Two data scenarios are investigated: one considering only the set of SNPs and another with the covariates age, sex, educational level and APOE4 genotype. The SNP rs2075650 and the APOE4 genotype are given as risk factors for AD, which is in line with the literature, and other four new SNPs are pointed, opening new hypotheses for in vivo analyses. These experiments demonstrate the potential of the new method for stable feature selection.

Cercospora beticola Sacc. is an ascomycete pathogen that causes Cercospora leaf spot in sugar beets (Beta vulgaris L.) and other related crops. The disease symptoms are characterized by circular lesions with an ash-gray centre and dark brown to purple margins on the leaves. It can lead to significant yield losses if not effectively managed. This study aims to characterize and assess the impact of bacterial biological control agents against C. beticola. Following a dual-culture screening, 18 bacteria exhibiting over 50% inhibition were selected, with six of them demonstrating more than 80% control. Identification of the bacteria was carried out through sequencing of the16S rDNA, revealing potentially 12 species belonging to six genera including Bacillus which was represented by four species. Additionally, the bacteria underwent in-depth characterization of biochemical and molecular properties, as well as plant-growth promotion. PCR detection of genes responsible for producing antifungal metabolites revealed that 83%, 78%, 89%, and 56% of the selected bacteria possessed the bacillomycin, iturin, fengycin, and surfactin encoding genes, respectively. Infrared spectroscopy analysis confirmed the presence of lipopeptide structure in the bacteria supernatant filtrate. Subsequently, the bacteria were assessed for their effect on sugar beet plants in controlled conditions. The bacteria exhibited notable capabilities, promoting growth in both roots and shoots, resulting in significant increases in root length, weight, and shoot length. The field experiment with four bacterial candidates demonstrated good performance against C. beticola compared to difenoconazole fungicide. These bacteria played a significant role in disease control, achieving a maximum efficacy of 77.42%, slightly below the 88.51% efficacy attained by difenoconazole. Additional field trials are necessary to verify the protective and growth-promoting effects of these candidates, whether applied individually, combined in consortia, or integrated with chemical inputs in sugar beet crop production.

Research Highlights: This study identified the cell cycle genes in birch that likely play important roles during plant growth and development. This analysis provides a basis for understanding the regulatory mechanism of various cell cycles in Betula pendula. Background and Objectives: The cell cycle factors not only influence cell cycle progression together, but also regulate accretion, division and differentiation of cells, and then regulate growth and development of plant. In this study, we identified the putative cell cycle genes in B. pendula genome, based on the annotated cell cycle genes in A. thaliana. It could serve as a foundation for further functional studies. Materials and Methods: The transcript abundance was determined for all the cell cycle genes in xylem, root, leaf and flower tissues using RNA-seq technology. Results: We identified 59 cell cycle gene models in the genome of B. pendula, 17 highly expression genes among them. These genes were BpCDKA.1, BpCDKB1.1, BpCDKB2.1, BpCKS1.2, BpCYCB1.1, BpCYCB1.2, BpCYCB2.1, BpCYCD3.1, BpCYCD3.5, BpDEL1, BpDpa2, BpE2Fa, BpE2Fb, BpKRP1, BpKRP2, BpRb1 and BpWEE1. Conclusions: We identified 17 core cell cycle genes in the genome of birch by combining phylogenetic analysis and tissue specific expression data.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), a novel evolutionarily divergent RNA virus etiological agent of COVID-19, is responsible for present devastating pandemic respiratory illness. To explore the genomic signatures, we comprehensively analyzed 2,492 complete and/or near-complete genome sequences of SARS-CoV-2 strains reported from across the globe to the GISAID database up to 30 March 2020. Genome-wide annotations revealed 1,407 nucleotide-level mutations at different positions throughout the entire genome of SARS-CoV-2. Moreover, nucleotide deletion analysis found nine deletions throughout the genome, including in polyprotein (n=6), ORF10 (n=1) and 3´-UTR (n=2). Evidence from the systematic gene-level mutational and protein profile analyses revealed a large number of amino acid (aa) substitutions (n=722), making the viral proteins heterogeneous. Notably, residues of receptor-binding domain (RBD) having crucial interactions with angiotensin-converting enzyme 2 (ACE2), and cross-reacting neutralizing antibody were found to be conserved among the analyzed SARS-CoV-2 strains, except for replacement of Lysine with Arginine at 378 position of the cryptic epitope of a Shanghai isolate, hCoV-19/Shanghai/SH0007/2020 (EPI_ISL_416320). Our method of genome annotation is a promising tool for monitoring and tracking the epidemic, the associated genetic variants, and their implications for the development of effective control and prophylaxis strategy.

Increasingly, genomics is being used for the prediction of specific traits and diseases (phenotypes) among humans. Wider availability of genomics data through multiple research projects (such as International HapMap Project1 and 1000 Genomes2) has been a catalyst in that direction. With the recent advances in machine learning and big data analysis, data computation resources and data models needed for genomics data analysis are readily available. However, the prediction of traits and diseases has its own challenges in terms of computational requirements and computational analysis, statistical analysis (example: confounding variables), and limited quality of data collection. Linear Mixed Models (LMM, a type of linear regression) is a common approach for Genome-wide Association Studies (GWAS) for the prediction of common traits among humans using genomics. This paper researches the existing LMM-based approaches for Genome-wide Association Studies (GWAS), describes the experiment performed on FaST-LMM approach from Microsoft Research, and then proposes an enhanced approach (called LMM-22) on how to address computational and statistical issues. LMM-22 focuses on the parallelization of LMM computations and execution of LMM-22 on General Purpose Graphics Processing Units (GPU) as against CPUs to accelerate the LMM approach for GWAS studies.

Zanthoxylum armatum is an forest tree with apomixis characteristic, and RWP-RKs are plant-specific transcription factors known to regulate nitrogen response and embryonic development. However, the genome-wide analysis of RWP-RK family genes in Z.armatum is unclear. In this study, 36 ZaRWP-RK transcription factors were identified in the genome of Z.armatum, among which 15 genes were belonged to the RKD subfamily and 21 belonged to the NLP subfamily. Duplication events of ZaRWP-RK genes were mainly segmental duplication, and synteny analysis revealed close phylogenetic relationship between Z.armatum and Arabidopsis. The analysis of cis-elements indicated that ZaRWP-RK genes may be involved in the regulation of embryonic development of Z.armatum by responding to plant hormones such as abscisic acid, auxin, and gibberellin. Real-time PCR showed that the expression levels of most ZaRWP-RK genes were significantly increased from flowers to young fruits. Protein-protein interaction network analysis further revealed the potential roles of the ZaRWP-RK proteins in apomixis. Collectively, this study provide a theoretical basis for future investigations into the ZaRWP-RK genes and their regulatory mechanisms in the apomixis process of Z.armatum.

Codon use bias is an important characteristic in the process of genetic information transmission of species. With the publication of the genome of C. panzhihuaensis, it is of great significance to investigate the codon use bias for understanding the genetic evolution of this species and for molecular breeding. In this study, Perl language and CodonW 1.4.2 software was used to systematically analyze the coding gene sequences of the C. panzhihuaensis genome, obtain the characteristics of codon use, and identify the sources of variation affecting codon use. The results showed that GC3s content, GCall content, and ENC value were 46.75%, 47.67%, and 52.8 respectively. As the first seed plant, the GC3 content of Cycas is similar to that of most gymnosperms and most dicotyledons, but considerably different from that of monocotyledons.RSCU value was greater than 1, among which 21 codons ended in U or A. Enc-plot analysis, PR2-plot neutral plot analysis, and correlation coefficient analysis of Axis1 with GC3s and CAI revealed that the codon bias was mainly affected by natural selection, but also by mutation pressure and other factors. Thirty-one optimal codons of Cycas genes were screened. This research can provide a reference for the phylogeny and genome codon evolution of Cycas.

Plant protein phosphatase 2Cs (PP2Cs) serve as negative regulators of protein kinase cascades activated in different processes and play important roles in plant development and abiotic-stress-mediated signaling pathways. In this study, a genome-wide study was conducted on the CqPP2C gene family. A total of putative 117 CqPP2C genes were identified. Comprehensive analyses of physical and chemical feature, chromosome localization and subcellular localization were conducted. According to phylogenetic analysis, CqPP2Cs were divided into 13 subfamilies. CqPP2Cs in the same subfamily had similar gene structure, conserved motif and all the CqPP2C proteins had the type 2C phosphatase domains. Gene duplication revealed that segmental duplication was the major driving force for CqPP2Cs expansion and all duplicated CqPP2Cs evolved from purifying selection. The expression of CqPP2Cs in various tissues under different abiotic stresses was analyzed using RNA-seq data. The results revealed that CqPP2C genes were involved in regulating development and stress responses of quinoa. Real-time quantitative reverse transcription PCR (qRT-PCR) analysis of six CqPP2C genes in subfamily A revealed that they were up-regulated or down-regulated under salt and drought treatments. Furthermore, the results of yeast two-hybrid assays revealed that subfamily A CqPP2Cs interacted not only with subclass III CqSnRK2s, but also with subclass II CqSnRK2s. Subfamily A CqPP2Cs could interact with CqSnRK2s in different combinations and intensities to respond to various biological processes and stresses. Overall, our results will be useful for understanding the functions of CqPP2C in regulating ABA signals and responding to abiotic stress.

The basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors widely distributed in eukaryotic organisms. In plants, they are not only involved in growth and development, defense and stress responses, and regulation of physiological processes, but also play a pivotal role in regulating secondary metabolism. To explore the function related to the bZIP gene family in Stevia rebaudiana Bertoni, we identified 105 SrbZIP genes at the genome-wide level and classified into 12 subfamilies using bioinformation methods. Analysis of physicochemical properties showed that the amino acid length readings ranged from 96 aa (SrbZIP13) to 704 aa (SrbZIP17), and the isoelectric point ranged from 4.64 (SrbZIP21) to 10.08 (SrbZIP101). The majority of these proteins (94.29%) were localized in the nucleus. Three main classes of cis-acting elements found in the SrbZIP promoter regions, including development-related elements (CAT-box), defense and stress-responsive elements (LTR/ARE/MBS), and phytohormone respon- sive elements (ABRE/TGACG-mitif/TCA-element). Through protein-protein interac- tion network of 105 SrbZIP proteins, SrbZIP proteins were mainly classified into four major categories: ABF2/ABF4/ABI5 (SrbZIP51/SrbZIP38/SrbZIP7), involved in phy- tohormone signaling, GBF1/GBF3/GBF4 (SrbZIP29/SrbZIP63/SrbZIP60) involved in environmental signaling, AREB3 (SrbZIP88), PAN (SrbZIP12), TGA1 (SrbZIP69), TGA4 (SrbZIP82), TGA7 (SrbZIP31), TGA9 (SrbZIP95), TGA10 (SrbZIP79) and HY5 (SrbZIP96) involved in cryptochrome signaling, and FD (SrbZIP23, SrbZIP24, SrbZIP47, SrbZIP67, SrbZIP68, SrbZIP71, SrbZIP72) promoted flowering. The transcriptomic data showed that SrbZIP genes differentially expressed in six S. rebaudiana cultivars (‘023’, ‘110’, ‘1188’, ‘11-14’, ‘GP’ and ‘GX’). Moreover, the expression levels of selected 15 SrbZIP genes in response to light, abiotic stress (low temperature, salt and drought), phytohormones (methyl jasmonate, gibberellic acid and salicylic acid) treatment and in different tissues were analyzed utilizing qRT-PCR. Six SrbZIP (SrbZIP54, SrbZIP63, SrbZIP32, SrbZIP45, SrbZIP60 and SrbZIP9) genes were further identified to be highly induced by factors affecting glycoside synthesis. Among them, three SrbZIP genes (SrbZIP54, SrbZIP63 and SrbZIP32) were predicted to be related to stress-responsive terpenoid synthesis in S. rebaudiana. The protein-protein interaction network expanded the potential functions of SrbZIP genes. This study firstly provided the comprehensive genome-wide report of the SrbZIP gene family, laying a foundation for further research on evolution, function and regulatory role of bZIP gene family in terpenoid synthesis in S. rebaudiana.

Entamoeba histolytica, the causative agent of amebiasis, is the third leading cause of death among parasitic diseases globally. Its life cycle includes encystation, which has been mostly studied in Entamoeba invadens, responsible of reptilian amebiasis. However, the molecular mechanisms underlying this process are not fully understood. Therefore, we focused on the identification and characterization of Myb proteins, that regulate the expression of encystation-related genes in various protozoan parasites. Through bioinformatic analysis, we identified 47 genes in E. invadens encoding MYB-domain-containing proteins. These were classified into single-repeat 1R (19), 2R-MYB proteins (27), and one 4R-MYB protein. The in silico analysis suggests that these proteins are multifunctional, participating in transcriptional regulation, chromatin remodeling, telomere maintenance, and splicing. Transcriptomic data analysis revealed expression signatures of eimyb genes, suggesting a potential orchestration in the regulation of early and late encystation-excystation genes. Furthermore, we identified 3135 probable target genes associated with reproduction, meiotic cell cycle, ubiquitin-dependent protein catabolism, and endosomal transport. In conclusion, our findings suggest that E. invadens Myb proteins regulate stage-specific proteins and a wide array of cellular processes. This study provides a foundation for further exploration of the molecular mechanisms governing encystation and unveils potential targets for therapeutic intervention in amebiasis.

Transcription factors (TFs) in the homeodomain-leucine zipper (HD-ZIP) family serve as essential regulators of plant development and control responses to environmental stimuli. To date, however, the HD-ZIP gene family in Prunus nana remains to be fully characterized. Accordingly, a genome-wide analysis of P. nana HD-ZIP family genes was performed using the most recent genomic data available, leading to the classification of 30 HD-ZIP TFs. These genes were annotated and subjected to systematic analyses of phylogenetic relationships, protein physicochemical properties, sequence-based structural characteristics, chromosomal distributions, and associated cis-acting regulatory elements. High levels of diversity were observed with respect to the gene structures for these PnaHD-Zip genes and the cis-regulatory elements found in the promoter regions upstream of these genes, suggesting that they play diverse roles in a range of biological contexts. These 30 PnaHD-Zip genes were further classified into four subgroups based on the results of phylogenetic, gene structure, and conserved motif analyses. Subsequent qPCR analyses indicated that PnaHD-Zip1 and PnaHD-Zip7 expression levels tended to increase was continuously inhibited in response to cold stress, suggesting that proteins in this HD-ZIP gene family may exhibit distinct responses to low-temperature stress exposure. Overall, these results offer a robust foundation for future studies seeking to explore the functional roles that HD-ZIP TFs play as regulators of cold stress tolerance in P. nana, in addition to offering more general insight regarding the regulatory functions and characteristics of these different HD-ZIP genes in P. nana.

A recent study claimed that genome-wide polygenic scores (GPSs) for five common diseases could identify individuals with risk equivalent to monogenic mutations. Receiver operator curve analyses were reported to have areas under the curve (AUCs) ranging from 0.63 for inflammatory bowel disease up to 0.81 for coronary artery disease (CAD) but these models also included age and sex, themselves strong predictors of risk. The GPS for CAD identified 8% of the population at threefold increased risk, which it was claimed was comparable to the excess risk from monogenic mutations. In the present study attempts were made to model the distribution of the GPS for CAD to match the information provided. These models were based on the reported distribution of prevalence by centile of GPS and on the distribution of GPS in controls and cases and were fitted to the reported results using linear approximations to the distributions and using simulations of a liability-threshold model. It was impossible to produce a compatible model in which the GPS produced an AUC as high as 0.81 and the most plausible estimate was that the true AUC was only 0.65. The reported distributions of the GPS in cases and controls overlap so much that they are not compatible with an AUC of 0.7 or higher. The AUC of the GPS for these diseases is modest. Furthermore, the literature robustly demonstrates that true CAD risk associated with monogenic mutations is much higher than the threefold increase which is predicted by the GPS. Together, these findings cast doubt on the clinical utility of the GPS.

Chromosome analysis (CA) and chromosomal microarray analysis (CMA) have been successfully used to diagnose genetic disorders. However, many conditions remain undiagnosed due to limitations in resolution (CA) and detection of only unbalanced events (CMA). Optical genome mapping (OGM) has the potential to address these limitations by capturing both structural variants (SVs) resulting in copy number changes and balanced rearrangements with high resolution. In this study, we investigated OGM’s concordance using 87 SVs previously identified by CA, CMA, or Southern blot. Overall, OGM was 98% concordant with only three discordant cases: (1) uncalled translocation with one breakpoint in a centromere; (2) uncalled duplication with breakpoints in the pseudoautosomal region 1; and (3) uncalled mosaic triplication originating from a marker chromosome. OGM provided diagnosis for three previously unsolved cases: (1) disruption of the SON gene due to a balanced reciprocal translocation; (2) disruption of NBEA gene due to an inverted insertion; (3) disruption of the TSC2 gene due to a mosaic deletion. We show that OGM is a valid method for detection of many types of SVs in a single assay and is highly concordant with legacy cytogenomic methods; however, it has limited SV detection capabilities in centromeric and pseudoautosomal regions.

The objective of this study were to conduct association mapping for drought tolerance at the seedling stage and yield-related traits. 60 cowpea accessions were used in the study. Single-nucleotide polymorphisms (SNPs) discovered through genotyping by sequencing (GBS) were used for genotyping. Association mapping was conducted using single-marker regression (SMR) in Q Gene, and general linear model (GLM) and mixed linear model (MLM) built in TASSEL. The population of the cowpea accessions were analysed using STRUCTURE 2.3.4 and the peak of delta K in the greenhouse showed seven population types, whereas the peak of delta K in the glasshouse indicated the presence of six population types. One SNP marker, 14083649|F|0-9 was associated with NP with a p value <0.001. Fifty SNP markers were associated with PWT at p <0.001. Four SNP markers, 14074781|F|0-16, 100047392|F|0-36, 14083801|F|0-28 and 100051488|F|0-49 were associated with AVSPD at p <0.001. SNP markers, 14074781|F|0-16, 14083801|F|0-28 and 100051488|F|0-49 were associated with PL at P <0.001. Five SNP markers, 100047392|F|0-36, 14083801|F|0-28, 100072738|F|0-34, 14076881|F|0-49 and 14076881|F|0-49 were associated with PWDTH at p <0.001. The 65 SNP markers identified can be used in cowpea molecular breeding to select for AVSPD, NP, PL, PWDTH, PWT, and RR through marker assisted selection (MAS).

Characterization of genetic diversity, population structure and linkage disequilibrium is prerequisite for proper management of breeding programs and conservation of genetic resources. In this study, 186 chickpea genotypes including advanced “Kabuli” breeding lines and Iranian landrace “Desi” chickpea genotypes were genotyped using DArTseq-Based SNP markers. Out of 3339 SNPs, 1152 markers with known chromosomal position were selected for genome diversity analysis. The number of mapped SNP markers varied from 52 (LG8) to 378 (LG4), with an average of 144 SNPs per linkage group. The chromosome size that covered by SNPs varied from 16236.36 kbp (LG8) to 67923.99 kbp (LG5), while LG4 showed higher number of SNPs, with an average of 6.56 SNPs per Mbp. Polymorphism information content (PIC) value of SNP markers ranged from 0.05 to 0.50, with an average of 0.32, while the markers on LG4, LG6 and LG8 showed higher mean PIC value than average. Un-weighted Neighbor Joining cluster analysis and Bayesian-based model population structure grouped chickpea genotypes into four distinct clusters. Principal component analysis (PCoA) and Discriminant Analysis of Principal Component (DAPC) results were consistent with that of the cluster and population structure analysis. Linkage disequilibrium (LD) was extensive and LD decay in chickpea germplasm was relatively low. A few markers showed r²≥0.8, while 2961 pairs of markers showed complete LD (r²=1) and a huge LD block was observed on LG4. High genetic diversity and low kinship value between pairs of genotypes suggesting the presence of a high genetic diversity among studied chickpea genotypes. This study also demonstrated the efficiency of DArTseq-based SNP genotyping for large scale genome analysis in chickpea. The genotypic markers provided in this study are useful for various association mapping studies when combined with phenotypic data of different traits such as seed yield, abiotic and biotic stresses and therefore can be efficiently used in breeding programs to improve chickpea.

(1) Background: Brassinosteroids (BRs) are an important hormone that participates in almost all stages of plant growth and development, and sterol dehydrogenase is a key enzyme that is involved in the biosynthesis of BRs. However, the steroid dehydrogenase gene family in Populus yunnanensis Dode (P. yunnanensis) has not been studied; (2) Methods: The PyDET2 gene family was identified and bioinformatics analyzed. Three genes were screened based on the RNA-seq of stem tip, the PyDET2e was further investigated by RT-qPCR (real-time quantitative polymerase chain reaction) and subcellular localization; (3) Results: We identified 14 DET2 family genes in P. yunnanensis that were classified into 4 groups according to sequence similarity and phylogenetic relationship, and 10 conserved protein motifs had been identified. The gene structure, chromosome distribution, collinearity, and codon preference of all PyDET2 genes identified in the genome of P. yunnanensis were analyzed. The codon of this family is biased to the end of A/U, and natural selection has a great influence on it. The expression level of PyDET2e in September was higher than in July, and it was significantly expressed in stems, stem tips, and leaves. The PyDET2e protein was localized in chloroplasts; (4) Conclusions: The PyDET2e is involved in the development of P. yunnanensis in the rapid growth stage. The systematic analysis provided a basis for the genome-wide identification of genes related to brassinolide biosynthesis processes in P. yunnanensis, and laid a foundation for the study on the mechanism of rapid growth of P. yunnanensis.

SARS-CoV2 (corona virus) has spread globally at an unprecedented rate; so far, increasing SARS-CoV2-infected individuals have been identified. Although the situation in China is improving and is currently under control, the outbreak in other countries and its pandemic management is only beginning to develop. Based on 154 SARS-CoV2 genome sequence analyses, we used receptor–ligand docking to identify one potential point mutation (V354F) on the spike structure which enhances spike binding to ACE2 receptors underlying potential super infection. Importantly, the V354F site on spike S1 had been identified in 5/10 infected French patients living in Paris, who sharing 100% identical SARS-CoV2 genomes. With Covid-19 cases increasing rapidly in France that could lead to a new explosion, we suggest that the French government should identify all potential super spreaders and treat them accordingly. In summary, our study provides on of the measures to avoid the potential second worldwide explosion of SARS-CoV2.

Mucosal melanomas (MM) are rare aggressive cancers in humans and one of the most common forms of oral cancers in dogs. Similar biological and histological features are shared between MM in both species making dogs a powerful model for comparative oncology studies of melanomas. Although exome sequencing recently identified recurrent coding mutations in canine MM, little is known about changes in non-coding gene expression and more particularly in canine long non-coding RNAs (lncRNAs), which are commonly dysregulated in human cancers. Here, we sampled a large cohort (n= 52) of canine normal/tumor oral MM from three predisposed breeds (poodles, Labrador retrievers and golden retrievers) and used deep transcriptome sequencing to identify more than 400 differentially expressed (DE) lncRNAs. We further prioritized candidate lncRNAs by comparative genomic analysis to pinpoint 26 dog-human conserved DE lncRNAs, including SOX21-AS, ZEB2-AS and CASC15 lncRNAs. Using unsupervised co-expression networks analysis with coding genes, we inferred potential functions of DE lncRNAs suggesting associations with cancer-related genes, cell cycle and carbohydrate metabolism GO terms. Finally, we exploited our multi-breed design to identify DE lncRNAs per breed. This study provides a unique transcriptomic resource for studying oral melanoma in dogs and highlights lncRNAs that may potentially be diagnostic or therapeutic targets for human and veterinary medicine.

Background: Thousands of genetic variants have been associated with hematological traits, though target genes remain unknown at most loci. Also, limited analyses have been conducted in African ancestry and Hispanic/Latino populations; hematological trait associated variants more common in these populations have likely been missed. Methods: To derive gene expression prediction models, we used ancestry-stratified datasets from the Multi-Ethnic Study of Atherosclerosis (MESA, including N=229 African American and N=381 Hispanic/Latino participants, monocytes) and the Depression Genes and Networks study (DGN, N = 922 European ancestry participants, whole blood). We then performed a transcriptome-wide association study (TWAS) for platelet count, hemoglobin, hematocrit, and white blood cell count in African (N = 27,955) and Hispanic/Latino (N = 28,324) ancestry participants. Results: Our results revealed 24 suggestive signals (p < 1×10^(-4)) that were conditionally distinct from known GWAS identified variants and successfully replicated these signals in European ancestry subjects from UK Biobank. We found modestly improved correlation of predicted and measured gene expression in an independent African American cohort (the Genetic Epidemiology Network of Arteriopathy (GENOA) study (N=802), lymphoblastoid cell lines) using the larger DGN reference panel; however, some genes were well predicted using MESA but not DGN. Conclusions: These analyses demonstrate the importance of performing TWAS and other genetic analyses across diverse populations and of balancing sample size and ancestry background matching when selecting a TWAS reference panel.

Sulfate transporters (SULTRs) are responsible for the uptake of the sulfate (SO42−) ions in the rhizosphere by the roots and their distribution in plant organs. In this study, SULTR family members in the genome of the two oilseed crops, Camelina sativa, and Brassica napus, were identified and characterized based on their sequence structure, duplication events, phylogenetic relationships, phosphorylation sites, and expression levels. Herein, 36 and 45 putative SULTR genes were recognized from the genome of C. sativa, and B. napus, respectively. SULTR proteins were predicted as basophilic proteins with low hydrophilicity in both studied species. According to phylogenetic relationships, we divided SULTRs into five groups, in which SULTRs 3 showed highest variation. Besides, several duplication events were observed between SULTRs. The first duplication event was predicted approximately five million years ago between three SULTRs 3.1 in C. sativa. Two subunits were indicated in the 3D structure of SULTRs that the active binding sites differed between C. sativa and B. napus. According to available RNA-seq data, SULTRs showed diverse expression in tissues and response to stimuli. SULTRs 3 showed an expression in all tissues. SULTRs 3.1 were more upregulated in response to abiotic stresses in C. sativa, while SULTRs 3.3, and SULTRs 2.1 showed an upregulation in B. napus. Furthermore, SULTRs 3 and SULTRs 4.1 showed an upregulation in response to biotic stresses in B. napus. Based on the distribution of cis-regulatory elements in the promoter region, we speculated that SULTRs might be controlled by phytohormones such as ABA, and MeJA. Therefore, it seems that SULTR genes in C. sativa have been more influenced by evolutionary processes and have acquired more diversity.

Glycans are involved in various metabolic processes via the functions of glycosyltransferases and glycoside hydrolases. Analysing the evolution of these enzymes is essential for improving the understanding of glycan metabolism and function. Based on our previous study of glycosyltransferases, we performed a genome-wide analysis of whole human glycoside hydrolases using the UniProt, BRENDA, CAZy, and KEGG databases. Using cluster analysis, 319 human glycoside hydrolases were classified into four clusters based on their similarity to enzymes conserved in chordates or metazoans (Class 1), metazoans (Class 2), metazoans and plants (Class 3), and eukaryotes (Class 4). The eukaryote and metazoan clusters included N- and O-glycoside hydrolases, respectively. The significant abundance of disordered regions within the most conserved cluster indicated a role for disordered regions in the evolution of glycoside hydrolases. These results suggest that the biological diversity of multicellular organisms is related to the acquisition of N- and O-linked glycans.

This paper is aimed at studying the temporal evolution of the surface displacements occurred over the past few years in the ocean-reclaimed platforms of the Shanghai megacity (China), which are mainly ascribable to consolidation processes of large dredger fills and alluvial deposits. With respect to previous analyses carried out over the same area, this work provides a joint multi-platform differential interferometry synthetic aperture radar (DInSAR) analysis, based on the application of the advanced Small BAseline Subset (SBAS) algorithm. This led us to retrieve long-term deformation time-series that are helpful for a better understanding of the on-going deformation phenomena. To this aim, we have exploited two sequences of SAR data collected by the ASAR/ENVISAT and by the COSMO-SkyMed (CSK) sensors, respectively, spanning the whole time period from 2007 to 2016. Unfortunately, the large time gap (of about three years) existing between the available ASAR/ENVISAT and CSK datasets gave rise to additional difficulties for their combination. Nevertheless, this problem has been faced by benefiting from the knowledge of a time-dependent model describing the temporal evolution of the expected deformations affecting the Shanghai ocean-reclaimed platforms.

Structural variations (SVs) play a key role in the pathogenicity of hematological malignancies. Optical genome mapping (OGM) is an emerging technology that enables genome-wide detection of all classes of SVs at a high resolution and sensitivity. Identification of cryptic SVs leading to gene disruption or predicted novel gene fusions could be important drivers for cancer development and/or portend a prognostic relevance, which could be used to modify the treatment plan. A cohort of 106 consented cases that had a successful OGM analysis performed were included in the study. Demographic, clinical, laboratory and treatment data were collected. Routine diagnostic and prognostic testing were done on the peripheral blood and bone marrow aspirate as indicated. Additional samples of peripheral blood and/or bone marrow were sent for OGM testing. OGM led to a change in risk stratification in 17/66 (25.75%) of patients with hematological malignancies, the majority of these patients (15/66, 22.72%), having their risk stratification upgraded with a resulting change in treatment of 14 patients. This study highlights the ability of OGM to detect rare, cryptic and clinically relevant variants that potentially impact disease diagnosis, risk stratification and actionable treatment targets.

Growth-regulating factor (GRF) is a plant-specific transcription factor family, which is involved in nearly all of the central developmental processes in plants. However, little is known about GRF family genes in cultivated sunflower. In this study, 17 GRF genes were identified and characterized from sunflower genome. Their gene structures, conserved motifs, chromosomal distributions and phylogenetic relationships were analyzed. The expression patterns of these genes were detected in various tissues of sunflower inbred line SK02R, which revealed that the 10 seed-specific GRF genes may play important roles during seed development in sunflower. Additionally, transcripts changes of the GRFs under two major abiotic stresses and phytohormones showed that most of the detected GRFs were reduced significantly by GA3 treatment, and other treatments(ABA, NaCl and PEG6000) differently regulated various sunflower growth-regulating factors at different time points. MiR396 target analysis indicated that there may exist a complicated homeostasis between miR396 and its targets GRF and WRKY transcription factor genes in cultivated sunflower. The phylogenetic and expression analyses of the GRF gene family in sunflower would be useful for further cloning and function exploration of the HaGRF genes.

Melatonin is a crucial regulator of plant growth and development as well as stress tolerance. However, we only have a limited understanding of the functions of endogenous melatonin. Tryp-tophan decarboxylase (TDC) serves as the initial rate-limiting enzyme in the melatonin synthesis pathway. Unfortunately, no cucumber TDC gene has been cloned and characterized. In this study, we identified two TDC genes (CsTDC1 and CsTDC2) in the cucumber genome. The subcellular localization analysis indicated CsTDC1 and CsTDC2 are predominantly localized in the cyto-plasm and plasma membrane. Tissue-specific expression analyses revealed that CsTDC1 and CsTDC2 are expressed in both vegetative and reproductive organs. Many cis-elements related to stress, hormone, and light responses as well as development were identified in the CsTDC pro-moter regions. Furthermore, the expression of CsTDC1 and CsTDC2 was strongly induced by treatments with various abiotic stresses and exogenous hormones. The transient overexpression of CsTDC1 and CsTDC2 in tobacco leaves resulted in increases in the TDC activity and melatonin content, which may be related to the enhanced ability of tobacco leaves to tolerate salt, drought, and low-temperature stresses. Notably, the overexpression of CsTDC2 had a more pronounced effect than the overexpression of CsTDC1. Accordingly, both CsTDC genes, but especially CsTDC2, may be important for regulating cucumber growth, development, and stress tolerance. The study findings provide a theoretical basis for future functional analyses of endogenous melatonin in cucumber.

MADS-box genes play a crucial role in fruit ripening, yet limited research has been conducted on mango. Based on the conserved domains of this gene family, 84 MADS-box genes were identified in the mango genome, including 22 type I and 62 type II MADS-box genes. Gene duplication analysis revealed that both tandem duplication and segmental replication significantly contributed to the expansion of MADS-box genes in the mango genome, with purifying selection playing a vital role in the segmental duplication events within the MiMADS gene family. Cis-acting element analysis demonstrated that most MiMADS genes were hormonally regulated and participated in the growth, development, and stress resistance of mango fruit. Moreover, through expression pattern analysis and phylogenetic tree construction, we identified six MiMADS genes belonging to the SEP1 subfamily and two belonging to the AG subfamily as potential candidates involved in mango ripening regulation. Notably, Mi08g17750 and Mi04g18430 from the SEP1 subfamily were identified as key regulators inhibiting mango fruit maturation; their interaction network was also analyzed. These findings provide a foundation for further investigation into the regulatory mechanisms underlying mango ripening.

Cotton leaf curl Kokhran virus (CLCuKoV) (genus, Begomovirus; family, Geminiviridae) is highly infectious, widespread and the most dangerous pathogen of cotton (Gossypium hirsutum L.) that is responsible to a serious disorder, cotton leaf curl disease (CLCuD). Begomoviruses are spread very efficiently by the whitefly Bemisia tabaci cryptic species, causing economic losses to cotton crop, all over the world. The ‘Lucknow’ strain of CLCuKoV has emerged as a divergent isolate that could cause CLCuD. The monopartite ssDNA genome of CLCuKoV-Lu (2.7 Kb) contains six open reading frames (ORFs) that was shown to encode four major proteins. RNA interference (RNAi)-based antiviral innate immunity is a sequence-specific biological phenomenon and a powerful tool to control plant viruses. The present study aims to determine cotton locus-derived microRNAs (ghr-miRNAs) that are identified for targeting the CLCuKoV-Lu ss-DNA-encoded mRNAs using a predictive approach that involves four computational algorithms, miRanda, RNA22, psRNATarget and RNAhybrid. Mature ghr-miRNA sequences (n=80) from allotetraploid upland cotton (2n = 4x = 52) were selected from the miRBase and were tested for alignment with the CLCuKoV-Lu genome. Among the 80 cotton locus-derived ghr-miRNAs evaluated, only one consensus cotton locus-derived ghr-miRNA (ghr-miR2950) was concluded to have effective ghr-miRNA target site at common nucleotide position 82 in the CLCuKoV-Lu genome respectively, using a stringent criterion, identified by all the algorithms used. The miRNA targeting is reliant on base pairing of miRNA-mRNA target pairings. Conservation of the hybridization binding site of the predicted ghr-miR2950 was validated using multiple sequence alignment within all the strains of CLCuKoV. We constructed a regulatory interaction network of miRNA–mRNA to identify novel targets. The efficacy of the predicted miRNAs against CLCuKoV-Lu was evaluated by RNAi-mediated targeted mRNA cleavage. The current investigated miRNA targets provide evidence for the development of CLCuD-resistant cotton plants.

Due to the global warming and dynamic changes in pathogenic virulence, leaf rust caused by Puccinia triticina has greatly expanded its epidermic region and become a severe threat to global wheat production. Genetic bases of wheat resistance to leaf rust mainly relies on the leaf rust resistance (Lr) gene or quantitative trait locus (QLr). Although these genetic loci have been insensitively studied during the last two decades, an updated overview of Lr/QLr in a genome-wide level is urgently needed. This review summarized recent progresses in genetic studies of wheat resistance to leaf rust. Wheat germplasms with great potentials in genetic improvement of resistance to leaf rust were highlighted. Key information about the genetic loci carrying Lr/QLr were summarized. A genome-wide chromosome distribution map for all the Lr/QLr was generated based on the released wheat reference genome. In conclusion, this review has provided valuable sources for both wheat breeders and researchers to understand the genetics of resistance to leaf rust in wheat.

Background: Telomere length (TL) shortening process is associated with several known environment and individual determinants. DNA methylation is the most studied epigenetic process and may be associated with TL. We investigated the associations between DNA methylation and TL in peripheral blood. Methods: Methylation wide association study was conducted in 47 women (37.1±8.8 years) with different nutritional status. Association between TL and DNA methylation levels were explored by univariate and multiple linear regression models, corrected by age and Body Mass Index. Corrections for multiple comparisons by Benjamini‐Hochberg test was also performed. WEBGestalt was used to identify pathways that are responsible for regulating TL. Results: We found negative correlations between TL and BMI (r = -0.641; p = 0.001), abdominal circumference (r = -0.622; p = 0.001) and fat mass (r = -0.656; p = 0.001). 44 CpGs sites were associated with TL, independent of age and BMI. The most of these sites were negatively correlated with TL. For the 7 remained sites, DNA hypomethylation were associated with shorter TL. These CpGs were related to nine different pathways, including thermogenesis, cancer, glutamatergic and serotonergic synapse. Conclusion: There is an epigenetic contribution in TL, independent of nutritional status and age. Genes related to TL are involved in important metabolic pathways.

Magnesium transporters (MGTs) play a prominent role in the absorption, transport and storage of magnesium in plant cells. In the present study, MGT gene family members were identified and characterized in two species of cucurbitaceae, including C. sativus and C. lanatus. 20, 19, and 20 MGT genes were recognized in C. lanatus, C. sativus, and C. melo, respectively. According to physicochemical properties, the members of each sub-class of MGTs in the species of cucurbitaceae showed the close relationship. Proteins from NIPA were identified as hydrophilic proteins with high stability. Based on phylogenetic analysis, MGT family members were classified into three groups, and NIPAs showed more diversity. Besides, duplication events were not identified between the MGT members in C. lanatus, and C. sativus. According to pocket analysis, residues such as L, V, S, I, and A were frequently observed in the binding sites of MGT proteins in both species. The prediction of post-translation modifications revealed that MSR2 proteins have high phosphorylation potential than other sub-classes in both studied plants. The expression profile of MGTs showed that MGTs are more expressed in root tissues. In addition, MGTs showed differential expression in response to abiotic/biotic stresses as well as hormone application and NIPAs were more induced in response to stimuli in watermelon. The results of this study, as the primary work of MGT gene family, can be used in programs related to cucurbitaceae breeding.

Calcineurin B-like proteins (CBL)-interacting protein kinases (CIPKs) by interacting with CBLs regulate developmental processes, hormone signalling transduction and mediate stress responses in plants. Although the genome of chickpea is available, information of CIPK gene family has been missing in chickpea. Here, a total of 22 CIPK encoding genes were identified in chickpea and characterized by in silico methods. We found a high structural conservation in chickpea CIPK family. Our analysis showed that chickpea CIPKs have evolved with dicots from common ancestors, and extensive gene duplication events have played an important role in evolution and expansion of CIPK family in chickpea. Most chickpea CIPK proteins localize in cytoplasm and nucleus. Promoter analysis revealed various cis-regulatory elements related to plant development, hormone signaling and abiotic stresses. Expression analysis indicated that CIPKs are significantly expressed in a spectrum of developmental stages, tissue/organs that hinted their important role in plant development. Several CIPK genes had specific and overlapping expressions in different abiotic stresses and seed development stages, suggesting the important role of CIPK family in abiotic stress signaling, and seed development in chickpea. Thus, this study provides the avenue for detailed functional characterization of CIPK family in chickpea and other legume crops.

The CDPK family genes play crucial roles in signal transduction pathways during plant develop-ment and stress response. In this study, we comprehensively analyzed the CDPK family genes in mustard (Brassica juncea L.), resulting in the identification of 101 genes (BjuCDPK1–101) located on chromosomes AA_Chr01 to BB_Chr08. RNA-seq analysis showed that most BjuCDPK genes were predominantly expressed in the root and flower bud, suggesting their organ-specific expression patterns. Furthermore, the expression levels of most BjuCDPK genes were significantly altered by cold stress. The promoter regions of BjuCDPK genes had abundant cis-acting elements related to hormone and stress. The interaction network analysis predicted that BjuCDPK proteins respond to abiotic stress probably via interaction with ABF1, ABF4, and/or DI19. Our findings provide valu-able information for further investigation of the cold stress adaption of mustard via the CDPK signaling pathway.

Alfalfa (Medicago sativa L.) is the most cultivated forage legume around the world. Under a variety of growing conditions, forage yield in alfalfa is stymied by biotic and abiotic stresses including heat, salt, drought, and disease. Given the sessile nature of plants, they use strategies such as differential gene expression to respond to environmental cues. Transcription factors control the expression of genes that contribute to or enable tolerance and survival during periods of stress. Basic-leucine zipper (bZIP) transcription factors have been demonstrated to play a critical role in regulating plant growth and development as well as mediate the responses to abiotic stress in several species, including Arabidopsis thaliana, Oryza sativa, Lotus japonicus and Medicago truncatula. However, there is little information about bZIP transcription factors in cultivated alfalfa. In the present study, 237 bZIP genes were identified in alfalfa from publicly available sequencing data. Multiple sequence alignments showed the presence of intact bZIP motifs in the identified sequences. Based on previous phylogenetic analyses in Arabidopsis thaliana, alfalfa bZIPs were similarly divided and fell into 10 groups. The physico-chemical properties, motif analysis and phylogenetic study of the alfalfa bZIPs revealed high specificity within groups. The differential expression of alfalfa bZIPs in a suite of tissues indicates that particular bZIP genes are specifically expressed at different developmental stages in alfalfa. Similarly, expression analysis in response to ABA, cold, drought and salt stresses, indicates that a subset of bZIP genes are also differentially expressed and likely play a role in abiotic stress signaling and/or tolerance. However, further functional characterization of bZIP transcription factors in alfalfa will help illuminate the role they play in stress tolerance mechanisms in legumes and facilitate the molecular breeding of stress tolerance in alfalfa.

Although the regulatory function of miRNAs and their targets have been characterized in model plants, a possible underlying role in the cotton response to herbivore infestation has not been determined. To investigate this, we performed small RNA and degradome sequencing between resistant and susceptible cotton cultivar following infestation with the generalist herbivore whitefly. In total, 260 miRNA families and 241 targets were identified. Quantitative-PCR analysis revealed that several miRNAs and their corresponding targets exhibited dynamic spatio-temporal expression patterns. Moreover, 17 miRNA precursors were generated from 29 long intergenic non-coding RNA (lincRNA) transcripts. Genome-wide analysis also led to the identification of 85 phased small interfering RNA (phasiRNA) loci. Among these, nine PHAS genes were triggered by miR167, miR390, miR482a, and two novel miRNAs, including those encoding a leucine-rich repeat (LRR) disease resistance protein, an auxin response factor (ARF) and MYB transcription factors. Through combined modeling and experimental data, we explored and expanded the miR390-tasiARF cascade during the cotton response to whitefly. Virus-induced gene silencing (VIGS) of ARF8 in whitefly-resistant cotton plants increased auxin and jasmonic acid (JA) accumulation, resulting in an increased tolerance to whitefly infestation. These results highlight the provides a useful transcriptomic resource for plant-herbivore interaction.

Nitrogen (N) is an important mineral nutrient for plant growth, as well as a limiting factor for crop yield, and how to improve the utilization efficiency of N fertilizer by plants is currently a research hotspot. This article uses bioinformatics methods to identify and analyze members of the glutamine synthetase (GS) and glutamate synthase (GOGAT) gene families in pecan. A total of 6 GS genes and 4 GOGAT genes were identified, and their physicochemical properties, gene structures, and homologous evolutionary relationships were analyzed. Analysis of tissue-specific expression of GS and GOGAT genes based on transcriptome data from pecan. The enzyme activities of GS and GOGAT and the gene expression were quantitatively analyzed under different N form ratios in pecan. According to the results, the promoter cis-acting elements of GS and GOGAT genes can be roughly divided into three types: light-responsive elements, hormone-responsive elements, and stress-responsive elements. The results of homologous evolution showed that there was no tandem duplication event for the two gene families, and GS and GOGAT have undergone purification during the evolutionary process. CiGS2s and CiFd-GOGATs were expressed mainly in leaves, and CiNADH-GOGATs were expressed mainly in fruits. The qPCR analysis results showed that T4 treatment significantly increased the expression levels of CiGS and CiGOGAT genes in the leaves. The enzyme activities of GS and GOGAT in pecan were significantly increased under T3, T4, and T5 treatments. In summary, a higher proportion of ammonium nitrogen (NH4+) in the nutrient solution was profit to pecan NH4+ assimilation. This study determined the appropriate nitrogen ratio for pecan, promoting a theoretical basis for reducing environmental pollution caused by nitrogen fertilizer and improving the nitrogen utilization efficiency of pecan. In summary, both CiGSs and CiGOGATs exhibit tissue specificity, and an ammonium-nitrate mixture with a higher proportion of NH4+ is more favorable for NH4+ assimilation in pecan. This study provides a reference basis for further understanding the functions of CiGSs and CiGOGATs in pecan, and offers a theoretical foundation for improving N use efficiency in pecan.

Magnesium (Mg) is a vital element, which involves in various key cellular processes in plants. Mg transporter (MGT) genes play an important role in magnesium distribution and ionic balance maintenance. Here, MGT family members were identified and characterized in three species of the plant family Malvaceae, including Theobroma cacao, Corchorus capsularis, and Gossypium hirsutum to improve our understanding of their structure, regulatory systems, functions, and possible interactions. We identified 18, 41, and 16 putative non-redundant MGT genes from the genome of T. cacao, G. hirsutum, and C. capsularis, respectively, which were clustered into five groups within maximum likelihood tree. Several segmental/tandem duplication events were determined between MGTs genes. It seems that MGTs were slowly evolving and have been evolved under a purifying selection. Analysis of gene promoter regions showed that MGTs have a high potential to respond to biotic/abiotic stresses as well as hormones. The expression patterns of MGT genes revealed their possible T. cacaorole in response to P. megakarya fungi in T. cacao, whereas MGT genes showed differential expression in various tissues and respond to several abiotic stresses, including cold, salt, drought, and heat stress in G. hirsutum. Besides, the co-expressions network of MGTs indicated that genes involved in auxin responsive, lipid metabolism, cell wall organization, and photoprotection can interact with MGTs.

Fibrillin (FBN) is a plastid lipid-associated protein found in photosynthetic organisms from cyanobacteria to plants. In this study, 10 CsaFBN genes were identified in genomic DNA sequences of cucumber (Chinese long and Gy14) through database searches using the conserved domain of FBN and the 14 FBN genes of Arabidopsis. Phylogenetic analysis of CsaFBN protein sequences showed that there was no counterpart of Arabidopsis and rice FBN5 in the cucumber genome. FBN5 is essential for growth in Arabidopsis and rice; its absence in cucumber may be because of incomplete genome sequences or that another FBN carries out its functions. Among the 10 CsaFBN genes, CsaFBN1 and CsaFBN9 were the most divergent in terms of nucleotide sequences. Most of the CsaFBN genes were expressed in the leaf, stem, and fruit. CsaFBN4 showed the highest mRNA expression levels in various tissues, followed by CsaFBN6, CsaFBN1, and CsaFBN9. High-light stress combined with low temperature decreased photosynthetic efficiency and highly induced transcript levels of CsaFBN1, CsaFBN6, and CsaFBN11, which decreased after 24 h treatment. Transcript levels of the other seven genes were changed only slightly. This result suggests that CsaFBN1, CsaFBN6, and CsaFBN11 may be involved in photoprotection under high-light conditions at low temperature.

Stalk rot caused by Fusarium fungi is one of the most widespread and devastating diseases of maize, and the introduction of resistant genotypes is one of the most effective strategies for controlling the disease. Breeding genotypes with genetically determined resistance will also allow less use of crop protection products. The aim of the research was to identify molecular markers and associated candidate genes determining maize plant resistance to Fusarium stalk rot. The plant material for the study consisted of 122 maize hybrids, which were sown in plots of 10 m2, in a complete randomized block design, in three repetitions, at two locations (Smolice: 51 42’ 58.904’’ N, 17 13’ 29.13’’ E and Kobierzyce: 50 58’ 19.411’’ N, 16 55’ 47.323’’ E). The analyzed genotypes were simultaneously subjected to next-generation sequencing using the Illumina platform. The results of the observation of the degree of infection and sequencing were used for association mapping, which ultimately resulted in the selection of 10 molecular markers important at both places. Among the identified markers, two SNP markers that are located inside candidate genes play an important role. Marker 4772836 is located inside the serine/threonine-protein kinase bsk3 gene, while marker 4765764 is located inside the histidine kinase 1 gene. Both genes can be associated with plant resistance to Fusarium stalk rot, and these genes can also be used in breeding programs to select resistant varieties.

Pv11, an insect cell line established from the midge Polypedilum vanderplanki, is capable of ametabolic desiccation tolerance, so-called anhydrobiosis. We previously discovered that heat shock factor 1 (HSF1) contributes to the acquisition of desiccation tolerance by Pv11 cells, but the mechanistic details have yet to be elucidated. Here, by analyzing the gene expression profiles of newly established HSF1-knockout and -rescue cell lines, we show that HSF1 has a genome-wide effect on gene regulation in Pv11. HSF1-knockout cells exhibit a reduced desiccation survival rate, but this is completely restored in HSF1-rescue cells. By comparing mRNA profiles of the two cell lines, we reveal that HSF1 induces anhydrobiosis-related genes, especially genes encoding late embryogenesis abundant proteins and thioredoxins, but represses a group of genes involved in basal cellular processes, thus promoting an ametabolic state in the cell. In addition, HSF1 binding motifs are enriched in the promoters of anhydrobiosis-related genes and we demonstrate binding of HSF1 to these promoters by ChIP-qPCR. Thus, HSF1 directly regulates the transcription of anhydrobiosis-related genes and consequently plays a pivotal role in the induction of anhydrobiotic ability in Pv11 cells.

Plant Cysteine Oxidase (PCO) is a plant O2-sensing enzyme of catalyzing oxidation of cysteine to Cys-sulfinic acid at the N-termini of target proteins. To better understand the Brassica napus PCO gene family, PCO genes in B. napus and related species were analyzed. In this study, 20, 7 and 8 PCO genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the PCOs were divided into five groups: PCO1, PCO2, PCO3, PCO4 and PCO5. Gene organization and motif distribution analysis suggested that PCO gene family was relatively conserved during evolution. Expression analysis showed that PCO genes were expressed in different tissues at different developmental stages. Interestingly, most of the Bna/Br/BoPCO5 members were expressed in leaves, roots, flowers, buds and siliques, suggesting an important role in both vegetative and reproductive development. Expression of BnaPCO was induced by various abiotic stress, especially waterlogging stress, which was consistent with the result of cis-element analysis. In this study, the PCO gene family of Brassicaceae was analyzed for the first time, which contributes to a comprehensive understanding of the origin and evolution of PCO genes in Brassicaceae and the function of BnaPCO in regulating plant responses to abiotic stresses.

Transcription factors (TFs) are the regulatory proteins that act as molecular switches in controlling stress responsive gene expression. Among them MYB transcription factor family is one of the largest TF family in plants, playing a significant role in plant growth, development, phytohormone signaling and stress-responsive processes. Pearl millet (Pennisetum glaucum L.) is one of the most important C4 crop plant of the arid and semi-arid regions of Africa and South-east Asia for sustaining food and fodder productions. To explore the evolutionary mechanism and functional diversity of the MYB family in pearl millet, we conducted a comprehensive genome-wide survey and identified 279 MYB TFs (PgMYB) in pearl millet and distributed unevenly across seven chromosomes of pearl millet. Phylogenetic analysis of identified PgMYBs classified them into 18 subgroups and members of the same group showed a similar gene structure and conserved motif/s pattern. Further, duplication events were identified in pearl millet that indicated towards evolutionary progression and expansion of the MYB family. Transcriptome data and relative expression analysis by qRT-PCR identified differentially expressed candidate PgMYBs (PgMYB2, PgMYB9, PgMYB88 and PgMYB151) under dehydration, salinity, heat and phytohormones (ABA, SA and MeJA) treatment. Taken together, this study provides valuable information for a prospective functional characterization of MYB family members of pearl millet and genetic improvement of crop plants.

Leucine-rich repeat receptor-like kinases (LRR-RLKs) have been reported to play important roles in plant growth, development and stress responses. However, no comprehensive analysis of this family has been performed in Gossypium, which are important economic crop and suffer various stresses in growth and development. Here we conducted a comprehensive analysis of LRR-RLK family in four Gossypium species (G. arboreum, G. barbadense, G. hirsutum and G. raimondii). A total of 1641 LRR-RLK genes were identified in the four Gossypium species involved in our study. Maximum-likelihood phylogenetic tree revealed that all the LRR-RLK genes were divided into 21 subgroups. Exon-intron organization structure of LRR-RLK genes kept relative conserved in subfamilies and between Arabidopsis and Gossypium. Subfamilies XI and XII were found dramatically expanded in Gossypium. Tandem duplication acted as an important mechanism in expansion of Gossypium LRR-RLK gene family. Function analysis suggested that plant hormone signaling and plant-pathogen interaction pathway were enriched in Gossypium LRR-RLK genes. Promoters analysis and expression profiles analysis revealed that Gossypium LRR-RLK genes were extensively regulated by TFs, phytohormone and various environmental stimuli, and play key roles in stress defense and diverse development processes. Our study provided valuable information for further function study of Gossypium LRR-RLK genes.

Multidisciplinary research efforts on potential COVID-19 vaccine and therapeutic candidates have increased since the pandemic outbreak of SARS-CoV-2 in 2019. This search has become imperative due to the increasing emergences and limited widely available medicines. The presence of bioactive anti-SARS-CoV-2 molecules was examined from various plant sources. Among them is a group of proteins called lectins that can bind carbohydrate moieties. In this article, we present ten novel, chitin-specific Hevein-like lectins that were derived from Selaginella moellendorffii v1.0's genome. The capacity of these lectin homologs to bind with the spike protein of SARS-CoV-2 was examined. Using the HDOCK server, 3D-modeled Hevein-domains were docked to the spike protein's receptor binding domain (RBD). The Smo446851, Smo125663, and Smo99732 interacted with Asn343-located complex N-glycan and RBD residues, respectively, with binding free energies of -17.5, -13.0, and -26.5 Kcal/mol. The normal-state analyses via torsional coordinate association for the Smo99732-RBD complex using iMODS is characterized by overall higher stability and minimum deformity than the other lectin complexes. The three lectins interacting with carbohydrates were docked against five individual mutations that frequently occur in major SARS-CoV-2 variants. These were in the spike protein’s receptor-binding motif (RBM), while Smo125663 and Smo99732 only interacted with the spike glycoprotein in a protein-protein manner. The precursors for the Hevein-like homologs underwent additional characterization and their expressional profile in different tissues was studied. These in-silico findings offered potential lectin candidates targeting key N-glycan sites crucial to the virus's virulence and infection.

Factors affecting the host-virus interaction must be understood for the effective application of bacteriophages to combat bacterial pathogens. Two novel E. coli phages, the T1-like siphophage LL5 and the rV5-like myophage LL12, were subjected to forward genetic screens against the Keio collection, a library of single non-essential gene deletions in E.coli str. BW25113. These genome-wide screens and subsequent experiments identified eight genes required for efficient propagation of phage LL5 and six genes required for propagation of LL12. The majority of the genes identified were involved in production of the phage receptors. E. coli mutants deficient in heptose II and the phosphoryl substituent of heptose I of the inner core lipopolysaccharide (LPS) were unable to propagate phage LL5, as were mutants deficient in the outer membrane protein TolC. Mutants lacking glucose I of the LPS outer core failed to propagate LL12. Two cytoplasmic chaperones, PpiB and SecB, were found to be required for efficient propagation of phage LL5 but not LL12. This approach may be useful for identifying phage receptors and required host factors in other phages, which would provide valuable information for their potential use as therapeutics and for phage engineering.

Background: Triple negative breast cancer (TNBC) is the most aggressive form of breast cancer with poor outcomes. The molecular basis of TNBC remains poorly understood. The objective of this study was to explore the relationship between obesity and TNBC in premenopausal and postmenopausal Caucasian women using whole genome transcription profiling. Methods: We compared gene expression levels of tumor samples drawn from normal weight, overweight and obese in pre and postmenopausal women diagnosed with TNBC. We performed hierarchical clustering to assess similarity in patterns of gene expression profiles, and conducted network and pathway analysis to identify molecular networks and biological pathways. Results: We discovered gene signatures distinguishing normal weight from obese, normal weight from overweight and overweight from obese individuals in both premenopausal and postmenopausal women. The analysis revealed molecular networks and biological pathways dysregulated in response to obesity. Among the discovered pathways included the unfolded protein response, endoplasmatic reticulum stress, B cell receptor and the autophagy signaling pathways in obese premenopausal women and the integrin, axonal guidance, ERK/MAPK and Glutathione biosynthesis signaling pathways obese postmenopausal women. Conclusions: The results suggest that both overweight and obesity are associated with TNBC, highlighting the need for conformation of these results in independent studies.

Background: Alfalfa, the most economically important forage legume worldwide, features modest genetic progress due to long selection cycles and large extent of non-additive genetic variance associated with its autotetraploid genome. Methods: To improve the efficiency of genomic selection in alfalfa, we explored the effects of genome parametrization (as autotetraploid and diploid dosage, plus allele ratios) and SNP marker subsetting (all available SNPs; only genic regions; only non-genic regions) on genomic regressions, together with various levels of filtering on reading depth and missing rates. We used Genotyping by Sequencing-generated data and focused on traits of different genetic complexity, i.e., dry biomass yield in moisture-favorable (FE) and drought stress (SE) environments, leaf size, and onset of flowering, assessed on 143 genotyped plants from a genetically-broad European reference population and their phenotyped half-sib progenies. Results: On average, allele ratio improved the predictive ability compared with other genome parametrizations (+7.9% vs. tetraploid dosage, +12.6% vs. diploid dosage), SNP subsetting offered an advantage compared with any tested subsettings (+3.7% vs. genic regions, +7.6% vs. non-genic regions). However, when focusing on specific traits, different combinations of genome parametrization and subsetting reached better performances. We also released Legpipe2, a SNP calling pipeline tailored for reduced representation (GBS, RAD) medium-sized genotyping experiments.

The largest gene families in plants were found to be Glutathione transferase (GSTs), which played significant roles in regulating plant growth, development, and stress response. Within the GSTs gene family, members were found to play a crucial role in the low-temperature response process of plants. A comprehensive study identified a total of 70 BraGSTs genes. Cluster analysis results demonstrated that the BraGSTs in Brassica rapa(B.rapa) could be categorized into eight sub-families and were unevenly distributed across ten chromosomes. The 39 BraGSTs genes were found to be organized into 15 tandem gene clusters, with the promoters containing multiple cis-elements associated with low-temperature response. Cold stress was observed to stimulate the expression of 15 genes, with the BraGSTF2 gene exhibiting the highest level of expression, sug-gesting its significant involvement in winter B.rapa's response to low-temperature stress. Sub-cellular localization analysis of the BraGSTF2 protein indicated its potential expression in both the cell membrane and nucleus. The analysis of stress resistance in BraGSTF2 transgenic Arabidopsis thaliana lines emon-strated that the over-expression of this gene resulted in significantly elevated levels of SOD, POD activity, and SP content compared to the wild type following exposure to low temperatures. These levels reached their peak after 24 hours of treatment. Conversely, the MDA content was lower in the transgenic plants compared to the wild type. Additionally, the survival rate of BraGSTF2 transgenic Arabidopsis thaliana was higher than that of the wild type, suggesting that the BraG-STF2 gene may play a crucial role in enhancing the cold stress tolerance of winter B. rapa. This study lay a foundation for further study on the role of BraGSTs gene in the molecular regulation of cold resistance in winter B.rapa.

The basic helix-loop-helix (bHLH) transcription factor family is the second-largest transcription factor in plants. Members of this family are involved in the process of growth and development, secondary metabolic biosynthesis, signal transduction, and plant resistance. Loropetalum chinense var. rubrum is a critical woody plant with higher ornamental and economic values, which used as ornamental architecture and traditional Chinese herbal medicine plant. However, the bHLH transcription factor in the L. chinense var. rubrum (Loropetalum chinense var. Rubrum) have yet to be systematically demonstrated, and their role in anthocyanin biosynthesis remains secret. Here, we identified 165 potential LcbHLHs genes by two methods, and they were unequally distributed on chromosomes 1 to 12 of the L. chinense genome. Based on a phylogenetic comparison with proteins from Arabidopsis, these bHLH proteins were categorized into 21 subfamilies. Most LcbHLHs in a given subfamily had similar gene structures and conserved motifs. The gene ontology annotation and cis-elements predicted that the LcbHLHs had many molecular functions and were involved in plant growth processes, including flavonoid and anthocyanin biosynthetic processes, etc. Transcriptomic analysis revealed different expression patterns among different tissues and cultivars in L. chinense. Many LcbHLHs were expressed in leaves, and only a few genes were highly expressed in flowers. Six LcbHLH genes differentially expressed between species and periods in color variation, which may play a key role in anthocyanin synthesis. These further functional predictions of LcbHLHs were established by quantitative real-time PCR analysis and protein-protein interaction network. This study identified the six potential bHLH genes related to anthocyanin in L. chinense and created a solid framework for future research on the function and evolution of bHLHs.

Middle East region experiences a high prevalence of hypovitaminosis D, yet most genetic studies on Vitamin D have focused on European populations. Furthermore, there is a lack of research on the genomic risk factors affecting the elderly population, who are more susceptible to health burden. We investigated the genetic determinants of 25-hydroxyvitamin D levels in elderly Lebanese individuals (n=199) through a whole exome-based genome-wide association study. We identified new loci with suggestive evidence of an association with Vitamin D levels, including rs141064014 in the MGAM gene (P-value of 4.40 × 10−06) and rs7036592 in PHF2 (P-value of 8.43 × 10−06). A meta-analysis of the Lebanese data and the largest European genome-wide association study confirmed consistence replication of numerous variants, including rs2725405 in SLC38A10 (P-value of 3.73 x 10-08). Despite the lower performance of European-derived polygenic risk scores compared to the European estimations, it still effectively predicted Vitamin D deficiency among elderly Lebanese individuals. Our findings provide novel insights into the genetic mechanisms of Vitamin D deficiency in Middle Eastern elderly populations, facilitating the development of personalized approaches for more effective management of hypovitaminosis D. Additionally, we demonstrated that whole exome-based genome-wide association study is an effective method for identifying genetic components associated with phenotypes.

For multi-sensor integrated systems, such as a mobile mapping system (MMS), data fusion at sensor-level, i.e., the 2D-3D registration between optical camera and LiDAR, is a prerequisite for higher level fusion and further applications. This paper proposes a line-based registration method for panoramic images and LiDAR point cloud collected by a MMS. We first introduce the system configuration and specification, including the coordinate systems of the MMS, the 3D LiDAR scanners, and the two panoramic camera models. We then establish the line-based transformation model for panoramic camera. Finally, the proposed registration method is evaluated for two types of camera models by visual inspection and quantitative comparison. The results demonstrate that the line-based registration method can significantly improve the alignment of the panoramic image and LiDAR datasets under either the ideal spherical or the rigorous panoramic camera model, though the latter is more reliable.

Rice, as a staple crop feeding billions, faces constant threats from various diseases jeopardizing global food security. Precise understanding of disease resistance mechanisms is crucial for developing resilient rice varieties. Traditional genetic mapping methods, such as QTL mapping, provide valuable insights into the genetic basis of diseases. However, the complex nature of rice diseases demands a holistic approach to gain accurate knowledge of it. Omics technologies, including genomics, transcriptomics, proteomics, and metabolomics, enable a comprehensive analysis of biological molecules, uncovering intricate molecular interactions within the rice plant. Integration of various mapping techniques using multi-omics data has revolutionized our understanding of rice disease resistance. By overlaying genetic maps with high-throughput omics datasets, researchers can pinpoint specific genes, proteins, or metabolites associated with disease resistance. This integration enhances the precision of disease-related biomarkers with better understanding of their functional roles in disease resistance. Improvement of rice breeding for disease resistance through this integration represents a significant stride in agricultural science because better understanding of the molecular intricacies and interactions underlying disease resistance architecture leads to more precise and efficient development of resilient and productive rice varieties. In this review, we explore how the integration of mapping and omics data can give a transformative impact on rice breeding for enhancing disease resistance.

The heavy-metal-associated (HMA) proteins are a class of PB1-Type ATPases related to the intracellular transport and detoxification of metals. However, due to a lack of information regarding the HMA gene family in the Cucurbitaceae family, a comprehensive genome-wide analysis of the HMA family was performed in ten Cucurbitaceae species: Citrullus amarus, Citrullus colocynthis, Citrullus lanatus, Citrullus mucosospermus, Cucumis melo, Cucumis sativus, Cucurbita maxima, Cucurbita moschata, Cucurbita pepo, and Legenaria siceraria. We identified 103 Cucurbit HMA proteins with various members, ranging from 8 (Legenaria siceraria) to 14 (Cucurbita pepo) across species. The phylogenetic and structural analysis confirmed that the Cucurbitaceae HMA protein family could be further classified into two major clades: Zn/Co/Cd/Pb and Cu/Ag. GO annotation-based subcellular localization analysis predicted that all HMA gene family members were localized on membranes. Moreover, the analysis of conserved motifs and gene structure (intron/exon) revealed the functional divergence between clades. Interspecies microsynteny analysis demonstrated that maximum orthologous genes were found between species of the Citrullus genera. Finally, nine candidates HMA genes were selected, and their expression analysis was carried out via qRT-PCR in root, leaf, flower, and fruit tissues of C. pepo under arsenic stress. The expression pattern of the CpeHMA genes showed a distinct pattern of expression in root and shoot tissues, with a remarkable expression of CpeHMA6 and CpeHMA3 genes from the Cu/Ag clade. Overall, this study provides insights into the functional analysis of the HMA gene family in Cucurbitaceae species and lays down the basic knowledge to explore the role and mechanism of the HMA gene family to cope with arsenic stress conditions.

Recognition, decoding and tracking of vortex patterns is of increasing importance in many fields reaching from astronomical observation of distant galaxies to turbulence phenomena in liquids or gases. Currently, coherent light beams with orbital angular momentum (OAM) are of particular interest for optical communication, metrology, micro-machining or particle manipulation. One common task is to identify characteristic spiral patterns in pixelated intensity maps at real-world signal-to-noise ratios. Recently, Fast Fourier Transform (FFT) was applied to the quantitative analysis of spiral interference patterns of OAM beams. Our study considers this approach in more detail. It is demonstrated that specific information can be extracted by combining conformal mapping by polar transform with determining one- or two-dimensional angular frequency spectra either in spatially fixed or scanning mode. Results of proof-of-principle experiments are presented.

We conducted a pilot study to analyze the differential methylation status of 20 primary acinar adenocarcinomas of the lungs. These adenocarcinomas had to be wild type in mutation analysis and had either high (TPS&gt;50%; n=10) or low (TPS&lt;1%; n=10) PD-L1 status to be integrated into our study. To examine the methylation of 866,895 specific sites, we utilized the Illumina Infinium EPIC bead chip array. Both hypermethylation and hypomethylation play significant roles in tumor development, progression, and metastasis. They also impact the formation of the tumor microenvironment, which plays a decisive role in tumor differentiation, epigenetics, dissemination, and immune evasion. The gained methylation patterns were correlated with PD-L1 expression. Our analysis has identified distinct methylation patterns in lung adenocarcinomas with high and low PD-L1 expression. After analyzing the correlation between methylation results of genes and promoters with their pathobiology, we found that tumors with high expression of PD-L1 tend to exhibit oncogenic effects through hypermethylation. On the other hand, tumors with low PD-L1 expression show loss of their suppressor functions through hypomethylation. The suppressor functions of hypermethylated genes and promoters are ineffective compared to simultaneously activated dominant oncogenic mechanisms. The tumor microenvironment supports tumor growth in both groups.

Tomato (Solanum lycopersicum) is one of the most cultivated vegetables in the world due to its consumption in a large variety of raw, cooked, or processed foods. Tomato breeding and productivity highly depend on the use of hybrid seeds and their higher yield, environmental adaption, and disease tolerance. However, the emasculation procedure during hybridization raises tomato seed production costs and labor expenses. Using male sterility is an effective way to reduce the cost of hybrid seeds and ensure cultivar purity. Recent developments in CRISPR genome editing technology enabled tomato breeders to investigate the male sterility genes and to develop male-sterile tomato lines. In the current study, the tomato Acotinase (SlACO) gene family was investigated via in-silico tools and functionally characterized with CRISPR/Cas9-mediated gene disruption. Genome-wide blast and HMM search represented two SlACO genes located on different tomato chromosomes. Both genes were estimated to have a segmental duplication in the tomato genome due to their identical motif and domain structure. One of these genes, SlACO2, showed a high expression profile in all generative cells of tomato. Therefore, the SlACO2 gene was targeted with two different gRNA/Cas9 construct to identify their functional role in tomato. The gene was mutated in a total of 6 genome-edited tomato lines, 2 of which were homozygous. Surprisingly, pollen viability was found to be extremely low in mutant plants compared to their wild-type (WT) counterparts. Likewise, the number of seeds per fruit also sharply decreased more than fivefold in mutant lines (10-12 seed) compared to that in WT (67 seed). The pollen shape, anther structures, and flower colors/shapes were not significantly varied between the mutant and WT tomatoes. The mutated lines were also subjected to salt and mannitol-mediated drought stress to test the effect of SlACO2 on abiotic stress tolerance. The results of the study indicated that mutant tomatoes have higher tolerance with significantly lower MDA content under stress conditions. This is the first CRISPR-mediated characterization of ACO genes on pollen viability, seed formation, and abiotic stress tolerance in tomato.

To further understand the role of phosphorylation in ISGF3- and STAT2/IRF9-mediated constitutive and long-term IFN-I-stimulated transcriptional responses, we performed RNA-Seq and ChIP-Seq, in combination with phosphorylation inhibition and anti-viral experiments. First, we identified a group of ISRE-containing ISGs that were commonly regulated in IFNα treated WT and STAT1-KO cells. Thus, in 2fTGH and Huh7.5 WT cells IFNα-inducible transcription and anti-viral activity relied on the recruitment of the ISGF3 components STAT1, STAT2 and IRF9 in a phosphorylation- and time-dependent manner. Likewise, in ST2-U3C and Huh-STAT1KO cells lacking STAT1, ISG expression correlated with DNA-binding of phosphorylated STAT2/IRF9. This pointed to a dominant role of classical ISGF3 and STAT2/IRF9, and not U-ISGF3 or U-STAT2/IRF9, in the regulation of early and prolonged ISG expression and viral protection, in WT and STAT1-KO cells. In addition, comparative experiments in U3C (STAT1-KO) cells overexpressing all ISGF3 components (ST1-ST2-IRF9-U3C), revealed a threshold-dependent role of U-ISFG3, and potentially U-STAT2/IRF9, in the regulation of constitutive and possibly long-term IFNα-treated ISG expression and anti-viral activity.

Short-range Internet of Things (IoT) sensor nodes operating at 2.4 GHz need to provide ubiquitous wireless sensor networks (WSNs) with energy-efficient, wide-range output power (POUT) and fully integrated on a single chip for wireless body area networks (WBANs) and wireless personal area networks (WPANs) using low-power Bluetooth (BLE) and Zigbee standards. The proposed fully integrated transmitter (TX) utilizes a digitally controllable current-mode class-D (CMCD) power amplifier (PA) with a second harmonic distortion (HD2) suppression to reduce VCO pulling in an integrated system while meeting harmonic limit regulations. The CMCD PA is divided into 7-bit slices that can be reconfigured between differential and single-ended topologies. Duty cycle distortion compensation is performed for HD2 suppression, and an HD2 rejection filter and a modified C-L-C low-pass filter (LPF) reduce HD2 further. Implemented in a 28 nm CMOS process, the TX achieves a wide POUT range of from 12.1 to -31 dBm and provides a maximum efficiency of 39.8% while consuming 41.1 mW at 12.1 dBm POUT. The calibrated HD2 level is -82.2 dBc at 9.93 dBm POUT, resulting in a transmitter figure of merit (TX_FoM) of -97.52 dB. Higher order harmonic levels remain below -41.2 dBm even at 12.1 dBm POUT, meeting regulatory requirements.

The emotion of disgust protects individuals against pathogens, and it has been found that during pregnancy, it is elevated. Physiological mechanisms discussed in relation to these changes include immune markers and progesterone levels. The aim of this study is to assess the association between steroids and disgust in pregnancy. In a sample of 179 pregnant women, we have analyzed blood serum steroid concentrations and measured disgust sensitivity by text-based questionnaires in the first and third trimester. We found positive correlations between disgust sensitivity and the levels of C19 steroids, including testosterone, but also its precursors in the Δ5 pathway (androstenediol, DHEA, and their sulfates) and the Δ4 pathway (androstenedione). Furthermore, positive correlations between disgust sensitivity and 5α/β-reduced C19 steroid metabolites were consistently present in both trimesters. In the first trimester, we found a positive association between disgust sensitivity and 17-hydroxypregnanolone as well as some estrogens. In the third trimester, we have observed positive associations between disgust sensitivity and cortisol as well as immunoprotective Δ5 C19 7α/β-hydroxy-steroids. Our study has thus shown that disgust sensitivity is positively correlated with immunomodulatory steroids, and in the third trimester, steroids may reflect potential maternal anxiety-related symptoms.

This literature review paper delves into the synergies between Climate-Smart Agriculture (CSA) and the effective management of food losses. With the escalating challenges posed by climate change and the pressing need to address food security concerns, the intersection of these two critical domains becomes paramount for sustainable agricultural practices. The review explores existing research on the implementation of climate-smart agricultural techniques and their impact on mitigating food losses throughout the agricultural value chain. The paper synthesizes findings related to climate-resilient crop varieties, precision farming, water-use efficiency, and sustainable soil management, highlighting their contributions to reducing post-harvest losses and enhancing overall food security. In addition to summarizing current knowledge, this paper identifies gaps in the literature and proposes avenues for future research. This review serves as a foundation for researchers, policymakers, and practitioners seeking to contribute to the advancement of climate-smart agriculture and the reduction of food losses, paving the way for a more sustainable and resilient global food system.

We developed and applied a method for characterizing bacterial promoters genome-wide by in vitro transcription coupled to transcriptome sequencing specific for native 5’-ends of transcripts. This method called ROSE (Run-Off transcription/RNA-SEquencing), only requires chromosomal DNA, ribonucleotides, RNA polymerase (RNAP) core enzyme, and a specific sigma factor, recognizing the corresponding promoters, which have to be analyzed. ROSE was performed on E. coli K-12 MG1655 genomic DNA using E. coli RNAP holoenzyme (including σ70) and yielded 3,226 transcription start sites, 2,167 of which were also identified in in vivo studies, and 598 were new. Many new promoters not yet identified by in vivo experiments might be repressed under the tested conditions. Complementary in vivo experiments with E. coli K-12 strain BW25113 and isogenic transcription factor gene knockout mutants of fis, fur, and hns were used to test this hypothesis. Comparative transcriptome analysis demonstrated that ROSE could identify bona fide promoters that were apparently repressed in vivo. In this sense, ROSE is well-suited as a bottom-up approach for characterizing transcriptional networks in bacteria and ideally complementary to top-down in vivo transcriptome studies.

Viroids are smallest pathogen that consist of non-capsidated, single-stranded non-coding RNA replicons and exploits host factors for their replication and propagation. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) is a serious threat, which spread rapidly within hop gardens. In this study, we employed comprehensive transcriptome analyses to dissect host-viroid interactions and identify gene expression changes associated with disease development in hop. Our analysis revealed that CBCVd-infection resulted in the massive modulation of activity of over 2000 genes. Expression of genes associated with plant immune responses (protein kinase and mitogen-activated protein kinase), hypersensitive responses, phytohormone signaling pathways, photosynthesis, pigment metabolism, protein metabolism, sugar metabolism and modification and others were altered, which could be attributed to systemic symptom development upon CBCVd-infection in hop. In addition, genes encoding RNA-dependent RNA polymerase, pathogenesis-related protein, chitinase as well as those related to basal defense responses were up-regulated. The expression levels of several genes identified from RNA sequencing analysis were confirmed by qRT-PCR. Our systematic comprehensive CBCVd-responsive transcriptome analysis provides a better understanding and insights into complex viroid-hop plant interaction. This information will assist further in the development of future measures for the prevention of CBCVd spread in hop fields.

Universal Stress Protein (USP) plays an important regulatory role in the response to abiotic stress. Most of the research related to USP so far has been conducted on plant models such as Arabidopsis, rice, and cotton. The potato is one of the four major food crops in the world. Therefore, it is of great significance to study the expression pattern of the potato StUSP gene family under abiotic stress conditions. In this study, a total of 108 USP genes were identified from the genome of the Atlantic potato, divided into four subgroups. And based on the genetic structure, the physical and chemical properties of proteins and other aspects of their biological characteristics are comprehensively analyzed. Collinear analysis showed that the homologous genes of StUSPs and four other representative species (Arabidopsis, rice, coyote tobacco, and tomato) were highly conserved. The cis-regulatory elements of the StUSPs promoter are involved in plant hormones, environmental stress, mechanical damage, and light response. RNA-seq analysis showed that under different abiotic stresses, there are differences in the expression patterns of members of each subgroup. WGCNA co-expression analysis of the central gene showed that the differential co-expression gene is mainly involved in the plant-pathogen response process, plant hormone signal transduction, and the biosynthesis process of secondary metabolites. Through qRT-PCR analysis, it was confirmed that StUSP13, StUSP14, StUSP15, and StUSP41 may be important candidate genes involved in the response to adversity stress in potatoes. The results of this study provide the basis for further research on the functional analysis of StUSPs in the potato response to adversity stress.

The SnRK gene family is a key regulator playing an important role in plant stress response by phosphorylating the target protein to regulate the signalling pathways. The function of SnRK gene family has been reported in many species but is limited to Triticum asetivum. In this study, SnRK gene family in the wheat genome was identified and its structural characteristics were described. One hundred forty-seven SnRK genes distributed across 21 chromosomes were identified in the Triticum aestivum genome and categorised into three subgroups (SnRK1/2/3) based on phylogenetic analyses and domain types. The gene intron-exon structure and protein-motif composition of SnRKs were similar within each subgroup but different amongst the groups. Gene duplication between the wheat, Arabidopsis, rice and barley genomes was also investigated in order to get insight into the evolutionary aspects of the TaSnRK family genes. SnRK genes showed differential expression patterns in leaves, roots, spike, and grains. Redundant stress-related cis-elements were also found in the promoters of 129 SnRK genes and their expression levels varied widely following drought, ABA and light regulated elements. In particular, TaSnRK2.11 had higher and increased expression under the abiotic stresses and can be a candidate gene for the abiotc stress tolerance. The findings will aid in the functional characterization of TaSnRK genes for further research.

Low temperature is a common biological abiotic stress in major cotton growing areas. Cold stress significantly affects the growth, yield and yield quality of cotton. Therefore, it is important to develop a more robust and cold stress tolerant cotton germplasms. Climate change and erratic climatic condition, plants have evolved various survival mechanisms, one of which induction of various stress responsive transcriptome factors, such as the C-repeat binding factor GthCBF4, which have been found to enhance cold tolerance in various plants. In this study detailed evaluation of the cotton C-repeat binding factor has been carried out. A total of29, 28, 25, 21, 30, 26 and 15 proteins encoded by the C-repeat binding factor were identified in G. herbaceum, G. arboreum, G. thurberi, G. raimondii, G. turneri, G. longicalyx and G. australe, respectively. Phylogeny evaluation revealed that the proteins were grouped into seven clades, with clade 1 and 6 being the largest. Moreover, majority of the proteins encoded by the genes were predicted to be located within. the nucleus, while some are distributed in other parts of the cell. Based on the transcriptome and RT-qPCR analysis, Gthu17439 (GthCBF4) was highly upregulated and was further validated through forward genetics. The Gthu17439 (GthCBF4) overexpressed plants showed a significantly tolerance to cold stress, with higher growth vigour compared to the wild types. The results showed that the Gthu17439 (GthCBF4) could be playing a significant role in enhancing cold stress tolerance in cotton and can be further exploited in developing a more cold stress tolerance cotton germplasm

This study aimed to determine shifts in microbial populations regarding richness and diversity from the daily use of a popular over-the-counter nasal spray. In addition, the finding of nasal commensal bacterial species that overlap with the oral microbiome may prove to be potential probiotics for the “gateway microbiomes.” Nasal swab samples were obtained before and after using the most popular over-the-counter (OTC) nasal spray in 10 participants aged 18–48. All participants were healthy volunteers with no significant medical histories. The participants were randomly assigned a number by randomizing software and consisted of five men and five women. The sampling consisted of placing a nasal swab atraumatically into the nasal cavity. The samples were preserved and sent to Northwestern University Sequencing Center for whole-genome deep sequencing. After 21 days of OTC nasal spray use twice daily, the participants returned for further nasal microbiome sampling. The microbial analysis included all bacteria, archaea, viruses, molds, and yeasts via deep sequencing for species analysis. The Northwestern University Sequencing Center utilized artificial intelligence analysis to determine shifts in species and strains following nasal spray use that resulted in changes in diversity and richness.

High hopes were associated with the new anti-angiogenic medicament, bevacizumab, which due to the possibility of binding to isoform A of vascular endothelial growth factor (VEGF), inhibited the creation of new blood vessels. However, not only blood vessels are responsible for the tumor cells' spread. During the process of tumor growth, lymphangiogenesis is mediated by other isoforms of the VEGF family, specifically VEGF-C and VEGF-D, which are independent of the bevacizumab action. Therefore, based on the mechanism of bevacizumab action and the processes of angio- and lymphangiogenesis, we formed three hypotheses: 1) if the lymph nodes in the primary ovarian cancers are metastatic, the outcome of bevacizumab treatment is worsened; 2) concerning the second-line treatment, bevacizumab will act in a weakened manner if the recurrence occurs in lymph nodes compared to local recurrences; 3) patients treated by bevacizumab are more likely to have recurrences in lymph nodes. These hypotheses raise the issue of the existing gap of knowledge, which concerns the effect of bevacizumab on the metastatic lymph nodes.

Global coverage and scale up of interventions to reduce newborn mortality remains low, though progress has been achieved in improving newborn survival in many low-income settings. An important factor in the success of newborn health interventions, and moving to scale, is appropriate design of community-based programs and strategies for local implementation. We report the results of formative research undertaken to inform the design of a newborn health intervention in Cambodia. Information was gathered on newborn care practices over a period of three months using multiple qualitative methods of data collection in the primary health facility and home setting. Analysis of the data indicated important gaps, both at home and facility level, between recommended newborn care practices and those typical in the study area. The results of this formative research have informed strategies for behavior change and improving referral of sick infants in the subsequent implementation study. Collection and dissemination of data on newborn care practices from settings such as these can contribute to efforts to advance survival, growth and development of newborns for intervention research, and for future newborn health programming.

Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as reference panel to detect associations with the phenotype at the gene level and were gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as reference panel to conduct tissue specific TWAS on cigarettes per day (CPD) over 13 brain tissues in two large cohorts: UK Biobank (UKBB; N=142,202) and the GWAS &amp; Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N=143,210), and then meta-analyzed the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues, partially homogeneous genes associated with CPD in cortex, cerebellum and hippocampus tissues, and lastly the tissue-specific genes associated with CPD in only few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.

The Work Ability Index (WAI) is the most widely used questionnaire for the self-assessment of working ability. Because of its different applications, shorter versions, and widespread use in healthcare activities, it is worthwhile assessing its characteristics. The WAI was distributed online among the workers of a healthcare company; the results were compared with data contained in the workers' personal health records and with absence registers. 340 out of 575 workers (59.1%) participated; 6.5% of them reported poor work ability. Exploratory factor analysis indicated that the one-factor version best described the characteristics of the WAI. The scores of the complete WAI, the shorter form without the list of diseases, and the minimal one-item version (WAS) had equal distribution and were significantly correlated. The WAI score was inversely related to age and significantly lower in women than in men, but it was higher in night workers than in colleagues due to the probable effect of selective factors. The WAI score was also correlated with absenteeism, but no differences were found between males and females in the average number of absences, suggesting that cultural or emotional factors influence the self-rating of the WAI. Workers tended to overreport illnesses in the online survey compared to data collected during occupational health checks. Musculoskeletal disorders were the most frequently reported illnesses (53%). Psychiatric illnesses affected 21% of workers and had the greatest impact on work ability. Multilevel ergonomic and human factor intervention seems to be needed to recover the working capacity of healthcare workers.