DNA Sequencing Methods Research Articles

Prevalence, risk factors, and BRAF mutation of colorectal sessile serrated lesions among Vietnamese patients.

Sessile serrated lesions (SSLs) are considered precancerous colorectal lesions that should be detected and removed to prevent colorectal cancer. Previous studies in Vietnam mainly investigated the adenoma pathway, with limited data on the serrated pathway. To evaluate the prevalence, risk factors, and BRAF mutations of SSLs in the Vietnamese population. This is a cross-sectional study conducted on patients with lower gastrointestinal symptoms who underwent colonoscopy at a tertiary hospital in Vietnam. SSLs were diagnosed on histopathology according to the 2019 World Health Organization classification. BRAF mutation analysis was performed using the Sanger DNA sequencing method. The multivariate logistic regression model was used to determine SSL-associated factors. There were 2489 patients, with a mean age of 52.1 ± 13.1 and a female-to-male ratio of 1:1.1. The prevalence of SSLs was 4.2% [95% confidence interval (CI): 3.5-5.1]. In the multivariate analysis, factors significantly associated with SSLs were age ≥ 40 [odds ratio (OR): 3.303; 95%CI: 1.607-6.790], male sex (OR: 2.032; 95%CI: 1.204-3.429), diabetes mellitus (OR: 2.721; 95%CI: 1.551-4.772), and hypertension (OR: 1.650, 95%CI: 1.045-2.605). The rate of BRAF mutations in SSLs was 35.5%. The prevalence of SSLs was 4.2%. BRAF mutations were present in one-third of SSLs. Significant risk factors for SSLs included age ≥ 40, male sex, diabetes mellitus, and hypertension.

Amplicon and Cas9-targeted nanopore sequencing of Varroa destructor at the onset of an outbreak in Australia

Varroa destructor is a major pest of European honey bees (Apis mellifera), causing significant economic and welfare impacts. Australia remained the last continent free from V. destructor. In June 2022, a detection of V. destructor was reported in sentinel colonies at the Port of Newcastle, Australia. Rapid and accurate identification of the species was critical for timely response and management. In this case study, two Nanopore DNA sequencing methods, PCR amplicon sequencing and Cas9-targeted sequencing, were used to rapidly diagnose the species and mitochondrial haplotype of Varroa mites in parasitized colonies. Nanopore PCR amplicon sequencing provided molecular identification of the species and halogroup determination within 24 hours based on a 458 bp amplicon of the mitochondrial Cytochrome c oxidase subunit I (COXI) gene. We also developed and applied a Cas9-targeted Nanopore sequencing technique that used eight guide RNAs to enrich for 5240 bp of the mitochondrial genome. This method delivered richer data for identification within the same timeframe. Our results underscore the efficacy of Nanopore amplicon sequencing and represent the first reported application of Cas9-targeted Nanopore sequencing within a biosecurity framework. These findings enhance the repertoire of diagnostic tools available for biosecurity applications.

Development of molecular identification methods for Dryophytes suweonensis and D. japonicus, and their hybrids.

As hybridization can reduce biodiversity or cause extinction, it is important to identify both purebred parental species and their hybrids prior to conserving them. The Suwon tree frog, Dryophytes suweonensis, is an endangered wildlife species in Korea that shares its habitat and often hybridizes with the Japanese tree frog, D. japonicus. In particular, D. suweonensis, D. japonicus, and their hybrids often have abnormal ovaries and gonads, which are known causes that could threaten their existence. We collected 57 individuals from six localities where D. suweonensis is known to be present. High-resolution melting curve (HRM) analysis of the mitochondrial 12S ribosomal RNA gene was performed to determine the maternal species. Thereafter, the DNA sequences of five nuclear genes (SIAH, TYR, POMC, RAG1, and C-MYC) were analyzed to determine their parental species and hybrid status. The HRM analysis showed that the melting temperature of D. suweonensis was in the range of 79.0-79.3°C, and that of D. japonicus was 77.7-78.0°C, which clearly distinguished the two tree frog species. DNA sequencing of the five nuclear genes revealed 37 single-nucleotide polymorphism (SNP) sites, and STRUCTURE analysis showed a two-group structure as the most likely grouping solution. No heterozygous position in the purebred parental sequences with Q values ≥ 0.995 were found, which clearly distinguished the two treefrog species from their hybrids; 11 individuals were found to be D. suweonensis, eight were found to be D. japonicus, and the remaining 38 individuals were found to be hybrids. Thus, it was possible to unambiguously identify the parental species and their hybrids using HRM analysis and DNA sequencing methods. This study provided fundamental information for D. suweonensis conservation and restoration research.

