Development of Breeder-Friendly KASP Markers for Low Concentration of Kunitz Trypsin Inhibitor in Soybean Seeds

Trypsin inhibitors (TI), a common anti-nutritional factor in soybean, prevent animals’ protein digestibility reducing animal growth performance. No commercial soybean cultivars with low or null concentration of TI are available. The availability of a high throughput genotyping assay will be beneficial to incorporate the low TI trait into elite breeding lines. The aim of this study is to develop and validate a breeder friendly Kompetitive Allele Specific PCR (KASP) assay linked to low Kunitz trypsin inhibitor (KTI) in soybean seeds. A total of 200 F_3:5 lines derived from PI 547656 (low KTI) X Glenn (normal KTI) were genotyped using the BARCSoySNP6K_v2 Beadchip. F_3:4 and F_3:5 lines were grown in Blacksburg and Orange, Virginia in three years, respectively, and were measured for KTI content using a quantitative HPLC method. We identified three SNP markers tightly linked to the major QTL associated to low KTI in the mapping population. Based on these SNPs, we developed and validated the KASP assays in a set of 93 diverse germplasm accessions. The marker Gm08_44814503 has 86% selection efficiency for the accessions with low KTI and could be used in marker assisted breeding to facilitate the incorporation of low KTI content in soybean seeds.     

A Sequence Variation in GmBADH2 Enhances Soybean Aroma and Is a Functional Marker for Improving Soybean Flavor

The vegetable soybean (Glycine max L. Merr.) plant is commonly consumed in Southeast Asian countries because of its nutritional value and desirable taste. A “pandan-like” aroma is an important value-added quality trait that is rarely found in commercial vegetable soybean varieties. In this study, three novel aromatic soybean cultivars with a fragrant volatile compound were isolated. We confirmed that the aroma of these cultivars is due to the potent volatile compound 2-acetyl-1-pyrroline (2AP) that was previously identified in soybean. A sequence comparison of GmBADH1/2 (encoding an aminoaldehyde dehydrogenase) between aromatic and non-aromatic soybean varieties revealed a mutation with 10 SNPs and an 11-nucleotide deletion in exon 1 of GmBADH2 in Quxian No. 1 and Xiangdou. Additionally, a 2-bp deletion was detected in exon 10 of GmBADH2 in ZK1754. The mutations resulted in a frame shift and the introduction of premature stop codons. Moreover, genetic analyses indicated that the aromatic trait in these three varieties was inherited according to a single recessive gene model. These results suggested that a mutated GmBADH2 may be responsible for the aroma of these three aromatic soybean cultivars. The expression and function of GmBADH2 in aromatic soybean seeds were confirmed by qRT-PCR and CRISPR/Cas9. A functional marker developed on the basis of the mutated GmBADH2 sequence in Quxian No. 1 and Xiangdou was validated in an F₂ population. A perfect association between the marker genotypes and aroma phenotypes implied that GmBADH2 is a major aroma-conferring gene. The results of this study are potentially useful for an in-depth analysis of the molecular basis of 2-AP formation in soybean and the marker-assisted breeding of aromatic vegetable soybean cultivars.     

DArTseq-Based High-Throughput SilicoDArT and SNP Markers Applied for Association Mapping of Genes Related to Maize Morphology

Today, agricultural productivity is essential to meet the needs of a growing population, and is also a key tool in coping with climate change. Innovative plant breeding technologies such as molecular markers, phenotyping, genotyping, the CRISPR/Cas method and next-generation sequencing can help agriculture meet the challenges of the 21st century more effectively. Therefore, the aim of the research was to identify single-nucleotide polymorphisms (SNPs) and SilicoDArT markers related to select morphological features determining the yield in maize. The plant material consisted of ninety-four inbred lines of maize of various origins. These lines were phenotyped under field conditions. A total of 14 morphological features was analyzed. The DArTseq method was chosen for genotyping because this technique reduces the complexity of the genome by restriction enzyme digestion. Subsequently, short fragment sequencing was used. The choice of a combination of restrictases allowed the isolation of highly informative low copy fragments of the genome. Thanks to this method, 90% of the obtained DArTseq markers are complementary to the unique sequences of the genome. All the observed features were normally distributed. Analysis of variance indicated that the main effect of lines was statistically significant (p < 0.001) for all 14 traits of study. Thanks to the DArTseq analysis with the use of next-generation sequencing (NGS) in the studied plant material, it was possible to identify 49,911 polymorphisms, of which 33,452 are SilicoDArT markers and the remaining 16,459 are SNP markers. Among those mentioned, two markers associated with four analyzed traits deserved special attention: SNP (4578734) and SilicoDArT (4778900). SNP marker 4578734 was associated with the following features: anthocyanin coloration of cob glumes, number of days from sowing to anthesis, number of days from sowing to silk emergence and anthocyanin coloration of internodes. SilicoDArT marker 4778900 was associated with the following features: number of days from sowing to anthesis, number of days from sowing to silk emergence, tassel: angle between the axis and lateral branches and plant height. Sequences with a length of 71 bp were used for physical mapping. The BLAST and EnsemblPlants databases were searched against the maize genome to identify the positions of both markers. Marker 4578734 was localized on chromosome 7, the closest gene was Zm00001d022467, approximately 55 Kb apart, encoding anthocyanidin 3-O-glucosyltransferase. Marker 4778900 was located on chromosome 7, at a distance of 45 Kb from the gene Zm00001d045261 encoding starch synthase I. The latter observation indicated that these flanking SilicoDArT and SNP markers were not in a state of linkage disequilibrium.     

