High-confidence Single Nucleotide Polymorphisms Research Articles

SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants.

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

Genetic diversity, population structure analysis and codon substitutions of Indicine Badri cattle using ddRAD sequencing.

The present study was carried out on 96 animals representing three distinct colour variants of Badri cattle to investigate the genetic diversity, population structure and substitution mutations in the genetic codons due to single nucleotide variations. The DNA samples of 96 Badri cows were genotyped using a double digestion restriction associated DNA (ddRAD) sequencing approach. A standardized bioinformatics pipeline was employed to identify single nucleotide polymorphisms (SNPs), initially detecting 7,168,552 SNPs through alignment with the Bos indicus reference genome assembly. Subsequent stringent quality filtration yielded 65,483 high-confidence SNPs for downstream analysis. Genetic diversity analysis of the Badri cattle population resulted in average values of 0.145, 0.088, and 0.091 for Shannon's diversity Index (I), Simpson's Diversity (h), and Simpson's Unbiased Diversity (uh), respectively. Genetic similarities between the black and brown, black and grey, and brown and grey Badri variantswere found to be 0.9972, 0.9980 and 0.9970, respectively. Tajima's D diversity value was observed to be significant and positive for 99.29% of high-confidence SNPs (65,483). STRUCTURE analysis showed admixture among the three Badri colour variants, suggesting a lack of genetic differentiation. Annotation of high-confidence SNPs regarding genetic codon changes indicated maximum substitutions in the GGC with GGT (22 occurrences), followed by AAC to AGC (20 occurrences), GAA to TAA (19 occurrences) and CAA to CAG (19 occurrences). The study concludes there are genetic similarities among colour variants, lack of rare alleles, balancing selection, sudden population contraction and genetic codon substitutions within the Badri cattle population. Insights derived from SNP data analysis hold potential significance for conservation initiatives and breed improvement programs for indicine cattle.

Genome-wide association study of kernel colour traits and mining of elite alleles from the major loci in maize

BackgroundMaize kernel colour is an important index for evaluating maize quality and value and mainly entails two natural pigments, carotenoids and anthocyanins. To analyse the genetic mechanism of maize kernel colour and mine single nucleotide polymorphisms (SNPs) related to kernel colour traits, an association panel including 244 superior maize inbred lines was used to measure and analyse the six traits related to kernel colour in two environments and was then combined with the about 3 million SNPs covering the whole maize genome in this study. Two models (Q + K, PCA + K) were used for genome-wide association analysis (GWAS) of kernel colour traits.ResultsWe identified 1029QTLs, and two SNPs contained in those QTLs were located in coding regions of Y1 and R1 respectively, two known genes that regulate kernel colour. Fourteen QTLs which contain 19 SNPs were within 200 kb interval of the genes involved in the regulation of kernel colour. 13 high-confidence SNPs repeatedly detected for specific traits, and AA genotypes of rs1_40605594 and rs5_2392770 were the most popular alleles appeared in inbred lines with higher levels. By searching the confident interval of the 13 high-confidence SNPs, a total of 95 candidate genes were identified.ConclusionsThe genetic loci and candidate genes of maize kernel colour provided in this study will be useful for uncovering the genetic mechanism of maize kernel colour, gene cloning in the future. Furthermore, the identified elite alleles can be used to molecular marker-assisted selection of kernel colour traits.

Genome-wide identification and annotation of SNPs and their mapping in candidate genes related to milk production and fertility traits in Badri cattle

