Library Preparation Research Articles

Detection of Amphibian Species based on Cytochrome b Metabarcoding using High-Throughput Sequencing Approach

Abstract Amphibians are essential in terrestrial and freshwater ecosystems, contributing to multiple ecological functions. The presence of amphibians in a particular area usually signifies a healthy ecosystem as it indicates that the habitat meets the requirements of their life cycles. However, the amphibian population have been declining due to human activities including pollution and habitat destruction. Researchers are using new methods to monitor amphibian diversity to address this decline. Conventional sampling methods were compared with new DNA metabarcoding techniques using advanced sequencing platforms. Therefore, this study aims to analyze the early detection of amphibian species using Cytochrome b (Cytb) gene with metabarcoding methods and next generation sequencing techniques. Thirty-one samples were used in this study and underwent processes such as extraction, amplification using Cytb gene, library preparation, and sequencing. The result showed that the total number of reads was 658,655, with reads in OTUs numbered 474,061. Out of 31 samples, only three species were successfully detected: Fejervarya limnocharis with the highest abundance (20.6%), followed by Polypedates leucomystax and Pulchrana baramica. The results indicate that Cytb has low capabilities in detecting amphibian species, which might be due to its highly conserved region and the lack of databases for this gene. For future studies, an enhanced DNA metabarcoding procedure should be applied, including better choice of genetic marker and bioinformatics pipelines to minimize error and to improve the accuracy and reliability of detecting amphibian’s species using NGS and DNA metabarcoding. This would enhance efforts to protect amphibian biodiversity.

Nanopore ultra-long sequencing and adaptive sampling spur plant complete telomere-to-telomere genome assembly

The pursuit of complete telomere-to-telomere (T2T) genome assembly in plants, challenged by genomic complexity, has been advanced by Oxford Nanopore Technologies (ONT), which offers ultra-long, real-time sequencing. Despite its promise, sequencing length and gap filling remain significant challenges. This study optimized DNA extraction and library preparation, achieving DNA lengths exceeding 485 kb; average N50 read lengths of 80.57 kb, reaching up to 440 kb; and maximum reads of 5.83 Mb. Importantly, we demonstrated that combining ultra-long sequencing and adaptive sampling can effectively fill gaps during assembly, evidenced by successfully filling the remaining gaps of a near-complete Arabidopsis genome assembly and resolving the sequence of an unknown telomeric region in watermelon genome. Collectively, our strategies improve the feasibility of complete T2T genomic assemblies across various plant species, enhancing genome-based research in diverse fields.

Разработка кастомных баркодов для секвенирования на платформе MGI

The MGI (MGI Tech Co. Ltd., China) next-generation sequencing platform, including the DNBSEQ-G50, -G400, and -T7 sequencers, is being actively adopted in research. Despite its widespread adoption, challenges persist in the form of limitations associated with the manufacturer's provided barcode set for library preparation. These limitations include constraints on the number of samples that can be concurrently sequenced, compatibility issues with barcodes from diverse or incomplete sets, and restrictions on the sample ratio. Purpose: to develop a universal method that allows sequencing of up to 252 samples simultaneously on a single sequencer lane, while eliminating barcode-related limitations. We proposed a “quad method” that provides 4 or 4n+2 equilibration of barcodes. This paper also delves into its comprehensive analysis, verification procedures, seamless integration into the sequencing process and validation of the method on the DNBSEQ G-400 platform. The quad method showed efficiency and reliability, allowing sequencing of up to 252 samples simultaneously without compromising data quality. The proposed method optimizes library preparation and improves the flexibility of sequencing on the MGI platform.

Automation of Liquid Handling for Long-Read Library Preparation

Automation of Liquid Handling for Long-Read Library Preparation

Introducing synthetic thermostable RNase inhibitors to single-cell RNA-seq

Single-cell RNA-sequencing (scRNAseq) is revolutionizing biomedicine, propelled by advances in methodology, ease of use, and cost reduction of library preparation. Over the past decade, there have been remarkable technical improvements in most aspects of single-cell transcriptomics. Yet, little to no progress has been made in advancing RNase inhibition despite maintained RNA integrity being critical during cell collection, storage, and cDNA library generation. Here, we demonstrate that a synthetic thermostable RNase inhibitor (SEQURNA) yields single-cell libraries of equal or superior quality compared to ubiquitously used protein-based recombinant RNase inhibitors (RRIs). Importantly, the synthetic RNase inhibitor provides additional unique improvements in reproducibility and throughput, enables new experimental workflows including retained RNase inhibition throughout heat cycles, and can reduce the need for dry-ice transports. In summary, replacing RRIs represents a substantial advancement in the field of single-cell transcriptomics.