Fine Mapping and Identification of a Candidate Gene of Downy Mildew Resistance,RPF2, in Spinach (Spinacia oleracea L.)

Downy mildew is a major threat to the economic value of spinach. The most effective approach to managing spinach downy mildew is breeding cultivars with resistance genes. The resistance allele RPF2 is effective against races 1–10 and 15 of Peronospora farinosa f. sp. Spinaciae (P. effusa) and is widely used as a resistance gene. However, the gene and the linked marker of RPF2 remain unclear, which limit its utilization. Herein, we located the RPF2 gene in a 0.61 Mb region using a BC₁ population derived from Sp39 (rr) and Sp62 (RR) cultivars via kompetitive allele specific PCR (KASP) markers. Within this region, only one R gene, Spo12821, was identified based on annotation information. The amino acid sequence analysis showed that there were large differences in the length of the LRR domain between the parents. Additionally, a molecular marker, RPF2-IN12821, was developed based on the sequence variation in the Spo12821, and the evaluation in the BC₁ population produced a 100% match with resistance/susceptibility. The finding of the study could be valuable for improving our understanding of the genetic basis of resistance against the downy mildew pathogen and breeding resistance lines in the future.     

Genome-Wide Identification and Characterisation of Wheat MATE Genes Reveals Their Roles in Aluminium Tolerance

The Multidrug and toxin efflux (MATE) gene family plays crucial roles in plant growth and development and response to adverse stresses. This work investigated the structural and evolutionary characteristics, expression profiling and potential functions involved in aluminium (Al) tolerance from a genome-wide level. In total, 211 wheat MATE genes were identified, which were classified into four subfamilies and unevenly distributed on chromosomes. Duplication analysis showed that fragments and tandem repeats played the main roles in the amplification of TaMATEs, and Type II functional disproportionation had a leading role in the differentiation of TaMATEs. TaMATEs had abundant Al resistance and environmental stress-related elements, and generally had a high expression level in roots and leaves and in response to Al stress. The 3D structure prediction by AlphaFold and molecular docking showed that six TaMATE proteins localised in the plasmalemma could combine with citrate via amino acids in the citrate exuding motif and other sites, and then transport citrate to soil to form citrate aluminium. Meanwhile, citrate aluminium formed in root cells might be transported to leaves by TaMATEs to deposit in vacuoles, thereby alleviating Al toxicity.     

Development of 101 Gene-based Single Nucleotide Polymorphism Markers in Sea Cucumber, Apostichopus japonicus

Single nucleotide polymorphisms (SNPs) are currently the marker of choice in a variety of genetic studies. Using the high resolution melting (HRM) genotyping approach, 101 gene-based SNP markers were developed for Apostichopus japonicus, a sea cucumber species with economic significance for the aquaculture industry in East Asian countries. HRM analysis revealed that all the loci showed polymorphisms when evaluated using 40 A. japonicus individuals collected from a natural population. The minor allele frequency ranged from 0.035 to 0.489. The observed and expected heterozygosities ranged from 0.050 to 0.833 and 0.073 to 0.907, respectively. Thirteen loci were found to depart significantly from Hardy-Weinberg equilibrium (HWE) after Bonferroni corrections. Significant linkage disequilibrium (LD) was detected in one pair of markers. These SNP markers are expected to be useful for future quantitative trait loci (QTL) analysis, and to facilitate marker-assisted selection (MAS) in A. japonicus.     

Unravelling the Contribution of the rs7041 and rs4588 Polymorphisms of the GC Gene and Serum VDBP Levels for Developing Metabolic Syndrome in the Mexican Population