This study was conducted in Badri cattle using a double digest restriction-site associated DNA sequencing approach. The study aimed to identify and annotate high confidence single nucleotide polymorphisms (SNPs) and their mapping in candidate genes related to production and fertility in dairy cattle. A total of 7,168,552 genome-wide SNPs were initially identified in Badri cattle by alignment with theBos indicusreference genome. After filtration of SNPs, 65,483 high confidence SNPs were retained and further used for downstream analysis. Annotation of high confidence SNPs revealed 99.197% SNPs had modifier impact, 0.326% SNPs were low impact, 0.036% were high impact, and 0.441% were moderate impact SNPs. Most SNPs in Badri cattle were found in intergenic, transcript and intronic regions. The candidate genes for milk production PRKCE, ABCG2, GHR, EPS8, CAST and NRXN1 were found to harbour maximum high confidence variants. Among candidate genes for fertility in cattle, ATP2B1, SOX5, WDR27, ARHGAP12, CACNA1D, ANKRD6, GRIA3, ZNF521 and CAST822 have maximum high confidence variants mapped in them. The SNPs found mapped in the candidate genes will be important genetic tools in the search for phenotype-modifying nucleotide changes and will aid in formulating relevant genetic improvement programmes for dairy cattle.

Genome-wide association study identifies 12 new genetic loci associated with growth traits in pigs

Growth traits are among the most important economic traits in pigs and are regulated by polygenes with complex regulatory mechanisms. As the major indicators of growth performance, the backfat thickness (BFT), loin eye area (LEA), and days to 100 kg (D100) traits are commonly used to the genetics improvement in pigs. However, the available genetic markers for these traits are limited. To uncover novel loci and candidate genes associated with growth performance, we collected the phenotypic information of BFT, LEA, and D100 in 1,186 pigs and genotyped all these individuals using the Neogen GGP porcine 80K BeadChip. We performed a genome-wide association study (GWAS) using 4 statistical models, including mixed linear models (MLM), fixed and random model circulating probability unification (FarmCPU), settlement of MLM under progressively exclusive relationships (SUPER), Bayesian-information and linkage-disequilibrium Iteratively nested keyway (Blink), and identified 5, 3, and 6 high-confidence single nucleotide polymorphisms (SNPs) associated with BFT, LEA, and D100, respectively. Variant annotation and quantitative trait locus (QTL) mapping analysis suggested that 6 genes (SKAP2, SATB1, PDE7B, PPP1R16B, WNT3, and WNT9B) were potentially associated with growth performance in pigs. Transcriptome analysis suggested that the expression of Src Kinase Associated Phosphoprotein 2 (SKAP2) was higher in prenatal muscles than in postnatal muscles, and the expression of Phosphodiesterase 7B (PDE7B) continuously increased during the prenatal stages and gradually decreased after birth, implying their potential roles in prenatal skeletal muscle development. Overall, this study provides new candidate loci and genes for the genetic improvement of pigs.

Genome-wide association study in quinoa reveals selection pattern typical for crops with a short breeding history.

Quinoa germplasm preserves useful and substantial genetic variation, yet it remains untapped due to a lack of implementation of modern breeding tools. We have integrated field and sequence data to characterize a large diversity panel of quinoa. Whole-genome sequencing of 310 accessions revealed 2.9 million polymorphic high confidence single nucleotide polymorphism (SNP) loci. Highland and Lowland quinoa were clustered into two main groups, with FST divergence of 0.36 and linkage disequilibrium (LD) decay of 6.5 and 49.8 kb, respectively. A genome-wide association study using multi-year phenotyping trials uncovered 600 SNPs stably associated with 17 traits. Two candidate genes are associated with thousand seed weight, and a resistance gene analog is associated with downy mildew resistance. We also identified pleiotropically acting loci for four agronomic traits important for adaptation. This work demonstrates the use of re-sequencing data of an orphan crop, which is partially domesticated to rapidly identify marker-trait association and provides the underpinning elements for genomics-enabled quinoa breeding.

Newly Developed MAGIC Population Allows Identification of Strong Associations and Candidate Genes for Anthocyanin Pigmentation in Eggplant.