Use of synthetic circular RNA spike-ins (SynCRS) for normalization of circular RNA sequencing data.

High-throughput RNA sequencing enables the quantification of transcript abundance and the identification of novel transcripts in biological samples. These include circular RNAs (circRNAs), a family of alternatively spliced RNA molecules that form a continuous loop. However, quantification and comparison of circRNAs between RNA sequencing libraries remain challenging due to confounding errors introduced during exonuclease digestion, library preparation and RNA sequencing itself. Here we describe a set of synthetic circRNA spike-ins-termed 'SynCRS'-that can be added directly to purified RNA samples before exonuclease digestion and library preparation. SynCRS, introduced either individually or in combinations of varying size and abundance, can be integrated into all next-generation sequencing workflows and, critically, facilitate the quantitative calibration of circRNA transcript abundance between samples, tissue types, species and laboratories. Our step-by-step protocol details the generation of SynCRS and guides users on the stoichiometry of SynCRS spike-in to RNA samples, followed by the bioinformatic steps required to facilitate quantitative comparisons of circRNAs between libraries. The laboratory steps to produce the SynCRS require an additional 3 d on top of the high throughput circRNA sequencing and bioinformatics. The protocol is suitable for users with basic experience in molecular biology and bioinformatics.

Updated annotation and meta-analysis of Brugia malayi transcriptomics data reveals consistent transcriptional profiles across time and space with some study-specific differences in adult female worm transcriptional profiles.

Genomics, transcriptomics, and proteomics have significantly advanced our understanding of obligately host-associated microbes, where interrogation of the biology is often limited by the complexity of the biological system and limited tools. This includes the causative agents of many neglected tropical diseases, including filarial nematodes. Therefore, numerous transcriptomics studies have been undertaken on filarial nematodes. Most of these transcriptomics studies focus on Brugia malayi, which causes lymphatic filariasis and is a laboratory model for human filarial disease. Here, we undertook a meta-analysis of the publicly available B. malayi transcriptomics data enabling the direct cross comparison of samples from almost a dozen studies. This reanalysis highlights the consistency of transcriptomics results across many different studies and experimental designs from across the globe for over a decade of research, across many different generations of a sequencing technology, library preparation protocols, and differential expression tools. Males and microfilariae across samples had similar expression profiles. However, female samples were clustered into two differential expression patterns that were significantly different from one another. Largely, we confirm previous results for all studies reanalyzed including tissue-specific gene expression and anti-Wolbachia doxycycline treatment of microfilaria. However, we did not detect previously reported differential expression upon in vitro or in vivo treatment with ivermectin, albendazole, and DEC, instead identifying a consistent lack of transcriptomic change upon exposure to these anthelminthic drugs. Updated annotation has been provided that denotes poorly supported genes including those overlapping rRNAs.

TIRTL-seq: Deep, quantitative, and affordable paired TCR repertoire sequencing.

