SNP Positions Research Articles

Unmutated RRAS2 emerges as a key oncogene in post-partum-associated triple negative breast cancer

BackgroundBreast cancer (BC) is the most common cancer in women, with triple negative BC (TNBC) accounting for 20% of cases. While early detection and targeted therapies have improved overall life expectancy, TNBC remains resistant to current treatments. Although parity reduces the lifetime risk of developing BC, pregnancy increases the risk of developing TNBC for years after childbirth. Although numerous gene mutations have been associated with BC, no single gene alteration has been identified as a universal driver. RRAS2 is a RAS-related GTPase rarely found mutated in cancer.MethodsConditional knock-in mice were generated to overexpress wild type human RRAS2 in mammary epithelial cells. A human sample cohort was analyzed by RT-qPCR to measure RRAS2 transcriptional expression and to determine the frequency of both a single-nucleotide polymorphism (SNP rs8570) in the 3’UTR region of RRAS2 and of genomic DNA amplification in tumoral and non-tumoral human BC samples.ResultsHere we show that overexpression of wild-type RRAS2 in mice is sufficient to develop TNBC in 100% of females in a pregnancy-dependent manner. In human BC, wild-type RRAS2 is overexpressed in 68% of tumors across grade, location, and molecular type, surpassing the prevalence of any previously implicated alteration. Still, RRAS2 overexpression is notably higher and more frequent in TNBC and young parous patients. The increased prevalence of the alternate C allele at the SNP position in tumor samples, along with frequent RRAS2 gene amplification in both tumors and blood of BC patients, suggests a cause-and-effect relationship between RRAS2 overexpression and breast cancer.ConclusionsHigher than normal expression of RRAS2 not bearing activating mutations is a key driver in the majority of breast cancers, especially those of the triple-negative type and those linked to pregnancy.

Measles outbreak in 2019: a warning for the post-COVID-19 pandemic period.

Introduction. At the end of 2019 and the year before, there was a significant spread of measles in the World Health Organization (WHO) European Region.Gap statement. Among the countries that reported, a measles outbreak was Bosnia and Herzegovina (BiH).Aim. To describe the measles outbreak in BiH (an entity of the Federation of BiH, FBiH) in 2019.Methodology. Confirmatory IgM serology, measles nucleic acid detection by real-time RT-PCR and virus genotyping were done in the WHO-accredited laboratory for measles and rubella at the Clinical Center of the University of Sarajevo, Unit for Clinical Microbiology. Genotype was determined in all measles-RNA-positive cases by sequence analysis of the 450 nt fragment coding the C-terminal of measles virus nucleoprotein (N).Results. From 1 January to 31 December 2019, 1332 measles cases were reported, with the peak observed in April 2019 (413/1332, 31.01 %). Sarajevo Canton had the highest incidence, number of cases and percentage (206.4; 868/1332; 65.17 %) of measles cases. Around four-fifths of infected persons were unvaccinated (1086/1332, 81.53 %), while 4.58 % of the patients (61/1332) were immunized with one dose of measles-containing vaccine. The highest proportion of cases was found in children 0-6 years of age (738/1332, 55.41 %). Measles IgM positivity was determined in 75.88 % (346/456), while virus RNA was detected in 82.46 % (47/57) of the swab samples. All measles virus sequences belonged to genotype B3. SNP (position 216: C=>T) was detected in 1 of the 40 sequences obtained during this outbreak.Conclusion. Due to suboptimal immunization coverage, BiH belongs to countries at a high risk for measles outbreaks. Post-COVID-19 (coronavirus disease 2019) pandemic, targeted and tailored strategies are required to ensure routine vaccination demand and acceptance and broad partner and stakeholder group participation.

Zinc finger knuckle genes are associated with tolerance to drought and dehydration in chickpea (Cicer arietinum L.).