Fish species mislabelling in China: evidence from a survey of frozen, dried, and surimi-based fish products marketed as cod-related species

‘Cod’-related species are among the most appreciated marine fish resources around the world, but are also prone to species mislabelling. In the present study, a total of 76 frozen, dried, and surimi-based fish products, sold as ‘Cod’ (59 products), ‘Atlantic authentic Cod’ (11 products), and ‘Authentic Cod’ (6 products), were collected in China. A species-specific LAMP (loop-mediated isothermal amplification) method was used to screen for the presence of Atlantic cod (Gadus morhua), Pacific cod (G. macrocephalus), Alaska pollock (G. chalcogrammus), Southern hake (Merluccius australis), which was cross-confirmed using real-time PCR and DNA sequencing methods. The results highlighted the greatest species diversity for ‘Cod’ products, and the identified species were from nine different families. It appears that the practice of assigning a specific type or category of species to the common name ‘Cod’ has not been widely advocated, and the misuse of this ambiguous common name has been a common practice for species adulteration, negatively impacting consumers’ rights and marine conservation. To rebuild consumers’ confidence, retail fish suppliers have differentiated their products by adding specific qualifiers in front of the common name ‘Cod’ on the label, such as ‘Authentic cod’ and ‘Atlantic authentic cod’. The endeavour is highly meaningful, since Gadus morhua was identified as the species for a significant majority of ‘Atlantic authentic cod’ and ‘Authentic cod’ products (64.7%, 11/17), with the remaining six products identified as Alaskan pollock (G. chalcogrammus), Pacific cod (G. macrocephalus) and North Pacific hake (Merluccius productus). Despite the positive effort to reverse species mislabelling from retail on-line fish suppliers, a standardized fish nomenclature stipulated by the responsible authorities remains crucial for enhancing transparency and continuing to reduce species mislabelling.

Genomic predictions for resistance to Aeromonas hydrophila in pacu (Piaractus mesopotamicus)

Infections caused by the bacterium Aeromonas hydrophila result in high mortality outbreaks in pacu (Piaractus mesopotamicus) production. An alternative to increase the resistance to A. hydrophila infection is the implementation of genomic selection programs, which can contribute to the sustainable development of the Brazilian aquaculture. In addition, low-density SNP panels, optimized by genotype imputation, represent a low-cost tool for implementing genomic selection on a commercial scale. This study aimed to evaluate the accuracy of genomic prediction using low (0.5 and 1 K) and medium (15 K imputed) density SNP panels; and compare breeding values accuracy with restriction-site associated DNA sequencing (18 K) and pedigree-based methods for individuals challenged for A. hydrophila. The genomic prediction analyses were conducted using genotypes from the medium-density SNP panel of pacu as reference and genotypes from a low-density SNP panel (0.5 and 1 K) using target genotype-by-sequencing technology (tGBS) for imputation. Approximately 15 K SNPs were successfully imputed with accuracy (r2) above 0.8. Time to death (TD) and binary survival (BS) against A. hydrophila infection were used as phenotypes for genomic prediction. Genomic values were estimated using the GBLUP model adopting different genotyping methods (SNP array and RAD-seq) and densities of SNPs (0.5 K, 1 K, and 15K_imputed for the array and 18 K for RAD-seq). Accuracies of genomic approaches using different methods and densities were compared to the pedigree method PBLUP to assess the benefits of incorporating genomic information. The accuracy of predicting breeding values using genomic methods outperformed the estimates obtained through PBLUP by 33%, 39%, and 67% for TD, and 36%, 39%, and 82% for BS, considering the densities of 0.5, 1, and 15K_imputed, respectively. The imputation strategy yielded better results in terms of genomic prediction accuracy for both traits. The RAD-seq genotyping method produced low-accuracy estimates (∼30% lower for both traits) compared to PBLUP, indicating that such method may not be suitable for genomic prediction of these traits. Genomic selection using genotype imputation is a viable and cost-effective approach to increase resistance to A. hydrophila infection in pacu and reduce genotyping costs without compromising genomic prediction accuracy.

