Habitat Transition in the Evolution of Bacteria and Archaea.

Related groups of microbes are widely distributed across Earth's habitats, implying numerous dispersal and adaptation events over evolutionary time. However, relatively little is known about the characteristics and mechanisms of these habitat transitions, particularly for populations that reside in animal microbiomes. Here, we review the literature concerning habitat transitions among a variety of bacterial and archaeal lineages, considering the frequency of migration events, potential environmental barriers, and mechanisms of adaptation to new physicochemical conditions, including the modification of protein inventories and other genomic characteristics. Cells dependent on microbial hosts, particularly bacteria from the Candidate Phyla Radiation, have undergone repeated habitat transitions from environmental sources into animal microbiomes. We compare their trajectories to those of both free-living cells-including the Melainabacteria, Elusimicrobia, and methanogenic archaea-and cellular endosymbionts and bacteriophages, which have made similar transitions. We conclude by highlighting major related topics that may be worthy of future study.

A scoutRNA Is Required for Some Type V CRISPR-Cas Systems.

CRISPR-Cas12c/d proteins share limited homology with Cas12a and Cas9 bacterial CRISPR RNA (crRNA)-guided nucleases used widely for genome editing and DNA detection. However, Cas12c (C2c3)- and Cas12d (CasY)-catalyzed DNA cleavage and genome editing activities have not been directly observed. We show here that a short-complementarity untranslated RNA (scoutRNA), together with crRNA, is required for Cas12d-catalyzed DNA cutting. The scoutRNA differs in secondary structure from previously described tracrRNAs used by CRISPR-Cas9 and some Cas12 enzymes, and in Cas12d-containing systems, scoutRNA includes a conserved five-nucleotide sequence that is essential for activity. In addition to supporting crRNA-directed DNA recognition, biochemical and cell-based experiments establish scoutRNA as an essential cofactor for Cas12c-catalyzed pre-crRNA maturation. These results define scoutRNA as a third type of transcript encoded by a subset of CRISPR-Cas genomic loci and explain how Cas12c/d systems avoid requirements for host factors including ribonuclease III for bacterial RNA-mediated adaptive immunity.

Artificial Intelligence and Machine Learning for Lead-to-Candidate Decision-Making and Beyond.

The use of artificial intelligence (AI) and machine learning (ML) in pharmaceutical research and development has to date focused on research: target identification; docking-, fragment-, and motif-based generation of compound libraries; modeling of synthesis feasibility; rank-ordering likely hits according to structural and chemometric similarity to compounds having known activity and affinity to the target(s); optimizing a smaller library for synthesis and high-throughput screening; and combining evidence from screening to support hit-to-lead decisions. Applying AI/ML methods to lead optimization and lead-to-candidate (L2C) decision-making has shown slower progress, especially regarding predicting absorption, distribution, metabolism, excretion, and toxicology properties. The present review surveys reasons why this is so, reports progress that has occurred in recent years, and summarizes some of the issues that remain. Effective AI/ML tools to derisk L2C and later phases of development are important to accelerate the pharmaceutical development process, ameliorate escalating development costs, and achieve greater success rates.

A kinetic model for solving a combination optimization problem in ab-initio Cryo-EM 3D reconstruction.

Cryo-Electron Microscopy (cryo-EM) is a widely used and effective method for determining the three-dimensional (3D) structure of biological molecules. For ab-initio Cryo-EM 3D reconstruction using single particle analysis (SPA), estimating the projection direction of the projection image is a crucial step. However, the existing SPA methods based on common lines are sensitive to noise. The error in common line detection will lead to a poor estimation of the projection directions and thus may greatly affect the final reconstruction results. To improve the reconstruction results, multiple candidate common lines are estimated for each pair of projection images. The key problem then becomes a combination optimization problem of selecting consistent common lines from multiple candidates. To solve the problem efficiently, a physics-inspired method based on a kinetic model is proposed in this work. More specifically, hypothetical attractive forces between each pair of candidate common lines are used to calculate a hypothetical torque exerted on each projection image in the 3D reconstruction space, and the rotation under the hypothetical torque is used to optimize the projection direction estimation of the projection image. This way, the consistent common lines along with the projection directions can be found directly without enumeration of all the combinations of the multiple candidate common lines. Compared with the traditional methods, the proposed method is shown to be able to produce more accurate 3D reconstruction results from high noise projection images. Besides the practical value, the proposed method also serves as a good reference for solving similar combinatorial optimization problems.

Identification of a second glycoform of the clinically prevalent O1 antigen from Klebsiella pneumoniae.

