Milletdb: a multi-omics database to accelerate the research of functional genomics and molecular breeding of millets.

Millets are a class of nutrient-rich coarse cereals with high resistance to abiotic stress; thus, they guarantee food security for people living in areas with extreme climatic conditions and provide stress-related genetic resources for other crops. However, no platform is available to provide a comprehensive and systematic multi-omics analysis for millets, which seriously hinders the mining of stress-related genes and the molecular breeding of millets. Here, a free, web-accessible, user-friendly millets multi-omics database platform (Milletdb, http://milletdb.novogene.com) has been developed. The Milletdb contains six millets and their one related species genomes, graph-based pan-genomics of pearl millet, and stress-related multi-omics data, which enable Milletdb to be the most complete millets multi-omics database available. We stored GWAS (genome-wide association study) results of 20 yield-related trait data obtained under three environmental conditions [field (no stress), early drought and late drought] for 2 years in the database, allowing users to identify stress-related genes that support yield improvement. Milletdb can simplify the functional genomics analysis of millets by providing users with 20 different tools (e.g., 'Gene mapping', 'Co-expression', 'KEGG/GO Enrichment' analysis, etc.). On the Milletdb platform, a gene PMA1G03779.1 was identified through 'GWAS', which has the potential to modulate yield and respond to different environmental stresses. Using the tools provided by Milletdb, we found that the stress-related PLATZs TFs (transcription factors) family expands in 87.5% of millet accessions and contributes to vegetative growth and abiotic stress responses. Milletdb can effectively serve researchers in the mining of key genes, genome editing and molecular breeding of millets.

Genomic analysis of Tibetan ground tits identifies molecular adaptations associated with cooperative breeding.

Cooperative breeding is a sophisticated altruistic social behavior that helps social animals to adapt to harsh environments. The Tibetan ground tit, Pseudopodoces humilis, is a high-altitude bird endemic to the Tibetan plateau. Recently, it has become an exciting system for studying the evolution of facultative cooperative breeding. To test for molecular adaptations associated with cooperative breeding, we resequenced the whole genome of ground tits from 6 wild populations that display remarkable variation in the frequency of cooperative breeding. Population structure analyses showed that the 6 populations were divided into 4 lineages, which is congruent with the major geographical distribution of the sampling sites. Using genome-wide selective sweep analysis, we identified putative positively selected genes (PSGs) in groups of tits that displayed high and low cooperative breeding rates. The total number of PSGs varied from 146 to 722 in high cooperative breeding rate populations, and from 272 to 752 in low cooperative breeding rate populations. Functional enrichment analysis of these PSGs identified several significantly enriched ontologies related to oxytocin signaling, estrogen signaling, and insulin secretion. PSGs involved in these functional ontologies suggest that molecular adaptations in hormonal regulation may have played important roles in shaping the evolution of cooperative breeding in the ground tit. Taken together, our study provides candidate genes and functional ontologies involved in molecular adaptations associated with cooperative breeding in Tibetan ground tits, and calls for a better understanding of the genetic roles in the evolution of cooperative breeding.

Comparative analyses of bat genomes identify distinct evolution of immunity in Old World fruit bats.

Bats have been identified as natural reservoir hosts of several zoonotic viruses, prompting suggestions that they have unique immunological adaptations. Among bats, Old World fruit bats (Pteropodidae) have been linked to multiple spillovers. To test for lineage-specific molecular adaptations in these bats, we developed a new assembly pipeline to generate a reference-quality genome of the fruit bat Cynopterus sphinx and used this in comparative analyses of 12 bat species, including six pteropodids. Our results reveal that immunity-related genes have higher evolutionary rates in pteropodids than in other bats. Several lineage-specific genetic changes were shared across pteropodids, including the loss of NLRP1, duplications of PGLYRP1 and C5AR2, and amino acid replacements in MyD88. We introduced MyD88 transgenes containing Pteropodidae-specific residues into bat and human cell lines and found evidence of dampened inflammatory responses. By uncovering distinct immune adaptations, our results could help explain why pteropodids are frequently identified as viral hosts.

Engineering sodium alginate-SiO2 composite beads for efficient removal of methylene blue from water.

The lack of sufficient active binding sites in commonly reported sodium alginate (SA)-based porous beads hampers their performances in adsorption of water contaminants. To address this problem, porous SA-SiO2 beads functionalized with poly(2-acrylamido-2-methylpropane sulfonic acid) (PAMPS) are reported in this work. Due to the porous properties and the existence of abundant sulfonate groups, the obtained composite material SA-SiO2-PAMPS shows excellent adsorption capacity toward cationic dye methylene blue (MB). The adsorption kinetic and adsorption isotherm studies reveal that the adsorption process fits closely to pseudo-second-order kinetic model and Langmuir isotherm model, respectively, suggesting the existence of chemical adsorption and monolayer adsorption behavior. The maximum adsorption capacity obtained from Langmuir model is found to be 427.36, 495.05, and 564.97 mg/g under 25, 35, and 45 °C, respectively. The calculated thermodynamic parameters indicate that MB adsorption on SA-SiO2-PAMPS is spontaneous and endothermic.

Post-synthetic modification (PSM) of MOFs with an ionic polymer for efficient adsorptive removal of methylene blue from water.

UiO-66-NH2 was functionalized with an ionic polymer poly(2-acrylamido-2-methylpropane sulfonic acid) (PAMPS) through a post-synthetic modification (PSM) strategy. Due to the excellent dispersibility in water and the existence of abundant active binding sites, the obtained UiO-66-PAMPS shows significantly improved adsorption capability toward methylene blue (MB) in aqueous solution.

Pangenomic analysis identifies structural variation associated with heat tolerance in pearl millet.

Pearl millet is an important cereal crop worldwide and shows superior heat tolerance. Here, we developed a graph-based pan-genome by assembling ten chromosomal genomes with one existing assembly adapted to different climates worldwide and captured 424,085 genomic structural variations (SVs). Comparative genomics and transcriptomics analyses revealed the expansion of the RWP-RK transcription factor family and the involvement of endoplasmic reticulum (ER)-related genes in heat tolerance. The overexpression of one RWP-RK gene led to enhanced plant heat tolerance and transactivated ER-related genes quickly, supporting the important roles of RWP-RK transcription factors and ER system in heat tolerance. Furthermore, we found that some SVs affected the gene expression associated with heat tolerance and SVs surrounding ER-related genes shaped adaptation to heat tolerance during domestication in the population. Our study provides a comprehensive genomic resource revealing insights into heat tolerance and laying a foundation for generating more robust crops under the changing climate.

Epirubicin may enhance the inhibition of hepatocellular carcinoma induced by iodine-125 seeds through downregulating WNT pathway.

AIM: To investigate the role of iodine-125 (125 I) combined with epirubicin (EPI) in inhibiting hepatocellular carcinoma (HCC) growth and promoting apoptosis. METHODS: Both in vivo and in vitro experiments were conducted. CCK-8 assay was performed to determine the cells viability after EPI treatment. HepG2 and SMMC7721 cells were treated with EPI or 125 I or in combination. Colony formation assays were performed to verify the antiproliferation effect. Annexin V-FITC/PI, terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assays, and western blotting were performed to analyze cellular apoptosis. Scratch wound healing assays and transwell assays were used to examine migration following different treatments. An isobaric tag for relative and absolute quantitation analysis was used to detect changes in protein expression after 125 I treatment, identifying the potential mediating protein cathelicidin (LL-37). LL-37 protein and WNT pathway-related proteins were detected by western blotting in SMMC7721 and HepG2 cells. Mice were treated with 125 I and EPI to evaluate whether EPI enhanced the antitumor effect of 125 I. RESULTS: EPI promoted 125 I-induced apoptosis and reduced the proliferation of HepG2 and SMMC7721 cells. EPI also prevented the migration of HepG2 and SMMC7721 cells. EPI combined with 125 I may interfere with the WNT signaling pathway by decreasing LL-37 to inhibit HCC development. The antitumor effects of EPI with 125 I were verified in mice. CONCLUSION: EPI combined with 125 I induced apoptosis and inhibited the proliferation of HCC cells by LL-37 downregulating the WNT pathway.

Genome-wide signatures of the geographic expansion and breeding of soybean.

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

Transarterial chemoembolization combined with metformin improves the prognosis of hepatocellular carcinoma patients with type 2 diabetes.

Objective: The study aims to investigate the effect of metformin on Hepatocellular carcinoma (HCC) patients with type 2 diabetes mellitus (T2DM) who received transarterial chemoembolization (TACE) for the first time. Methods: From January 2016 to December 2019, T2DM patients diagnosed with HCC in Shandong Cancer Hospital and treated with TACE were included in this retrospective study. Overall survival (OS) and Progression-free survival (PFS) were compared between patients treated with metformin and other antidiabetics. Univariate and multivariate Cox regression models were used to evaluate the independent risk factors associated with OS and PFS. And sub-analysis was performed to investigate whether metformin could give a survival advantage in each Barcelona Clinic Liver Cancer (BCLC) stage of HCC. Propensity score matched (PSM) analyses based on patient and tumor characteristics were also conducted. Results: A total of 123 HCC patients with T2DM underwent TACE, of which 50 (40.65%) received treatment with metformin. For the whole cohort, the median OS (42 vs 32 months, p=0.054) and PFS (12 vs 7 months, P=0.0016) were longer in the metformin group than that in the non-metformin group. Multi-analysis revealed that BCLC stage, BMI (Body Mass Index), and metformin use were independent predictors of OS. Metformin use was independently associated with recurrence. After PSM, 39 matched pairs were identified. The use of metformin was associated with a numerically longer m OS (43 vs 35 months, P=0.183) than the use of other anti-diabetics. And the difference in median PFS (13 vs 7 months, p=0.018) between the metformin group and non-metformin group remained significant. Conclusion: The combination of transarterial chemoembolization and metformin may be associated with better OS and PFS in HCC patients with T2DM.

Genomic insights into the evolutionary history and diversification of bulb traits in garlic.

BACKGROUND: Garlic is an entirely sterile crop with important value as a vegetable, condiment, and medicine. However, the evolutionary history of garlic remains largely unknown. RESULTS: Here we report a comprehensive map of garlic genomic variation, consisting of amazingly 129.4 million variations. Evolutionary analysis indicates that the garlic population diverged at least 100,000 years ago, and the two groups cultivated in China were domesticated from two independent routes. Consequently, 15.0 and 17.5% of genes underwent an expression change in two cultivated groups, causing a reshaping of their transcriptomic architecture. Furthermore, we find independent domestication leads to few overlaps of deleterious substitutions in these two groups due to separate accumulation and selection-based removal. By analysis of selective sweeps, genome-wide trait associations and associated transcriptomic analysis, we uncover differential selections for the bulb traits in these two garlic groups during their domestication. CONCLUSIONS: This study provides valuable resources for garlic genomics-based breeding, and comprehensive insights into the evolutionary history of this clonal-propagated crop.

Multi-omics provides new insights into the domestication and improvement of dark jute (Corchorus olitorius).

Jute (Corchorus sp.) is the most important bast fiber crop worldwide; however, the mechanisms underlying domestication and improvement remain largely unknown. We performed multi-omics analysis by integrating de novo sequencing, resequencing, and transcriptomic and epigenetic sequencing to clarify the domestication and improvement of dark jute Corchorus olitorius. We demonstrated that dark jute underwent early domestication and a relatively moderate genetic bottleneck during improvement breeding. A genome-wide association study of 11 important agronomic traits identified abundant candidate loci. We characterized the selective sweeps in the two breeding stages of jute, prominently, soil salinity differences played an important role in environmental adaptation during domestication, and the strongly selected genes for improvement had an increased frequency of favorable haplotypes. Furthermore, we speculated that an encoding auxin/indole-3-acetic acid protein COS07g_00652 could enhance the flexibility and strength of the stem to improve fiber yield. Our study not only provides valuable genetic resources for future fiber breeding in jute, but also is of great significance for reviewing the genetic basis of early crop breeding.

Effects of Chlorella extracts on growth of Capsicum annuum L. seedlings.

The long-term application of chemical fertilizers has caused to the farmland soil compaction, water pollution, and reduced the quality of vegetable to some extent. So, its become a trend in agriculture to find new bio-fertilizers. Chlorella extract is rich in amino acids, peptides, nucleic acids, growth hormones, potassium, calcium, magnesium, iron, zinc ions, vitamin E, B1, B2, C, B6, folic acid, free biotin and chlorophyll. Chlorella extract can promote biological growth, mainly by stimulating the speed of cell division, thereby accelerating the proliferation rate of cells and playing a role in promoting plant growth. Whether Chlorella extract can be used to improve the growth of pepper (Capsicum annuum), needs to be verified. In current study, a pepper variety 'Chao Tian Jiao' was used as experiment material, by determining the changes of the related characteristics after spraying the seedlings with Chlorella extract, and its effect on growth of Capsicum annuum plants was investigated. The results showed that the Chlorella extract significantly increased plant height of pepper seedlings (treatment: 32.2 ± 0.3 cm; control: 24.2 ± 0.2 cm), stem diameter (treatment: 0.57 ± 0.02 cm; control: 0.41 ± 0.03 cm) and leaf area (treatment: 189.6 ± 3.2 cm2; control: 145.8 ± 2.5 cm2). Particularly, the pepper seedlings treated with Chlorella extract, developed the root system in better way, significantly increased the chlorophyll a, and the activities of SOD, POD and CAT enzymes were also improved significantly. Based on our results, we can speculate that it is possible to improve the growth of Capsicum annuum seedlings and reduce the application of chemical fertilizers in pepper production by using Chlorella extract.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Comparative 3D genome architecture in vertebrates.

BACKGROUND: The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. RESULTS: This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. CONCLUSIONS: Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding.

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

The role of bile acid subtypes in the diagnosis of cholangiocarcinoma.

OBJECTIVES: To find the potential biomarkers of cholangiocarcinoma, form a biomarker package, evaluate its efficiency, and validate it. METHODS: R software was used to analyze the differential expression of mRNAs between cholangiocarcinoma and adjacent nontumorous tissues, obtained from The Cancer Genome Atlas (TCGA), and enrich the KEGG pathway. Metabo-Profile Inc. performed the comprehensive bile acid profiling and quantitation. The training set concluded 20 cholangiocarcinoma and 20 nontumorous volunteers. Receiver operating characteristic (ROC) curve and accompanying area under the curve (AUC) was calculated. The top four bile acids formed a new biomarker package. The validation set included 15 cholangiocarcinoma and 15 nontumorous, and the sensitivity and specificity of the new biomarker package were tested. RESULTS: Gene expression of 36 cholangiocarcinoma and nine adjacent nontumorous tissues was obtained in January 2020. Totally 9887 differential genes were eligible (logFC ≥ 1 or ≤ -1, P < 0.05, and adjust P < 0.01). GO analysis showed that 20 KEGG pathways were enriched, including primary bile acid biosynthesis and bile secretion. Comprehensive bile acid profiling and quantitation showed 15 differential bile acid types, and the ROC-AUC was between 0.953 and 0.750. HDCA, isoLCA, bCDCA, and DCA were selected to form a biomarker package. The Logit (p = cholangiocarcinoma) = 7.898 - 3.70*(1isoLCA) - 0.444*(bCDCA) + 0.415*(HDCA) + 0.041*DCA. Its ROC-AUC was 0.944. In the validation set, the sensitivity was 0.933 and the specificity was 0.867. CONCLUSION: Bile acid types package was efficient to distinguish nontumorous population and cholangiocarcinoma. The difference might be associated to the downregulation of primary bile acid biosynthesis and bile secretion pathway of cholangiocarcinoma.

Association of female reproductive tract microbiota with egg production in layer chickens.

BACKGROUND: The microbiota of the female reproductive tract is increasingly recognized as playing fundamental roles in animal reproduction. To explore the relative contribution of reproductive tract microbiomes to egg production in chickens, we investigated the microbiota in multiple reproductive and digestive tract sites from 128 female layer (egg-producing) chickens in comparable environments. RESULTS: We identified substantial differences between the diversity, composition, and predicted function of site-associated microbiota. Differences in reproductive tract microbiota were more strongly associated with egg production than those in the digestive tract. We identified 4 reproductive tract microbial species, Bacteroides fragilis, Bacteroides salanitronis, Bacteroides barnesiae, and Clostridium leptum, that were related to immune function and potentially contribute to enhanced egg production. CONCLUSIONS: These findings provide insights into the diverse microbiota characteristics of reproductive and digestive tracts and may help in designing strategies for controlling and manipulating chicken reproductive tract microbiota to improve egg production.

High-quality genome assembly and resequencing of modern cotton cultivars provide resources for crop improvement.

Cotton produces natural fiber for the textile industry. The genetic effects of genomic structural variations underlying agronomic traits remain unclear. Here, we generate two high-quality genomes of Gossypium hirsutum cv. NDM8 and Gossypium barbadense acc. Pima90, and identify large-scale structural variations in the two species and 1,081 G. hirsutum accessions. The density of structural variations is higher in the D-subgenome than in the A-subgenome, indicating that the D-subgenome undergoes stronger selection during species formation and variety development. Many structural variations in genes and/or regulatory regions potentially influencing agronomic traits were discovered. Of 446 significantly associated structural variations, those for fiber quality and Verticillium wilt resistance are located mainly in the D-subgenome and those for yield mainly in the A-subgenome. Our research provides insight into the role of structural variations in genotype-to-phenotype relationships and their potential utility in crop improvement.

Building a Chinese pan-genome of 486 individuals.

Pan-genome sequence analysis of human population ancestry is critical for expanding and better defining human genome sequence diversity. However, the amount of genetic variation still missing from current human reference sequences is still unknown. Here, we used 486 deep-sequenced Han Chinese genomes to identify 276 Mbp of DNA sequences that, to our knowledge, are absent in the current human reference. We classified these sequences into individual-specific and common sequences, and propose that the common sequence size is uncapped with a growing population. The 46.646 Mbp common sequences obtained from the 486 individuals improved the accuracy of variant calling and mapping rate when added to the reference genome. We also analyzed the genomic positions of these common sequences and found that they came from genomic regions characterized by high mutation rate and low pathogenicity. Our study authenticates the Chinese pan-genome as representative of DNA sequences specific to the Han Chinese population missing from the GRCh38 reference genome and establishes the newly defined common sequences as candidates to supplement the current human reference.

Vulture Genomes Reveal Molecular Adaptations Underlying Obligate Scavenging and Low Levels of Genetic Diversity.

Obligate scavenging on the dead and decaying animal matter is a rare dietary specialization that in extant vertebrates is restricted to vultures. These birds perform essential ecological services, yet many vulture species have undergone recent steep population declines and are now endangered. To test for molecular adaptations underlying obligate scavenging in vultures, and to assess whether genomic features might have contributed to their population declines, we generated high-quality genomes of the Himalayan and bearded vultures, representing both independent origins of scavenging within the Accipitridae, alongside a sister taxon, the upland buzzard. By comparing our data to published sequences from other birds, we show that the evolution of obligate scavenging in vultures has been accompanied by widespread positive selection acting on genes underlying gastric acid production, and immunity. Moreover, we find evidence of parallel molecular evolution, with amino acid replacements shared among divergent lineages of these scavengers. Our genome-wide screens also reveal that both the Himalayan and bearded vultures exhibit low levels of genetic diversity, equating to around a half of the mean genetic diversity of other bird genomes examined. However, demographic reconstructions indicate that population declines began at around the Last Glacial Maximum, predating the well-documented dramatic declines of the past three decades. Taken together, our genomic analyses imply that vultures harbor unique adaptations for processing carrion, but that modern populations are genetically depauperate and thus especially vulnerable to further genetic erosion through anthropogenic activities.