A search for bacteria identified from cerebrospinal fluid shunt infections in previous surgical events.

Shunt infections are a common complication when treating hydrocephalus by cerebrospinal fluid (CSF) shunt placement. The source of infecting pathogens is not well understood. One hypothesis, which we explored here, is that microorganisms persist chronically in the host long before a symptomatic infection occurs and may be detectable in surgical events preceding infection. A cohort of 13 patients was selected, for which CSF samples were available from an infection episode and from a previous surgery event, which was either an initial shunt placement or a revision. Microbiota were analyzed both directly from CSF and from isolates cultured from CSF on aerobic and anaerobic media. The detection and identification of bacteria was done with high throughput DNA sequencing methods and mass spectrometry. The presence of bacteria was confirmed in 4 infection samples, of which 2 were after initial placement and 2 after revision surgery. Taxonomic identification was consistent with clinical microbiology laboratory results. Bacteria were not detected in any of the CSF samples collected at the time of the previous surgical events. While our findings do not provide direct evidence for long-term persistence of pathogens, they suggest the need for consideration of additional source material, such as biofilm and environmental swabs, and/or the use of more sensitive and specific analytical methods.

Plasma TNF-α Elevation in Biologic Naive Rheumatoid Arthritis Patients Belonging to a Population with New Mutations in TLR4 and CYP51A1 genes without Association with Disease-Related Antibodies Levels.

In a system biology-based study, we previously reported that IL-6 and IL6R -specific m-RNA levels were elevated in leukocytes of patients with Rheumatoid arthritis (RA). Here, the association of toll-like receptor4 (TLR4) rs 141534085 and cytochrome P450 family 51 subfamily A member 1(CYP51A1) rs6 with tumor necrosis factor-α (TNF- α), rheumatoid factor (RF)- and Anti- cyclic citrullinated peptide (anti-CCP) antibody -positivity was investigated in almost the same subjects. Forty-six patients and 48 normal subjects were recruited in this study. The blood leucocytes TLR4 rs 141534085 and CYP51A1 rs6 -comprising DNA sequences were amplified by using tetra-primer amplification refractory mutation system polymerase chain reaction (T-ARMS-PCR) technique and the PCR products were checked by Sanger DNA sequencing method. ELISA method was used to determine plasma levels of TNF- α, anti-CCP antibody and RF positivity of plasma was evaluated through a latex agglutination test. The TNF- α level was significantly higher in the patient group than control subjects (p= 0.001). Moreover, we were not able to find any correlation between TNF-α levels and RF as well as anti-CCP antibodies when we used the K2/ Fisher's exact test and Pearson test respectively. Our DNA sequencing data revealed the following new mutations in TLR4 rs141534085 comprising regions: A>T in position 1050, T>A in position 1052, and C>A in position 1054; and for CYP51A1 rs6 encompassing region, the new mutations were; G>A in position 21680, the T nucleotide was inserted in position 21762 and the G nucleotide was inserted in position 21763, G>T in position 21764. The data of this study showed that both TLR4 rs141534085 and CYP51A1 rs6 related DNA regions should be considered as hotspot areas in RA pathogenicity. Moreover, these data indicated that, TNF- α did not alter the production of anti-CCP and RF pathogenic antibodies in patients with long-term RA.

The polymorphism of the ovine insulin like growth factor-2 (IGF2) gene and their associations with growth related traits in Tibetan sheep.