Carbapenemase and extended ß-lactamase-producing Klebsiella pneumoniae isolates represent a major health threat, stimulating increasing interest in immunotherapeutic approaches for combating Klebsiella infections. Lipopolysaccharide O antigen polysaccharides offer viable targets for immunotherapeutic development, and several studies have described protection with O-specific antibodies in animal models of infection. O1 antigen is produced by almost half of clinical Klebsiella isolates. The O1 polysaccharide backbone structure is known, but monoclonal antibodies raised against the O1 antigen showed varying reactivity against different isolates that could not be explained by the known structure. Reinvestigation of the structure by NMR spectroscopy revealed the presence of the reported polysaccharide backbone (glycoform O1a), as well as a previously unknown O1b glycoform composed of the O1a backbone modified with a terminal pyruvate group. The activity of the responsible pyruvyltransferase (WbbZ) was confirmed by western immunoblotting and in vitro chemoenzymatic synthesis of the O1b terminus. Bioinformatic data indicate that almost all O1 isolates possess genes required to produce both glycoforms. We describe the presence of O1ab-biosynthesis genes in other bacterial species and report a functional O1 locus on a bacteriophage genome. Homologs of wbbZ are widespread in genetic loci for the assembly of unrelated glycostructures in bacteria and yeast. In K. pneumoniae, simultaneous production of both O1 glycoforms is enabled by the lack of specificity of the ABC transporter that exports the nascent glycan, and the data reported here provide mechanistic understanding of the capacity for evolution of antigenic diversity within an important class of biomolecules produced by many bacteria.

African swine fever virus I73R is a critical virulence-related gene: A potential target for attenuation.

African swine fever virus (ASFV) is a large, double-stranded DNA virus that causes a fatal disease in pigs, posing a threat to the global pig industry. Whereas some ASFV proteins have been found to play important roles in ASFV-host interaction, the functional roles of many proteins are still largely unknown. In this study, we identified I73R, an early viral gene in the replication cycle of ASFV, as a key virulence factor. Our findings demonstrate that pI73R suppresses the host innate immune response by broadly inhibiting the synthesis of host proteins, including antiviral proteins. Crystallization and structural characterization results suggest that pI73R is a nucleic-acid-binding protein containing a Zα domain. It localizes in the nucleus and inhibits host protein synthesis by suppressing the nuclear export of cellular messenger RNA (mRNAs). While pI73R promotes viral replication, the deletion of the gene showed that it is a nonessential gene for virus replication. In vivo safety and immunogenicity evaluation results demonstrate that the deletion mutant ASFV-GZΔI73R is completely nonpathogenic and provides effective protection to pigs against wild-type ASFV. These results reveal I73R as a virulence-related gene critical for ASFV pathogenesis and suggest that it is a potential target for virus attenuation. Accordingly, the deletion mutant ASFV-GZΔI73R can be a potent live-attenuated vaccine candidate.

PROSE: phenotype-specific network signatures from individual proteomic samples.

Proteomic studies characterize the protein composition of complex biological samples. Despite recent advancements in mass spectrometry instrumentation and computational tools, low proteome coverage and interpretability remains a challenge. To address this, we developed Proteome Support Vector Enrichment (PROSE), a fast, scalable and lightweight pipeline for scoring proteins based on orthogonal gene co-expression network matrices. PROSE utilizes simple protein lists as input, generating a standard enrichment score for all proteins, including undetected ones. In our benchmark with 7 other candidate prioritization techniques, PROSE shows high accuracy in missing protein prediction, with scores correlating strongly to corresponding gene expression data. As a further proof-of-concept, we applied PROSE to a reanalysis of the Cancer Cell Line Encyclopedia proteomics dataset, where it captures key phenotypic features, including gene dependency. We lastly demonstrated its applicability on a breast cancer clinical dataset, showing clustering by annotated molecular subtype and identification of putative drivers of triple-negative breast cancer. PROSE is available as a user-friendly Python module from https://github.com/bwbio/PROSE.

Early-stage use of hearing aids preserves auditory cortical structure in children with sensorineural hearing loss.

Hearing is critical to spoken language, cognitive, and social development. Little is known about how early auditory experiences impact the brain structure of children with bilateral sensorineural hearing loss. This study examined the influence of hearing aid use and residual hearing on the auditory cortex of children with severe to profound congenital sensorineural hearing loss. We evaluated cortical preservation in 103 young pediatric cochlear implant candidates (55 females and 48 males) by comparing their multivoxel pattern similarity of auditory cortical structure with that of 78 age-matched children with typical hearing. The results demonstrated that early-stage hearing aid use preserved the auditory cortex of children with bilateral congenital sensorineural hearing loss. Children with less residual hearing experienced a more pronounced advantage from hearing aid use. However, this beneficial effect gradually diminished after 17 months of hearing aid use. These findings support timely fitting of hearing aids in conjunction with early implantation to take advantage of neural preservation to maximize auditory and spoken language development.