α/β T cells are key players in adaptive immunity. The specificity of T cells is determined by the sequences of the hypervariable T cell receptor (TCR) α and β chains. Although bulk TCR sequencing offers a cost-effective approach for in-depth TCR repertoire profiling, it does not provide chain pairings, which are essential for determining T cell specificity. In contrast, single-cell TCR sequencing technologies produce paired chain data, but are limited in throughput to thousands of cells and are cost-prohibitive for cohort-scale studies. Here, we present TIRTL-seq (Throughput-Intensive Rapid TCR Library sequencing), a novel approach that generates ready-to-sequence TCR libraries from live cells in less than 7 hours. The protocol is optimized for use with non-contact liquid handlers in an automation-friendly 384-well plate format. Reaction volume miniaturization reduces library preparation costs to <$0.50 per well. The core principle of TIRTL-seq is the parallel generation of hundreds of libraries providing multiple biological replicates from a single sample that allows precise inference of both frequencies of individual clones and TCR chain pairings from well-occurrence patterns. We demonstrate scalability of our approach up to 1 million unique paired αβTCR clonotypes corresponding to over 30 million T cells per sample at a cost of less than $2000. For a sample of 10 million cells the cost is ~$200. We benchmarked TIRTL-seq against state-of-the-art 5'RACE bulk TCR-seq and 10x Genomics Chromium technologies on longitudinal samples. We show that TIRTL-seq is able to quantitatively identify expanding and contracting clonotypes between timepoints while providing accurate TCR chain pairings, including distinct temporal dynamics of SARS-CoV-2-specific and EBV-specific CD8+ T cell responses after infection. While clonal expansion was followed by sharp contraction for SARS-CoV-2 specific TCRs, EBV-specific TCRs remained stable once established. The sequences of both α and β TCR chains are essential for determining T cell specificity. As the field moves towards greater applications in diagnostics and immunotherapy that rely on TCR specificity, we anticipate that our scalable paired TCR sequencing methodology will be instrumental for collecting large paired-chain datasets and ultimately extracting therapeutically relevant information from the TCR repertoire.

Benchmarking sample pooling for epigenomics of natural populations.

DNA methylation (DNAm) is a mechanism for rapid acclimation to environmental conditions. In natural systems, small effect sizes relative to noise necessitates large sampling efforts to detect differences. Large numbers of individually sequenced libraries are costly. Pooling DNA prior to library preparation may be an efficient way to reduce costs and increase sample size, yet there are to date no recommendations in ecological epigenetics research. We test whether pooled and individual libraries yield comparable DNAm signals in a natural system exposed to different pollution levels by generating whole-epigenome data from two invasive molluscs (Corbicula fluminea, Dreissena polymorpha) collected from polluted and unpolluted localities (Italy). DNA of the same individuals were used for pooled and individual epigenomic libraries and sequenced with equivalent resources per individual. We found that pooling effectively captures similar genome-wide and global methylation signals as individual libraries, highlighting that pooled libraries are representative of the global population signal. However, pooled libraries yielded orders of magnitude more data than individual libraries, which was a consequence of higher coverage. We would therefore recommend aiming for a high initial coverage of individual libraries (15×) in future studies. Consequently, we detected many more differentially methylated regions (DMRs) with the pooled libraries and a significantly lower statistical power for regions from individual libraries. Computationally pooled data from the individual libraries produced fewer DMRs and the overlap with wet-lab pooled DMRs was relatively low. We discuss possible causes for discrepancies, list benefits and drawbacks of pooling, and provide recommendations for future epigenomic studies.

PSVI-24 Gene expression profiles of granulosa and peripheral white blood cells from fertile and sub-fertile beef heifers

Abstract Beef heifer fertility issues contribute to a major economic loss in the cow-calf production industry. Therefore, identifying beef heifers with superior genetic potential for improved fertility would increase profitability. This study aimed to identify differences in the transcriptome profiles from granulosa and peripheral white blood cells (pWBCs) of beef heifers with varying reproductive potential. For this, Angus-Simmental crossbred heifers were subjected to an estrus synchronization and fixed-time artificial insemination (AI) protocol (7-D CO-Synch + CIDR) followed by exposure for 60-d to a fertile bull. Depending on the presence or absence of conceptus 120 d post-AI, heifers were classified as fertile (pregnant by AI) or sub-fertile (non-pregnant by AI or bull-breeding). Pregnancies were terminated, and all animals in both groups (fertile, n = 8; and sub-fertile, n = 5) were cycling when the blood and ovaries were collected from each heifer. Total RNA was extracted from the pWBCs and granulosa cells and subjected to library preparation and sequencing on the Nova-Seq platform. The read counts were obtained after data quality control using FastQC v0.11.9 and MultiQC v1.12 and alignment to the Ensemble’s ARS UCD1.2 Bos taurus genome reference using STAR aligner v2.7.5. The filtered data were subjected to differential expression analysis using DESeq2. We identified 1,061 and 72 significantly differentially expressed genes (DEGs) with P-values ≤ 0.05 and absolute (log2 fold change) ≥ 0.5 from pWBC and granulosa, respectively. Notably, 12 targets including 9 protein coding genes (PLCL1, DNER, GNAS, CDH3, PER1, ITGA2B, CXCL12, ENSBTAG00000048613, ENSBTAG00000051519), bta-mir-2887-1, 5-8S-rRNA and U5 were found as DEGs in both the tissues. Based on a differential co-expression analysis using PCIT, we identified GNAS and DNER as hub genes in pWBC and the granulosa cells of the sub-fertile heifer group. The 12 shared genes were over-represented for pathways such as NF kappa B and chemokine signaling, regulation of actin cytoskeleton, and platelet activation. Some of the identified genes have been previously associated with fertility, while others are novel. A detailed understanding of the underlying biological mechanisms of the top genes and a follow-up study with a larger sample size at different time points could validate the candidates identified in this study for their role as potential therapeutic targets.

