SNP Array Research Articles

Decoding the fish genome opens a new era in important trait research and molecular breeding in China.

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Genomic profiles and clinical presentation of chordoma

Chordoma is a rare bone cancer with variable clinical outcomes. Here, we recruited 184 sporadic chordoma patients from the US and Canada and collected their clinical and treatment data. The average age at diagnosis was 45.5 years (Range 5–78) and the chordoma site distribution was 49.2% clivus, 26.2% spinal, and 24.0% sacral. Most patients (97.5%) received surgery as the primary treatment, among whom 85.3% also received additional treatment. Except for the most prevalent cancers like prostate, lung, breast, and skin cancer, there was no discernible enrichment for any specific cancer type among patients or their family members. Among a subset of patients (N = 70) with tumor materials, we conducted omics analyses and obtained targeted panel sequencing and SNP array genotyping data for 51 and 49 patients, respectively. The most recurrent somatic driver mutations included PIK3CA (12%), followed by chromatin remodeling genes PBRM1 and SETD2. Amplification of the 6q27 region, containing the chordoma susceptibility gene TBXT, was detected in eight patients (16.3%). Clival patients appeared to be less likely to carry driver gene mutations, chromosome arm level deletion events (e.g., 5p, 5p, and 9p), or 6q27 amplification compared to sacral patients. After adjusting for age, sex, tumor site, and additional treatment, patients with somatic deletions of 14q (OR = 13.73, 95% CI 1.96–96.02, P = 0.008) and 18p (OR = 13.68, 95% CI 1.77–105.89, P = 0.012) were more likely to have persistent chordoma. The study highlights genomic heterogeneity in chordoma, potentially linked to location and clinical progression.

Interplay between polygenic risk score and solar insolation: Implication for systemic lupus erythematosus diagnosis and pathogenesis

ObjectivesThis research elucidates the correlation between solar radiation insolation, polygenic risk score (PRS), and systemic lupus erythematosus (SLE) diagnosis, utilizing genomic, environmental, and clinical data. MethodsWe included 1,800 SLE participants and 1,800 controls from the Taiwan Precision Medicine Initiative, genotyped via the Affymetrix Genome-Wide TWB 2.0 SNP Array. The study employed a SLE-PRS tailored for individuals of Taiwanese ancestry, comprising 27 single nucleotide polymorphisms (SNPs). QGIS computed solar radiation insolation from participants' residences. We employed logistic regression to investigate the associations between SLE-PRS, solar insolation susceptibility, and SLE. Additive and multiplicative interactions were utilized to assess the interactions between solar insolation and SLE-PRS regarding the risk of SLE. ResultsSLE patients showed decreased solar insolation (p < 0.001). The highest decile of SLE-PRS exhibited a statistically significant lower solar insolation 1, 3, 6, and 12 months prior to diagnosis as compared to the lowest decile. Specifically, there were significant differences observed at 1 and 12 months (p = 0.025 and p = 0.004, respectively). It suggests that higher SLE-PRS correlated with reduced solar insolation tolerance. We observed an increase in SLE risk across ascending SLE-PRS percentiles exclusively in the high solar insolation group, not in the low solar insolation group. However, the interaction effect of SLE-PRS and solar insolation on SLE risk is not statistically significant. Compared to the lowest decile, the highest SLE-PRS decile showed a 10.98-fold increase in SLE risk (95 % CI, 3.773–31.952, p < 0.001). High SLE-PRS scores in conjunction with high solar insolation contribute to SLE incidence. ConclusionsOur study unveils the intertwined nature of UV insolation and polygenic risks in SLE. Future studies should explore the preventative potential of robust solar radiation protection for high-risk individuals before the disease onset.

High‐density single nucleotide polymorphism‐based linkage mapping reveals novel quantitative trait loci for rice (Oryza sativa) sheath blight resistance from wild rice (Oryza rufipogon)‐derived introgression line