Metabolic syndrome (MetS) is a multifactorial disorder integrated by a constellation of cardiovascular risk factors. The genetic and environmental determinants of MetS are not fully elucidated. This study investigated the association of two common single nucleotide polymorphisms (SNPs) on GC, rs7041 and rs4588, derived haplotypes, and serum vitamin D binding protein (VDBP) levels with the susceptibility to suffer MetS in Mexican adults. We included 1924 individuals; clinical and biochemical data were obtained through standard methods. Genotyping was performed through predesigned TaqMan assays. Logistic regression models were used to assess the associations of interest. Prevalence of MetS was 52.9% in the whole population, being more frequent in women. We observed that some association results differed between sexes. The GG genotype of the rs7041 was associated with increased odds of MetS in women. For the rs4588, the CA genotype had a protective effect against MetS in women. The haplotype GC2 was associated with reduced odds for MetS and some of its components in women. Our data suggest that VDBP serum levels were influenced by genotypes/haplotypes and this interplay seems to influence the risk of MetS. Our data provide reliable evidence regarding the association of GC polymorphisms with MetS risk in Mexican women.     

Development and Characterization of New Single Nucleotide Polymorphism Markers from Expressed Sequence Tags in Common Carp (Cyprinus carpio)

The common carp (Cyprinus carpio) is an important aquaculture fish worldwide but only limited single nucleotide polymorphism (SNP) markers are characterized from expressed sequence tags (ESTs) in this species. In this study, 1487 putative SNPs were bioinformatically mined from 14,066 online ESTs mainly from the European common carp, with the occurrence rate of about one SNP every 173 bp. One hundred and twenty-one of these SNPs were selected for validation using PCR fragment sequencing, and 48 out of 81 primers could amplify the expected fragments in the Chinese common carp genome. Only 26 (21.5%) putative SNPs were validated, however, 508 new SNPs and 68 indels were identified. The ratios of transitions to transversions were 1.77 for exon SNPs and 1.05 for intron SNPs. All the 23 SNPs selected for population tests were polymorphic, with the observed heterozygosity (Ho) ranging from 0.053 to 0.526 (mean 0.262), polymorphism information content (PIC) from 0.095 to 0.357 (mean 0.246), and 21 SNPs were in Hardy-Weinberg equilibrium. These results suggest that different common carp populations with geographic isolation have significant genetic variation at the SNP level, and these new EST-SNP markers are readily available for genetics and breeding studies in common carp.     

Using Staphylococcus aureus Cas9 to Expand the Scope of Potential Gene Targets for Genome Editing in Soybean

The CRISPR (clustered regularly interspaced short palindromic repeat)/Cas (CRISPR-associated) is a revolutionary genome editing technology that has been used to achieve site-specific gene knock-out, large fragment deletion, or base editing in many plant species including soybean (Glycine max). The Streptococcus pyogenes Cas9 (SpCas9) is widely used in plants at present, although there are some reports describing the application of CRISPR/Cpf1 in soybean. Therefore, the selection range of PAM (protospacer adjacent motif) sequences for soybean is currently limited to 5′-NGG-3′ (SpCas9) or 5′-TTTN-3′ (Cpf1), which in turn limits the number of genes that can be mutated. Another Cas9 enzyme from Staphylococcus aureus (SaCas9) recognizes the PAM sequence 5′-NNGRRT-3′ (where R represents A or G), which can provide a wider range of potential target sequences. In this study, we developed a CRISPR/SaCas9 system and used this tool to specifically induce targeted mutations at five target sites in the GmFT2a (Glyma.16G150700) and GmFT5a (Glyma.16G044100) genes in soybean hairy roots. We demonstrated that this tool can recognize the PAM sequences 5′-AAGGGT-3′, 5′-GGGGAT-3′, 5′-TTGAAT-3′, and 5′-TAGGGT-3′ in soybean, and it achieved mutation rates ranging from 34.5% to 73.3%. Our results show that we have established a highly efficient CRISPR/SaCas9 tool that is as suitable as SpCas9 for genome editing in soybean, and it will be useful for expanding the range of target sequences for genome editing.     

Identification and Validation of Quantitative Trait Loci Mapping for Spike-Layer Uniformity in Wheat