A retrospective whole genome sequencing of SARS-CoV-2 isolate from respiratory sample of an individual from Ota – Ogun State, Nigeria

Severe acute respiratory syndrome coronavirus – 2 (SARS-CoV-2) which was responsible for the COVID-19 pandemic has now been considered an endemic virus, that will continually produce sporadic outbreaks in different communities around the globe. Although, there are several surveillance studies on SARS-CoV-2 globally, including Nigeria, there is still an important need to understand the uniqueness of the strains of the virus that was circulating immediately after the pandemic, to add to the global database of information that will aid vaccine production and future preparedness against another SARS - related CoV pandemic.Towards the end of the pandemic, between February – August 2022, SARS-CoV-2 was detected in a surveillance project in Ota – Ogun State Nigeria. We carried out a retrospective whole genome sequencing (WGS) analysis of SARS-CoV-2 in the previously reported positive sample, at Inqaba biotech in South Africa. RNA extraction was carried out using the Quick – RNA viral kit (Zymo) and library preparation was done by NEBNext ARTIC SARS-CoV-2 FS Library Prep kit (Illumina) according to the manufacturer's instruction. The WGS was carried out on the NextSeq500 platform by Illumina. The fastq file containing the unassembled raw sequence reads were submitted to NCBI SARS-CoV-2 resources, and a BioProject accession number PRJNA1076330 was issued. The DRAGEN Targeted Microbial, GISAID-CoVsurver mutation, BLAST and MAFFT applications were used for analysis.The spike (s) gene of the study sequence possessed seven mutations (G142D, A163V, V213G, D614G, H655Y, N679K, P681H) and shared 99.4 % identity with those of the Wuhan reference sequence WIV04. The non-structural proteins (NSP) – 7,8,9,10 and 16 shared 100 % identity with bat/Yunnan/RaTG13 (a SARS- related CoV found in bats) sequence on GISAID. Although, the study sequence was obtained from an asymptomatic individual, the S mutations observed are known to be related to virulence, antigenic shift, enhanced transmissibility and host change. Thus, upon successful exploitation of this asymptomatic variant, it may possibly be transformed to a potential candidate for SARS-CoV-2 vaccine in future.

OPTIMIZING ULTRASONIC gDNA FRAGMENTATION FOR NGS SEQUENCING

The fragmentation of genomic DNA (gDNA) is a crucial step in many Illumina library preparation protocols, including the TSO 500 assay, aimed at achieving fragment sizes ranging from 90 to 250 base pairs (bp) for next-generation sequencing (NGS). However, challenges arise in practice due to variations in fragmentation times, which are influenced by different extraction methods and sample sources. To optimize ultrasonic gDNA fragmentation across various extraction methods for NGS sequencing applications. Shearing of gDNA was performed using the Covaris S220 system, employing microtubes containing a total volume of 130 microliters of diluted Tris EDTA (pH 8.0) and 100 ng of gDNA, following the manufacturer's protocol. gDNA samples were extracted using three different kits: AllPrep DNA/RNA FFPE (Qiagen) from FFPE blocks, ReliaPrep FFPE gDNA Miniprep system (Promega) from glass slides, and Gentra Puregene (Qiagen) from cell lines. The shearing times varied depending on the extraction method and sample source. gDNA extracted using the AllPrep DNA/RNA FFPE kit from FFPE blocks required shearing times ranging from 300 to 400 seconds. In contrast, gDNA extracted using the ReliaPrep FFPE gDNA Miniprep system from glass slides exhibited shearing times of approximately 480 to 500 seconds. Lastly, gDNA extracted using the Gentra Puregene kit from cell lines showed shearing times ranging from 150 to 180 seconds. Our study demonstrates the influence of gDNA extraction methods and sample sources on ultrasonic fragmentation efficiency for NGS sequencing. The observed differences in shearing times highlight the need for method optimization to achieve consistent and optimal fragment sizes across diverse sample types. These findings contribute to the development of standardized protocols for gDNA fragmentation, facilitating reliable and reproducible NGS data generation in various research and clinical settings. Further investigations into the underlying factors affecting fragmentation efficiency are warranted to advance our understanding and improve NGS library preparation methodologies.