AbstractRice sheath blight is one of the most important fungal disease‐causing yield losses up to 50%. An Oryza rufipogon‐derived introgression line (Pusa 1908‐13‐12‐5) identified to be highly resistant to sheath blight was used for developing a RIL population by crossing with a highly susceptible rice cultivar, Pusa Basmati 1 (PB 1). The RIL population was genotyped with 80K Rice Pangenome (RPAN) single nucleotide polymorphism (SNP) array and phenotyped for disease response under artificial inoculation for 2 consecutive years. Inclusive composite interval mapping could identify six quantitative trait loci (QTLs) for two traits, namely, relative lesion height and relative yield loss per plant. Among these QTLs, four were major QTLs and three are novel QTLs. Three major QTLs were co‐localised within 1.4 Mb region on chromosome 1. In silico analysis revealed the presence of 12 candidate genes with potential role in disease resistance. Kompetitive Allele Specific PCR (KASP™) assay was developed for the linked SNPs and used for validation in an independent BC1F2:3 population (PB 1121/Pusa 1908‐13‐12‐5). One SNP (AX‐182186537) on chromosome 1 was validated, which offers scope for its utilisation in marker‐assisted breeding.

Whole paternal uniparental disomy of chromosome 4 with a novel homozygous IDUA splicing variant, c.159-9T>A, in a Chinese patient with mucopolysaccharidosis type I.

Mucopolysaccharidosis type I (MPS-I) is a rare autosomal recessive genetic lysosomal storage disorder that is caused by pathogenic variants of the α-L-iduronidase (IDUA) gene. This study aimed to identify the genetic causes of MPS-I in a Chinese patient and construct a minigene of IDUA to analyze its variants upon splicing. Whole-exome sequencing (WES) and Sanger sequencing were used to confirm the potential causative variants. Single-nucleotide polymorphism (SNP) array was subsequently performed to confirm uniparental disomy (UPD). Minigene assay was performed to analyze the effect on splicing of mRNA. We meanwhile explored the conservative analysis and protein homology simulation. A novel homozygous splicing mutation of IDUA, c.159-9T>A, was identified in an individual presenting with overlapping features of MPS-I. Interestingly, only the father and sisters, but not the mother, carried the variant in a heterozygous state. WES and SNP array analyses validated paternal UPD on chromosome 4. Minigene splicing revealed two aberrant splicing events: exon 2 skipping and intron 1 retention. Moreover, the specific structure of the mutant protein obviously changed according to the results of the homologous model. This study describes a rare autosomal recessive disorder with paternal UPD of chromosome 4 leading to the homozygosity of the IDUA splicing variant in patients with MPS-I for the first time. This study expands the variant spectrum of IDUA and provides insights into the splicing system, facilitating its enhanced diagnosis and treatment.

Structural Variants in the SMC1A Gene Associated With Near-Haploidy in Undifferentiated Pleomorphic Sarcomas.

Near-haploidization, that is, loss of one copy of most chromosomes, is a relatively rare phenomenon in most tumors, but is enriched among certain soft tissue sarcomas, including undifferentiated pleomorphic sarcoma (UPS). Presumably, near-haploidization can arise through many mechanisms. This study aimed to identify gene rearrangements that could cause near-haploidization. We here present two UPS in which near-haploidization was an early event, identified through single nucleotide polymorphism (SNP) array analysis. One of the cases was studied further using whole genome and transcriptome sequencing, as well as cytogenetic and molecular cytogenetic methods. Both tumors had chromosomal rearrangements in the form of copy number shifts/structural variants affecting the SMC1A gene. These findings suggest that cohesin defects could contribute to mitotic errors resulting in massive loss of chromosomes. SMC1A encodes one of the components of the cohesin multiprotein complex, which is critical for proper alignment of the sister chromatids during S-phase and separation to opposite spindle poles. Further studies should explore the role of cohesin defects in near-haploidization in other sarcomas and to clarify its role in tumor development.

Application effectiveness of "Youxin-1" in genetic diversity and structure analysis of local chickens.