Spike-layer uniformity (SLU), the consistency of the spike distribution in the vertical space, is an important trait. It directly affects the yield potential and appearance. Revealing the genetic basis of SLU will provide new insights into wheat improvement. To map the SLU-related quantitative trait loci (QTL), 300 recombinant inbred lines (RILs) that were derived from a cross between H461 and Chinese Spring were used in this study. The RILs and parents were tested in fields from two continuous years from two different pilots. Phenotypic analysis showed that H461 was more consistent in the vertical spatial distribution of the spike layer than in Chinese Spring. Based on inclusive composite interval mapping, four QTL were identified for SLU. There were two major QTL on chromosomes 2BL and 2DL and two minor QTL on chromosomes 1BS and 2BL that were identified. The additive effects of QSlu.sicau-1B, Qslu.sicau-2B-2, and QSlu.sicau-2D were all from the parent, H461. The major QTL, QSlu.sicau-2B-2 and QSlu.sicau-2D, were detected in each of the conducted trials. Based on the best linear unbiased prediction values, the two loci explained 23.97% and 15.98% of the phenotypic variation, respectively. Compared with previous studies, the two major loci were potentially novel and the two minor loci were overlapped. Based on the kompetitive allele-specific PCR (KASP) marker, the genetic effects for QSlu.sicau-2B-2 were validated in an additional RIL population. The genetic effects ranged from 26.65% to 32.56%, with an average value of 30.40%. In addition, QSlu.sicau-2B-2 showed a significant (p < 0.01) and positive influence on the spike length, spikelet number, and thousand kernel weight. The identified QTL and the developed KASP marker will be helpful for fine-mapping these loci, finally contributing to wheat breeding programs in a marker-assisted selection way.     

The Influence of PSCA Gene Variation on Its Expression and Gastric Adenocarcinoma Susceptibility in the Northwest Chinese Population

Gastric adenocarcinoma (GAC) imposes a considerable health burden around the world. Gene variation in prostate stem cell antigen gene (PSCA) has been identified to be associated with GAC risk, while the results showed regional variation. To explore the influence of PSCA gene variation on its expression and GAC risk in the Northwest Chinese population, four single nucleotide polymorphisms (SNPs) of PSCA were genotyped in 476 GAC cases and 481 controls using MassARRAY system. Two SNPs of rs2294008 (C>T) and rs2976392 (G>A) were identified to be associated with GAC risk. rs2294008, rs2976392 and rs10216533 made up two statistically significant haplotypes (Hap-CGG and Hap-TAG). Additionally, PSCA expression was analyzed by quantitative real time PCR, immunohistochemistry and tissue microarray. The results showed that PSCA expression was decreased in GAC tissues compared with adjacent normal tissues. For normal tissues, PSCA expression was higher with Hap-TA than that with Hap-CG. For GAC tissues, the differentiation degree of Hap-TA was higher than that of Hap-CG. The expression distribution of PSCA in multiple human organs showed disparity. These results suggest that PSCA gene variation has a potential effect on its expression and GAC risk in the Northwest Chinese population.     

Ontogenic Expression Pattern and Genetic Polymorphisms of the Fatty Acid Transport Protein 4 (FATP4) Gene in Chinese Chicken Populations

In the current research, the polymorphism of FATP4 gene was analyzed in Erlang Mountainous chickens. A total of nine genetic variants were identified by FATP4 gene sequencing analysis across the chicken samples. Significant associations (p < 0.05) were observed for two SNPs (g.5608778C>T and g.5608814G>A in exon 6) with certain carcass traits (such as live weight, carcass weight, eviscerated weight) in S01 and S05 populations, respectively. Meanwhile, in S05 population, haplotype 3 (T-G) and haplotype 2 (C-A) were associated with higher and lower partial carcass traits such as live weight, carcass weight, eviscerated weight and semi-eviscerated weight, respectively. Moreover, we investigated the expression profile of this gene during ontogenesis in Mountainous black-boned chicken. Quantitative real-time PCR analysis showed that FATP4 mRNA had the highest expression level in small intestine tissue over all other tissues examined. The FATP4 mRNA levels presented remarkable developmental changes with age in the various tissues. These results suggested that the FATP4 gene might play an important role in controlling chicken carcass traits.     

Expression of Specific Alleles of Zinc-Finger Transcription Factors,HvSAP8 and HvSAP16, and Corresponding SNP Markers,Are Associated with Drought Tolerance in Barley Populations

Two genes, HvSAP8 and HvSAP16, encoding Zinc-finger proteins, were identified earlier as active in barley plants. Based on bioinformatics and sequencing analysis, six SNPs were found in the promoter regions of HvSAP8 and one in HvSAP16, among parents of two barley segregating populations, Granal × Baisheshek and Natali × Auksiniai-2. ASQ and Amplifluor markers were developed for HvSAP8 and HvSAP16, one SNP in each gene, and in each of two populations, showing simple Mendelian segregation. Plants of F₆ selected breeding lines and parents were evaluated in a soil-based drought screen, revealing differential expression of HvSAP8 and HvSAP16 corresponding with the stress. After almost doubling expression during the early stages of stress, HvSAP8 returned to pre-stress level or was strongly down-regulated in plants with Granal or Baisheshek genotypes, respectively. For HvSAP16 under drought conditions, a high expression level was followed by either a return to original levels or strong down-regulation in plants with Natali or Auksiniai-2 genotypes, respectively. Grain yield in the same breeding lines and parents grown under moderate drought was strongly associated with their HvSAP8 and HvSAP16 genotypes. Additionally, Granal and Natali genotypes with specific alleles at HvSAP8 and HvSAP16 were associated with improved performance under drought via higher 1000 grain weight and more shoots per plant, respectively.     

