DNA Sequence Variation Research Articles

Green algae (Viridiplantae) in sediments from three lakes on Vega Island, Antarctica, assessed using DNA metabarcoding.

Vega Island is located off the eastern tip of the Antarctic Peninsula (Maritime Antarctica), in the Weddell Sea. In this study, we used metabarcoding to investigate green algal DNA sequence diversity present in sediments from three lakes on Vega Island (Esmeralda, Copépodo, and Pan Negro Lakes). Total DNA was extracted and the internal transcribed spacer 2 region of the nuclear ribosomal DNA was used as a DNA barcode for molecular identification. Green algae were represented by sequences representing 78 taxa belonging to Phylum Chlorophyta, of which 32% have not previously been recorded from Antarctica. Sediment from Pan Negro Lake generated the highest number of DNA reads (11,205), followed by Esmeralda (9085) and Copépodo (1595) Lakes. Esmeralda Lake was the richest in terms of number of taxa (59), with Copépodo and Pan Negro Lakes having 30 taxa each. Bray-Curtis dissimilarity among lakes was high (~ 0.80). The Order Chlamydomonadales (Chlorophyceae) gave the highest contribution in terms of numbers of taxa and DNA reads in all lakes. The most abundant taxon was Chlorococcum microstigmatum. The study confirms the utility of DNA metabarcoding in assessing potential green algal diversity in Antarctic lakes, generating new Antarctic records.

Validation of lipid-related therapeutic targets for coronary heart disease prevention using human genetics

Drug target Mendelian randomization (MR) studies use DNA sequence variants in or near a gene encoding a drug target, that alter the target’s expression or function, as a tool to anticipate the effect of drug action on the same target. Here we apply MR to prioritize drug targets for their causal relevance for coronary heart disease (CHD). The targets are further prioritized using independent replication, co-localization, protein expression profiles and data from the British National Formulary and clinicaltrials.gov. Out of the 341 drug targets identified through their association with blood lipids (HDL-C, LDL-C and triglycerides), we robustly prioritize 30 targets that might elicit beneficial effects in the prevention or treatment of CHD, including NPC1L1 and PCSK9, the targets of drugs used in CHD prevention. We discuss how this approach can be generalized to other targets, disease biomarkers and endpoints to help prioritize and validate targets during the drug development process.

Mitochondrial DNA sequence variation and risk of meningioma.

Risk factors for meningioma include female gender, African American race, high body mass index (BMI), and exposure to ionizing radiation. Although genome-wide association studies (GWAS) have identified two nuclear genome risk loci for meningioma (rs12770228 and rs2686876), the relation between mitochondrial DNA (mtDNA) sequence variants and meningioma is unknown. We examined the association of 42 common germline mtDNA variants (minor allele frequency ≥ 5%), haplogroups, and genes with meningioma in 1080 controls and 478 meningioma cases from a case-control study conducted at medical centers in the southeastern United States. Associations were examined separately for meningioma overall and by WHO grade (n = 409 grade I and n = 69 grade II/III). Overall, meningioma was significantly associated with being female (OR 2.85; 95% CI 2.21-3.69), self-reported African American race (OR 2.38, 95% CI 1.41-3.99), and being overweight (OR 1.48; 95% CI 1.11-1.97) or obese (OR 1.70; 95% CI 1.25-2.31). The variant m.16362T > C (rs62581341) in the mitochondrial control region was positively associated with grade II/III meningiomas (OR 2.33; 95% CI 1.14-4.77), but not grade I tumors (OR 0.99; 95% CI 0.64-1.53). Haplogroup L, a marker for African ancestry, was associated with meningioma overall (OR 2.92; 95% CI 1.01-8.44). However, after stratifying by self-reported race, this association was only apparent among the few self-reported Caucasians with this haplogroup (OR 6.35; 95% CI 1.56-25.9). No other mtDNA variant, haplogroup, or gene was associated with meningioma. Common mtDNA variants and major mtDNA haplogroups do not appear to have associations with the odds of developing meningioma.

Beyond GWAS: from simple associations to functional insights.