Multi-parent advanced generation inter-cross (MAGIC) populations facilitate the genetic dissection of complex quantitative traits in plants and are valuable breeding materials. We report the development of the first eggplant MAGIC population (S3 Magic EGGplant InCanum, S3MEGGIC; 8-way), constituted by the 420 S3 individuals developed from the intercrossing of seven cultivated eggplant (Solanum melongena) and one wild relative (S. incanum) parents. The S3MEGGIC recombinant population was genotyped with the eggplant 5k probes SPET platform and phenotyped for anthocyanin presence in vegetative plant tissues (PA) and fruit epidermis (FA), and for the light-insensitive anthocyanic pigmentation under the calyx (PUC). The 7,724 filtered high-confidence single-nucleotide polymorphisms (SNPs) confirmed a low residual heterozygosity (6.87%), a lack of genetic structure in the S3MEGGIC population, and no differentiation among subpopulations carrying a cultivated or wild cytoplasm. Inference of haplotype blocks of the nuclear genome revealed an unbalanced representation of the founder genomes, suggesting a cryptic selection in favour or against specific parental genomes. Genome-wide association study (GWAS) analysis for PA, FA, and PUC detected strong associations with two myeloblastosis (MYB) genes similar to MYB113 involved in the anthocyanin biosynthesis pathway, and with a COP1 gene which encodes for a photo-regulatory protein and may be responsible for the PUC trait. Evidence was found of a duplication of an ancestral MYB113 gene with a translocation from chromosome 10 to chromosome 1 compared with the tomato genome. Parental genotypes for the three genes were in agreement with the identification of the candidate genes performed in the S3MEGGIC population. Our new eggplant MAGIC population is the largest recombinant population in eggplant and is a powerful tool for eggplant genetics and breeding studies.

Revealing the coexistence of differentiation and communication in an endemic hare, Lepus yarkandensis (Mammalia, Leporidae) using specific-length amplified fragment sequencing

BackgroundThe Yarkand hare (Lepus yarkandensis Günther, 1875) is endemic to oasis and desert areas around the Tarim Basin in the Xinjiang Uyghur Autonomous Region of northwest China; however, genome-wide information for this species remains limited. Moreover, the genetic variation, genetic structure, and phylogenetic relationships of Yarkand hare from the plateau mountain regions have not been reported. Thus, we used specific-length amplified fragment sequencing (SLAF-seq) technology to evaluate the genetic diversity of 76 Yarkand hares from seven geographic populations in the northern and southwestern parts of the Tarim Basin to investigate single-nucleotide polymorphism (SNP) marker-based population differentiation and evolutionary processes. Selective sweep analysis was conducted to identify genetic differences between populations.ResultsUsing SLAF-seq, a total of 1,835,504 SNPs were initially obtained, of which 308,942 high-confidence SNPs were selected for further analysis. Yarkand hares exhibited a relatively high degree of genetic diversity at the SNP level. Based on pairwise FST estimates, the north and southwest groups showed a moderate level of genetic differentiation. Phylogenetic tree and population structure analyses demonstrated evident systematic phylogeographical structure patterns consistent with the geographical distribution of the hares. Hierarchical analysis of molecular variation further indicated that genetic variation was mainly observed within populations. Low to moderate genetic differentiation also occurred among populations despite a common genomic background, likely due to geographical barriers, genetic drift, and differential selection pressure of distinct environments. Nevertheless, the observed lineage-mixing pattern, as indicated by the evolutionary tree, principal component analysis, population structure, and TreeMix analyses, suggests a certain degree of gene flow between the north and southwest groups. This may be related to the migration of hares to high-altitude water sources southwest of the basin during glacial climatic oscillations, as well as river re-diffusion and oasis restoration in the basin following the glacial period. We also identified candidate genes, and their associated gene ontology terms and pathways, related to the adaptation of Yarkand hares to different environmental habitats.ConclusionsThe identified genome-wide SNPs, genetic diversity, and population structure of Yarkand hares expand our understanding of the genetic background of this endemic species and provide valuable insights into its environmental adaptation, allowing for further exploration of the underlying mechanisms.

Integration of association and computational methods reveals functional variants of LEPR gene for abdominal fat content in chickens