Type IV pili trigger episymbiotic association of Saccharibacteria with its bacterial host.

Recent characterization of the obligate episymbiont Saccharibacteria (TM7) belonging to the candidate phyla radiation (CPR) has expanded the extent of microbial diversity. However, the episymbiotic lifestyle of TM7 is still underexploited due to the deficiency of cultivated representatives. Here, we describe gene-targeted TM7 cultivation guided by repurposing epicPCR (emulsion, paired isolation, and concatenation PCR) to capture in situ TM7‒host associations. Using this method, we obtained a novel Saccharibacteria isolate TM7i and its host Leucobacter aridicollis J1 from Cicadae Periostracum, the castoff shell of cicada. Genomic analyses and microscopic characterizations revealed that TM7i could bind to J1 through twitching-like motility mediated by type IV pili (T4P). We further showed that the inhibition of T4P extrusion suppressed the motility and host adherence of TM7i, resulting in its reduced growth. However, the inactivation of T4P had little effect on the growth of TM7i that had already adhered to J1, suggesting the essential role of T4P in host recognition by TM7i. By capturing CPR‒host association and elaborating the T4P-dependent episymbiotic association mechanism, our studies shed light on the distinct yet widespread lifestyle of CPR bacteria.

Plancitoxin-1 mediates extracellular trap evasion by the parasitic helminth Trichinella spiralis.

BACKGROUND: Trichinella spiralis (T. spiralis) is a parasitic helminth that causes a globally prevalent neglected zoonotic disease, and worms at different developmental stages (muscle larvae, adult worms, newborn larvae) induce immune attack at different infection sites, causing serious harm to host health. Several innate immune cells release extracellular traps (ETs) to entrap and kill most pathogens that invade the body. In response, some unicellular pathogens have evolved a strategy to escape capture by ETs through the secretion of nucleases, but few related studies have investigated multicellular helminths. RESULTS: In the present study, we observed that ETs from neutrophils capture adult worms of T. spiralis, while ETs from macrophages trap muscle larvae and newborn larvae, and ETs had a killing effect on parasites in vitro. To defend against this immune attack, T. spiralis secretes plancitoxin-1, a DNase II-like protein, to degrade ETs and escape capture, which is essential for the survival of T. spiralis in the host. CONCLUSIONS: In summary, these findings demonstrate that T. spiralis escapes ET-mediated capture by secreting deoxyribonuclease as a potential conserved immune evasion mechanism, and plancitoxin-1 could be used as a potential vaccine candidate.

The mrp-3 gene is involved in haem efflux and detoxification in a blood-feeding nematode.

BACKGROUND: Haem is essential but toxic for metazoan organisms. Auxotrophic nematodes can acquire sufficient haem from the environment or their hosts in the meanwhile eliminate or detoxify excessive haem through tightly controlled machinery. In previous work, we reported a role of the unique transporter protein HRG-1 in the haem acquisition and homeostasis of parasitic nematodes. However, little is known about the haem efflux and detoxification via ABC transporters, particularly the multiple drug resistance proteins (MRPs). RESULTS: Here, we further elucidate that a member of the mrp family (mrp-3) is involved in haem efflux and detoxification in a blood-feeding model gastrointestinal parasite, Haemonchus contortus. This gene is haem-responsive and dominantly expressed in the intestine and inner membrane of the hypodermis of this parasite. RNA interference of mrp-3 resulted in a disturbance of genes (e.g. hrg-1, hrg-2 and gst-1) that are known to be involved in haem homeostasis and an increased formation of haemozoin in the treated larvae and lethality in vitro, particularly when exposed to exogenous haem. Notably, the nuclear hormone receptor NHR-14 appears to be associated the regulation of mrp-3 expression for haem homeostasis and detoxification. Gene knockdown of nhr-14 and/or mrp-3 increases the sensitivity of treated larvae to exogenous haem and consequently a high death rate (> 80%). CONCLUSIONS: These findings demonstrate that MRP-3 and the associated molecules are essential for haematophagous nematodes, suggesting novel intervention targets for these pathogens in humans and animals.

Identification of candidate genes associated with peanut pod length by combined analysis of QTL-seq and RNA-seq.