CLONAL HEMATOPOIESIS AND ITS ROLE IN THE DEVELOPMENT OF CARDIOVASCULAR DISEASES

Clonal hematopoiesis (CH) is the increase in somatic clones carrying mutations in hematopoietic stem cells. This biological condition can lead to cancer and is implicated in the development of atherosclerosis and cardiovascular diseases. CH of indeterminate potential (CHIP) refers to somatic mutations in genes that are associated with the development of Coronary heart disease (CHD) and leukemia. Now, with next-generation sequencing, we can detect mutations even if they are present in only a small number of cells in the tissue sample being studied. Consequently, numerous studies have shown that mutations in the DNMT3A, TET2, ASXL1, and JAK2 genes found in the blood are associated with coronary heart disease. The purpose of our study is to detect a genetic variant of genes in Kazakh patients with atherosclerosis, which will further allow us to study the role of CHIP in increasing the risk of developing coronary microvascular dysfunction and CHD. We perform high-throughput whole-exome sequencing in patients with atherosclerosis. The object of the study is samples delivered from the clinic “National Scientific Cardiac Surgery Center” in Astana. The study is carried out in accordance with the principles of the Declaration of Helsinki of the World Medical Association. Genomic DNA extraction was performed on 245 samples using the Illustra Blood Kit (Cytiva, USA). Qualitative and quantitative assessment of the extracted genomic DNA was conducted using the NanoDrop 2000 (Thermo Scientific) and Qubit 2.0 (Thermo Fisher Scientific). Subsequently, DNA library preparation was carried out using the Illumina DNA Prep with Exome 2.5 Enrichment kit. Quality control of the resulting DNA libraries was conducted on a 2100 Bioanalyzer system using the Agilent DNA 1000 Kit and High Sensitivity DNA kit (Agilent Technologies). Following assessment of the DNA library concentration on Qubit 2.0, they were loaded onto the high-performance sequencer NovaSeq 6000 (Illumina, USA) following the manufacturer’s protocol. We conducted human whole-exome sequencing on 245 samples to identify somatic clones in the peripheral blood of patients with atherosclerosis. Patients were divided into three groups – low risk, middle, and high risk of atherosclerotic complications. Bioinformatic analysis of sequence data were performed. We created list of genes associated with CH to deep analysis and to interpreter the data in patient groups. Subsequently, DNA extraction from the remaining samples is performed, followed by sequencing on the NovaSeq 6000 platform. Analysis of the resulting sequenced data is currently underway. Further research and monitoring of CHIP will enable us to understand the mechanisms underlying the development of cardiovascular diseases and to devise a comprehensive management plan for patients with this pathology. Supported by Committee of Science of the Ministry of Science and Higher Education of Republic of Kazakhstan (AP14869903), (BR21881970) and Nazarbayev University CRP (211123CRP1608).

The Comparative Full-Length Genome Characterization of African Swine Fever Virus Detected in Thailand.