Chickpea (Cicer arietinum L.) is a very important food legume and needs improved drought tolerance for higher seed production in dry environments. The aim of this study was to determine diversity and genetic polymorphism in zinc finger knuckle genes with CCHC domains and their functional analysis for practical improvement of chickpea breeding. Two CaZF-CCHC genes, Ca04468 and Ca07571, were identified as potentially important candidates associated with plant responses to drought and dehydration. To study these genes, various methods were used including Sanger sequencing, DArT (Diversity array technology) and molecular markers for plant genotyping, gene expression analysis using RT-qPCR, and associations with seed-related traits in chickpea plants grown in field trials. These genes were studied for genetic polymorphism among a set of chickpea accessions, and one SNP was selected for further study from four identified SNPs between the promoter regions of each of the two genes. Molecular markers were developed for the SNP and verified using the ASQ and CAPS methods. Genotyping of parents and selected breeding lines from two hybrid populations, and SNP positions on chromosomes with haplotype identification, were confirmed using DArT microarray analysis. Differential expression profiles were identified in the parents and the hybrid populations under gradual drought and rapid dehydration. The SNP-based genotypes were differentially associated with seed weight per plant but not with 100 seed weight. The two developed and verified SNP molecular markers for both genes, Ca04468 and Ca07571, respectively, could be used for marker-assisted selection in novel chickpea cultivars with improved tolerance to drought and dehydration.

An SNP based genotyping assay for genes associated with drought tolerance in bread wheat.

SnRK2 plays vital role in responding to adverse abiotic stimuli. The applicability of TaSnRK2.4 and TaSnRK2.9 was investigated to leverage the potential of these genes in indigenous wheat breeding programs. Genetic diversity was assessed using pre-existing markers for TaSnRK2.4 and TaSnRK2.9. Furthermore, new markers were also developed to enhance their broader applicability. KASP markers were designed for TaSnRK2.4, while CAPS-based markers were tailored for TaSnRK2.9. Analysis revealed lack of polymorphism in TaSnRK2.4 among Pakistani wheat germplasm under study. To validate this finding, available gel-based markers for TaSnRK2.4 were employed, producing consistent results and offering limited potential for application in marker-assisted wheat breeding with Pakistani wheat material. For TaSnRK2.9-5A, CAPS2.9-5A-1 and CAPS2.9-5A-2 markers were designed to target SNP positions at 308 nt and 1700 nt revealing four distinct haplotypes. Association analysis highlighted the significance of Hap-5A-1 of TaSnRK2.9-5A, which exhibited association with an increased number of productive tillers (NPT), grains per spike (GPS), and reduced plant height (PH) under well-watered (WW) conditions. Moreover, it showed positive influence on NPT under WW conditions, GPS under water-limited (WL) conditions, and PH under both WW and WL conditions. High selection intensity observed for Hap-5A-1 underscores the valuable role it has played in Pakistani wheat breeding programs. Gene expression studies of TaSnRK2.9-5A revealed the involvement of this gene in response to PEG, NaCl, low temperature and ABA treatments. These findings propose that TaSnRK2.9 can be effectively employed for improving wheat through marker-assisted selection in wheat breeding efforts.

An artificial-restriction fragment length polymorphism (A-RFLP) method for genotyping intronic SNP rs7587633 C/T of ATG16L1

SNPs could either cause a disorder or directly alter the efficacy of a particular treatment and act as biological markers. The SNP rs7587633 C/T present in the intronic region of the ATG16L1 gene has been studied for its role in psoriasis vulgaris and Palmoplantar pustulosis. To genotype rs7587633 C/T using PCR-RFLP no restriction site is present for any of the restriction enzymes at the SNP position. To develop an artificial-RFLP method for genotyping rs7587633 C/T, the forward primer was designed in such a way that it resulted in the creation of an EcoRI restriction site in the amplified product which could further be digested with EcoRI to find the genotype of the individual. The newly developed A-RFLP method was applied to genotype the SNP rs7587633 C/T in DNA samples of 100 healthy control individuals. The allelic and genotypic frequencies of the SNPs were 0.80(C), 0.20(T) and 65%(CC), 31%(CT) and 4%(TT), respectively. In conclusion, we developed an A-RFLP method to genotype the SNP rs7587633 C/T which is not present in any of the natural restriction sites and this method could be applied to genotype this SNP in various populations/diseases to find its role.

Development of a next generation SNP genotyping array for wheat.