China's local chicken breeds are rich in resources, and have formed different germplasm characteristics in the process of long-term selection and evolution. Scientific assessment of population genetic diversity and identification of inter-breed genetic structure are of great value to the protection and innovative utilization of local chicken breed resource. In order to evaluate the application effectiveness of 23K SNP chip "Youxin-1" in the analysis of genetic diversity and genetic structure of local chickens, we used RADseq to identify genomic genetic variation of 21 local chicken breeds and developed 23K chip "Youxin-1". The genetic statistics of each variety were calculated based on two sets of SNP data, and correlation, fitting and phylogenetic analysis were carried out to evaluate the application effectiveness of the chip. The results showed that the observed heterozygosity (Ho), polymorphism information content (PIC), inbred coefficient (FROH) and genetic differentiation coefficient (Fst) calculated based on the two SNP data sets were basically consistent in the 21 local chicken breeds. The genetic diversity of Langya chicken (LA), Piao chicken (PJ) and Wenchang chicken (WC) was relatively rich. The genetic diversity of Bian chickens (BJ), Langshan chickens (LS), Gushi chickens (GS), Dongxiang blue-eggshell chickens (DX) and Beijing fatty chickens (BY) was relatively poor, and the correlation coefficients of Ho, PIC, FROH and average Fst in the two groups were 0.794, 0.901, 0.926 and 0.984, respectively, all reaching extremely significant levels (P<0.01) with a high degree of fit (P<0.001) and R2 were 0.644, 0.827, 0.916 and 0.927. For the two sets of SNP data, the evolutionary tree constructed by neighbor-joining (NJ) method and maximum likelihood (ML) method was reasonable, and the 21 local chicken breeds were generally divided into six categories, which was consistent with the formation history and geographical distribution of the varieties. The 23K chip also realized reasonable clustering of the five new varieties without individual deviation. There are some differences in the estimation of genetic statistics using SNP with different densities, and data standardization is needed. 23K chip has good efficacy in the analysis of genetic diversity and structure of local chickens.

To be or not to be tetraploid-the impact of marker ploidy on genomic prediction and GWAS of potato.

Cultivated potato, Solanum tuberosum L., is considered an autotetraploid with 12 chromosomes with four homologous phases. However, recent evidence found that, due to frequent large phase deletions in the genome, gene ploidy is not constant across the genome. The elite cultivar "Otava" was found to have an average gene copy number of 3.2 across all loci. Breeding programs for elite potato cultivars rely increasingly on genomic prediction tools for selection breeding and elucidation of quantitative trait loci underpinning trait genetic variance. These are typically based on anonymous single nucleotide polymorphism (SNP) markers, which are usually called from, for example, SNP array or sequencing data using a tetraploid model. In this study, we analyzed the impact of using whole genome markers genotyped as either tetraploid or observed allele frequencies from genotype-by-sequencing data on single-trait additive genomic best linear unbiased prediction (GBLUP) genomic prediction (GP) models and single-marker regression genome-wide association studies of potato to evaluate the implications of capturing varying ploidy on the statistical models employed in genomic breeding. A panel of 762 offspring of a diallel cross of 18 parents of elite breeding material was used for modeling. These were genotyped by sequencing and phenotyped for five key performance traits: chipping quality, length/width ratio, senescence, dry matter content, and yield. We also estimated the read coverage required to confidently discriminate between a heterozygous triploid and tetraploid state from simulated data. It was found that using a tetraploid model neither impaired nor improved genomic predictions compared to using the observed allele frequencies that account for true marker ploidy. In genome-wide associations studies (GWAS), very minor variations of both signal amplitude and number of SNPs supporting both minor and major quantitative trait loci (QTLs) were observed between the two data sets. However, all major QTLs were reproducible using both data sets.

Development and identification of a novel wheat-Thinopyrum ponticum disomic substitution line DS5Ag(5D) with new genes conferring resistance to powdery mildew and leaf rust

BackgroundPowdery mildew (caused by Blumeria graminis f. sp. tritici (Bgt)) and leaf rust (caused by Puccinia triticina (Pt)) are prevalent diseases in wheat (Triticum aestivum L.) production. Thinopyrum ponticum (2n = 10x = 70, EeEeEbEbExExStStStSt) contains genes that confer high levels of resistance to these diseases.ResultsAn elite wheat-Th. ponticum disomic substitution line, DS5Ag(5D), was developed in the Bainong Aikang 58 (AK58) background. The line was assessed using genomic in situ hybridization (GISH), oligo-nucleotide probe multiplex (ONPM) fluorescence in situ hybridization (FISH), and molecular markers. Twenty eight chromosome-specific molecular markers were identified for the alien chromosome, and 22 of them were co-dominant. Additionally, SNP markers from the wheat 660 K SNP chip were utilized to confirm chromosome identification and they provide molecular tools for tagging the chromosome in concern. The substitution line demonstrated high levels of resistance to powdery mildew throughout its growth period and to leaf rust at the adult stage. Based on the resistance evaluation of five F5 populations between the substitution lines and wheat genotypes with different levels of sensitivity to the two diseases. Results showed that the resistance genes located on 5Ag confered stable resistance against both diseases across different backgrounds. Resistance spectrum analysis combined with diagnostic marker detection of known resistance genes of Th. ponticum revealed that 5Ag contained two novel genes, Pm5Ag and Lr5Ag, which conferred resistance to powdery mildew and leaf rust, respectively.ConclusionsIn this study, a novel wheat-Th. ponticum disomic substitution line DS5Ag(5D) was successfully developed. The Th. ponticum chromosome 5Ag contain new resistance genes for powdery mildew and leaf rust. Chromosomic—specific molecular markers were generated and they can be used to track the 5Ag chromosome fragments. Consequently, this study provides new elite germplasm resources and molecular markers to facilitate the breeding of wheat varieties that is resistant to powdery mildew and leaf rust.