Development of 20 Microsatellite Markers for Solenocera crassicornis and Their Cross-Species Application in Solenocera melantho

Twenty microsatellite markers were isolated and characterized for Solenocera crassicornis from a (GT)13-enriched genomic library. Their polymorphisms were investigated using 44 wild individuals from the South Yellow Sea. Our investigation revealed that all the markers were polymorphic. The number of alleles per locus varied from 6 to 19 with an average of 12.35. The observed and expected heterozygosities ranged from 0.400 to 0.977 and from 0.609 to 0.940, with averages of 0.788 and 0.859, respectively. Four loci significantly deviated from Hardy-Weinberg equilibrium after Bonferroni's correction. Cross-species amplification was also conducted in Solenocera melantho collected from the East China Sea. The result showed that 14 loci could be amplified from Solenocera melantho DNAs. These polymorphic markers would be useful for assessment of genetic variation and population structure of S. crassicornis and S. melantho.     

A Simple Strategy for Development of Single Nucleotide Polymorphisms from Non-Model Species and Its Application in Panax

Single nucleotide polymorphisms (SNPs) are widely employed in the studies of population genetics, molecular breeding and conservation genetics. In this study, we explored a simple route to develop SNPs from non-model species based on screening the library of single copy nuclear genes (SCNGs). Through application of this strategy in Panax, we identified 160 and 171 SNPs from P. quinquefolium and P. ginseng, respectively. Our results demonstrated that both P. ginseng and P. quinquefolium possessed a high level of nucleotide diversity. The number of haplotype per locus ranged from 1 to 12 for P. ginseng and from 1 to 9 for P. quinquefolium, respectively. The nucleotide diversity of total sites (π_T) varied between 0.000 and 0.023 for P. ginseng and 0.000 and 0.035 for P. quinquefolium, respectively. These findings suggested that this approach is well suited for SNP discovery in non-model organisms and is easily employed in standard genetics laboratory studies.     

Candidate Genes in Testing Strategies for Linkage Analysis and Bioinformatic Sorting of Whole Genome Sequencing Data in Three Small Japanese Families with Idiopathic Superior Oblique Muscle Palsy

Idiopathic superior oblique muscle palsy is a major type of paralytic, non-comitant strabismus and presents vertical and cyclo-torsional deviation of one eye against the other eye, with a large vertical fusion range and abnormal head posture such as head tilt. Genetic background is considered to play a role in its development, as patients with idiopathic superior oblique muscle palsy have varying degrees of muscle hypoplasia and, rarely, the complete absence of the muscle, that is, aplasia. In this study, whole genome sequencing was performed, and single nucleotide variations and short insertions/deletions (SNVs/InDels) were annotated in two patients each in three small families (six patients in total) with idiopathic superior oblique muscle palsy, in addition to three normal individuals in one family. At first, linkage analysis was carried out in the three families and SNVs/InDels in chromosomal loci with negative LOD scores were excluded. Next, SNVs/InDels shared by the six patients, but not by the three normal individuals, were chosen. SNVs/InDels were further narrowed down by choosing low-frequency (<1%) or non-registered SNVs/InDels in four databases for the Japanese population, and then by choosing SNVs/InDels with functional influence, leading to one candidate gene, SSTR5-AS1 in chromosome 16. The six patients were heterozygous for 13-nucleotide deletion in SSTR5-AS1, except for one homozygous patient, while the three normal individuals were wild type. Targeted polymerase chain reaction (PCR) and direct sequencing of PCR products confirmed the 13-nucleotide deletion in SSTR5-AS1. In the face of newly-registered SSTR5-AS1 13-nucleotide deletion at a higher frequency in a latest released database for the Japanese population, the skipping of low-frequency and non-registration sorting still resulted in only 13 candidate genes including SSTR5-AS1 as common variants. The skipping of linkage analysis also led to the same set of 13 candidate genes. Different testing strategies that consisted of linkage analysis and simple unintentional bioinformatics could reach candidate genes in three small families with idiopathic superior oblique muscle palsy.