Each human, when born, has slightly different DNA sequences, which make each of us unique. The variations in DNA sequences are called genetic variants. The primary aim of genome-wide association study (GWAS) is to detect associations between genetic variants and human phenotypes. Since GWAS focuses on germ-line variants, there is no reverse causation. Therefore, GWAS is one of the few tools that can assess the causality of human diseases. In the past 10years, many large-scale GWAS have been conducted. Although the primary outputs of GWAS are just a series of statistics, its downstream analyses provided many insights beyond simple associations: the causal mechanisms for autoimmune diseases and shared etiology between diseases. Moreover, GWAS downstream analyses generated scores potentially helpful in predicting clinical outcomes of each patient. This review focuses on GWAS for autoimmune diseases and introduces significant achievements of its downstream analyses. We also provide future directions that potentially overcome current limitations. We restrict our discussion to common autoimmune diseases (e.g., rheumatoid arthritis) since rare Mendelian diseases possess distinct genetic etiologies and are not tested by GWAS.

Beyond Genes: Germline Disruption in the Etiology of Autism Spectrum Disorders

Investigations into the etiology of autism spectrum disorders have been largely confined to two realms: variations in DNA sequence and somatic developmental exposures. Here we suggest a third route—disruption of the germline epigenome induced by exogenous toxicants during a parent’s gamete development. Similar to cases of germline mutation, these molecular perturbations may produce dysregulated transcription of brain-related genes during fetal and early development, resulting in abnormal neurobehavioral phenotypes in offspring. Many types of exposures may have these impacts, and here we discuss examples of anesthetic gases, tobacco components, synthetic steroids, and valproic acid. Alterations in parental germline could help explain some unsolved phenomena of autism, including increased prevalence, missing heritability, skewed sex ratio, and heterogeneity of neurobiology and behavior.

Annotation of 1350 Common Genetic Variants of the 19 ALDH Multigene Family from Global Human Genome Aggregation Database (gnomAD).