African swine fever virus (ASFV) has been responsible for the globally devastating epidemics in wild and domesticated pigs. Of the 24 identified ASFV genotypes, genotype II is the primary cause for the pandemic occurring in Europe and Asia since its emergence in Georgia in 2007. The current study aimed to characterize the full-length genomic pattern of the ASFV strain from Thailand, TH1_22/CR (Accession No. PP915735), which was then compared with genomic diversity across other Asian isolates using Georgia 2007/1 (Accession No. FR682468) as the reference. Viral DNA was isolated from the pig spleen sample following library preparation and paired-end sequencing using the MiSeq Illumina platform. The sequenced TH1_22/CR isolate spanned 189,395 nucleotides encoding 193 open reading frames (ORFs), exhibiting maximum nucleotide similarity (99.99%) with Georgian (Georgia 2007/1) and Chinese (Wuhan 2019-1 and China HLJ) isolates. Based on phylogenetic analysis, the TH1_22/CR isolate (Accession No. PP915735) was characterized as genotype II, serogroup 8, and IGR-II due to the presence of three tandem repeat sequences (TRSs). Genetic variations including SNPs and single and polynucleotide indels were identified in TH1_22/CR in agreement with other Asian isolates. For comprehensive analysis, the genome was divided into four regions (I-IV) based on gene location. Overall, the TH1_22/CR isolate demonstrated eight SNPs and indels in its genome. Two unique SNPs were reported in the coding regions of the TH1_22/CR isolate, out of which, a C-591-T substitution was seen in MGF 360-4L and a C-297-T was found in A238L, and four unique SNPs were reported in non-coding regions (NCRs). Furthermore, a 29 bp deletion was observed in the IGR between MGF 110-13La and MGF 110-13Lb, as well as 52 bp deletion in the ASFV G ACD 00350 gene. This comparative analysis establishes the foundational information for future studies on the diversity and phylogeography of this regionally significant genetic sub-group of ASFV.

Laser-Assisted Microdissection and High-Throughput RNA Sequencing of the Arabidopsis Gynoecium Medial and Lateral Domains.

For obtaining insights into gene networks during plant reproductive development, having transcriptomes of specific cells from developmental stages as starting points is very useful. During development, there is a balance between cell proliferation and differentiation, and many cell and tissue types are formed. While there is a wealth of transcriptome data available, it is mostly at the organ level and not at specific cell or tissue type level. Therefore, methods to isolate specific cell and tissue types are needed. One method is fluorescent activated cell sorting (FACS), but it has limitations such as requiring marker lines and protoplasting. Recently, single-cell/nuclei isolation methods have been developed; however, a minimum amount of genetic information (marker genes) is needed to annotate/predict the resulting cell clusters in these experiments. Another technique that has been known for some time is laser-assisted microdissection (LAM), where specific cells are microdissected and collected using a laser mounted on a microscope platform. This technique has advantages over the others because no fluorescent marker lines must be made, no marker genes must be known, and no protoplasting must be done. The LAM technique consists in tissue fixation, tissue embedding and sectioning using a microtome, microdissection and collection of the cells of interest on the microscope, and finally RNA extraction, library preparation, and RNA sequencing. In this protocol, we implement the use of normal slides instead of the membrane slides commonly used for LAM. We applied this protocol to obtain the transcriptomes of specific tissues during the development of the gynoecium of Arabidopsis. Key features • Laser-assisted microdissection (LAM) allows the isolation of specific cells or tissues. • Normal slides can be used for LAM. • It allows the identification of the transcriptional profiles of specific tissues of the Arabidopsis gynoecium.

A low-cost dpMIG-seq method for elucidating complex inheritance in polysomic crops: a case study in tetraploid blueberry.

Next-generation sequencing (NGS) library construction often requires high-quality DNA extraction, precise adjustment of DNA concentration, and restriction enzyme digestion to reduce genome complexity, which results in increased time and cost in sample preparation and processing. To address these challenges, a PCR-based method for rapid NGS library preparation, named dpMIG-seq, has been developed and proven effective for high-throughput genotyping. However, the application of dpMIG-seq has been limited to diploid and polyploid species with disomic inheritance. In this study, we obtained genome-wide single nucleotide polymorphism (SNP) markers for tetraploid blueberry to evaluate genotyping and downstream analysis outcomes. Comparison of genotyping qualities inferred across samples with different DNA concentrations and multiple bioinformatics approaches revealed high accuracy and reproducibility of dpMIG-seq-based genotyping, with Pearson's correlation coefficients between replicates in the range of 0.91 to 0.98. Furthermore, we demonstrated that dpMIG-seq enables accurate genotyping of samples with low DNA concentrations. Subsequently, we applied dpMIG-seq to a tetraploid F1 population to examine the inheritance probability of parental alleles. Pairing configuration analysis supported the random meiotic pairing of homologous chromosomes on a genome-wide level. On the other hand, preferential pairing was observed on chr-11, suggesting that there may be an exception to the random pairing. Genotypic data suggested quadrivalent formation within the population, although the frequency of quadrivalent formation varied by chromosome and cultivar. Collectively, the results confirmed applicability of dpMIG-seq for allele dosage genotyping and are expected to catalyze the adoption of this cost-effective and rapid genotyping technology in polyploid studies.