Review on camel genetic diversity: ecological and economic perspectives.

Camels, known as the "Ship of the Desert," play a vital role in the ecosystems and economies of arid and semi-arid regions. They provide meat, milk, transportation, and other essential services, and their resilience to harsh environments makes them invaluable. Despite their similarities, camel breeds exhibit notable differences in size, color, and structure, with over 40million camels worldwide. This number is projected to increase, underscoring their growing significance. Economically, camels are crucial for food production, tourism, and trade, with camel racing being particularly significant in Arab countries. Their unique physiological traits, such as low disease susceptibility and efficient water conservation, further enhance their value. Camel products, especially meat and milk, offer substantial nutritional and therapeutic benefits, contributing to their high demand. Genetic diversity studies have advanced our understanding of camels' adaptation to extreme environments. Functional genomics and whole-genome sequencing have identified genes responsible for these adaptations, aiding breeding programs and conservation efforts. High-throughput sequencing has revealed genetic markers linked to traits like milk production and disease resistance. The development of SNP chips has revolutionized genetic studies by providing a cost-effective alternative to whole-genome sequencing. These tools facilitate large-scale genotyping, essential for conserving genetic diversity and improving breeding strategies. To prevent the depletion of camel genetic diversity, it is crucial to streamline in situ and ex situ conservation efforts to maintain their ecological and economic value. A comprehensive approach to camel conservation and genetic preservation, involving advanced genomic technologies, reproductive biotechniques, and sustainable management practices, will ensure their continued contribution to human societies.

Development and validation of a 40 K liquid SNP array for the mud crab (Scylla paramamosain)

Genomic selection (GS) and genome-wide association studies (GWAS) have been widely implemented to accelerate the genetic improvement of important economic traits. Developing an efficient and reliable high-throughput genotyping tool is the premise and basis of implementing GS and GWAS on the mud crab (Scylla paramamosain). In this study, leveraging the whole-genome resequencing data of 146 mud crabs from the southeast coast of China, a total of 16,714,943 high-quality SNPs were discovered. Based on polymorphic information content and evenly distribution on 49 chromosomes, a total of 40,352 SNPs was selected to develop a 40 K liquid SNP array (named “Xiexin No. 1”) using the genotyping by target sequencing (GBTS) technology. The genotyping performance test showed that this SNP array had a high site polymorphism (over 92.70%) and individual calling rate (97.02%). Furthermore, the mean of SNP concordance, genotype concordance, non-reference sensitivity, and non-reference discrepancy between “Xiexin No. 1” SNP array and resequencing were identified to be 0.986, 0.954, 0.942, and 0.097, respectively. Using this SNP array data, the genomic heritabilities were evaluated as 0.887 ± 0.065, 0.643 ± 0.108, 0.599 ± 0.111, and 0.624 ± 0.111 for body weight (BW), carapace length (CL), carapace length (CW), and body height (BH) in a cultured population (233 crabs). Furthermore, five significant SNPs associated with growth traits were identified by GWAS. After annotation, one candidate gene (XP_027221997.1) associated with BW and one candidate gene (XP_027210446.1) associated with CW and BH were identified. Most importantly, the prediction accuracy of GBLUP for BW, BH, CL, and CW were 0.535 ± 0.064, 0.544 ± 0.101, 0.574 ± 0.087, and 0.585 ± 0.083, with the unbiasedness values of 0.988 ± 0.22, 1.015 ± 0.274, 1.011 ± 0.196 and 1.037 ± 0.235. In conclusion, the “Xiexin No. 1” SNP array provides an efficient and reliable platform for high-quality genotyping, which would greatly promote the genetic studies and thus accelerate the genetic improvement of mud crabs.