In the current study, the role of the ovine IGF2 as a potential candidate gene was investigated as though marker-assisted selection in Chinese Tibetan sheep. The Sanger DNA sequencing method explored five single nucleotide polymorphisms (SNPs) in 5'UTR of the ovine IGF2 gene (C15640T, G15801A, G15870A, C15982G and G15991A) in Chinese Tibetan sheep. The frequencies of four SNPs were within the Hardy-Weinberg Equilibrium (chi-square test) except C15982G. The statistical analysis indicated that the C15640T and G15801A were significantly associated with body height, body length, chest circumference, and body weight (P < 0.05 or P < 0.01). Furthermore, C15982G variant exhibited significant correlation with the body weight (P < 0.01). These findings suggests that the promoter variants of IGF2 gene could be used as a candidate gene through marker-assisted selection for the body weight and body measurement traits in Tibetan sheep breeding program.

From Dormant Collections to Repositories for the Study of Habitat Changes: The Importance of Herbaria in Modern Life Sciences.

In recent decades, the advent of new technologies for massive and automatized digitization, together with the availability of new methods for DNA sequencing, strongly increased the interest and relevance of herbarium collections for the study of plant biodiversity and evolution. These new approaches prompted new projects aimed at the creation of a large dataset of molecular and phenological data. This review discusses new challenges and opportunities for herbaria in the context of the numerous national projects that are currently ongoing, prompting the study of herbarium specimens for the understanding of biodiversity loss and habitat shifts as a consequence of climate changes and habitat destruction due to human activities. With regard to this, the National Biodiversity Future Center (active in Italy since 2022) started a large-scale digitization project of the Herbarium Centrale Italicum in Florence (Italy), which is the most important Italian botanical collection, consisting of more than 4 million samples at present.

Long-term decitabine/retinoic acid maintenance treatment in an elderly sAML patient with high-risk genetics

Elderly patients with AML ineligible for induction have a dismal prognosis; hence disease stabilization is a primary treatment goal. This case of a 75-year-old patient with secondary AML receiving the combination of decitabine and ATRA (within the DECIDER trial, NCT00867672) demonstrates an above-average survival. The therapy administered over 52 cycles led to complete molecular and hematological remission and resulted in 5.3 years overall survival. Clonal evolution of the leukemic clone could be demonstrated using DNA sequencing methods. According to the literature, this case constitutes the longest continued HMA exposure in an elderly AML patient ineligible for standard chemotherapy.

Genomic Divergence and the Evolution of Ecotypes in Bottlenose Dolphins (Genus Tursiops).

Climatic changes have caused major environmental restructuring throughout the world's oceans. Marine organisms have responded to novel conditions through various biological systems, including genomic adaptation. Growing accessibility of next-generation DNA sequencing methods to study nonmodel species has recently allowed genomic changes underlying environmental adaptations to be investigated. This study used double-digest restriction-site associated DNA (ddRAD) sequence data to investigate the genomic basis of ecotype formation across currently recognized species and subspecies of bottlenose dolphins (genus Tursiops) in the Southern Hemisphere. Subspecies-level genomic divergence was confirmed between the offshore common bottlenose dolphin (T. truncatus truncatus) and the inshore Lahille's bottlenose dolphin (T. t. gephyreus) from the southwestern Atlantic Ocean (SWAO). Similarly, subspecies-level divergence is suggested between inshore (eastern Australia) Indo-Pacific bottlenose dolphin (T. aduncus) and the proposed Burrunan dolphin (T. australis) from southern Australia. Inshore bottlenose dolphin lineages generally had lower genomic diversity than offshore lineages, a pattern particularly evident for T. t. gephyreus, which showed exceptionally low diversity. Genomic regions associated with cardiovascular, musculoskeletal, and energy production systems appear to have undergone repeated adaptive evolution in inshore lineages across the Southern Hemisphere. We hypothesize that comparable selective pressures in the inshore environment drove similar adaptive responses in each lineage, supporting parallel evolution of inshore bottlenose dolphins. With climate change altering marine ecosystems worldwide, it is crucial to gain an understanding of the adaptive capacity of local species and populations. Our study provides insights into key adaptive pathways that may be important for the long-term survival of cetaceans and other organisms in a changing marine environment.

Genome-wide screening for selection signatures in native and cosmopolitan cattle breeds reared in Türkiye.