Leptin receptor (LEPR) plays a vital role in obesity in humans and animals. The objective of this study is to assess LEPR functional variants for chicken adipose deposition by integration of association and in-silico analysis using a unique chicken population, the Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF). Five online bioinformatics tools were used to predict the functionality of the single nucleotide polymorphisms (SNPs) in coding region. Further, the possible structure–function relationship of high confidence SNPs was determined by bioinformatics analyses, including the conservation and stability analysis based on amino acid residues, prediction of protein ligand-binding sites, and the superposition of protein tertiary structure. Meanwhile, we analyzed the association between abdominal fat traits and 20 polymorphisms of chicken LEPR gene. The integrated results showed that rs731962924 (N867I) and rs13684622 (C1002R) could lead to striking changes in the structure and function of proteins, of which rs13684622 (C1002R) was significantly associated with abdominal fat weight (AFW, P=0.0413) and abdominal fat percentage (AFP, P=0.0260) in chickens. Therefore, we are of the opinion that rs13684622 (C1002R) may be an essential functional SNP affecting chicken abdominal fat deposition, and potentially applied to improvement of broiler abdominal fat in molecular marker-assisted selection (MAS) program. Additionally, the coupling of association with computer electronic predictive analysis provides a new avenue to identify important molecular markers for breeders.

Genome-wide association study and transcriptome analysis discover new genes for bacterial leaf blight resistance in rice (Oryza sativa L.)

BackgroundRice (Oryza sativa) bacterial leaf blight (BLB), caused by the hemibiotrophic Xanthomonas oryzae pv. oryzae (Xoo), is one of the most devastating diseases affecting the production of rice worldwide. The development and use of resistant rice varieties or genes is currently the most effective strategy to control BLB.ResultsHere, we used 259 rice accessions, which are genotyped with 2 888 332 high-confidence single nucleotide polymorphisms (SNPs). Combining resistance variation data of 259 rice lines for two Xoo races observed in 2 years, we conducted a genome-wide association study (GWAS) to identify quantitative trait loci (QTL) conferring plant resistance against BLB. The expression levels of genes, which contains in GWAS results were also identified between the resistant and susceptible rice lines by transcriptome analysis at four time points after pathogen inoculation. From that 109 candidate resistance genes showing significant differential expression between resistant and susceptible rice lines were uncovered. Furthermore, the haplotype block structure analysis predicted 58 candidate genes for BLB resistance based on Chr. 7_707158 with a minimum P-value (–log 10 P = 9.72). Among them, two NLR protein-encoding genes, LOC_Os07g02560 and LOC_Os07g02570, exhibited significantly high expression in the resistant line, but had low expression in the susceptible line of rice.ConclusionsTogether, our results reveal novel BLB resistance gene resources, and provide important genetic basis for BLB resistance breeding of rice crops.

DdRAD sequencing: an emerging technology added to the biosecurity toolbox for tracing the origin of brown marmorated stink bug, Halyomorpha halys (Hemiptera: Pentatomidae)