Asian Screening Array and Next-Generation Sequencing Based Panels Applied to Preimplantation Genetic Testing for Monogenic Disorders Preclinical Workup in 294 Families: A Retrospective Analysis.

Currently, the most commonly used methods for linkage analysis of pre-implantation genetic testing for monogenic disorders (PGT-M) are next generation sequencing (NGS) and SNP array. We aim to investigate whether the application efficacy of Asian screening array (ASA) in PGT-M preclinical workup for the Chinese population is superior to NGS based single nucleotide polymorphism (SNP) panels. We conducted a retrospective analysis by reviewing 294 couples from a single center over the past 4years and compared the detection results between NGS-based SNP panels and ASA. Using the numbers of informative SNPs upstream and downstream flanking of variants, we assessed the detection efficiency of both methods in monogenic diseases, chromosomal microdeletion syndrome and males with de novo variants, among other scenarios. Results indicate that ASA offers a greater number of informative SNPs compared with NGS-based SNP panels. Additionally, data analysis for ASA is generally more straightforward and may require less computational resources. While ASA can address most PGT-M challenges, we have also identified certain genes in previous tests that are not suitable for PGT-M using ASA. The application of ASA in PGT-M preclinical workup for Chinese populations has good practical value as it can perform linkage analysis for most genetic variants. However, for certain variants, NGS or other testing methods, such as mutated allele revealed by sequencing with aneuploidy and linkage analysis (MARSALA), may still be necessary for completion.

Detection and evaluation of parameters influencing the identification of heterozygous-enriched regions in Holstein cattle based on SNP chip or whole-genome sequence data

BackgroundA heterozygous-enriched region (HER) is a genomic region with high variability generated by factors such as balancing selection, introgression, and admixture processes. In this study, we evaluated the genomic background of HERs and the impact of different parameters (i.e., minimum number of SNPs in a HER, maximum distance between two consecutive SNPs, minimum length of a HER, maximum number of homozygous allowed in a HER) and scenarios [i.e., different SNP panel densities and whole-genome sequence (WGS)] on the detection of HERs. We also compared HERs characterized in Holstein cattle with those identified in Angus, Jersey, and Norwegian Red cattle using WGS data.ResultsThe parameters used for the identification of HERs significantly impact their detection. The maximum distance between two consecutive SNPs did not impact HERs detection as the same average of HERs (269.31 ± 787.00) was observed across scenarios. However, the minimum number of markers, maximum homozygous markers allowed inside a HER, and the minimum length size impacted HERs detection. For the minimum length size, the 10 Kb scenario showed the highest average number of HERs (1,364.69 ± 1,483.64). The number of HERs decreased as the minimum number of markers increased (621.31 ± 1,271.83 to 6.08 ± 21.94), and an opposite pattern was observed for the maximum homozygous markers allowed inside a HER (54.47 ± 195.51 to 494.89 ± 1,169.35). Forty-five HER islands located in 23 chromosomes with high Tajima’s D values and differential among the observed and estimated heterozygosity were detected in all evaluated scenarios, indicating their ability to potentially detect regions under balancing selection. In total, 3,440 markers and 28 genes previously related to fertility (e.g., TP63, ZSCAN23, NEK5, ARHGAP44), immunity (e.g., TP63, IGC, ARHGAP44), residual feed intake (e.g., MAYO9A), stress sensitivity (e.g., SERPINA6), and milk fat percentage (e.g., NOL4) were identified. When comparing HER islands among breeds, there were substantial overlaps between Holstein with Angus (95.3%), Jersey (94.3%), and Norwegian Red cattle (97.1%), indicating conserved HER across taurine breeds.ConclusionsThe detection of HERs varied according to the parameters used, but some HERs were consistently identified across all scenarios. Heterozygous genotypes observed across generations and breeds appear to be conserved in HERs. The results presented could serve as a guide for defining HERs detection parameters and further investigating their biological roles in future studies.

Genetic Analysis of Stripe Rust Resistance in the Chinese Wheat Cultivar Luomai 163.