Pod length (PL) is one of the major traits determining pod size and yield of peanut. Discovering the quantitative trait loci (QTL) and identifying candidate genes associated with PL are essential for breeding high-yield peanut. In this study, quantitative trait loci sequencing (QTL-seq) was performed using the F2 population constructed by a short-pod variety Tifrunner (Tif) and a long-pod line Lps, and a 0.77 Mb genomic region on chromosome 07 was identified as the candidate region for PL. Then, the candidate region was narrowed to a 265.93 kb region by traditional QTL approach. RNA-seq analysis showed that there were four differentially expressed genes (DEGs) in the candidate region, among which Arahy.PF2L6F (AhCDC48) and Arahy.P4LK2T (AhTAA1) were speculated to be PL-related candidate genes. These results were informative for the elucidation of the underlying regulatory mechanism in peanut pod length and would facilitate further identification of valuable target genes.

Genome-wide association analysis reveals the genetic basis of thermal tolerance in dwarf surf clam Mulinia lateralis.

Recently, elevated seawater temperatures have resulted numerous adverse effects, including significant mortality among bivalves. The dwarf surf clam, Mulinia lateralis, is considered a valuable model species for bivalve research due to its rapid growth and short generation time. The successful cultivation in laboratory setting throughout its entire life cycle makes it an ideal candidate for exploring the potential mechanisms underlying bivalve responses to thermal stress. In this study, a total of 600 clams were subjected to a 17-day thermal stress experiment at a temperature of 30 °C which is the semi-lethal temperature for this species. Ninety individuals who perished initially were classified as heat-sensitive populations (HSP), while 89 individuals who survived the experiment were classified as heat-tolerant populations (HTP). Subsequently, 179 individuals were then sequenced, and 21,292 single nucleotide polymorphisms (SNPs) were genotyped for downstream analysis. The heritability estimate for survival status was found to be 0.375 ± 0.127 suggesting a genetic basis for thermal tolerance trait. Furthermore, a genome-wide association study (GWAS) identified three SNPs and 10 candidate genes associated with thermal tolerance trait in M. lateralis. These candidate genes were involved in the ETHR/EHF signaling pathway and played pivotal role in signal sensory, cell adhesion, oxidative stress, DNA damage repair, etc. Additionally, qPCR results indicated that, excluding MGAT4A, ZAN, and RFC1 genes, all others exhibited significantly higher expression in the HTP (p < 0.05), underscoring the critical involvement of the ETHR/EHF signaling pathway in M. lateralis' thermal tolerance. These results unveil the presence of standing genetic variations associated with thermal tolerance in M. lateralis, highlighting the regulatory role of the ETHR/EHF signaling pathway in the bivalve's response to thermal stress, which contribute to comprehension of the genetic basis of thermal tolerance in bivalves.

Systems genetics and bioinformatics analyses using ESR1-correlated genes identify potential candidates underlying female bone development.

Estrogen receptor α (ESR1) is involved in E2 signaling and plays a major role in postmenopausal bone loss. However, the molecular network underlying ESR1 has not been explored. We used systems genetics and bioinformatics to identify important genes associated with Esr1 in postmenopausal bone loss. We identified ~2300 Esr1-coexpressed genes in female BXD bone femur, functional analysis of which revealed 'osteoblast signaling' as the most enriched pathway. PPI network led to the identification of 25 'female bone candidates'. The gene-regulatory analysis revealed RUNX2 as a key TF. ANKRD1 and RUNX2 were significantly different between osteoporosis patients and healthy controls. Sp7, Col1a1 and Pth1r correlated with multiple femur bone phenotypes in BXD mice. miR-3121-3p targeted Csf1, Ankrd1, Sp7 and Runx2. ß-estradiol treatment markedly increased the expression of these candidates in mouse osteoblast. Our study revealed that Esr1-correlated genes Ankrd1, Runx2, Csf1 and Sp7 may play important roles in female bone development.

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility.

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences.

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs. IMPORTANCE: Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Epitope screening and vaccine molecule design of PRRSV GP3 and GP5 protein based on immunoinformatics.

