Correction: The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

[This corrects the article DOI: 10.1371/journal.pgen.1009726.].

The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

Selective breeding for desirable traits in strictly controlled populations has generated an extraordinary diversity in canine morphology and behaviour, but has also led to loss of genetic variation and random entrapment of disease alleles. As a consequence, specific diseases are now prevalent in certain breeds, but whether the recent breeding practice led to an overall increase in genetic load remains unclear. Here we generate whole genome sequencing (WGS) data from 20 dogs per breed from eight breeds and document a ~10% rise in the number of derived alleles per genome at evolutionarily conserved sites in the heavily bottlenecked cavalier King Charles spaniel breed (cKCs) relative to in most breeds studied here. Our finding represents the first clear indication of a relative increase in levels of deleterious genetic variation in a specific breed, arguing that recent breeding practices probably were associated with an accumulation of genetic load in dogs. We then use the WGS data to identify candidate risk alleles for the most common cause for veterinary care in cKCs-the heart disease myxomatous mitral valve disease (MMVD). We verify a potential link to MMVD for candidate variants near the heart specific NEBL gene in a dachshund population and show that two of the NEBL candidate variants have regulatory potential in heart-derived cell lines and are associated with reduced NEBL isoform nebulette expression in papillary muscle (but not in mitral valve, nor in left ventricular wall). Alleles linked to reduced nebulette expression may hence predispose cKCs and other breeds to MMVD via loss of papillary muscle integrity.

Parental legacy, demography, and admixture influenced the evolution of the two subgenomes of the tetraploid Capsella bursa-pastoris (Brassicaceae).

Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species nor how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe, and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (CbpCg) and the other one from the selfing and genetically depauperate Capsella orientalis (CbpCo). For each subgenome, we assessed its relationship with the diploid relatives, temporal changes of effective population size (Ne), signatures of positive and negative selection, and gene expression patterns. In all three regions, Ne of the two subgenomes decreased gradually over time and the CbpCo subgenome accumulated more deleterious changes than CbpCg. There were signs of widespread admixture between C. bursa-pastoris and its diploid relatives. The two subgenomes were impacted differentially depending on geographic region suggesting either strong interploidy gene flow or multiple origins of C. bursa-pastoris. Selective sweeps were more common on the CbpCg subgenome in Europe and the Middle East, and on the CbpCo subgenome in Asia. In contrast, differences in expression were limited with the CbpCg subgenome slightly more expressed than CbpCo in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy but their evolutionary trajectory strongly varied across geographic regions.

Reanalysis of genome sequences of tomato accessions and its wild relatives: development of Tomato Genomic Variation (TGV) database integrating SNPs and INDELs polymorphisms.

MOTIVATION: Facilitated by technological advances and expeditious decrease in the sequencing costs, whole-genome sequencing is increasingly implemented to uncover variations in cultivars/accessions of many crop plants. In tomato (Solanum lycopersicum), the availability of the genome sequence, followed by the resequencing of tomato cultivars and its wild relatives, has provided a prodigious resource for the improvement of traits. A high-quality genome resequencing of 84 tomato accessions and wild relatives generated a dataset that can be used as a resource to identify agronomically important alleles across the genome. Converting this dataset into a searchable database, including information about the influence of single-nucleotide polymorphisms (SNPs) on protein function, provides valuable information about the genetic variations. The database will assist in searching for functional variants of a gene for introgression into tomato cultivars. RESULTS: A recent release of better-quality tomato genome reference assembly SL3.0, and new annotation ITAG3.2 of SL3.0, dropped 3857 genes, added 4900 novel genes and updated 20 766 genes. Using the above version, we remapped the data from the tomato lines resequenced under the '100 tomato genome resequencing project' on new tomato genome assembly SL3.0 and made an online searchable Tomato Genomic Variations (TGVs) database. The TGV contains information about SNPs and insertion/deletion events and expands it by functional annotation of variants with new ITAG3.2 using SIFT4G software. This database with search function assists in inferring the influence of SNPs on the function of a target gene. This database can be used for selecting SNPs, which can be potentially deployed for improving tomato traits. AVAILABILITY AND IMPLEMENTATION: TGV is freely available at http://psd.uohyd.ac.in/tgv.

ChronQC: a quality control monitoring system for clinical next generation sequencing.