High-throughput genotyping arrays have provided a cost-effective, reliable and interoperable system for genotyping hexaploid wheat and its relatives. Existing, highly cited arrays including our 35K Wheat Breeder's array and the Illumina 90K array were designed based on a limited amount of varietal sequence diversity and with imperfect knowledge of SNP positions. Recent progress in wheat sequencing has given us access to a vast pool of SNP diversity, whilst technological improvements have allowed us to fit significantly more probes onto a 384-well format Axiom array than previously possible. Here we describe a novel Axiom genotyping array, the 'Triticum aestivum Next Generation' array (TaNG), largely derived from whole genome skim sequencing of 204 elite wheat lines and 111 wheat landraces taken from the Watkins 'Core Collection'. We used a novel haplotype optimization approach to select SNPs with the highest combined varietal discrimination and a design iteration step to test and replace SNPs which failed to convert to reliable markers. The final design with 43 372 SNPs contains a combination of haplotype-optimized novel SNPs and legacy cross-platform markers. We show that this design has an improved distribution of SNPs compared to previous arrays and can be used to generate genetic maps with a significantly higher number of distinct bins than our previous array. We also demonstrate the improved performance of TaNGv1.1 for Genome-wide association studies (GWAS) and its utility for Copy Number Variation (CNV) analysis. The array is commercially available with supporting marker annotations and initial genotyping results freely available.

Genome-wide association studies of salinity tolerance in local aman rice.