Via long-term natural and artificial selection pressure, homozygosity may extend across the genome, leaving genomic patterns called selection signatures. This study is the first attempt to assess genome-wide selection signatures in six native Turkish and two cosmopolitan cattle breeds by 211.119 bi-allelic SNPs recovered using the double digest restriction associated DNA sequencing method. The integrated haplotype score (iHS) statistic was utilised to reveal selection signatures within populations, whereas the cross-population extended haplotype homozygosity (XP-EHH) and fixation index (FST ) approaches were preferred to reveal differently fixed genomic regions between native Turkish and cosmopolitan cattle breeds. Selection signatures in 142 genomic regions containing 305 genes were detected within eight cattle breeds by iHS statistics. The XP-EHH and FST approaches revealed that 197 and 114 SNPs were under selection pressure, respectively, which overlapped with 144 and 190 genes, respectively. A total of 18 genes were detected by at least two approaches. Six genes related to disease resistance (TTP2), meat yield (DIAPH3 and METTL21C), meat quality (ZNF24 and ZNF397) and first calving interval (ZSCAN30) turned out to be differently fixed between native Turkish and cosmopolitan cattle breeds, as they were identified by both XP-EHH and FST approaches. In addition, the iHS approach revealed that eight genes associated with visual modality (LSGN), olfaction (MOXD2, OR4C1F and OR4C1F), and immune response (TRBV3-1 and CLDN10) were under selection pressure in both native and cosmopolitan cattle breeds. Owing to their being significantly related to survival traits, these regions may have played a key role in cattle genome evolution. Future studies utilising denser genetic data are required to obtain deeper knowledge on effects of natural and artificial selection in Anatolian cattle breeds. © 2023 Stichting International Foundation for Animal Genetics.

Biomolecular mycorrhizal diversity in Azadirechta excelsa (Jack) Jacobs roots in habitats with different altitude

Azadirechta excelsa is a member of the Meliaceae family with high economic value. This tree is widely planted in South Sumatra, especially in Bengkulu Province. The mutualism symbiosis between several Meliaceae plants and mycorrhiza plants has been studied. However, the diversity of mycorrhiza in the roots of A. excelsa is not well known. We conducted a biomolecular analysis of the roots of this tree species to find out in more detail the types of mycorrhiza associated with the roots and investigated habitat differences in altitude (highland and lowland) that affect mycorrhizal diversity. We used the DNA sequencing method from the roots of A. excelsa and amplified it through PCR using ITS primers forward ITS5.8S and reverse ITS4. We took 14 root samples of E. excelsa; 6 trees from the highland (959 m asl) and eight trees from the lowland (63m asl). The average age of the tree is eight years. The results of the sequencing identification from NCBI using the BLAST method found ten types of mycorrhiza from 4 families associated with A. excelsa roots. They are seven species from the Mycenaceae family (Mycena sp, Filoboletus manipularis, M. amicta, M. pura, M. rosea, M. citrinomarginata, and Favolaschia manipularis), one species from the Psathyrellaceae family (Psathyrella sp), one species from the Tricholomataceae family (Tricholomataceae sp), and one species from the Agaricaceae family (Agaricales sp). All of those types of mycorrhiza belong to the Agaricales order, the Agaricomycetes class from the Basidiomycota phylum. There was no specific type of mycorrhiza that occupied the roots of the tree. Thus, the difference in altitude did not affect the type of mycorrhiza on A. excelsa roots.

EP20.43: The rapid long read DNA sequencing method identified bacterial species causing silent intra‐amniotic infection

EP20.43: The rapid long read DNA sequencing method identified bacterial species causing silent intra‐amniotic infection

Blast from the past: a study of decades-old fungal cultures resolves a long-standing tree disease mystery

A root disease in plantations of Pinus radiata and Pinus pinaster, where trees died in distinct patches, was present in the Western Cape province of South Africa during the 1970s and 1980s. Phytophthora cinnamomi was initially believed to be the cause, but the disease was later ascribed to the insect-associated fungus Leptographium serpens, a fungal species residing in the Ophiostomatales. Doubt regarding the cause of the disease was raised in a later study due to the fact that most Leptographium spp., particularly those that colonise ray parenchyma tissues, which is the case for L. serpens, are not typically primary disease agents. In this study, cultures of an unidentified sterile fungus collected from the dying trees were revived and identified using DNA sequencing methods, which were not available when the disease was first studied. These cultures were identified as the pyrophillic pathogen Rhizina undulata, well-known to cause patch death of conifers in South Africa and elsewhere in the world. While the patches of dying trees no longer exist and the disease cannot be newly studied, it is most likely that the tree death originally thought to be caused by L. serpens was due primarily to R. undulata. The study provides a vivid example of the value of preserving cultures of fungi for later study and the power of modern techniques to identify fungal pathogens.