Genome-wide discovery of selection signatures in four Anatolian sheep breeds revealed by ddRADseq

High-density genomic data analyzed by accurate statistical methods are of potential to enlighten past breeding practices such as selection by unraveling fixed regions. In this study, four native Turkish sheep breeds (80 samples) were genotyped via 296.097 single nucleotide polymorphisms (SNPs) detected by double-digest restriction site-associated DNA (ddRADseq) library preparation combined with the Illumina HiSeq X Ten instrument in order to identify genes under selection pressure. A total of 32, 136, 133, and 119 protein-coding genes were detected under selection pressure by runs of homozygosity (ROH), integrated haplotype score (iHS), the ratio of extended haplotype homozygosity (Rsb), and fixation index (FST) approaches, respectively. Of these, a total of 129 genes were identified by at least two statistical models which overlapped with a total of 52 quantitative trait loci (QTL)-associated SNPs, known to be related to fiber diameter, milk content, body weight, carcass traits, some blood parameters, and entropion. A total of six genes under selection pressure were validated by three statistical approaches five of which are of potential to be integrated into animal breeding since they were associated with wool fiber diameter (ZNF208B), behaviors related to neurocognitive development (CBX1 and NFE2L1), adaptation to high-altitude (SDK1), and anxiety causing internal stress (GSG1L). The sixth gene (COPZ1) turned out to play an important role in coping with different types of cancer in mammals. In particular, ROH analysis uncovered significant findings that the Güney Karaman (GKR) had experienced different selection practices than the Akkaraman (AKR) breed. Moreover, some genes specifically under selection in the GKR breed turned out to be associated with olfaction (OR6K6, OR6N1, OR6N2, and OR4C16), survival during the gestation period (PRR15L), and heat stress (CDK5RAP9). The results of this study imply that GKR may become genetically different from the AKR breed at the genome level due to most probably experiencing different adaptation processes occurring in raised climatic conditions. These differences should be conserved to face future challenges, while other native Turkish sheep breeds could be monitored via genome-wide high-density SNP data to obtain deeper knowledge about the effects of natural selection.

Investigation of the Relationship Between DNA Mismatch Repair Genes and Microsatellite Instability in Solid Tumors

Aim: In this study, we aimed to detect the MSI status and somatic mutations in MMR genes (MSH2, MSH6, MLH1, PMS2) of a total of 55 solid tumors diagnosed with colorectal, endometrium and ovarian cancer by NGS method and to reveal the relationship between them. Material and Method: DNA isolation was performed by taking 10-micron sections from paraffin-embedded tissue samples of 55 patients diagnosed with kolorektal, endometrium ve ovarian solid tümörlerinin and Kapa NGS DNA extraction kit was used for sequence analysis. The purity and concentration of the DNA obtained was measured by Qubit fluoremeter, and NadPrep DNA Universal Library Preparation Kit was used for high quality library preparation. Bioinformatics analyses were performed on the Genomize Seq platform. The MSI value was analysed by Roche Navify Mutation Caller software and percentage MSI values were determined using the MSIsensor2 pipeline for secondary analysis of NGS. Results: All ovarian tumors were in the MSI-Stable category. The average MSI value was 2.01. One sample had an MSI of zero. In addition, no mutations in the MSH2 and MLH1 genes were detected in any of the ovarian tumors. 3 of 4 endometrial tumors were in the MSI-Stable category, and 1 tumor was in the MSI-low (MSI value: 13.2) category. No variants were detected in the MSH2 and PMS2 genes in endometrial tumors. 2 of the 41 colorectal tumors (Case4, Case16) were in the MSI-High category. No variants were detected in the MMR genes in 19 tumors. Conclusion: Although frame-shift and stop-gain mutations were detected in 23 tumors that would cause protein deficiency in MMR genes, MSI-H was not detected, as expected, except for two colorectal tumors. Therefore, the results of our study emphasise the need to define new predictive biomarkers clinically within the framework of algorithms to predict response to immunotherapy and determine prognosis.