Stripe or yellow rust (YR) caused by Puccinia striiformis tritici (Pst) is an important foliar disease affecting wheat production globally. Resistant varieties are the most economically and environmentally effective way to manage this disease. The common winter wheat (Triticum aestivum L.) cultivar Luomai 163 exhibited resistance to Pst races CYR32 and CYR33 at the seedling stage and showed a high level adult plant resistance in the field. To understand the genetic basis of YR resistance in this cultivar, 142 F5 recombinant inbred lines (RILs) derived from cross Apav#1 × LM163 and both parents were genotyped with the 16K SNP array and bulked segregant analysis sequencing (BSA-Seq). The analysis detected a major gene, YrLM163, at the seedling stage associated with the 1BL.1RS translocation. Additionally, three genes for resistance at the adult plant stage were detected on chromosome arms 1BL (Lr46/Yr29/Pm39/Sr58), 6BS and 6BL in Luomai 163, whereas Apav#1 contributed resistance at a QTL on 2BL. These QTL explained YR disease severity variations ranging from 6.9 to 54.8%. KASP markers KASP-2BL, KASP-6BS and KASP-6BL for three novel loci QYr.hzau-2BL, QYr.hzau-6BS and QYr.hzau-6BL were developed and validated. QYr.hzau-1BL, QYr.hzau-2BL and QYr.hzau-6BS showed varying degrees of resistance to YR when present individually or in combination based on genotype and phenotype analysis of a panel of 570 wheat accessions. Six RILs combining resistance alleles of all QTL, showing higher resistance to YR in the field than Luomai 163 with disease severities of 10.7-16.0%, are important germplasm resources for breeding programs to develop YR resistant wheat varieties with good agronomic traits.

Imputation strategies for low-coverage whole-genome sequencing data and their effects on genomic prediction and genome-wide association studies in pigs

The uncertainty resulting from missing genotypes in low-coverage whole-genome sequencing (LCWGS) data complicates genotype imputation. The aim of this study is to find out an optimal strategy for accurately imputing LCWGS data and assess its effectiveness for genomic prediction (GP) and genome-wide association study (GWAS) on economically important traits of Large White pigs. The LCWGS data of 1 423 Large White pigs were imputed using three different strategies: (1) using the high-coverage whole-genome sequencing (HCWGS) of 30 key progenitors as the reference panel (Ref_LG); (2) mixing HCWGS of key progenitors with LCWGS (Mix_HLG) and (3) self-imputation in LCWGS (Within_LG). Additionally, to compare the imputation effects of LCWGS, we also imputed SNP chip data of 1 423 Large White pigs to the whole-genome sequencing level using the reference panel consisting of key progenitors (Ref_SNP). To evaluate effects of the imputed sequencing data, we compared the accuracies of GP and statistical power of GWAS for four reproductive traits based on the chip data, sequencing data imputed from chip data and LCWGS data using an optimal strategy. The average imputation accuracies of the Within_LG, Ref_LG and Mix_HLG were 0.9893, 0.9899 and 0.9875, respectively, which were higher than that of the Ref_SNP (0.8522). Using the imputed sequencing data from LCWGS with the Ref_LG imputation strategy, the accuracies of GP for four traits improved by approximately 0.31–1.04% compared to the chip data, and by 0.7–1.05% compared to the imputed sequencing data from chip data. Furthermore, by using the sequence data imputed from LCWGS with the Ref_LG, 18 candidate genes were identified to be associated with the four reproductive traits of interest in Large White pigs: total number of piglets born - EPC2, MBD5, ORC4 and ACVR2A; number of piglets born healthy - IKBKE; total litter weight of piglets born alive - HSPA13 and CPA1; gestation length - GTF2H5, ITGAV, NFE2L2, CALCRL, ITGA4, STAT1, HOXD10, MSTN, COL5A2 and STAT4. With the exception of EPC2, ORC4, ACVR2A and MSTN, others represent novel candidates. Our findings can provide a reference for the application of LCWGS data in livestock and poultry.

Accurate genotype imputation from low-coverage whole-genome sequencing data of rainbow trout.