Porcine reproductive and respiratory syndrome (PRRS) is a respiratory disease in pigs that causes severe economic losses. Currently, live PRRSV vaccines are commonly used but fail to prevent PRRS outbreaks and reinfection. Inactivated PRRSV vaccines have poor immunogenicity, making PRRSV a significant threat to swine health globally. Therefore, there is an urgent need to develop an effective PRRSV vaccine. This study used immunoinformatics to predict, screen, design and construct a candidate vaccine that fused B-cell epitopes, CTL- and HTL-dominant protective epitopes of PRRSV strain's GP3 and GP5 proteins. The study identified 12 B-cell epitopes, 6 CTL epitopes and 5 HTL epitopes of GP3 and GP5 proteins. The candidate vaccine was constructed with 50S ribosomal protein L7/L1 molecular adjuvant, which has antigenicity, solubility, stability, non-allergenicity and a high affinity for its target receptor, TLR-3. The C-ImmSim immunostimulation results showed significant increases in cellular and humoral responses (B cells and T cells) and production of TGF-ß, IL-2, IL-10, IFN-γ and IL-12. The constructed vaccine was stable and immunogenic, and it can effectively induce strong T-cell and B-cell immune responses against PRRSV. Therefore, it is a promising candidate vaccine for controlling and preventing PRRSV outbreaks.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection.

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

Immunogenicity, Reactogenicity, and Safety of AS01E-adjuvanted RSV Prefusion F Protein-based Candidate Vaccine (RSVPreF3 OA) When Co-administered With a Seasonal Quadrivalent Influenza Vaccine in Older Adults: Results of a Phase 3, Open-Label, Randomized Controlled Trial.

BACKGROUND: Co-administration of vaccines against respiratory syncytial virus (RSV) and influenza can be considered given their overlapping seasonality, and may increase vaccine uptake and compliance. In this phase 3, open-label, randomized study, we evaluated the immunogenicity, reactogenicity, and safety of the AS01E-adjuvanted RSV prefusion F protein-based candidate vaccine (RSVPreF3 OA) when co-administered with a seasonal quadrivalent influenza vaccine (FLU-QIV) in older adults. METHODS: Participants aged ≥60 years (randomized 1:1) received either RSVPreF3 OA and FLU-QIV simultaneously on day 1 (Co-Ad group) or FLU-QIV on day 1 followed by RSVPreF3 OA on day 31 (sequential administration [SA] group). The co-primary objectives were to demonstrate noninferiority of RSVPreF3 OA in terms of RSV-A neutralization geometric mean titer (GMT) ratio and FLU-QIV in terms of hemagglutination inhibition GMT ratio for each FLU-QIV strain, when co-administered versus when administered alone at 1-month post-vaccination. Noninferiority was demonstrated if the upper limit of the 95% confidence interview of the group GMT ratio (SA/Co-Ad) was ≤1.5. Secondary descriptive objectives comprised additional immunogenicity assessments, reactogenicity, and safety. RESULTS: Of the 885 participants who received one dose of the study vaccines, 837 were included in the per protocol set. Demographic and baseline characteristics were balanced between the groups. Both co-primary objectives were met for both vaccines. Reported adverse events in both groups were mild-to-moderate, with a low frequency of grade 3 events. CONCLUSIONS: Data from this study demonstrate that RSVPreF3 OA can be co-administered with FLU-QIV. Co-administration is well tolerated, with an acceptable safety profile. CLINICALTRIALS.GOV REGISTRATION: NCT04841577.

Identification of QTLs and candidate genes for water-soluble protein content in soybean seeds.

Soybean represents a vital source of premium plant-based proteins for human nutrition. Importantly, the level of water-soluble protein (WSP) is crucial for determining the overall quality and nutritional value of such crops. Enhancing WSP levels in soybean plants is a high-priority goal in crop improvement. This study aimed to elucidate the genetic basis of WSP content in soybean seeds by identifying quantitative trait loci (QTLs) and set the foundation for subsequent gene cloning and functional analysis. Using 180 F10 recombinant inbred lines generated by crossing the high-protein soybean cultivar JiDou 12 with the wild variety Ye 9, our researcher team mapped the QTLs influencing protein levels, integrating Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and gene expression profiling to identify candidate genes. During the 2020 and 2022 growing seasons, a standard bell-shaped distribution of protein content trait data was observed in these soybean lines. Eight QTLs affecting protein content were found across eight chromosomes, with LOD scores ranging from 2.59 to 7.30, explaining 4.15-11.74% of the phenotypic variance. Notably, two QTLs were newly discovered, one with a elite allele at qWSPC-15 from Ye 9. The major QTL, qWSPC-19, on chromosome 19 was stable across conditions and contained genes involved in nitrogen metabolism, amino acid biosynthesis, and signaling. Two genes from this QTL, Glyma.19G185700 and Glyma.19G186000, exhibited distinct expression patterns at maturity, highlighting the influence of these genes on protein content. This research revealed eight QTLs for WSP content in soybean seeds and proposed a gene for the key QTL qWSPC-19, laying groundwork for gene isolation and enhanced soybean breeding through the use of molecular markers. These insights are instrumental for developing protein-rich soybean cultivars.