DNA Sequence Research Articles

Sampling technique to pool genetic materials of microorganism communities in blue-swimming crab processing plant industry

The crab and seafood processing industry must fulfill standard requirements for sanitation, hygiene, and good manufacturing methods to ensure the safety of the products and free from foodborne bacteria. However, equipment and processing unit surfaces are challenging to clean optimally, which can cause persistent bacteria to emerge. Eliminating persistent bacteria is the latest challenge in the fish processing industry for optimal disinfection, preventing cross-contamination, and controlling foodborne outbreaks. Microbiological testing in industry has been limited to selective culture-media techniques; thus, a rapid, sensitive, accurate, and routine applicable analytical method is urgently needed. The significant reduction in the costs of high-throughput sequencing technologies supports the possibility of routine applications in the industry. This study aimed to determine the profile of the microbial community on the surface of the production room and blue-swimming crab processing unit equipment using short-read metagenomic techniques. The analysis included the stages of sampling, bacterial incubation, bacterial DNA isolation, sequencing, and bioinformatics analysis. The first important step to increase the possibility of routine adoption in the seafood industry is to reduce the cost, complexity, and time required to complete the analysis. Therefore, in this protocol, we generate a scalable, flexible, cost-effective, and auditable workflow.• Collection of bacterial samples by swabbing the surface of the equipment using a sterile cotton swab and sterile cloth, which is easy to apply and follow in the blue-swimming crab processing plant industry.• Effective and efficient sample-pooling is an important step in identifying bacterial communities by metagenomic analysis.

Exotic myxozoan parasites establish complex life cycles in farm pond aquaculture

Myxozoans are obligate parasites with complex life cycles, typically infecting fish and annelids. Here, we examined annelids from fish farm pond sediments in the Beit Shean Valley, in the Syrian-African Rift Valley, Israel, for myxozoan infections. We examined 1486 oligochaetes, and found 74 (5 %) were infected with actinospore stages. We used mitochondrial 16S sequencing to infer identity of 25 infected annelids as species of Potamothrix, Psammoryctides, Tubifex and Dero. We identified 7 myxozoan types from collective groups Neoactinomyxum and Sphaeractinomyxon, and characterized them by small subunit ribosomal DNA sequencing. The Neoactinomyxum type was genetically most similar (∼93 %) to cyprinid fish-infecting Myxobolus spp. The six Sphaeractinomyxon types were genetically similar (93–100 %) to Mugilid-infecting Myxobolus spp.; with one being the previously unknown actinospore stage of a myxospore that infects mullet from aquaculture from the Israeli coast of the Mediterranean Sea. As the farm pond system is artificial and geographically isolated from the Mediterranean, the presence of at least seven myxozoans in their annelid hosts demonstrates introduction and establishment of these parasites in a novel, brackish environment.

The development of non-natural type nucleoside to stabilize triplex DNA formation against CG and TA inversion site.

Based on the sequence-specific recognition of target duplex DNA by triplex-forming oligonucleotides (TFOs) at the major groove side, the antigene strategy has been exploited as a gene-targeting tool with considerable attention. Triplex DNA is formed via the specific base triplets by the Hoogsteen or reverse Hoogsteen hydrogen bond interaction between TFOs and the homo-purine strand from the target duplex DNA, leading to the established sequence-specificity. However, the presence of inversion sites, which are known as non-natural nucleosides that can form satisfactory interactions with 2-deoxythymidine (dT) and 2-deoxycytidine (dC) in TA and CG base pairs in the target homo-purine DNA sequences, drastically restricts the formation of classically stable base triplets and even the triplex DNA. Therefore, the design of non-natural type nucleosides, which can effectively recognize CG or/and TA inversion sites with satisfactory selectivity, should be of great significance to expanding the triplex-forming sequence. Here, this review mainly provides a comprehensive review of the current development of novel non-natural nucleosides to recognize CG or/and TA inversion sites in triplex DNA formation against double-strand DNA (dsDNA).