BackgroundBrown marmorated stink bug (BMSB), Halyomorpha halys (Hemiptera: Pentatomidae) is native to East Asia but has invaded many countries in the world. BMSB is a polyphagous insect pest and causes significant economic losses to agriculture worldwide. Knowledge on the genetic diversity among BMSB populations is scarce but is essential to understand the patterns of colonization and invasion history of local populations. Efforts have been made to assess the genetic diversity of BMSB using partial mitochondrial DNA sequences but genetic divergence on mitochondria is not high enough to precisely accurately identify and distinguish various BMSB populations. Therefore, in this study, we applied a ddRAD (double digest restriction-site associated DNA) sequencing approach to ascertain the genetic diversity of BMSB populations collected from 12 countries (2 native and 10 invaded) across four continents with the ultimate aim to trace the origin of BMSBs intercepted during border inspections and post-border surveillance.ResultA total of 1775 high confidence single nucleotide polymorphisms (SNPs) were identified from ddRAD sequencing data collected from 389 adult BMSB individuals. Principal component analysis (PCA) of the identified SNPs indicated the existence of two main distinct genetic clusters representing individuals sampled from regions where BMSB is native to, China and Japan, respectively, and one broad cluster comprised individuals sampled from countries which have been invaded by BMSB. The population genetic structure analysis further discriminated the genetic diversity among the BMSB populations at a higher resolution and distinguished them into five potential genetic clusters.ConclusionThe study revealed hidden genetic diversity among the studied BMSB populations across the continents. The BMSB populations from Japan were genetically distant from the other studied populations. Similarly, the BMSB populations from China were also genetically differentiated from the Japanese and other populations. Further genetic structure analysis revealed the presence of at least three genetic clusters of BMSB in the invaded countries, possibly originating via multiple invasions. Furthermore, this study has produced novel set of SNP markers to enhance the knowledge of genetic diversity among BMSB populations and demonstrates the potential to trace the origin of BMSB individuals for future invasion events.