Fluorometric Detection of SARS-CoV-2 Single-Nucleotide Variant L452R Using Ligation-Based Isothermal Gene Amplification.

Since the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant was first discovered, several variants showing different infectivity and immune responses have emerged globally. As the conventional method, whole-genome sequencing following polymerase chain reaction (PCR) is currently used for diagnosis of SARS-CoV-2 mutations. However, these conventional PCR-based direct DNA sequencing methods are time-consuming, complicated, and require expensive DNA sequencing modules. Here, we developed a fluorometric method for the accurate detection of a single missense mutation of U to G in the spike (S) gene that changes leucine to arginine (L452R) in SARS-CoV-2 genomic RNA. Our method for the detection of single-nucleotide mutations (SNM) in the viral RNA genome includes RNA sequence-dependent DNA ligation and tandem isothermal gene amplification methods, such as strand displacement amplification (SDA) and rolling circle amplification (RCA) generating G-quadruplex (GQ). In the presence of SNM in the viral RNA, ligation of both ends of the probe DNAs occurs between 5'-phosphorylated hairpin DNA and linear probe DNA that can discriminate a single base mismatch. The ligated DNAs were then extended to generate long-stem hairpin DNAs that are subjected to the first isothermal gene amplification (SDA). SDA produces multitudes of short ssDNA from the long-stem hairpin DNAs, which then serve as primers by annealing to circular padlock DNA for the second isothermal gene amplification (RCA). RCA produces a long stretch of ssDNA containing GQ structures. Thioflavin T (ThT) is then intercalated into GQ and emits green fluorescence, which allows the fluorometric identification of SARS-CoV-2 variants. This fluorometric analysis sensitively distinguished SNM in the L452R variant of SARS-CoV-2 RNA as low as 10 pM within 2 h. Hence, this fluorometric detection method using ligation-assisted tandem isothermal gene amplification can be applied for the diagnosis of SARS-CoV-2 SNM variants with high accuracy and sensitivity, without the need for cumbersome whole-genome DNA sequencing.

Potential value of high-throughput single-cell DNA sequencing of Juvenile myelomonocytic leukemia: report of two cases

Juvenile myelomonocytic leukemia (JMML) is a rare myeloproliferative disease of early childhood that develops due to mutations in the genes of the RAS-signaling pathway. Next-generation high throughput sequencing (NGS) enables identification of various secondary molecular genetic events that can facilitate JMML progression and transformation into secondary acute myeloid leukemia (sAML). The methods of single-cell DNA sequencing (scDNA-seq) enable overcoming limitations of bulk NGS and exploring genetic heterogeneity at the level of individual cells, which can help in a better understanding of the mechanisms leading to JMML progression and provide an opportunity to evaluate the response of leukemia to therapy. In the present work, we applied a two-step droplet microfluidics approach to detect DNA alterations among thousands of single cells and to analyze clonal dynamics in two JMML patients with sAML transformation before and after hematopoietic stem cell transplantation (HSCT). At the time of diagnosis both of our patients harbored only “canonical” mutations in the RAS signaling pathway genes detected by targeted DNA sequencing. Analysis of samples from the time of transformation JMML to sAML revealed additional genetic events that are potential drivers for disease progression in both patients. ScDNA-seq was able to measure of chimerism level and detect a residual tumor clone in the second patient after HSCT (sensitivity of less than 0.1% tumor cells). The data obtained demonstrate the value of scDNA-seq to assess the clonal evolution of JMML to sAML, response to therapy and engraftment monitoring.

Exploration of pathogenic microorganism within the small intestine of necrotizing enterocolitis.