Efficacy of Cefiderocol, a Novel Siderophore Cephalosporin, against Multidrug Resistant Acinetobacter baumannii Clinical Isolates in Japan.

Multidrug-resistant Acinetobacter baumannii (MDRAB) is an important pathogen that causes nosocomial infections and is resistant to almost all antibiotics, including carbapenems. Cefiderocol is a novel siderophore cephalosporin active against a broad spectrum of gram-negative bacteria. However, the susceptibility of MDRAB to cefiderocol has not yet been reported in Japan. In this study, we measured the minimum inhibitory concentrations (MICs) of antibiotics including cefiderocol against MDRAB clinical isolates collected during a nosocomial outbreak between 2009 and 2010 at the Teikyo University Hospital in Japan. We found that all 10 MDRAB clinical isolates tested were susceptible to cefiderocol, with an MIC range of 0.12 to 1 μg/mL. All the isolates also exhibited resistance to ampicillin-sulbactam and an intermediate resistance to colistin, whereas nine of them were susceptible to tigecycline. DNA sequencing revealed that all strains harbored an OXA-51-like carbapenemase, a major cause of carbapenem resistance in A. baumannii in Japan. In conclusion, this study showed that the cefiderocol susceptibility of MDRAB clinical isolates in Japan was equivalent to that to colistin or tigecycline, and thus cefiderocol is a potential treatment option for MDRAB infections.

Identifying widespread and recurrent variants of genetic parts to improve annotation of engineered DNA sequences.

Engineered plasmids have been workhorses of recombinant DNA technology for nearly half a century. Plasmids are used to clone DNA sequences encoding new genetic parts and to reprogram cells by combining these parts in new ways. Historically, many genetic parts on plasmids were copied and reused without routinely checking their DNA sequences. With the widespread use of high-throughput DNA sequencing technologies, we now know that plasmids often contain variants of common genetic parts that differ slightly from their canonical sequences. Because the exact provenance of a genetic part on a particular plasmid is usually unknown, it is difficult to determine whether these differences arose due to mutations during plasmid construction and propagation or due to intentional editing by researchers. In either case, it is important to understand how the sequence changes alter the properties of the genetic part. We analyzed the sequences of over 50,000 engineered plasmids using depositor metadata and a metric inspired by the natural language processing field. We detected 217 uncatalogued genetic part variants that were especially widespread or were likely the result of convergent evolution or engineering. Several of these uncatalogued variants are known mutants of plasmid origins of replication or antibiotic resistance genes that are missing from current annotation databases. However, most are uncharacterized, and 3/5 of the plasmids we analyzed contained at least one of the uncatalogued variants. Our results include a list of genetic parts to prioritize for refining engineered plasmid annotation pipelines, highlight widespread variants of parts that warrant further investigation to see whether they have altered characteristics, and suggest cases where unintentional evolution of plasmid parts may be affecting the reliability and reproducibility of science.

Oesophagostomum stephanostomum causing parasitic granulomas in wild chimpanzees (Pan troglodytes verus) of Taï National Park, Côte d'Ivoire.

Nematodes belonging to the genus Oesophagostomum frequently infect wild chimpanzees (Pan troglodytes) across widely separated field sites. Nodular lesions (granulomas) containing Oesophagostomum are commonly seen in the abdomen of infected chimpanzees post-mortem. At Taï National Park, Côte d'Ivoire, previous studies have identified larvae of a variety of Oesophagostomum spp. in wild chimpanzee stool, based on sequencing of larval DNA, and nodular lesions associated with Oesophagostomum, identified morphologically to the genus level but not sequenced. Here we present three recent cases of parasitic granulomas found post-mortem in chimpanzees at Taï. We complement descriptions of gross pathology, histopathology and parasitology with PCR and sequencing of DNA isolated from the parasitic nodules and from adult worms found inside the nodules. In all three cases, we identify Oesophagostomum stephanostomum as the causative agent. The sequences from this study were identical to the only other published sequences from nodules in nonhuman primates-those from the wild chimpanzees of Gombe, Tanzania.