O51 EXPLORING GENETIC BASIS OF RESISTANT HYPERTENSION IN EUROPEAN AND AFRICAN POPULATIONS

Background: Approximately 15% of hypertensive (HTN) patients experience resistant hypertension (RHTN) despite the use of ≥ 3 antihypertensive agents (diuretic and blockers of calcium channels and the renin-angiotensin system) at maximum tolerated doses. Intriguingly, genetic variation and alterations in gene content may contribute to RHTN developing. Aim: To comprehend the genetic determinants of RHTN across diverse ethnicities, this study has compared the gene content and presence of single nucleotide variants (SNVs) in RHTN versus controlled-HTN and normotensive participants within European (EU) and South African (SA) populations. Methods: Blood-isolated DNA was collected from 6 participants from EU and 15 from SA. Patients were categorized into 3 groups: RHTN (SA; n=5 and EU; n=6), controlled-HTN (n=5), and normotensive (n=5). DNA was used for Whole Exome Sequencing (WES) analysis. Library preparation and exome capture were done using the MGIEasy Universal DNA library prep and exome capture v5 probe set. Mapping was done using GATK, and the variant calling and annotation were performed using with the SIFT tool and R software. Results: PCA analysis and minimum spanning network showed significant genetic diversity between EU and SA populations. WES identified 16 common genes, with 207 and 220 genes associated to RHTN in EU and SA, respectively. These genes were involved in cardiac electrical conduction (K+/Ca++ channels) and lipid storage (HDL deficiency). In addition, SNV within pharmacogenes in SA (CYP2D6, CYP3A5, APOE) and EU (ADD1) cohorts alongside genes linked to disease susceptibility (RyR1 and PPP1CA) were also identified. RyR1 and PPP1CA are both implicated in HTN development through their roles in calcium signaling and vascular smooth muscle cell function. Noteworthy, two SNVs, rs121918006 and rs61732908, within RyR1 were previously linked to HTN. Conclusions: The large genetic differences observed between EU and SA subjects highlight the diverse genetic landscape surrounding RHTN. Notably, specific genes and their related SNV that associated to cardiac calcium channels and lipid storage could show potential targets for personalized therapeutic intervention. Scientific Session 15: Exploring Novel Therapeutic Approaches to Kidney Disease

Visible/near-infrared luminescence and concentration effects of Pr3+-doped Sr2Al2GeO7 downconversion phosphors

The Pr3+ ion has been widely doped into various materials as a red and near-infrared (NIR) emitting center for applications in lighting and solar spectrum downconversion. Herein, the preparation of a new library of Pr3+-doped Sr2Al2GeO7 phosphors was proved by powder x-ray diffraction patterns and Rietveld refinements and characterized by a scanning electron microscope with energy-dispersive x-ray spectrometry. The Sr2Al2GeO7:Pr3+ sample strongly absorbs blue photons over 420–500 nm and yields intense visible emissions with dominant peaks around 490 nm from the Pr3+ 3P0 → 3H4 transition, as well as robust NIR emission bands over 800–1200 nm. In addition to the typical transitions of 1D2 → 3F2 at 880 nm, 1G4 → 3H4 at 1000 nm, and 1D2 → 3F3,4 at 1070 nm, the distinguishable NIR emission at 929 nm was demonstrated from the 3P0 → 1G4 transition via static and dynamic spectroscopic analysis. Most interestingly, for the 3P0 blue-excited state, a considerably elevated concentration of about 10%Pr3+ was optimal for the visible/NIR emissions, in stark contrast to the diluted optimal 1%Pr3+ for the 1D2 state. The relevant cross-relaxation from the 3P0 and 1D2 states between Pr3+ was comprehensively treated by theoretical speculations and experimental results. Such concentrated Pr3+ blue activators would significantly facilitate the blue-to-NIR downconversion through a desired two-step sequential transition from the 3P0 initial state to the 1G4 intermediate level for quantum efficiency exceeding unity. The current results would consolidate the basis of concentrated Pr3+ donors to promote the novel Pr3+/Yb3+ codoping downconversion for greatly increasing Si solar cell efficiency.