The present study aimed to identify and characterize new sources of salt tolerance among 94 rice varieties from varied geographic origins. The genotypes were divided into five groups based on their morphological characteristics at both vegetative and reproductive stages using salinity scores from the Standard Evaluation System (SES). The experiment was designed as per CRD (Completely Randomized Design) with 2 sets of salinity treatments for 8 dS/meter and 12 dS/meter, respectively compared with one non-salinized control set. Using a Soil Plant Analysis Development (SPAD) meter, assessments of the apparent chlorophyll content (greenness) of the genotypes were done to comprehend the mechanism underlying their salt tolerance. To evaluate molecular genetic diversity, a panel of 1 K RiCA SNP markers was employed. Utilizing TASSEL 5.0 software, 598 filtered SNPs were used for molecular analysis. Whole-genome association studies (GWAS) were also used to investigate panicle number per plant (pn, tiller number per plant (till), SPAD value (spad), sterility (percent) (str), plant height (ph) and panicle length (pl. It is noteworthy that these characteristics oversee conveying the visible signs of salt damage in rice. Based on genotype data, diversity analysis divided the germplasm groups into four distinct clusters (I, II, III and IV). For the traits studied, thirteen significant marker-trait associations were discovered. According to the phenotypic screening, seven germplasm genotypes namely Koijuri, Asha, Kajal, Kaliboro, Hanumanjata, Akundi and Dular, are highly tolerant to salinity stress. The greenness of these genotypes was found to be more stable over time, indicating that these genotypes are more resistant to stress. Regarding their tolerance levels, the GWAS analysis produced comparable results, supporting that salinity-tolerant genotypes having minor alleles in significant SNP positions showed more greenness during the stress period. The Manhattan plot demonstrated that at the designated significant SNP position, the highly tolerant genotypes shared common alleles. These genotypes could therefore be seen as important genomic resources for accelerating the development and release of rice varieties that are tolerant to salinity.

Genome-wide association study (GWAS) reveals key loci associated with rhizoctonia aerial blight resistance in soybean

Rhizoctonia Aerial Blight (RAB) is caused by Rhizoctonia solani is one of the most worldwide devastating diseases in soybean. In present study, a GWAS was conducted for rhizoctonia aerial blight (RAB) resistance by using 269 diverse soybean genotypes for two consecutive year 2020 & 2021. GWAS analysis revealed two most significant SNPs, located very close each other on chromosome 12 (3193687 & 3193885) and 16 (6074131 & 6076330). Some candidate genes such as related to leucine-rich repeats containing (LRRNT; Glyma.16g064200), protein ubiquitination (F-box-like) and amino acid transmembrane transport (Trp_Tyr_perm; Glyma.16g062500 & Glyma.16g062600) were identified near the significant SNP positions. Our findings may be relevant for future studies and provide insights into the genetic basis for RAB resistance in soybean.

Genome-wide association studies on charcoal rot resistance in soybean (Glycine max, L.)

A genome-wide association study (GWAS) was conducted for Charcoal rot resistance by using 267 soybean germplasm lines. GWAS analysis revealed two most significant SNPs at chromosome 3 (S3_8215775 and S3_8309181) for multiple traits associated with adult plant resistance. Candidate genes related to calcium ion signaling, response to nematode infestation, abscisic acid stimulus and defense response to bacteria were identified within the flanking regions of the significant SNP positions. Among all, one most important gene (Glyma.03G057100) involved in jasmonic acid biosynthesis was identified. Two SNPs each on chromosome 8 (S8_22467783 and S8_22467802) and chromosome 9 (S9_2923863 and S9_2935518), associated with seedlingresistance have also been identified.

Simultaneous genome-wide haplotyping and copy number detection enabling PGT-A and PGT-M testing in equine in vitro produced embryos

The number of equine embryos being produced increases yearly (37,094 +10% in 2021; newsletter IETS December 2022). The reasons include the growing popularity of in vitro embryo production (IVP), which increases the genetic yield from high quality mares, maximizes the use of expensive straws of semen from high quality stallions and enables IVP to be performed outside the reproductive season. On the 4th October 2022, a record price of €95,000 was paid for an embryo. However, despite the embryos’ promising pedigrees, it remains unknown if transfer of such embryos will result in a healthy foal with the aspired phenotypical characteristics. This is largely because in horses, preimplantation genetic testing for aneuploidies and monogenic diseases or traits (PGT-A and PGT-M) remains underdeveloped especially when compared to humans. In addition, equine IVP has been reported to double early pregnancy losses, possibly because of an increased incidence of genetic errors. Here, we explored the possibility of concurrent genome-wide PGT-A and PGT-M in horses and assessed the incidence of chromosomal aberrations in IVP arrested or blastocyst stage embryos. We used GGP Equine SNP arrays (Illumina) with over 70,000 SNP positions to analyze whole-genome amplified DNA (Repli-G Single cell kit, Qiagen) from 14 trophectoderm biopsies and corresponding whole blastocysts and 26 individual cells from six embryos arrested at cleavage stage following IVP. Haplarithmisis (Zamani–Esteki et al. Am. J. Hum.Genet. 2015; 96:894-912) was applied to simultaneously determine the genome-wide copy number and haplotype. Five biopsies and seven single cells failed to give a result. From six arrested embryos, five had chromosomal errors in most of their cells, including (combinations of) triploidy and one or more monosomies and trisomies of meiotic and / or mitotic origin. Two biopsies showed a loss of maternal chromosome 28 and 31 respectively, which was confirmed in the corresponding blastocysts, suggesting meiotic errors. In one of those biopsies, additional complex aneuploidies not present in the blastocyst were found, suggesting a mitotic origin. Thus, aneuploidy and ploidy errors are highly prevalent in IVP arrested embryos compared to blastocysts. Like other mammals, these errors may contribute to embryonic arrest before or after clinical pregnancy detection. As a proof of principle, we determined the inherited haplotypes for four common disease-associated genes with high carrier frequency in different horse breeds (GBE1, PLOD1, B3GALNT2, MUTYH), and for one color coat-associated gene (STX17) in one biopsy-embryo combination, confirming the potential of generic genome-wide PGT-M concurrent with PGT-A testing. Similar to humans, we postulate that embryo selection will improve embryo transfer success rates in horses and enable selection for traits and against diseases thus accelerating breeding programs.

Abstract 3124: Identification of sub-clonal loss of heterozygosity (LOH) in glioblastomas analyzed with spatial transcriptomics

Abstract We propose a method to leverage spatial transcriptomics and bioinformatic analysis to identify sub-clonal loss of heterozygosity that could otherwise be missed by standard bulk sequencing in glioblastoma tumors. Glioblastomas are highly invasive and aggressive brain tumors with a low five-year survival rate of less than 10%. Current treatment strategies are not sufficient for long term disease outcome, and innovative approaches to analyze glioblastoma are imperative for future improved options. Many factors contribute to difficulty of treatment, though one critical consideration is the diffuse nature and known intratumoral heterogeneity of this cancer. It is important to determine LOH heterogeneity, as this genomic alteration is associated with copy number variation and specific classifications of gliomas. Loss of heterozygosity is an irreversible deletion of DNA that is a common alteration in cancer. Application of spatial transcriptomics can increase our understanding of genomic complexity across a tissue, and unsupervised graph clustering can computationally separate spatial regions with similar gene expression patterns. Our method investigates spatial clusters individually at known heterozygous SNP sites to identify evidence for LOH with Bayesian inference and a hidden Markov model. Our analysis centers on the balance of reference (A) and non-reference (B) spatial transcriptomics read counts at SNPs of interest. B allele frequency is the number of reads aligned to the alternative allele divided by the total number of reads counted at the SNP. At a normal two-copy state at a heterozygous site, B allele frequency is ideally 50 percent. LOH can be identified as this value deviates toward 0 or 100, and copy number alterations impact this frequency. Our method utilizes Bayes factor values, or likelihood ratios between two models, using allele counts to provide support for either LOH or a heterozygous event. For our analysis, spatial transcriptomics BAMs are processed with a graph clustering algorithm and split into cluster-specific BAMs. We calculate allele count coverage for each cluster at pre-selected heterozygous SNP positions identified through germline exome sequencing. Bayes factor values are then calculated at each SNP, resulting in evidence for LOH or a heterozygous event. Chromosomal regions are segmented by a hidden Markov approach. Cumulative metrics for each segment are evaluated to determine a “state” label of heterozygous, LOH, or undefined. We tested our method on three WHO grade 4 IDH wild-type EGFRvIII positive fresh frozen glioblastoma tissue samples. Sub-clonal heterogeneity was identified in two samples and confirmed by bulk exome sequencing, bulk RNA sequencing, and analysis of initial clinical reports. Citation Format: Michelle G. Webb, Frances Chow, Carmel McCullough, Bohan Zhang, Kyle Hurth, John Carpten, Gabriel Zada, David W. Craig. Identification of sub-clonal loss of heterozygosity (LOH) in glioblastomas analyzed with spatial transcriptomics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3124.

QTL-Seq Approach Identified Pi63 Conferring Blast Resistance at the Seedling and Tillering Stages of Thai Indigenous Rice Variety “Phaladum”

Rice blast (BL) caused by Magnaporthe oryzae is a fungal disease causing significant yield losses in rice production worldwide. To overcome the breakdown of resistance by the rapid adaptation of pathogens, identifying resistance (R) genes or QTLs in indigenous rice, which harbors the R genes that co-evolved with the local pathogen race, is necessary. In this study, a recombinant inbred line (RIL) population derived from a cross between RD6 and Phaladum (PLD) was used to map quantitative trait loci (QTL) for BL resistance through a QTL-seq approach. A single QTL (qBLchr4) associated with BL resistance at the seedling and maximum tillering stages was mapped on the long arm of chromosome 4. Five genes, LOC_Os04g0616600, LOC_Os04g0617900 (OsGLP4-1), LOC_Os04g0619600 (OsRLCK161), LOC_Os04g0620800 (Pi63), and LOC_Os04g0621500, were considered the candidate genes representing qBLchr4. Subsequently, the Kompetitive Allele-Specific PCR (KASP) markers specific for the SNP variant and position of each gene were designed for validation in the mapping population. These markers showed the high phenotypic variance explained (PVE) values in all testing methods and/or environments, signifying the major effect of qBLchr4. Among these markers, the Pi63-KASP marker explained the highest and most stable phenotypic variation across all testing methods and/or environments, with 84.18%, 80.34%, and 23.43% in the upland short row (USR) method, Sila environment, and Mueang environment, respectively. Therefore, Pi63 was suggested to be the strongest candidate gene. These results represent the potential utility of future BL resistance breeding and/or pyramiding using marker-assisted selection (MAS).

Genomic Mutations of BK Polyomavirus in Patients after Kidney Transplantation: A Cross-Sectional Study in Vietnam

Objectives: The purpose of this study was to identify the SNP sites and determine the BKV genotype circulating in kidney-transplant Vietnamese recipients based on the VP1 gene region. Methods: 344 samples were collected from post-kidney-transplant recipients at the 103 Vietnam Military Hospital to investigate the number of BKV infections. Positive samples with a sufficient virus concentration were analyzed by nested PCR in the VP1 region, sequencing detected genotyping and single-nucleotide polymorphism. Results: BKV infection was determined in 214 patients (62.2%), of whom 11 (5.1%) were diagnosed with BKV-associated nephropathy. Among the 90 BKV-I strains sequenced, 89 (98.88%) were strains of I/b-1 and 1 (1.12%) was strain I/b-2. The 60 BKV-IV strains had a greater diversity of subgroups, including 40% IV/a-1, 1.66% IV/a-2, 56.68% IV/c-1, and 1.16% IV/c-2. Additionally, of 11 cases diagnosed with BKVN, seven belonged to subgroup I/b-1 (63.6%) and four to subgroup IV/c-1 (36.4%). Moreover, 22 specific SNPs that were genotype I or IV were determined in this Vietnamese population. Specifically, at position 1745, for the Vietnamese BKV-IV strains, the SNP position (A→G) appeared in 57/60 samples (95%). This causes transformation of the amino acid N→S. This SNP site can enable detection of genotype IV in Vietnam. It represents a unique evolution pattern and mutation that has not been found in other international strains. Conclusion: The BKV-I genotype was more common than BKV-IV; however, mutations that occur on the VP1 typing region of BKV-IV strains were more frequent than in BKV-I strains.

Genome-Wide Genomic and Functional Association Study for Workability and Calving Traits in Holstein Cattle.

Simple SummaryTraits related to the calving abilities and workability of cows are essential both for the dairy economy and for the welfare of cows. We analysed the direct and maternal components of calving ease, direct and maternal components of stillbirth, milking speed, and temperament by looking for SNPs, gene ontologies, and metabolic pathways significantly associated with the additive genetic variability of those traits. At the genetic level, 150 significant SNPs were identified whereas at the functional level we found 45 significant gene ontology terms and 5 metabolic pathways.The goal of our study was to identify the SNPs, metabolic pathways (KEGG), and gene ontology (GO) terms significantly associated with calving and workability traits in dairy cattle. We analysed direct (DCE) and maternal (MCE) calving ease, direct (DSB) and maternal (MSB) stillbirth, milking speed (MSP), and temperament (TEM) based on a Holstein-Friesian dairy cattle population consisting of 35,203 individuals. The number of animals, depending on the trait, ranged from 22,301 bulls for TEM to 30,603 for DCE. We estimated the SNP effects (based on 46,216 polymorphisms from Illumina BovineSNP50 BeadChip Version 2) using a multi-SNP mixed model. The SNP positions were mapped to genes and the GO terms/KEGG pathways of the corresponding genes were assigned. The estimation of the GO term/KEGG pathway effects was based on a mixed model using the SNP effects as dependent variables. The number of significant SNPs comprised 59 for DCE, 25 for DSB and MSP, 17 for MCE and MSB, and 7 for TEM. Significant KEGG pathways were found for MSB (2), TEM (2), and MSP (1) and 11 GO terms were significant for MSP, 10 for DCE, 8 for DSB and TEM, 5 for MCE, and 3 for MSB. From the perspective of a better understanding of the genomic background of the phenotypes, traits with low heritabilities suggest that the focus should be moved from single genes to the metabolic pathways or gene ontologies significant for the phenotype.

Whole Exome-Sequencing of Pooled Genomic DNA Samples to Detect Quantitative Trait Loci in Esotropia and Exotropia of Strabismus in Japanese.

Background: Esotropia and exotropia are two major phenotypes of comitant strabismus. It remains controversial whether esotropia and exotropia would share common genetic backgrounds. In this study, we used a quantitative trait locus (QTL)-sequencing pipeline for diploid plants to screen for susceptibility loci of strabismus in whole exome sequencing of pooled genomic DNAs of individuals. Methods: Pooled genomic DNA (2.5 ng each) of 20 individuals in three groups, Japanese patients with esotropia and exotropia, and normal members in the families, was sequenced twice after exome capture, and the first and second sets of data in each group were combined to increase the read depth. The SNP index, as the ratio of variant genotype reads to all reads, and Δ(SNP index) values, as the difference of SNP index between two groups, were calculated by sliding window analysis with a 4 Mb window size and 10 kb slide size. The rows of 200 “N”s were inserted as a putative 200-b spacer between every adjoining locus to depict Δ(SNP index) plots on each chromosome. SNP positions with depth < 20 as well as SNP positions with SNP index of <0.3 were excluded. Results: After the exclusion of SNPs, 12,242 SNPs in esotropia/normal group and 12,108 SNPs in exotropia/normal group remained. The patterns of the Δ(SNP index) plots on each chromosome appeared different between esotropia/normal group and exotropia/normal group. When the consecutive groups of SNPs on each chromosome were set at three patterns: SNPs in each cytogenetic band, 50 consecutive sliding SNPs, and SNPs in 4 Mb window size with 10 kb slide size, p values (Wilcoxon signed rank test) and Q values (false discovery rate) in a few loci as Manhattan plots showed significant differences in comparison between the Δ(SNP index) in the esotropia/normal group and exotropia/normal group. Conclusions: The pooled DNA sequencing and QTL mapping approach for plants could provide overview of genetic background on each chromosome and would suggest different genetic backgrounds for two major phenotypes of comitant strabismus, esotropia and exotropia.

A linkage disequilibrium-based approach to position unmapped SNPs in crop species

BackgroundHigh-density SNP arrays are now available for a wide range of crop species. Despite the development of many tools for generating genetic maps, the genome position of many SNPs from these arrays is unknown. Here we propose a linkage disequilibrium (LD)-based algorithm to allocate unassigned SNPs to chromosome regions from sparse genetic maps. This algorithm was tested on sugarcane, wheat, and barley data sets. We calculated the algorithm’s efficiency by masking SNPs with known locations, then assigning their position to the map with the algorithm, and finally comparing the assigned and true positions.ResultsIn the 20-fold cross-validation, the mean proportion of masked mapped SNPs that were placed by the algorithm to a chromosome was 89.53, 94.25, and 97.23% for sugarcane, wheat, and barley, respectively. Of the markers that were placed in the genome, 98.73, 96.45 and 98.53% of the SNPs were positioned on the correct chromosome. The mean correlations between known and new estimated SNP positions were 0.97, 0.98, and 0.97 for sugarcane, wheat, and barley. The LD-based algorithm was used to assign 5920 out of 21,251 unpositioned markers to the current Q208 sugarcane genetic map, representing the highest density genetic map for this species to date.ConclusionsOur LD-based approach can be used to accurately assign unpositioned SNPs to existing genetic maps, improving genome-wide association studies and genomic prediction in crop species with fragmented and incomplete genome assemblies. This approach will facilitate genomic-assisted breeding for many orphan crops that lack genetic and genomic resources.

Genome Wide Association Study of Rice (Oryza sativa L.) during Heading Stage under a High Temperature

At the reproductive development stage of rice (Oryza sativa L.), temperature stress can decrease spikelet fertility, ultimately resulting in a yield loss. In this study, a total of 98 rice varieties were used in genome-wide association study (GWAS) to understand spikelet fertility under a high temperature (SFHT). GWAS results revealed that two lead SNPs were significantly associated with SFHT. Candidate genes located within ± 250 kb of the corresponding SNP position were discovered, resulting in a total of 21 candidate genes on chromosome 10 and 18 candidate genes on chromosome 11. Based on previously reported function and haplotype analysis, Os10g0177200 (EF-HAND 2domain containing protein) as one candidate gene showed significant differences among groups of haplotypes. This candidate gene will be further evaluated for its function to determine whether it is useful for improving molecular breeding studies and developing new high temperature tolerant rice varieties.

Genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan.

Tularemia is a highly dangerous zoonotic infection due to the bacteria Francisella tularensis. Low genetic diversity promoted the use of polymorphic tandem repeats (MLVA) as first-line assay for genetic description. Whole genome sequencing (WGS) is becoming increasingly accessible, opening the perspective of a time when WGS might become the universal genotyping assay. The main goal of this study was to describe F. tularensis strains circulating in Kazakhstan based on WGS data and develop a MLVA assay compatible with in vitro and in silico analysis. In vitro MLVA genotyping and WGS were performed for the vaccine strain and for 38 strains isolated in Kazakhstan from natural water bodies, ticks, rodents, carnivores, and from one migratory bird, an Isabellina wheatear captured in a rodent burrow. The two genotyping approaches were congruent and allowed to attribute all strains to two F. tularensis holarctica lineages, B.4 and B.12. The seven tandem repeats polymorphic in the investigated strain collection could be typed in a single multiplex PCR assay. Identical MLVA genotypes were produced by in vitro and in silico analysis, demonstrating full compatibility between the two approaches. The strains from Kazakhstan were compared to all publicly available WGS data of worldwide origin by whole genome SNP (wgSNP) analysis. Genotypes differing at a single SNP position were collected within a time interval of more than fifty years, from locations separated from each other by more than one thousand kilometers, supporting a role for migratory birds in the worldwide spread of the bacteria.

Patterns in Genotype Composition of Indian Isolates of the Bombyx mori Nucleopolyhedrovirus and Bombyx mori Bidensovirus.

The mulberry silkworm, Bombyx mori (L.), is a model organism of lepidopteran insects with high economic importance. The viral diseases of the silkworm caused by Bombyx mori nucleopolyhedrovirus (BmNPV) and Bombyx mori bidensovirus (BmBDV) inflict huge economic losses and significantly impact the sericulture industry of India and other countries. To understand the distribution of Indian isolates of the BmNPV and to investigate their genetic composition, an in-depth population structure analysis was conducted using comprehensive and newly developed genomic analysis methods. The seven new Indian BmNPV isolates from Anantapur, Dehradun, Ghumarwin, Jammu, Kashmir, Mysore and Salem grouped in the BmNPV clade, and are most closely related to Autographa californica multiple nucleopolyhedrovirus and Rachiplusia ou multiple nucleopolyhedrovirus on the basis of gene sequencing and phylogenetic analyses of the partial polh, lef-8 and lef-9 gene fragments. The whole genome sequencing of three Indian BmNPV isolates from Mysore (-My), Jammu (-Ja) and Dehradun (-De) was conducted, and intra-isolate genetic variability was analyzed on the basis of variable SNP positions and the frequencies of alternative nucleotides. The results revealed that the BmNPV-De and BmNPV-Ja isolates are highly similar in their genotypic composition, whereas the population structure of BmNPV-My appeared rather pure and homogenous, with almost no or few genetic variations. The BmNPV-De and BmNPV-Ja samples further contained a significant amount of BmBDV belonging to the Bidnaviridae family. We elucidated the genotype composition within Indian BmNPV and BmBDV isolates, and the results presented have broad implications for our understanding of the genetic diversity and evolution of BmNPV and co-occurring BmBDV isolates.

DyPS: Dynamic, Private and Secure GWAS

Abstract Genome-Wide Association Studies (GWAS) identify the genomic variations that are statistically associated with a particular phenotype (e.g., a disease). The confidence in GWAS results increases with the number of genomes analyzed, which encourages federated computations where biocenters would periodically share the genomes they have sequenced. However, for economical and legal reasons, this collaboration will only happen if biocenters cannot learn each others’ data. In addition, GWAS releases should not jeopardize the privacy of the individuals whose genomes are used. We introduce DyPS, a novel framework to conduct dynamic privacy-preserving federated GWAS. DyPS leverages a Trusted Execution Environment to secure dynamic GWAS computations. Moreover, DyPS uses a scaling mechanism to speed up the releases of GWAS results according to the evolving number of genomes used in the study, even if individuals retract their participation consent. Lastly, DyPS also tolerates up to all-but-one colluding biocenters without privacy leaks. We implemented and extensively evaluated DyPS through several scenarios involving more than 6 million simulated genomes and up to 35,000 real genomes. Our evaluation shows that DyPS updates test statistics with a reasonable additional request processing delay (11% longer) compared to an approach that would update them with minimal delay but would lead to 8% of the genomes not being protected. In addition, DyPS can result in the same amount of aggregate statistics as a static release (i.e., at the end of the study), but can produce up to 2.6 times more statistics information during earlier dynamic releases. Besides, we show that DyPS can support a larger number of genomes and SNP positions without any significant performance penalty.