Epitope binning for multiple antibodies simultaneously using mammalian cell display and DNA sequencing

Epitope binning, an approach for grouping antibodies based on epitope similarities, is a critical step in antibody drug discovery. However, conventional methods are complex, involving individual antibody production. Here, we established Epitope Binning-seq, an epitope binning platform for simultaneously analyzing multiple antibodies. In this system, epitope similarity between the query antibodies (qAbs) displayed on antigen-expressing cells and a fluorescently labeled reference antibody (rAb) targeting a desired epitope is analyzed by flow cytometry. The qAbs with epitope similar to the rAb can be identified by next-generation sequencing analysis of fluorescence-negative cells. Sensitivity and reliability of this system are confirmed using rAbs, pertuzumab and trastuzumab, which target human epidermal growth factor receptor 2. Epitope Binning-seq enables simultaneous epitope evaluation of 14 qAbs at various abundances in libraries, grouping them into respective epitope bins. This versatile platform is applicable to diverse antibodies and antigens, potentially expediting the identification of clinically useful antibodies.

Enhancing the epidemiological surveillance of SARS-CoV-2 using Sanger sequencing to identify circulating variants and recombinants.

Since the emergence of SARS-CoV-2 in December 2019, more than 12,000 mutations in the virus have been identified. These could cause changes in viral characteristics and directly impact global public health. The emergence of variants is a great concern due to the chance of increased transmissibility and infectivity. Sequencing for surveillance and monitoring circulating strains is extremely necessary as the early identification of new variants allows public health agencies to make faster and more effective decisions to contain the spread of the virus. In the present study, we identified circulating variants in samples collected in Belo Horizonte, Brazil, and detected a recombinant lineage using the Sanger method. The identification of lineages was done through gene amplification of SARS-CoV-2 by Reverse Transcription-Polymerase Chain Reaction (RT-PCR). By using these specific fragments, we were able to differentiate one variant of interest and five circulating variants of concern. We were also able to detect recombinants. Randomly selected samples were sequenced by either Sanger or Next Generation Sequencing (NGS). Our findings validate the effectiveness of Sanger sequencing as a powerful tool for monitoring variants. It is easy to perform and allows the analysis of a larger number of samples in countries that cannot afford NGS.

Formamide denaturation of double-stranded DNA for fluorescence in situ hybridization (FISH) distorts nanoscale chromatin structure.

As imaging techniques rapidly evolve to probe nanoscale genome organization at higher resolution, it is critical to consider how the reagents and procedures involved in sample preparation affect chromatin at the relevant length scales. Here, we investigate the effects of fluorescent labeling of DNA sequences within chromatin using the gold standard technique of three-dimensional fluorescence in situ hybridization (3D FISH). The chemical reagents involved in the 3D FISH protocol, specifically formamide, cause significant alterations to the sub-200 nm (sub-Mbp) chromatin structure. Alternatively, two labeling methods that do not rely on formamide denaturation, resolution after single-strand exonuclease resection (RASER)-FISH and clustered regularly interspaced short palindromic repeats (CRISPR)-Sirius, had minimal impact on the three-dimensional organization of chromatin. We present a polymer physics-based analysis of these protocols with guidelines for their interpretation when assessing chromatin structure using currently available techniques.

Development of label-free electrochemical OMP-DNA probe biosensor as a highly sensitive system to detect of citrus huanglongbing