With the rapid and significant cost reduction of next-generation sequencing, low-coverage whole-genome sequencing (lcWGS), followed by genotype imputation, is becoming a cost-effective alternative to single-nucleotide polymorphism (SNP)-array genotyping. The objectives of this study were 2-fold: (1) construct a haplotype reference panel for genotype imputation from lcWGS data in rainbow trout (Oncorhynchus mykiss); and (2) evaluate the concordance between imputed genotypes and SNP-array genotypes in 2 breeding populations. Medium-coverage (12×) whole-genome sequences were obtained from a total of 410 fish representing 5 breeding populations with various spawning dates. The short-read sequences were mapped to the rainbow trout reference genome, and genetic variants were identified using GATK. After data filtering, 20,434,612 biallelic SNPs were retained. The reference panel was phased with SHAPEIT5 and was used as a reference to impute genotypes from lcWGS data employing GLIMPSE2. A total of 90 fish from the Troutlodge November breeding population were sequenced with an average coverage of 1.3×, and these fish were also genotyped with the Axiom 57K rainbow trout SNP array. The concordance between array-based genotypes and imputed genotypes was 99.1%. After downsampling the coverage to 0.5×, 0.2×, and 0.1×, the concordance between array-based genotypes and imputed genotypes was 98.7, 97.8, and 96.7%, respectively. In the USDA odd-year breeding population, the concordance between array-based genotypes and imputed genotypes was 97.8% for 109 fish downsampled to 0.5× coverage. Therefore, the reference haplotype panel reported in this study can be used to accurately impute genotypes from lcWGS data in rainbow trout breeding populations.

Genome-wide association study identifies loci and candidate genes for RVA parameters in wheat (Triticum aestivum L.).

The gelatinization and retrogradation characteristics of wheat starch affect the eating quality of Chinese-style food. Rapid Visco Analyzer (RVA) parameters have been widely used as important indicators to evaluate and improve the quality of wheat starch. However, the genetic basis of RVA parameters remains to be further explored. In the present study, a natural population was genotyped using 90K single nucleotide polymorphism (SNP) arrays, and the RVA parameters of this population grown in five environments were evaluated. The results showed that 22,068 high-quality SNP markers were identified and distributed unequally on the chromosomes. According to the genetic distance, 214 wheat materials were divided into four groups. Except for the pasting temperature (PTT), six parameters followed a normal distribution. Based on the general linear model, 969 significant association SNPs were detected by genome-wide association studies (GWAS), and chromosomes 7A and 2B had the most associated SNPs. Breakdown viscosity (BV) was associated with the most SNPs (n = 238), followed by PTT (n = 186), peak viscosity (PV; n = 156), trough viscosity (TV; n = 127), and final viscosity (FV; n = 126). According to the average linkage disequilibrium (LD), 33 stable quantitative trait loci (QTLs) were identified for single parameters in multiple environments, of which 12 were associated with BV, followed by peak time (PT; n = 8) and PTT (n = 7). On the other hand, 67 pleiotropic QTLs were identified for multiple parameters. Three candidate genes-TasbeIIa, TasbeI, and TassIIa-were screened for phenotyping analysis. The grain width and the weight of the TasbeIIa and TaSSIIa knockout (KO) lines were significantly lower than those of the TasbeI KO lines and the control (CK). The KO lines had smaller endosperm cells, smaller A-type starch granules, and higher amylose content. The TasbeI KO lines showed normal RVA curves, while the TasbeIIa KO lines showed flat curves. However, the TaSSIIa lines failed to paste under the RVA temperatures. Conclusively, the SNPs/QTLs significantly associated with the RVA parameters and genetic resources with novel haplotypes could be used to improve the quality of wheat starch.

The addition of Psathyrostachys Huashanica Keng 6Ns large segment chromosomes has positive impact on stripe rust resistance and plant spikelet number of common wheat

BackgroundDeveloping novel germplasm by using wheat wild related species is an effective way to rebuild the wheat resource bank. The Psathyrostachys huashanica Keng (P. huashanica, 2n = 2x = 14, NsNs) is regarded as a superior species to improve wheat breeding because of its multi-resistance, early maturation and numerous tiller traits. Introducing genetic components of P. huashanica into the common wheat background is the most important step in achieving the effective use. Therefore, the cytogenetic characterization and influence of the introgressed P. huashanica large segment chromosomes in the wheat background is necessary to be explored.ResultsIn this study, we characterized a novel derived line, named D88-2a, a progeny of the former characterized wheat-P. huashanica partial amphiploid line H8911 (2n = 7x = 49, AABBDDNs). Cytological identification showed that the chromosomal composition of D88-2a was 2n = 44 = 22II, indicating the addition of exogenous chromosomes. Genomic in situ hybridization demonstrated that the supernumerary chromosomes were a pair of homologues from the P. huashanica and could be stably inherited in the common wheat background. Molecular markers and 15 K SNP array indicated that the additional chromosomes were derived from the sixth homoeologous group (i.e., 6Ns) of P. huashanica. Based on the distribution of the heterozygous single-nucleotide polymorphism sites and fluorescence in situ hybridization karyotype of each chromosome, this pair of additional chromosomes was confirmed as P. huashanica 6Ns large segment chromosomes, which contained the entire short arm and the proximal centromere portion of the long arm. In terms of the agronomic traits, the addition line D88-2a exhibited enhanced stripe rust resistance, improved spike characteristics and increased protein content than its wheat parent line 7182.ConclusionsThe new wheat germplasm D88-2a is a novel cytogenetically stable wheat-P. huashanica 6Ns large segment addition line, and the introgressed large segment alien chromosome has positive impact on plant spikelet number and stripe rust resistance. Thus, this germplasm can be used for genetic improvement of cultivated wheat and the study of functional alien chromosome segment.