Discovery of single nucleotide polymorphisms in bone morphogenetic protein (

Context Two native goat breeds from Kerala, Malabari and Attappady Black, differ significantly in prolificacy (i.e. no. of kids born/kidding). Prolificacy is an important economic trait and the subject of genetic research showing that bone morphogenetic protein (BMP) genes have a significant effect. Double digest restriction-site associated DNA sequencing (ddRADseq) is a highly efficient and low cost technology for high density discovery of single nucleotide polymorphisms (SNPs), which could serve as predictive markers for animal breeding programs. Aims The study was aimed at discovering SNPs in BMP genes that affect prolificacy, using ddRADseq followed by validation of selected SNP. Methods Blood DNA samples of 10 highly prolific Malabari and 10 less prolific Attappady Black goats were pooled by group and subjected to ddRADseq. SNPs observed in BMP genes were investigated and compared between groups. A validation study was done for the c.614–32789C>T variant in 100 Malabari and 50 Attappady Black goats by using PCR-RFLP. Key results In total, 6333 variants were identified by ddRADseq. Three variants were identified in BMP genes, which included two intronic variants c.614–32789C>T and c.490+6793T>C, in genes BMP6 and BMP5 and a downstream gene variant near the BMPR1B gene. According to ddRADseq data, variants in BMP5 and BMP6 differed in allelic distribution between Malabari and Attappady Black goats. For c.490+6793T>C in BMP5, the CC genotype was predominant in the highly prolific Malabari whereas TC was present in the Attappady Black group. The variant c.614–32789C>T in BMP6 was genotyped as TC in Malabari and CC in Attappady Black goats by ddRADseq. This variant was predicted to have an effect on splicing, according to the tool SplicePort. On the basis of bioinformatics analysis and the role of BMP6 gene in follicular dynamics, the variant in BMP6 was selected for further validation studies. All three genotypes were identified by PCR-RFLP; the C allele was the rare allele in the population with an allele frequency of 0.36. Presence of both alleles C and T and the three genotypes CC, TC and TT in this larger population substantiated the robustness of ddRADseq technique. Conclusions The technique discovered high confidence SNPs, which could be used for further validation and association studies to develop markers for selection of animals and for genetic improvement of this complex trait. Implications Techniques such as ddRADseq provide a large number of SNPs, and investigation of those polymorphisms found across the genome will help to identify new loci affecting traits of interest. This, in turn, will aid in exploring genetically complex traits in a faster and cheaper manner.

Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii.

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii. We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures.

Genome-Wide Association Study Reveals the Genetic Architecture of Stripe Rust Resistance at the Adult Plant Stage in Chinese Endemic Wheat.

Chinese endemic wheat, comprising Tibetan semi-wild wheat (Triticum aestivum ssp. tibetanum), Yunnan hulled wheat (T. aestivum ssp. yunnanense), and Xinjiang rice wheat (T. petropavlovskyi), are genetically and morphologically unique. To examine the adult plant resistance to stripe rust among Chinese endemic wheat germplasms, a panel of 213 accessions was inoculated with mixed virulent races of wheat stripe rust (Puccinia striiformis f. sp. tritici) in four different field environments. Four traits associated with stripe rust resistance, infection type, final disease severity, disease index, and area under the disease progress curve, were used to evaluate the accessions. The phenotypic datasets were used for 55K single-nucleotide polymorphism (SNP) array-based genome-wide association studies to identify effective resistance loci. Eighty-nine accessions with stable resistance were identified in at least three of the four environments by phenotypic evaluation. Eleven markers located on chromosomes 1A, 2B, 5A, 5D, 7B, and 7D by the genome-wide association studies analysis showed significant associations with at least two resistance-associated traits in two of the environments. These loci, corresponding to seven genomic regions based on linkage disequilibrium decay distance, explained 9.3 to 26.0% of the total phenotypic variation. Five quantitative trait loci (QTLs) on chromosomes 1A, 2B, 7B, and 7D overlapped or were in close proximity to previously reported QTLs based on the consensus and physical maps using the reference sequence of bread wheat (IWGSC RefSeq v1.0). The other two QTLs were potential novel QTLs given their physical positions. Haplotype variants of QTL QYr.sicau-2BS showed subspecies-specific inheritance of the stripe rust resistance locus. Resistant loci among Chinese endemic wheat germplasms could be introduced into common wheat cultivars, and the high-confidence SNP markers will aid in marker-assisted selection in breeding for stripe rust disease resistance.

Genome-wide analysis of expression quantitative trait loci (eQTLs) reveals the regulatory architecture of gene expression variation in the storage roots of sweet potato

Dissecting the genetic regulation of gene expression is critical for understanding phenotypic variation and species evolution. However, our understanding of the transcriptional variability in sweet potato remains limited. Here, we analyzed two publicly available datasets to explore the landscape of transcriptomic variations and its genetic basis in the storage roots of sweet potato. The comprehensive analysis identified a total of 724,438 high-confidence single nucleotide polymorphisms (SNPs) and 26,026 expressed genes. Expression quantitative trait locus (eQTL) analysis revealed 4408 eQTLs regulating the expression of 3646 genes, including 2261 local eQTLs and 2147 distant eQTLs. Two distant eQTL hotspots were found with target genes significantly enriched in specific functional classifications. By combining the information from regulatory network analyses, eQTLs and association mapping, we found that IbMYB1-2 acts as a master regulator and is the major gene responsible for the activation of anthocyanin biosynthesis in the storage roots of sweet potato. Our study provides the first insight into the genetic architecture of genome-wide expression variation in sweet potato and can be used to investigate the potential effects of genetic variants on key agronomic traits in sweet potato.

Genome-wide association study for starch content and constitution in sorghum (Sorghum bicolor (L.) Moench)

Starch is the most important component in endosperm of sorghum grain. Usually, two types of starch are present: amylose (AM) and amylopectin (AP). The levels of AM and AP contents play a significant role in the appearance, structure, and quality of sorghum grains and in marketing applications. In the present study, a panel of 634 sorghum (Sorghum bicolor (L.) Moench) accessions were evaluated for starch, AM, and AP contents of grain, which included a mini core collection of 242 accessions from the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) in India, and 252 landraces and 140 cultivars from China. The average starch content was 67.64% and the average AM and AP contents were 20.19 and 79.81%, respectively. We developed a total of 260 000 high-confidence single nucleotide polymorphism (SNP) markers in the panel of 634 accessions of S. bicolor using specific locus amplified fragment sequencing (SLAF-seq). We performed genome-wide association studies (GWAS) of starch, AM, and AM/AP of grain and SNP markers based on a mixed linear model (MLM). In total, 70 significant association signals were detected for starch, AM, and AM/AP ratio of grain with P<4.452×10–7, of which 10 SNPs were identified with significant starch, 51 SNPs were associated with AM, and nine SNPs were associated with the AM/AP ratio. The Gene Ontology (GO) analysis identified 12 candidate genes at five QTLs associated with starch metabolism within the 200-kb intervals, located on chromosomes 1, 5, 6, and 9. Of these genes, Sobic.006G036500.1 encodes peptidyl-prolyl cis-trans-isomerase CYP38 responsible for hexose monophosphate shunt (HMS) and Sobic.009G071800 encodes 6-phospho-fructokinase (PFK), which is involved in the embden-meyerhof pathway (EMP). Kompetitive allele specific PCR (KASP) markers were developed to validate the GWAS results. The C allele is correlated with a high starch content, while the T allele is linked with a low level of starch content, and provides reliable haplotypes for MAS in sorghum quality improvement.