The fabrication of the first label-free electrochemical DNA probe biosensor for highly sensitive detection of Candidatus Liberibacter asiaticus (CLas), as the causal agent of citrus huanglongbing disease, is conducted here. An OMP probe was designed based on the hybridization with its target-specific sequence in the outer membrane protein (OMP) gene of CLas. The characterization of the steps of biosensor fabrication and hybridization process between the immobilized OMP-DNA probe and the target ssDNA oligonucleotides (OMP-complementary and three mismatches OMP or OMP-mutation) was monitored using cyclic voltammetry and electrochemical impedance spectroscopy based on increasing or decreasing in the electron transfer in [Fe (CN)6]3−/4− on the modified gold electrode surface. The biosensor sensitivity indicated that the peak currents were linear over ranges from 20 to 100 nM for OMP-complementary with the detection limit of 0.026 nM (S/N = 3). The absence of any cross-interference with other biological DNA sequences confirmed a high selectivity of fabricated biosensor. Likewise, it showed good specificity in discriminating the mutation oligonucleotides from complementary target DNAs. The functional performance of optimized biosensor was achieved via the hybridization of OMP-DNA probe with extracted DNA from citrus plant infected with CLas. Therefore, fabricated biosensor indicates promise for sensitivity and early detection of citrus huanglongbing disease.

Interpretable deep learning reveals the role of an E-box motif in suppressing somatic hypermutation of AGCT motifs within human immunoglobulin variable regions

IntroductionSomatic hypermutation (SHM) of immunoglobulin variable (V) regions by activation induced deaminase (AID) is essential for robust, long-term humoral immunity against pathogen and vaccine antigens. AID mutates cytosines preferentially within WRCH motifs (where W=A or T, R=A or G and H=A, C or T). However, it has been consistently observed that the mutability of WRCH motifs varies substantially, with large variations in mutation frequency even between multiple occurrences of the same motif within a single V region. This has led to the notion that the immediate sequence context of WRCH motifs contributes to mutability. Recent studies have highlighted the potential role of local DNA sequence features in promoting mutagenesis of AGCT, a commonly mutated WRCH motif. Intriguingly, AGCT motifs closer to 5’ ends of V regions, within the framework 1 (FW1) sub-region1, mutate less frequently, suggesting an SHM-suppressing sequence context.MethodsHere, we systematically examined the basis of AGCT positional biases in human SHM datasets with DeepSHM, a machine-learning model designed to predict SHM patterns. This was combined with integrated gradients, an interpretability method, to interrogate the basis of DeepSHM predictions.ResultsDeepSHM predicted the observed positional differences in mutation frequencies at AGCT motifs with high accuracy. For the conserved, lowly mutating AGCT motifs in FW1, integrated gradients predicted a large negative contribution of 5’C and 3’G flanking residues, suggesting that a CAGCTG context in this location was suppressive for SHM. CAGCTG is the recognition motif for E-box transcription factors, including E2A, which has been implicated in SHM. Indeed, we found a strong, inverse relationship between E-box motif fidelity and mutation frequency. Moreover, E2A was found to associate with the V region locale in two human B cell lines. Finally, analysis of human SHM datasets revealed that naturally occurring mutations in the 3’G flanking residues, which effectively ablate the E-box motif, were associated with a significantly increased rate of AGCT mutation.DiscussionOur results suggest an antagonistic relationship between mutation frequency and the binding of E-box factors like E2A at specific AGCT motif contexts and, therefore, highlight a new, suppressive mechanism regulating local SHM patterns in human V regions.

Migratory Tumor Cells Cooperate with Cancer Associated Fibroblasts in Hormone Receptor-Positive and HER2-Negative Breast Cancer

Hormone receptor-positive and HER2-negative breast cancer (HR+/HER2-BC) is the most common type with a favorable prognosis under endocrine therapy. However, it still demonstrates unpredictable progression and recurrences influenced by high tumoral diversity and microenvironmental status. To address these heterogeneous molecular characteristics of HR+/HER2-BC, we aimed to simultaneously characterize its transcriptomic landscape and genetic architecture at the same resolution. Using advanced single-cell RNA and DNA sequencing techniques together, we defined four distinct tumor subtypes. Notably, the migratory tumor subtype was closely linked to genomic alterations of EGFR, related to the tumor-promoting behavior of IL6-positive inflammatory tumor-associated fibroblast, and contributing to poor prognosis. Our study comprehensively utilizes integrated analysis to uncover the complex dynamics of this breast cancer subtype, highlighting the pivotal role of the migratory tumor subtype in influencing surrounding cells. This sheds light on potential therapeutic targets by offering enhanced insights for HR+/HER2-BC treatment.