Necrotizing enterocolitis (NEC) is the most common severe gastrointestinal emergency in neonates. We designed this study to identify the pathogenic microorganisms of NEC in the microbiota of the small intestine of neonates. Using the 16S ribosomal DNA (rDNA) sequencing method, we compared and analyzed the structure and diversity of microbiotas in the intestinal feces of different groups of neonates: patients undergoing jejunostomy to treat NEC (NP group), neonates undergoing jejunostomy to treat other conditions (NN group), and neonates with NEC undergoing conservative treatment (NC group). We took intestinal feces and saliva samples from patients at different time points. The beta diversities of the NP, NN, and NC groups were all similar. When comparing the beta diversities between different time points in the NP group, we found similar beta diversities at time points E1 to E3 but significant differences between the E2-E3 and E4 time points: the abundances of Klebsiella and Enterococcus (Proteobacteria) were higher at the E1-E3 time points; the abundance of Escherichia-Shigella (Proteobacteria) increased at the E2 time point, and the abundance of Klebsiella decreased significantly, whereas that of Streptococcus increased significantly at the E4 time point. Our results suggest that the pathological changes of intestinal necrosis in the small intestine of infants with NEC are not directly caused by excessive proliferation of pathogenic bacteria in the small intestine. The sources of microbiota in the small intestine of neonates, especially in premature infants, may be affected by multiple factors.

Scaling-up RADseq methods for large datasets of non-invasive samples: Lessons for library construction and data preprocessing.

Genetic non-invasive sampling (gNIS) is a critical tool for population genetics studies, supporting conservation efforts while imposing minimal impacts on wildlife. However, gNIS often presents variable levels of DNA degradation and non-endogenous contamination, which can incur considerable processing costs. Furthermore, the use of restriction-site-associated DNA sequencing methods (RADseq) for assessing thousands of genetic markers introduces the challenge of obtaining large sets of shared loci with similar coverage across multiple individuals. Here, we present an approach to handling large-scale gNIS-based datasets using data from the spotted hyena population inhabiting the Ngorongoro Crater in Tanzania. We generated 3RADseq data for more than a thousand individuals, mostly from faecal mucus samples collected non-invasively and varying in DNA degradation and contamination level. Using small-scale sequencing, we screened samples for endogenous DNA content, removed highly contaminated samples, confirmed overlap fragment length between libraries, and balanced individual representation in a sequencing pool. We evaluated the impact of (1) DNA degradation and contamination of non-invasive samples, (2) PCR duplicates and (3) different SNP filters on genotype accuracy based on Mendelian error estimated for parent-offspring trio datasets. Our results showed that when balanced for sequencing depth, contaminated samples presented similar genotype error rates to those of non-contaminated samples. We also showed that PCR duplicates and different SNP filters impact genotype accuracy. In summary, we showed the potential of using gNIS for large-scale genetic monitoring based on SNPs and demonstrated how to improve control over library preparation by using a weighted re-pooling strategy that considers the endogenous DNA content.

Identification of a novel variant c.761C>T on ABO*B.01 gene in ABO glycosyltransferases associated with Bweak phenotype.

ABO antigens are produced from H antigen by the activity of glycosyltransferase enzyme encoded by the ABO gene. Variants in the ABO gene can produce a weak ABO phenotype. In this study, we identify a novel ABO*BW allele and investigate the underlying mechanism leading to the Bweak phenotype. The ABO phenotype and genotype of the sample were determined using serological and direct DNA sequencing methods. We assessed the impact of the novel variant by three-dimensional modelling to predict protein stability changes (ΔΔG), and carried out an in vitro expression assay. The total glycosyltransferase transfer capacity in the supernatant of transfected cells was also examined. Serological analysis confirmed the Bweak phenotype in the subject, and gene sequencing identified a novel variant c.761C>T (p.A254V) on the ABO*B.01 allele, resulting in a BW-var/O.01.02 genotype. In silico analysis suggested that the p.A254V variant on the B allele may reduce the stability of glycosyltransferase B (GTB), as indicated by the ΔΔG values. In vitro expression studies showed that the variant p.A254V impaired H to B antigen conversion, although it did not affect the expression of GTB. We identified a novel BW allele and demonstrated that the variant c.761C>T (p.A254V) can cause the Bweak phenotype by reducing the stability of GTB.