Whole-genome resequencing of Cucurbita pepo morphotypes to discover genomic variants associated with morphology and horticulturally valuable traits

Cucurbita pepo contains two cultivated subspecies, each of which encompasses four fruit-shape morphotypes (cultivar groups). The Pumpkin, Vegetable Marrow, Cocozelle, and Zucchini Groups are of subsp. pepo and the Acorn, Crookneck, Scallop, and Straightneck Groups are of subsp. ovifera. Recently, a de novo assembly of the C. pepo subsp. pepo Zucchini genome was published, providing insights into its evolution. To expand our knowledge of evolutionary processes within C. pepo and to identify variants associated with particular morphotypes, we performed whole-genome resequencing of seven of these eight C. pepo morphotypes. We report for the first time whole-genome resequencing of the four subsp. pepo (Pumpkin, Vegetable Marrow, Cocozelle, green Zucchini, and yellow Zucchini) morphotypes and three of the subsp. ovifera (Acorn, Crookneck, and Scallop) morphotypes. A high-depth resequencing approach was followed, using the BGISEQ-500 platform that enables the identification of rare variants, with an average of 33.5X. Approximately 94.5% of the clean reads were mapped against the reference Zucchini genome. In total, 3,823,977 high confidence single-nucleotide polymorphisms (SNPs) were identified. Within each accession, SNPs varied from 636,918 in green Zucchini to 2,656,513 in Crookneck, and were distributed homogeneously along the chromosomes. Clear differences between subspecies pepo and ovifera in genetic variation and linkage disequilibrium are highlighted. In fact, comparison between subspecies pepo and ovifera indicated 5710 genes (22.5%) with Fst > 0.80 and 1059 genes (4.1%) with Fst = 1.00 as potential candidate genes that were fixed during the independent evolution and domestication of the two subspecies. Linkage disequilibrium was greater in subsp. ovifera than in subsp. pepo, perhaps reflective of the earlier differentiation of morphotypes within subsp. ovifera. Some morphotype-specific genes have been localized. Our results offer new clues that may provide an improved understanding of the underlying genomic regions involved in the independent evolution and domestication of the two subspecies. Comparisons among SNPs unique to particular subspecies or morphotypes may provide candidate genes responsible for traits of high economic importance.

Genome Sequencing Verifies Relapsed Infection of Helicobacter cinaedi.