Elucidating epigenetic mechanisms governing odontogenic differentiation in dental pulp stem cells: an in-depth exploration

Epigenetics refers to the mechanisms such as DNA methylation and histone modification that influence gene expression without altering the DNA sequence. These epigenetic modifications can regulate gene transcription, splicing, and stability, thereby impacting cell differentiation, development, and disease occurrence. The formation of dentin is intrinsically linked to the odontogenic differentiation of dental pulp stem cells (DPSCs), which are recognized as the optimal cell source for dentin-pulp regeneration due to their varied odontogenic potential, strong proliferative and angiogenic characteristics, and ready accessibility Numerous studies have demonstrated the critical role of epigenetic regulation in DPSCs differentiation into specific cell types. This review thus provides a comprehensive review of the mechanisms by which epigenetic regulation controls the odontogenesis fate of DPSCs.

Uncovering further diversity of Ochoterenatrema Caballero, 1943 (Digenea: Lecithodendriidae) in South American bats

Ochoterenatrema Caballero, 1943 is a genus of lecithodendriid digeneans that prior to this study included 8 species parasitic in bats in the Western Hemisphere. Species of Ochoterenatrema possess a unique morphological feature in form of the pseudogonotyl on the sinistral side of the ventral sucker. In this study, we describe 2 new species of Ochoterenatrema from bats in Ecuador. The new species are readily differentiated from their congeners by a combination of morphological characters, including the distribution of vitelline follicles, length of oesophagus, sucker ratio and the body shape, among other features. We have generated partial nuclear 28S rDNA and mitochondrial cox1 gene DNA sequences from both new species. The newly obtained sequences were used to differentiate among species and study the phylogenetic interrelationships among Ochoterenatrema spp. The internal topology of the clade was weakly supported, although the cox1 tree was much better resolved than the 28S tree. Comparison of sequences revealed 0–1.2% interspecific divergence in 28S and 3.3–20.5% interspecific divergence in cox1 among Ochoterenatrema spp. The new findings demonstrate that bats in South America likely harbor multiple additional undescribed species of Ochoterenatrema. More extensive sampling from broader geographic and host ranges, especially in North America, should allow for a better understanding of the evolution of host associations and morphological traits of this lineage of lecithodendriid digeneans.

An efficient exact algorithm for planted motif search on large DNA sequence datasets.

DNA motif is the pattern shared by similar fragments in DNA sequences, which plays a key role in regulating gene expression, and DNA motif discovery has become a key research topic. Exact planted (l,d)-motif search (PMS) is one of the motif discovery approaches, which aims to find from t sequences all the (l,d)-motifs that are motifs of l length appearing in at least qt sequences with at most d mismatches. The existing exact PMS algorithms are only suitable for small datasets of DNA sequences. The development of high-throughput sequencing technology generates vast amount of DNA sequence data, which brings challenges to solving exact PMS problems efficiently. Therefore, we propose an efficient exact PMS algorithm called PMmotif for large datasetsof DNA sequences, after analyzing the time complexity of the existing exact PMS algorithms. PMmotif finds (l,d) -motifs with strategy by searching the branches on the pattern tree that may contain (l,d) -motifs. It is verified by experiments that the running time ratio of the existing excellentPMS algorithmstoPMmotif isbetween14.83and 58.94. In addition, for the first time, PMmotif can solve the (15,5) and(17,6) challenge problem instances on large DNA sequence datasets within 24 hours.

Size-dependent enhancement of gene expression by Plasmodium 5'UTR introns.