Summary: ChronQC is a quality control (QC) tracking system for clinical implementation of next-generation sequencing (NGS). ChronQC generates time series plots for various QC metrics to allow comparison of current runs to historical runs. ChronQC has multiple features for tracking QC data including Westgard rules for clinical validity, laboratory-defined thresholds and historical observations within a specified time period. Users can record their notes and corrective actions directly onto the plots for long-term recordkeeping. ChronQC facilitates regular monitoring of clinical NGS to enable adherence to high quality clinical standards. Availability and implementation: ChronQC is freely available on GitHub (https://github.com/nilesh-tawari/ChronQC), Docker (https://hub.docker.com/r/nileshtawari/chronqc/) and the Python Package Index. ChronQC is implemented in Python and runs on all common operating systems (Windows, Linux and Mac OS X). Contact: tawari.nilesh@gmail.com or pauline.c.ng@gmail.com. Supplementary information: Supplementary data are available at Bioinformatics online.

In silico, in vitro and in vivo characterization of host-associated Latilactobacillus curvatus strains for potential probiotic applications in farmed Atlantic salmon (Salmo salar).

Salmon aquaculture is the fastest growing animal protein production system in the world; however, intensive farming leads to poor weight gain, stress, and disease outbreaks. Probiotics offer the potential to enhance growth performance and feed efficiency in Atlantic salmon, as well as immunostimulate fish against common pathogens, benefitting farmers and consumers with more efficient production. Here, we isolated and identified 900 native microbial isolates including 18 Lactobacilli from the farmed salmon intestines. Based on whole-genome sequencing and phylogenetic analysis, the Lactobacillus candidates belonged to Latilactobacillus curvatus (L. curvatus) species and formed two distinct phylogenetic groups. Using bioinformatics and in vitro analyses, we selected two candidates L. curvatus ATCC PTA-127116 and L. curvatus ATCC PTA-127117, which showed desirable safety and probiotic properties. The two L. curvatus candidates were evaluated for safety and efficacy (higher final weight) in Atlantic salmon alongside spore-forming Bacilli isolated from salmon, poultry, and swine. All the tested candidates were safe to salmon with no adverse effects. While we did not see efficacy in any Bacillus supplemented groups, compared to untreated group, the group administered with the two L. curvatus strains consortium in feed for seven weeks in freshwater showed indicators of improvement in final body weight by 4.2%. Similarly, the two L. curvatus candidates were also evaluated for safety and efficacy in Atlantic salmon in saltwater; the group administered with the two L. curvatus strains consortium in feed for 11 weeks showed indicators of improvement in final body weight by 4.7%. Comprehensive metabolomics analyses in the presence of different prebiotics and/or additives identified galactooligosaccharide as a potential prebiotic to enhance the efficacy of two L. curvatus candidates. All together, these data provide comprehensive genomic, phenotypic and metabolomic evidence of safety and desirable probiotic properties as well as indicators of in vivo efficacy of two novel endogenous L. curvatus candidates for potential probiotic applications in Atlantic salmon. The in vivo findings need to be confirmed in larger performance studies, including field trials.

In silico, in vitro and in vivo safety evaluation of Limosilactobacillus reuteri strains ATCC PTA-126787 & ATCC PTA-126788 for potential probiotic applications.

The last two decades have witnessed a tremendous growth in probiotics and in the numbers of publications on their potential health benefits. Owing to their distinguishing beneficial effects and long history of safe use, species belonging to the Lactobacillus genus are among the most widely used probiotic species in human food and dietary supplements and are finding increased use in animal feed. Here, we isolated, identified, and evaluated the safety of two novel Limosilactobacillus reuteri (L. reuteri) isolates, ATCC PTA-126787 & ATCC PTA-126788. More specifically, we sequenced the genomes of these two L. reuteri strains using the PacBio sequencing platform. Using a combination of biochemical and genetic methods, we identified the two strains as belonging to L. reuteri species. Detailed in silico analyses showed that the two strains do not encode for any known genetic sequences of concern for human or animal health. In vitro assays confirmed that the strains are susceptible to clinically relevant antibiotics and do not produce potentially harmful by-products such as biogenic amines. In vitro bile and acid tolerance studies demonstrated that the two strains have similar survival profiles as the commercial L. reuteri probiotic strain DSM 17938. Most importantly, daily administration of the two probiotic strains to broiler chickens in drinking water for 26 days did not induce any adverse effect, clinical disease, or histopathological lesions, supporting the safety of the strains in an in vivo avian model. All together, these data provide in silico, in vitro and in vivo evidence of the safety of the two novel candidates for potential probiotic applications in humans as well as animals.

Predictive models for nucleoside bisubstrate analogs as inhibitors of siderophore biosynthesis in Mycobacterium tuberculosis: pharmacophore mapping and chemometric QSAR study.

Inhibitors of aryl acid adenylating enzymes (AAAE), known as MbtA, involved in siderophore biosynthesis in Mycobacterium tuberculosis, are being explored as potential antitubercular agents. In this article, we report the development of a robust pharmacophore model and investigation of structure-activity relationship of several nucleoside bisubstrate analogs reported as MbtA inhibitors. The developed pharmacophore model revealed the importance of two hydrogen bond donors and one hydrogen bond acceptor features. Furthermore, it was found that an aromatic ring at the distal part of molecule away from the two aromatic rings of adenyl moiety is a critical requirement for the tight binding of inhibitor. The generated pharmacophore-based alignment was used to derive a predictive atom-based 3D-QSAR model for training set (r (2) = 0.97, SD = 0.23, F = 310.6, N = 48) and test set (Q (2) = 0.71, RMSE = 0.65, Pearson-R = 0.85, N = 15). Structure-activity relationship investigation further revealed that bulky substitutions at the C-6 position of adenyl moiety is detrimental to activity, while hydrophobic substitutions can be tolerated at C-2 position. Taken together, the PLS-generated QSAR regression cubes along with developed pharmacophore model provide a qualitative picture of the active site and can be used as a powerful tool for the rational modification of bisubstrate inhibitors of MbtA in search of better antitubercular agents. Furthermore, a three-class classification chemometric QSAR model was developed using molecular descriptors for the prediction of whole-cell activity which could be used in the predictive layer for screening of compounds before synthesis.

Multi-Omics Characterization of Host-Derived Bacillus spp. Probiotics for Improved Growth Performance in Poultry.

Microbial feed ingredients or probiotics have been used widely in the poultry industry to improve production efficiency. Spore-forming Bacillus spp. offer advantages over traditional probiotic strains as Bacillus spores are resilient to high temperature, acidic pH, and desiccation. This results in increased strain viability during manufacturing and feed-pelleting processes, extended product shelf-life, and increased stability within the animal's gastrointestinal tract. Despite numerous reports on the use of Bacillus spores as feed additives, detailed characterizations of Bacillus probiotic strains are typically not published. Insufficient characterizations can lead to misidentification of probiotic strains in product labels, and the potential application of strains carrying virulence factors, toxins, antibiotic resistance, or toxic metabolites. Hence, it is critical to characterize in detail the genomic and phenotypic properties of these strains to screen out undesirable properties and to tie individual traits to clinical outcomes and possible mechanisms. Here, we report a screening workflow and comprehensive multi-omics characterization of Bacillus spp. for use in broiler chickens. Host-derived Bacillus strains were isolated and screened for desirable probiotic properties. The phenotypic, genomic and metabolomic analyses of three probiotic candidates, two Bacillus amyloliquefaciens (Ba ATCC PTA126784 and ATCC PTA126785), and a Bacillus subtilis (Bs ATCC PTA126786), showed that all three strains had promising probiotic traits and safety profiles. Inclusion of Ba ATCC PTA12684 (Ba-PTA84) in the feed of broiler chickens resulted in improved growth performance, as shown by a significantly improved feed conversion ratio (3.3%), increased of European Broiler Index (6.2%), and increased average daily gain (ADG) (3.5%). Comparison of the cecal microbiomes from Ba PTA84-treated and control animals suggested minimal differences in microbiome structure, indicating that the observed growth promotion presumably was not mediated by modulation of cecal microbiome.

Singapore Undiagnosed Disease Program: Genomic Analysis aids Diagnosis and Clinical Management.

OBJECTIVE: Use next-generation sequencing (NGS) technology to improve our diagnostic yield in patients with suspected genetic disorders in the Asian setting. DESIGN: A diagnostic study conducted between 2014 and 2019 (and ongoing) under the Singapore Undiagnosed Disease Program. Date of last analysis was 1 July 2019. SETTING: Inpatient and outpatient genetics service at two large academic centres in Singapore. PATIENTS: Inclusion criteria: patients suspected of genetic disorders, based on abnormal antenatal ultrasound, multiple congenital anomalies and developmental delay. EXCLUSION CRITERIA: patients with known genetic disorders, either after clinical assessment or investigations (such as karyotype or chromosomal microarray). INTERVENTIONS: Use of NGS technology-whole exome sequencing (WES) or whole genome sequencing (WGS). MAIN OUTCOME MEASURES: (1) Diagnostic yield by sequencing type, (2) diagnostic yield by phenotypical categories, (3) reduction in time to diagnosis and (4) change in clinical outcomes and management. RESULTS: We demonstrate a 37.8% diagnostic yield for WES (n=172) and a 33.3% yield for WGS (n=24). The yield was higher when sequencing was conducted on trios (40.2%), as well as for certain phenotypes (neuromuscular, 54%, and skeletal dysplasia, 50%). In addition to aiding genetic counselling in 100% of the families, a positive result led to a change in treatment in 27% of patients. CONCLUSION: Genomic sequencing is an effective method for diagnosing rare disease or previous 'undiagnosed' disease. The clinical utility of WES/WGS is seen in the shortened time to diagnosis and the discovery of novel variants. Additionally, reaching a diagnosis significantly impacts families and leads to alteration in management of these patients.

Pharmacophore mapping and electronic feature analysis for a series of nitroaromatic compounds with antitubercular activity.

A five point pharmacophore was generated using PHASE for a series of nitroaromatic compounds and their congeners as antitubercular agents. The generated pharmacophore yielded significant 3D-QSAR model with r(2) of 0.890 for a training set of 92 molecules. The model also showed excellent predictive power with correlation coefficient Q(2) of 0.857 for a test set of 31 compounds. The pharmacophore indicated that presence of a nitro group, a piperazine moiety, one aromatic ring feature and two acceptor features are necessary for potent antitubercular activity. The pharmacophore was supported by electronic property analysis using density functional theory (DFT) at B3LYP/3-21*G level. Molecular electrostatic profile of the compounds was consistent with the generated pharmacophore model, particularly appearance of localized negative potential regions near both the oxygen atoms of nitro group extending laterally to the isoxazole ring system/amide bond in the most active compounds. Calculated data further revealed that all active compounds have smaller LUMO energies located over the nitro group, furan ring, and isoxazole ring/amide bond attached to it. Higher negative values of LUMO energies concentrated over the nitro group are indicative of the electron acceptor capacity of the compounds, suggesting that these compounds are prodrugs and must be activated by TB-nitroreductase. The results obtained from this study should aid in efficient design and development of nitroaromatic compounds as antitubercular agents.

Design, synthesis, and biological evaluation of 4-(5-nitrofuran-2-yl)prop-2-en-1-one derivatives as potent antitubercular agents.

Based on stereoelectronic feature analysis using density functional theory (DFT) at B3LYP/3-21∗G level, a series of 4-(5-nitrofuran-2-yl)prop-2-en-1-one derivatives with low LUMO energies (<-0.10eV); concentrated over the nitro group, furan moiety and α,ß-unsaturated carbonyl bridge were envisaged as potential antitubercular agents. The target compounds were prepared by condensation of 5-nitro-2-furaldehyde with various ketones under acidic condition. The compounds were evaluated for antitubercular activity against Mycobacterium tuberculosis H37Rv and their cytotoxicity in VERO cell line. Several synthesized compounds showed good antitubercular activity of <5µM along with low cytotoxicity. In particular, compound ((E)-3-(5-nitrofuran-2-yl)-1-(4-(piperidin-1-yl)phenyl)prop-2-en-1-one) (3v) was found to be very potent (MIC: 0.19µM) with good selectivity index (MIC(90)/CC(50): >1800). Thus, this study shows the potential of stereoelectronic property analysis in developing improved nitroaromatics as antitubercular agents.

Novel molecular hybrids of cinnamic acids and guanylhydrazones as potential antitubercular agents.

In an attempt to identify potential new agents active against tuberculosis, 20 novel phenylacrylamide derivatives incorporating cinnamic acids and guanylhydrazones were synthesized using microwave assisted synthesis. Activity of the synthesized compounds was evaluated using resazurin microtitre plate assay (REMA) against Mycobacterium tuberculosis H37Rv. Based on empirical structure-activity relationship data it was observed that both steric and electronic parameters play major role in the activity of this series of compounds. Compound 7s (2E)-N-((-2-(3,4-dimethoxybenzylidene) hydrazinyl) (imino) methyl)-3-(4-methoxyphenyl) acrylamide showed MIC of 6.49microM along with good safety profile of >50-fold in VERO cell line. Thus, this compound could act as a potential lead for further antitubercular studies.

Design, synthesis, biological evaluation and computational investigation of novel inhibitors of dihydrofolate reductase of opportunistic pathogens.

The present work deals with design, synthesis and biological evaluation of novel, diverse compounds as potential inhibitors of dihydrofolate reductase (DHFR) from opportunistic microorganisms; Pneumocystis carinii (pc), Toxoplasma gondii (tg) and Mycobacterium avium (ma). A set of 14 structurally diverse compounds were designed with varying key pharmacophoric features of DHFR inhibitors, bulky distal substitutions and different bridges joining the distal part and 2,4-diaminopyrimidine nucleus. The designed compounds were synthesized and evaluated in enzyme assay against pc, tg and ma DHFR. The rat liver (rl) DHFR was used as mammalian standard. As the next logical step of the project, flexible molecular docking studies were carried out to predict the binding modes of these compounds in pcDHFR active site and the obtained docked poses were post processed using MM-GBSA protocol for prediction of relative binding affinity. The predicted binding modes were able to rationalize the experimental results in most cases. Of particular interest, both the docking scores and MM-GBSA predicted Delta G(bind) were able to distinguish between the active and low active compounds. Furthermore, good correlation coefficient of 0.797 was obtained between the IC(50) values and MM-GBSA predicted Delta G(bind). Taken together, the current work provides not only a novel scaffold for further optimization of DHFR inhibitors but also an understanding of the specific interactions of inhibitors with DHFR and structural modifications that improve selectivity.

Synthesis and biological evaluation of biguanide and dihydrotriazine derivatives as potential inhibitors of dihydrofolate reductase of opportunistic microorganisms.

Twenty-one biguanide and dihydrotriazine derivatives were synthesized and evaluated as inhibitors of dihydrofolate reductase (DHFR) from opportunistic microorganisms: Pneumocystis carinii (pc), Toxoplasma gondii (tg), Mycobacterium avium (ma), and rat liver (rl). The most potent compound in the series was B2-07 with 12 nM activity against tgDHFR. The most striking observation was that B2-07 showed similar potency to trimetrexate, approximately 233-fold improved potency over trimethoprim and approximately 7-fold increased selectivity as compared to trimetrexate against tgDHFR. Molecular docking studies in the developed homology model of tgDHFR rationalized the observed potency of B2-07. This molecule can act as a good lead for further design of molecules with better selectivity and improved potency.

3-(5-Nitrofuran-2-yl)prop-2-en-1-one Derivatives, with Potent Antituberculosis Activity, Inhibit A Novel Therapeutic Target, Arylamine N-acetyltransferase, in Mycobacteria.

In this study, the inhibitory potential of 3-(5-nitrofuran-2-yl)prop-2-en-1-one derivatives was evaluated against a panel of bacteria, as well as mammalian cell lines to determine their therapeutic index. In addition, we investigated the mechanism of antibiotic action of the derivatives to identify their therapeutic target. We discovered compound 2 to be an extremely potent inhibitor of Mycobacterium tuberculosis H37Rv growth (MIC: 0.031 mg/L) in vitro, performing better than the currently used first-line antituberculosis drugs such as isoniazid, rifampicin, ethambutol, and pretomanid in vitro. Furthermore, compound 3 was equipotent to pretomanid against a multidrug-resistant M. tuberculosis clinical isolate. The derivatives were selective and bactericidal towards slow-growing mycobacteria. They showed low cytotoxicity towards murine RAW 264.7 and human THP-1 cell lines, with high selectivity indices. Compound 1 effectively eliminated the intracellular mycobacteria in a mycobacteria-infected macrophage model. The derivatives were assessed for their potential to inhibit mycobacterial arylamine N-acetyltransferase (NAT) and were identified as good inhibitors of recombinant mycobacterial NAT, a novel target essential for the intracellular survival of M. tuberculosis. This study provided hits for designing new potent and selective antituberculosis leads, having mycobacterial NAT inhibition as their possible endogenous mechanisms of action.

The genetic basis of inbreeding depression in potato.

Inbreeding depression confers reduced fitness among the offspring of genetic relatives. As a clonally propagated crop, potato (Solanum tuberosum L.) suffers from severe inbreeding depression; however, the genetic basis of inbreeding depression in potato is largely unknown. To gain insight into inbreeding depression in potato, we evaluated the mutation burden in 151 diploid potatoes and obtained 344,831 predicted deleterious substitutions. The deleterious mutations in potato are enriched in the pericentromeric regions and are line specific. Using three F2 populations, we identified 15 genomic regions with severe segregation distortions due to selection at the gametic and zygotic stages. Most of the deleterious recessive alleles affecting survival and growth vigor were located in regions with high recombination rates. One of these deleterious alleles is derived from a rare mutation that disrupts a gene required for embryo development. This study provides the basis for genome design of potato inbred lines.

High throughput sequencing reveals Drosophila suzukii responses to insecticides.

Global climate change and acquired resistance to insecticides are threats to world food security. Drosophila suzukii, a devastating invasive pest in many parts of the world, causes substantial economic losses to fruit production industries, forcing farmers to apply broad-spectrum insecticides frequently. This could lead to the development of insecticide resistance. We determined the Lethal Concentration 50 (median lethal concentration, LC50 ) values of zeta-cypermethrin, spinosad, and malathion insecticides against D. suzukii colonies established from Clarke and Pierce county Georgia, United States. The LC50 values were 3 fold higher in the Pierce county population for all insecticide treatments. We then used RNA sequencing to analyze the responses of Pierce and Clarke population flies surviving a LC50 treatment of the 3 insecticides. We identified a high number of differentially expressed genes that are likely involved in detoxification and reduced cuticular penetration, especially in the Pierce population, with extensive overlap in differentially expressed genes between the 3 insecticide treatments. Finally, we predicted fewer nonsynonymous single nucleotide variants having deleterious effects on protein function among detoxification, insecticide target, and cuticular protein encoding genes in Pierce flies. Thus a combination of increased gene expression and fewer deleterious single nucleotide variants highlights molecular mechanisms underlying the higher LC50 values for Pierce population flies.

Population genomic analyses of the chocolate tree, Theobroma cacao L., provide insights into its domestication process.

Domestication has had a strong impact on the development of modern societies. We sequenced 200 genomes of the chocolate plant Theobroma cacao L. to show for the first time to our knowledge that a single population, the Criollo population, underwent strong domestication ~3600 years ago (95% CI: 2481-13,806 years ago). We also show that during the process of domestication, there was strong selection for genes involved in the metabolism of the colored protectants anthocyanins and the stimulant theobromine, as well as disease resistance genes. Our analyses show that domesticated populations of T. cacao (Criollo) maintain a higher proportion of high-frequency deleterious mutations. We also show for the first time the negative consequences of the increased accumulation of deleterious mutations during domestication on the fitness of individuals (significant reduction in kilograms of beans per hectare per year as Criollo ancestry increases, as estimated from a GLM, P = 0.000425).

Whole-Genome Resequencing Identifies the Molecular Genetic Cause for the Absence of a Gy5 Glycinin Protein in Soybean PI 603408.

During ongoing proteomic analysis of the soybean (Glycine max (L.) Merr) germplasm collection, PI 603408 was identified as a landrace whose seeds lack accumulation of one of the major seed storage glycinin protein subunits. Whole genomic resequencing was used to identify a two-base deletion affecting glycinin 5 The newly discovered deletion was confirmed to be causative through immunological, genetic, and proteomic analysis, and no significant differences in total seed protein content were found to be due to the glycinin 5 loss-of-function mutation per se In addition to focused studies on this one specific glycinin subunit-encoding gene, a total of 1,858,185 nucleotide variants were identified, of which 39,344 were predicted to affect protein coding regions. In order to semiautomate analysis of a large number of soybean gene variants, a new SIFT 4G (Sorting Intolerant From Tolerated 4 Genomes) database was designed to predict the impact of nonsynonymous single nucleotide soybean gene variants, potentially enabling more rapid analysis of soybean resequencing data in the future.