BackgroundRecurrent infections of Helicobacter cinaedi are often reported, and long-term antimicrobial treatment is empirically recommended to prevent such infections. However, there have been no studies examining whether recurrent infections are relapses of former infections or reinfections with different clones.MethodsA 69-year-old woman presented with recurrent H cinaedi bacteremia-associated cellulitis after a 51-day interval. We isolated 10 colonies from the blood cultures obtained during each of the 2 episodes and subjected them to whole-genome sequencing (WGS). High-confidence single-nucleotide polymorphisms (SNPs) were identified by an assembly based method. Heterogeneous SNP sites were identified by read mapping. The susceptibility of a representative isolate to 14 antimicrobials was also examined.ResultsWhole-genome sequence analysis revealed only 6 SNP sites among the 20 isolates at the whole-genome level. Based on the 6 SNPs, 5 within-host variants (referred to as genotypes) were identified. All 5 genotypes were detected in the first infection; however, only 2 genotypes were detected in the second infection. Although the H cinaedi clone showed a higher minimum inhibitory concentration to fluoroquinolones and macrolides and responsible mutations were identified, none of the 6 SNPs appeared related to additional resistance.ConclusionsThe second infection analyzed here was a relapse of the first infection. A certain level of within-host genomic heterogeneity of the H cinaedi clone was already present in the first infection. Our results suggest the importance of longer treatment courses to eradicate H cinaedi for preventing the relapse of its infection.

Enhancement of Glen Moy x Latham raspberry linkage map using GbS to further understand control of developmental processes leading to fruit ripening

BackgroundThe changing climate is altering timing of key fruit ripening processes and increasing the occurrence of fruit defects. To improve our understanding of the genetic control of raspberry fruit development an enhanced genetic linkage map was developed and used to examine ripening phenotypic data.ResultsIn this study we developed an enhanced genetic linkage map for the raspberry cvs. Glen Moy x Latham reference mapping population using genotyping by sequencing (GbS). Alignment to a newly sequenced draft reference genome of red raspberry, cultivar (cv.) Glen Moy, identified 8019 single nucleotide polymorphisms (SNPs). After stringent filtering to take account of read coverage over all the progeny individuals, association with a single chromosome, heterozygosity and marker regression mapping, 2348 high confidence SNPs were retained and integrated with an existing raspberry genetic map. The linkage map contained many more SNPs segregating in Latham than in Glen Moy. This caused difficulties in quantitative trait loci (QTL) mapping with standard software and a novel analysis based on a hidden Markov model was used to improve the mapping. QTL mapping using the newly generated dense genetic map not only corroborated previously identified genetic locations but also provided additional genetic elements controlling fruit ripening in raspberry.ConclusionThe high-density GbS map located the QTL peaks more precisely than in earlier studies, aligned the QTLs with Glen Moy genome scaffolds, narrowed the range of potential candidate genes to these regions that can be utilised in other populations or in gene expression studies to confirm their role and increased the repertoire of markers available to understand the genetic control of fruit ripening traits.

Diversity and Genome Analysis of Australian and Global Oilseed Brassica napus L. Germplasm Using Transcriptomics and Whole Genome Re-sequencing.

Intensive breeding of Brassica napus has resulted in relatively low diversity, such that B. napus would benefit from germplasm improvement schemes that sustain diversity. As such, samples representative of global germplasm pools need to be assessed for existing population structure, diversity and linkage disequilibrium (LD). Complexity reduction genotyping-by-sequencing (GBS) methods, including GBS-transcriptomics (GBS-t), enable cost-effective screening of a large number of samples, while whole genome re-sequencing (WGR) delivers the ability to generate large numbers of unbiased genomic single nucleotide polymorphisms (SNPs), and identify structural variants (SVs). Furthermore, the development of genomic tools based on whole genomes representative of global oilseed diversity and orientated by the reference genome has substantial industry relevance and will be highly beneficial for canola breeding. As recent studies have focused on European and Chinese varieties, a global diversity panel as well as a substantial number of Australian spring types were included in this study. Focusing on industry relevance, 633 varieties were initially genotyped using GBS-t to examine population structure using 61,037 SNPs. Subsequently, 149 samples representative of global diversity were selected for WGR and both data sets used for a side-by-side evaluation of diversity and LD. The WGR data was further used to develop genomic resources consisting of a list of 4,029,750 high-confidence SNPs annotated using SnpEff, and SVs in the form of 10,976 deletions and 2,556 insertions. These resources form the basis of a reliable and repeatable system allowing greater integration between canola genomics studies, with a strong focus on breeding germplasm and industry applicability.