Eukaryotic genes contain introns that are removed by the spliceosomal machinery during mRNA maturation. Introns impose a huge energetic burden on a cell; therefore, they must play an essential role in maintaining genome stability and/or regulating gene expression. Many genes (> 50%) in Plasmodium parasites contain predicted introns, including introns in 5' and 3' untranslated regions (UTR). However, the roles of UTR introns in the gene expression of malaria parasites remain unknown. In this study, an episomal dual-luciferase assay was developed to evaluate gene expression driven by promoters with or without a 5'UTR intron from four Plasmodium yoelii genes. To investigate the effect of the 5'UTR intron on endogenous gene expression, the pytctp gene was tagged with 3xHA at the N-terminal of the coding region, and parasites with or without the 5'UTR intron were generated using the CRISPR/Cas9 system. We showed that promoters with 5'UTR introns had higher activities in driving gene expression than those without 5'UTR introns. The results were confirmed in recombinant parasites expressing an HA-tagged gene (pytctp) driven by promoter with or without 5'UTR intron. The enhancement of gene expression was intron size dependent, but not the DNA sequence, e.g. the longer the intron, the higher levels of expression. Similar results were observed when a promoter from one strain of P. yoelii was introduced into different parasite strains. Finally, the 5'UTR introns were alternatively spliced in different parasite development stages, suggesting an active mechanism employed by the parasites to regulate gene expression in various developmental stages. Plasmodium 5'UTR introns enhance gene expression in a size-dependent manner; the presence of alternatively spliced mRNAs in different parasite developmental stages suggests that alternative slicing of 5'UTR introns is one of the key mechanisms in regulating parasite gene expression and differentiation.

Genomic signatures of past megafrugivore‐mediated dispersal in Malagasy palms

Abstract Seed dispersal affects gene flow and hence genetic differentiation of plant populations. During the Late Quaternary, most fruit‐eating and seed‐dispersing megafauna went extinct, but whether these animals have left signatures in the population genetics of their food plants, particularly those with large, ‘megafaunal’ fruits (i.e. >4 cm—megafruits), remains unclear. Here, we assessed the population history, genetic differentiation and recent migration among populations of four animal‐dispersed palm (Arecaceae) species with large (Borassus madagascariensis), medium‐sized (Hyphaene coriacea, Bismarckia nobilis) and small (Chrysalidocarpus madagascariensis) fruits on Madagascar. We integrated double‐digest restriction‐site‐associated DNA sequencing (ddRAD) of 167 individuals from 25 populations with (past) distribution ranges for extinct (e.g., giant lemurs and elephant birds) and extant seed‐dispersing animals, landscape and human impact data, and applied linear mixed‐effects models to explore the drivers of genetic variation in Malagasy palms. Palm populations that shared more megafrugivore species in the past had lower genetic differentiation than populations that shared fewer megafrugivore species. This suggests that megafrugivore‐mediated seed dispersal in the past may have led to frequent gene flow among populations. In comparison, extant frugivore diversity only decreased genetic differentiation in the small‐fruited palm. Furthermore, genetic differentiation of all palm species decreased with landscape connectivity (i.e. environmental suitability, forest cover and river density) and human impact (i.e. road density), while migration rates of the small‐fruit palm increased with road density. Synthesis. Our results suggest that the legacy of megafrugivores regularly achieving long dispersal distances is still reflected in the population genetics of palms that were formerly dispersed by such animals. Furthermore, low genetic differentiation was possibly maintained after the megafauna extinctions through alternative dispersal (e.g. human‐ or river‐mediated), long generation times and long lifespans of these megafruit palms. Our study illustrates how species interactions that happened >1000 years ago can leave imprints in their population genetics.

Genomic and T cell repertoire biomarkers associated with malignant mesothelioma survival.