Human aldehyde dehydrogenase (ALDH) is a multigene family with 19 functional members encoding a class of diverse but important enzymes for detoxification or biotransformation of different endogenous and exogenous aldehyde substrates. Genetic mutations in the ALDH genes can cause the accumulation of toxic aldehydes and abnormal carbonyl metabolism and serious human pathologies. However, the physiological functions and substrate specificity of many ALDH genes are still unknown. Although many genetic variants of the ALDH gene family exist in human populations, their phenotype or clinical consequences have not been determined. Using the most comprehensive global human Genome Aggregation Database, gnomAD, we annotated here 1350 common variants in the 19 ALDH genes. These 1350 common variants represent all known genetic polymorphisms with a variant allele frequency of ≥0.1% (or an expected occurrence of ≥1 carrier per 500 individuals) in any of the seven major ethnic groups recorded by gnomAD. We detailed 13 types of DNA sequence variants, their genomic positions, SNP ID numbers, and allele frequencies among the seven major ethnic groups worldwide for each of the 19 ALDH genes. For the 313 missense variants identified in the gnomAD, we used two software algorithms, Polymorphism Phenotyping (PolyPhen) and Sorting Intolerant From Tolerant (SIFT), to predict the consequences of the variants on the structure and function of the enzyme. Finally, gene constraint analysis was used to predict how well genetic mutations were tolerated by selection forces for each of the ALDH genes in humans. Based on the ratio of observed and expected variant numbers in gnomAD, the three ALDH1A gene members, ALDH1A1, ALDH1A2, and ALDH1A3, appeared to have the lowest tolerance for loss-of-function mutations as compared to the other ALDH genes (# observed/# expected ratio 0.15–0.26). These analyses suggest that the ALDH1A1, ALDH1A2, and ALDH1A3 enzymes may serve a more essential function as compared with the other ALDH enzymes; functional loss mutations are much less common in healthy human populations than expected. This informatic analysis may assist the research community in determining the physiological function of ALDH isozymes and associate common variants with clinical phenotypes.

Importin α2 association with chromatin: Direct DNA binding via a novel DNA-binding domain.

The nuclear transport of proteins is important for facilitating appropriate nuclear functions. The importin α family proteins play key roles in nuclear transport as transport receptors for copious nuclear proteins. Additionally, these proteins possess other functions, including chromatin association and gene regulation. However, these nontransport functions of importin α are not yet fully understood, especially their molecular-level mechanisms and consequences for functioning with chromatin. Here, we report the novel molecular characteristics of importin α binding to diverse DNA sequences in chromatin. We newly identified and characterized a DNA-binding domain-the Nucleic Acid Associating Trolley pole domain (NAAT domain)-in the N-terminal region of importin α within the conventional importin β binding (IBB) domain that is necessary for nuclear transport of cargo proteins. Furthermore, we found that the DNA binding of importin α synergistically coupled the recruitment of its cargo protein to DNA. This is the first study to delineate the interaction between importin α and chromatin DNA via the NAAT domain, indicating the bifunctionality of the importin α N-terminal region for nuclear transport and chromatin association.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation.

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

Isolation and molecular characterization of clinical and environmental dematiaceous fungi and relatives from Iran.

Background and Purpose: The frequency and genetic diversity of black fungi in environmental and clinical settings have not been fully studied in Iran. This study aimed to identify and evaluate intra- and inter-species DNA sequence variation and also understand the phylogenetic relationships of melanized fungi and relatives isolated from different geographical regions of Iran. Materials and Methods: In total, 111 clinical and environmental strains of dematiaceous fungi were isolated, and their internal transcribed spacer ribosomal DNA (rDNA) regions were sequenced and analyzed. Results: An inter-species nucleotide sequence diversity rate of 1 to 464 nucleotides was observed between the species. Intra-species differences were found in the strains of Alternaria alternata, Cladosporium cladosporioides, Alternaria tenuissima, Curvularia spicifera, Aureobasidium pullulans, Curvularia hawaiiensis, Neoscytalidium dimidiatum, Alternaria terricola, Alternaria chlamydospora, Didymella glomerata, and Drechslera dematioidea by 0–59, 0–22, 0–4, 0–4, 0–3, 0–2, 0–2, 0–2, 0–2, 0–1, and 0–1 nucleotide, respectively. Conclusion: The internal transcribed spacer rDNA is useful for the discrimination of several taxa of dematiaceous fungi. However, a better understanding of the taxonomy of species of Alternaria requires a larger rDNA region or a library of other gene sequences.

High levels of mitochondrial DNA sequence divergence among isolated populations of Fontigens (Truncatelloidea: Emmericiidae) in eastern USA

ABSTRACT Populations of the minute emmericiid gastropod Fontigens inhabit isolated springs, caves and groundwaters widely scattered across eastern and midwestern USA. Nine species are currently recognized, with a tenth species, Fontigens cryptica, of uncertain status. Here, we survey intrapopulation, interpopulation and interspecific sequence variation of cytochrome c oxidase subunit I among six of the nine species (13 populations, 87 individuals), documenting up to 14.4% sequence divergence among conspecific populations. Interspecific divergence ranges from 8.7% to 21.9%. We also include in our analysis a single putative individual of F. cryptica, recently sampled from Kentucky, evaluating whether the population it represents may be genetically distinct. Fit within a larger phylogenetic context, the 14.4% sequence divergence demonstrated between this individual and its nearest neighbour, F. bottimeri sampled from Washington, DC, supports the specific status of F. cryptica.

Small vs. large eggs: comparative population connectivity and demographic history along a depth gradient in deep-sea crangonid Argis shrimps

Abstract The crangonid shrimps Argis hozawai, A. lar and A. toyamaensis, co-distributed in the Sea of Japan, exhibit intriguing differences in geographical and bathymetric distributions and in reproductive biology. Argis hozawai (150–250 m depth) and A. lar (200–300 m) are broadly distributed in the north-western Pacific Ocean and spawn relatively large numbers of small eggs, whereas A. toyamaensis (250–2000 m) is distributed in the Sea of Japan and spawns a small number of large eggs. We examined the relationship between egg size and dispersal patterns in the deep sea by comparing genetic population structures using mitochondrial DNA sequence variation. We found little or no genetic divergence within the Sea of Japan for A. hozawai and A. lar, whereas there was a slight but significantly higher genetic differentiation in A. toyamaensis. This suggests that A. toyamaensis has lower dispersal ability than A. hozawai and A. lar, and therefore might maximize larval survival through larger size at hatching, with either direct or abbreviated larval development, to adapt to the deep-sea environment in the Sea of Japan. We also detected the effects of drastic environmental changes during the Pleistocene glacial periods on their demographic processes in the Sea of Japan.

DNA sequence and community structure diversity of multi-year soil fungi in Grape of Xinjiang

This study is designed to understand the community structure and diversity of fungi in the rhizosphere soil of grape. As the sample for this study, the rhizosphere soil of Crimson seedless grape with different planting years was collected from Shihezi in Xinjiang to carry out high-throughput sequencing, by which the complete sequence of soil fungi DNA was identified, and accordingly, the richness and diversity index of fungi were determined. The results showed that the dominant phyla of fungi in the grape rhizosphere soil with different planting years were Ascomycota and Basidiomycota, and the dominant classes of fungi were Sordariomycetes and Dothideomycetes. Soil organic matter, total potassium, total nitrogen and available phosphorus were the main soil fertility factors affecting the abundance and diversity of soil fungal communities, among which soil organic matter had the most significant influence. In addition, the fungal diversity and richness were highest in the middle layer (20–35 cm) of the grape rhizosphere soil with 12 planting years and lowest in the lower layer (35–50 cm) of the grape rhizosphere soil with 5 planting years. Linear discriminant analysis suggested that there were more biomarkers in the vineyard rhizosphere soil with 10 planting years, which meant there were more fungal communities with significant difference in the soil, especially in the middle layer (20–35). The results of this study can provide data reference and theoretical basis for improving vineyard soil quality, evaluating soil microecological effects and improving ecological environment of vineyard soil.

Epigenetics for Crop Improvement in Times of Global Change.

Simple SummaryResearch on plant epigenetics aims to understand how endogenous, biotic, and abiotic factors regulate plant development and growth independent of changes in the genome sequence. Often, the epigenetic changes are heritable across generations and modulate plant growth and crop tolerance, particularly in response to environmental stimuli. To take advantage of epigenetic adaptation, recent research has focused on implementing targeted epigenetic diversity to engineer plants that harbour advantageous traits for optimal crop production. Epigenetics has the potential to provide a powerful toolbox for crop breeders; however, most mechanistic studies are based on information from model plant species due to the challenges that arise when working with crops. Here, we summarise the contribution of epigenetics to optimising crop adaptation in response to climate change and overview potential future applications as well as challenges.Epigenetics has emerged as an important research field for crop improvement under the on-going climatic changes. Heritable epigenetic changes can arise independently of DNA sequence alterations and have been associated with altered gene expression and transmitted phenotypic variation. By modulating plant development and physiological responses to environmental conditions, epigenetic diversity—naturally, genetically, chemically, or environmentally induced—can help optimise crop traits in an era challenged by global climate change. Beyond DNA sequence variation, the epigenetic modifications may contribute to breeding by providing useful markers and allowing the use of epigenome diversity to predict plant performance and increase final crop production. Given the difficulties in transferring the knowledge of the epigenetic mechanisms from model plants to crops, various strategies have emerged. Among those strategies are modelling frameworks dedicated to predicting epigenetically controlled-adaptive traits, the use of epigenetics for in vitro regeneration to accelerate crop breeding, and changes of specific epigenetic marks that modulate gene expression of traits of interest. The key challenge that agriculture faces in the 21st century is to increase crop production by speeding up the breeding of resilient crop species. Therefore, epigenetics provides fundamental molecular information with potential direct applications in crop enhancement, tolerance, and adaptation within the context of climate change.

Immunomodulation—A Molecular Solution to Treating Patients with Severe Bladder Pain Syndrome?

BackgroundPatients with bladder pain syndrome experience debilitating pain and extreme frequency of urination. Numerous therapeutic approaches have been tested, but as the molecular basis of disease has remained unclear, specific therapies are not available. ObjectiveRecently, a systematic gene deletion strategy identified interleukin-1 (IL-1) hyperactivation as a cause of severe cystitis in a murine model. Treatment with an IL-1 receptor antagonist (IL-1RA) restored health in genetically susceptible mice, linking IL-1–dependent inflammation to pain and pathology in the bladder mucosa. The study objective was to investigate whether IL-1RA treatment might be beneficial in patients with bladder pain syndrome. Design, setting, and participantsPatients diagnosed with bladder pain syndrome were invited to participate and subjected to daily IL-1RA injections for 1 wk, followed by a treatment break. Patients with other urological disorders accompanied by pain were included as controls. Outcome measurements and statistical analysisWhen symptoms returned, treatment was resumed and responding patients were maintained on treatment long term, with individualized dosing regimens. Symptom scores were recorded and molecular effects were quantified by neuropeptide and gene expression analysis. DNA samples were subjected to exome genotyping. Results and limitationsIL-1RA treatment reduced bladder pain and the frequency of urination in 13/17 patients (p < 0.001). Substance P levels in urine were lowered, and responders returned to a more normal lifestyle. Neuroinflammatory-dependent and IL-1–dependent gene networks were inhibited, as well as regulators of innate immunity. Genotyping revealed disease-associated IL1R1, NLRP3, and IL1RN DNA sequence variants in the responders. Controls did not benefit from IL-1RA treatment, except for one patent with cystitis cystica. ConclusionsIn this clinical study, IL-1RA treatment is proposed to reduce chronic bladder pain, immediately and in the long term. Despite the limited number of study patients, the potent acute effect and lasting symptom relief indicate that this therapeutic approach may be worth exploring in controlled clinical trials. Patient summaryTreatment with an interleukin-1 (IL-1) receptor antagonist is proposed for treating bladder pain syndrome, as it can result in symptom relief and increase quality of life. Reduced neuroinflammation and IL-1 signaling provided molecular evidence of the treatment effects.

Machine learning-enabled phenotyping for GWAS and TWAS of WUE traits in 869 field-grown sorghum accessions.

Sorghum (Sorghum bicolor) is a model C4 crop made experimentally tractable by extensive genomic and genetic resources. Biomass sorghum is studied as a feedstock for biofuel and forage. Mechanistic modeling suggests that reducing stomatal conductance (gs) could improve sorghum intrinsic water use efficiency (iWUE) and biomass production. Phenotyping to discover genotype-to-phenotype associations remains a bottleneck in understanding the mechanistic basis for natural variation in gs and iWUE. This study addressed multiple methodological limitations. Optical tomography and a machine learning tool were combined to measure stomatal density (SD). This was combined with rapid measurements of leaf photosynthetic gas exchange and specific leaf area (SLA). These traits were the subject of genome-wide association study and transcriptome-wide association study across 869 field-grown biomass sorghum accessions. The ratio of intracellular to ambient CO2 was genetically correlated with SD, SLA, gs, and biomass production. Plasticity in SD and SLA was interrelated with each other and with productivity across wet and dry growing seasons. Moderate-to-high heritability of traits studied across the large mapping population validated associations between DNA sequence variation or RNA transcript abundance and trait variation. A total of 394 unique genes underpinning variation in WUE-related traits are described with higher confidence because they were identified in multiple independent tests. This list was enriched in genes whose Arabidopsis (Arabidopsis thaliana) putative orthologs have functions related to stomatal or leaf development and leaf gas exchange, as well as genes with nonsynonymous/missense variants. These advances in methodology and knowledge will facilitate improving C4 crop WUE.

Production and characterization of a human antisperm monoclonal antibody against CD52g for topical contraception in women.

BackgroundApproximately 40% of human pregnancies are unintended, indicating a need for more acceptable effective contraception methods. New antibody production systems make it possible to manufacture reagent-grade human monoclonal antibodies (mAbs) for clinical use. We used the Nicotiana platform to produce a human antisperm mAb and tested its efficacy for on-demand topical contraception.MethodsHeavy and light chain variable region DNA sequences of a human IgM antisperm antibody derived from an infertile woman were inserted with human IgG1 constant region sequences into an agrobacterium and transfected into Nicotiana benthamiana. The product, an IgG1 mAb [“Human Contraception Antibody” (HCA)], was purified on Protein A columns, and QC was performed using the LabChip GXII Touch protein characterization system and SEC-HPLC. HCA was tested for antigen specificity by immunofluorescence and western blot assays, antisperm activity by sperm agglutination and complement dependent sperm immobilization assays, and safety in a human vaginal tissue (EpiVaginal™) model.FindingsHCA was obtained at concentrations ranging from 0.4 to 4 mg/ml and consisted of > 90% IgG monomers. The mAb specifically reacted with a glycan epitope on CD52g, a glycoprotein produced in the male reproductive tract and found in abundance on sperm. HCA potently agglutinated sperm under a variety of relevant physiological conditions at concentrations ≥ 6.25 µg/ml, and mediated complement-dependent sperm immobilization at concentrations ≥ 1 µg/ml. HCA and its immune complexes did not induce inflammation in EpiVaginal™ tissue.InterpretationHCA, an IgG1 mAb with potent sperm agglutination and immobilization activity and a good safety profile, is a promising candidate for female contraception.FundingThis research was supported by grants R01 HD095630 and P50HD096957 from the National Institutes of Health.

Incidental detection of acquired variants in germline genetic and genomic testing: a points to consider statement of the American College of Medical Genetics and Genomics (ACMG)

With recent advances in DNA sequencing technology, it is now possible to begin to appreciate the full scope of DNA variation that arises over the course of an individual’s lifetime.1,2 Our understanding of how the human genome changes over time and in response to external exposures is growing with the improved availability of next-generation sequencing (NGS) based testing, including exome/genome sequencing of large patient cohorts. Clinical laboratories employing NGS-based methodologies can detect many types of DNA sequence variation including those that are present at a reduced variant allele fraction (VAF) (i.e., less than the 50% expected for a heterozygous germline finding).

Digital-Droplet PCR to Detect Indels Mutations in Genetically Modified Anopheline Mosquito Populations

Recent advances in mosquito genomics and genetic engineering technologies have fostered a need for quick and efficient methods for detecting targeted DNA sequence variation on a large scale. Specifically, detecting insertions and deletions (indels) at gene-edited sites generated by CRISPR guide RNA (gRNA)/Cas9-mediated non-homologous end-joining (NHEJ) is important for assessing the fidelity of the mutagenesis and the frequency of unintended changes. We describe here a protocol for digital-droplet PCR (ddPCR) that is well-suited for high-throughput NHEJ analysis. While this method does not produce data that identifies individual sequence variation, it provides a quantitative estimate of the sequence variation within a population. Additionally, with appropriate resources, this protocol can be implemented in a field-site laboratory setting more easily than next-generation or Sanger sequencing. ddPCR also has a faster turn-around time for results than either of those methods, which allows a more quick and complete analysis of genetic variation in wild populations during field trials of genetically-engineered organisms.

Digital-Droplet PCR to Detect Indels Mutations in Genetically Modified Anopheline Mosquito Populations.

Recent advances in mosquito genomics and genetic engineering technologies have fostered a need for quick and efficient methods for detecting targeted DNA sequence variation on a large scale. Specifically, detecting insertions and deletions (indels) at gene-edited sites generated by CRISPR guide RNA (gRNA)/Cas9-mediated non-homologous end-joining (NHEJ) is important for assessing the fidelity of the mutagenesis and the frequency of unintended changes. We describe here a protocol for digital-droplet PCR (ddPCR) that is well-suited for high-throughput NHEJ analysis. While this method does not produce data that identifies individual sequence variation, it provides a quantitative estimate of the sequence variation within a population. Additionally, with appropriate resources, this protocol can be implemented in a field-site laboratory setting more easily than next-generation or Sanger sequencing. ddPCR also has a faster turn-around time for results than either of those methods, which allows a more quick and complete analysis of genetic variation in wild populations during field trials of genetically-engineered organisms.

TnpAREP and REP sequences dissemination in bacterial genomes: REP recognition determinants.

REP, diverse palindromic DNA sequences found at high copy number in many bacterial genomes, have been attributed important roles in cell physiology but their dissemination mechanisms are poorly understood. They might represent non-autonomous transposable elements mobilizable by TnpAREP, the first prokaryotic domesticated transposase associated with REP. TnpAREP, fundamentally different from classical transposases, are members of the HuH superfamily and closely related to the transposases of the IS200/IS605 family. We previously showed that Escherichia coli TnpAREP processes cognate single stranded REP in vitro and that this activity requires the integrity of the REP structure, in particular imperfect palindromes interrupted by a bulge and preceded by a conserved DNA motif. A second group of REPs rather carry perfect palindromes, raising questions about how the latter are recognized by their cognate TnpAREP. To get insight into the importance of REP structural and sequence determinants in these two groups, we developed an in vitro activity assay coupled to a mutational analysis for three different TnpAREP/REP duos via a SELEX approach. We also tackled the question of how the cleavage site is selected. This study revealed that two TnpAREP groups have co-evolved with their cognate REPs and use different strategies to recognize their REP substrates.