Functional ex vivo DNA fibre assay to measure replication dynamics in breast cancer tissue.

Replication stress (RS) is a key trait of cancer cells, and a potential actionable target in cancer treatment. Accurate methods to measure RS in tumour samples are currently lacking. DNA fibre analysis has been used as a common technique to measure RS in cell lines. Here, we investigated DNA fibre analysis on fresh breast cancer specimens and correlated DNA replication kinetics to known RS markers and genomic alterations. Fresh, treatment-naïve primary breast cancer samples (n = 74) were subjected to ex vivo DNA fibre analysis to measure DNA replication kinetics. Tumour cell proliferation was confirmed by EdU incorporation and cytokeratin AE1/AE3 (CK) staining. The RS markers phospho-S33-RPA and γH2AX and the RS-inducing proto-oncogenes Cyclin E1 and c-Myc were analysed by immunohistochemistry. Copy number variations (CNVs) were assessed from genome-wide single nucleotide polymorphism (SNP) arrays. We found that the majority of proliferating (EdU-positive) cells in each sample were CK-positive and therefore considered to be tumour cells. DNA fibre lengths varied largely in most tumour samples. The median DNA fibre length showed a significant inverse correlation with pRPA expression (r = -0.29, p = 0.033) but was not correlated with Cyclin E1 or c-Myc expression and global CNVs in this study. Nuclear Cyclin E1 expression showed a positive correlation with pRPA levels (r = 0.481, p < 0.0001), while cytoplasmic Cyclin E1 expression exhibited an inverse association with pRPA expression (r = -0.353, p = 0.002) and a positive association with global CNVs (r = 0.318, p = 0.016). In conclusion, DNA fibre analysis performed with fresh primary breast cancer samples is feasible. Fibre lengths were associated with pRPA expression. Cyclin E1 expression was associated with pRPA and the percentage of CNVs. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Mapping of Aegilops speltoides derived leaf rust and stripe rust resistance genes using 35K SNP array

Wheat is an essential food commodity cultivated throughout the world. However, this crop faces continuous threats from fungal pathogens, leaf rust (LR) and stripe rust (YR). To continue feeding the growing population, these major destructors of wheat must be effectively countered by enhancing the genetic diversity of cultivated germplasm. In this study, an introgression line with hexaploid background (ILsp3603) carrying resistance against Pt pathotypes 77−5 (121R63-1), 77−9 (121R60-1) and Pst pathotypes 46S119 (46E159), 110S119 (110E159), 238S119 (238E159) was developed from donor wheat wild progenitor, Aegilops speltoides acc pau 3603. To understand the genetic basis of resistance and map these genes (named Lrsp3603 and Yrsp3603), inheritance studies were carried out in F6 and F7 mapping population, developed by crossing ILsp3603 with LR and YR susceptible cultivar WL711, which revealed a monogenic (single gene) inheritance pattern for each of these traits. Bulk segregant analysis combined with 35 K Axiom SNP array genotyping mapped both genes as separate entities on the short arm of chromosome 6B. A genetic linkage map, comprising five markers, 1 SNP, 1 PLUG and three gene based SSRs, covered a genetic distance of 12.65 cM. Lrsp3603 was flanked by markers Tag-SSR14 (located proximally at 2.42 cM) and SNP AX-94542331 (at 3.28 cM) while Yrsp3603 was mapped at one end closest to AX-94542331 at 6.62 cM distance. Functional annotation of Lrsp3603 target region (∼ 1 Mbp) revealed 10 gene IDs associated with disease resistance mechanisms including three encoding typical R gene domains.