Malignant mesothelioma (MM) is an exceedingly rare tumor with poor prognosis due to the limited availability of effective treatment. Immunotherapy has emerged as a novel treatment approach for MM, but less than 40% of the patients benefit from it. Thus, it is necessary to identify accurate and effective biomarkers that can predict the overall survival (OS) and immunotherapy efficacy for MM. DNA sequencing was used to identify the genomic landscape based on the data from 86 Chinese patients. T cell receptor (TCR) sequencing was used to characterize MM TCR repertoires of 28 patients between October 2016 and April 2023. Patients with TP53, NF2, or CDKN2A variants at the genomic level, as well as those exhibiting lower Shannon index (<6.637), lower evenness (<0.028), or higher clonality (≥0.194) according to baseline tumor tissue TCR indexes, demonstrated poorer OS. Furthermore, patients with TP53, CDKN2A, or CDKN2B variants and those with a lower evenness (<0.030) in baseline tumor tissue showed worse immunotherapy efficacy. The present study is the first to identify five special TCR Vβ-Jβ rearrangements associated with MM immunotherapy efficacy. The present study reported the largest-scale genomic landscape and TCR repertoire of MM in Chinese patients and identified genomic and TCR biomarkers for the prognosis and immunotherapy efficacy in MM. The study results might provide new insights for prospective MM trials using specific genes, TCR indexes, and TCR clones as biomarkers and offer a reference for future antitumor drugs based on TCR-specific clones.

SEDA 2024 update: enhancing the SEquence DAtaset builder for seamless integration into automated data analysis pipelines.

The initial version of SEDA assists life science researchers without programming skills with the preparation of DNA and protein sequence FASTA files for multiple bioinformatics applications. However, the initial version of SEDA lacks a command-line interface for more advanced users and does not allow the creation of automated analysis pipelines. The present paper discusses the updates of the new SEDA release, including the addition of a complete command-line interface, new functionalities like gene annotation, a framework for automated pipelines, and improved integration in Linux environments. SEDA is an open-source Java application and can be installed using the different distributions available ( https://www.sing-group.org/seda/download.html ) as well as through a Docker image ( https://hub.docker.com/r/pegi3s/seda ). It is released under a GPL-3.0 license, and its source code is publicly accessible on GitHub ( https://github.com/sing-group/seda ). The software version at the time of submission is archived at Zenodo (version v1.6.0, http://doi.org/10.5281/zenodo.10201605 ).

Exploring mitochondrial heteroplasmy in neonates: implications for growth patterns and overweight in the first years of life.

Mitochondrial heteroplasmy reflects genetic diversity within individuals due to the presence of varying mitochondrial DNA (mtDNA) sequences, possibly affecting mitochondrial function and energy production in cells. Rapid growth during early childhood is a critical development with long-term implications for health and well-being. In this study, we investigated if cord blood mtDNA heteroplasmy is associated with rapid growth at 6 and 12 months and overweight in childhood at 4-6 years. This study included 200 mother-child pairs of the ENVIRONAGE birth cohort. Whole mitochondrial genome sequencing was performed to determine mtDNA heteroplasmy levels (in variant allele frequency; VAF) in cord blood. Rapid growth was defined for each child as the difference between WHO-SD scores of predicted weight at either 6 or 12 months and birth weight. Logistic regression models were used to determine the association of mitochondrial heteroplasmy with rapid growth and childhood overweight. Determinants of relevant cord blood mitochondrial heteroplasmies were identified using multiple linear regression models. One % increase in VAF of cord blood MT-D-Loop16362T > C heteroplasmy was associated with rapid growth at 6 months (OR = 1.03; 95% CI: 1.01-1.05; p = 0.001) and 12 months (OR = 1.02; 95% CI: 1.00-1.03; p = 0.02). Furthermore, this variant was associated with childhood overweight at 4-6 years (OR = 1.01; 95% CI 1.00-1.02; p = 0.05). Additionally, rapid growth at 6 months (OR = 3.00; 95% CI: 1.49-6.14; p = 0.002) and 12 months (OR = 4.05; 95% CI: 2.06-8.49; p < 0.001) was also associated with childhood overweight at 4-6 years. Furthermore, we identified maternal age, pre-pregnancy BMI, maternal education, parity, and gestational age as determinants of cord blood MT-D-Loop16362T > C heteroplasmy. Our findings, based on mitochondrial DNA genotyping, offer insights into the molecular machinery leading to rapid growth in early life, potentially explaining a working mechanism of the development toward childhood overweight.