Intragenomic Variation Research Articles

Empirical evidence for concerted evolution in the 18S rDNA region of the planktonic diatom genus Chaetoceros

Concerted evolution is a process of homogenisation of repetitive sequences within a genome through unequal crossing over and gene conversion. This homogenisation is never fully achieved because mutations always create new variants. Classically, concerted evolution has been detected as “noise” in electropherograms and these variants have been characterised through cloning and sequencing of subsamples of amplified products. However, this approach limits the number of detectable variants and provides no information about the abundance of each variant. In this study, we investigated concerted evolution by using environmental time-series metabarcoding data, single strain high-throughput sequencing (HTS) and a collection of Sanger reference barcode sequences. We used six species of the marine planktonic diatom genus Chaetoceros as study system. Abundance plots obtained from environmental metabarcoding and single strain HTS showed the presence of a haplotype far more abundant than all the others (the “dominant” haplotype) and identical to the reference sequences of that species obtained with Sanger sequencing. This distribution fitted best with Zipf’s law among the rank abundance/ dominance models tested. Furthermore, in each strain 99% of reads showed a similarity of 99% with the dominant haplotype, confirming the efficiency of the homogenisation mechanism of concerted evolution. We also demonstrated that minor haplotypes found in the environmental samples are not only technical artefacts, but mostly intragenomic variation generated by incomplete homogenisation. Finally, we showed that concerted evolution can be visualised inferring phylogenetic networks from environmental data. In conclusion, our study provides an important contribution to the understanding of concerted evolution and to the interpretation of DNA barcoding and metabarcoding data based on multigene family markers.

Population genomics and history of speciation reveal fishery management gaps in two related redfish species (Sebastes mentellaandSebastes fasciatus).

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter‐ and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

Molecular forms of Anopheles subpictus and Anopheles sundaicus in the Indian subcontinent

BackgroundAnopheles subpictus and Anopheles sundaicus are closely related species, each comprising several sibling species. Ambiguities exist in the classification of these two nominal species and the specific status of members of these species complexes. Identifying fixed molecular forms and mapping their spatial distribution will help in resolving the taxonomic ambiguities and understanding their relative epidemiological significance.MethodsDNA sequencing of Internal Transcribed Spacer-2 (ITS2), 28S-rDNA (D1-to-D3 domains) and cytochrome oxidase-II (COII) of morphologically identified specimens of two nominal species, An. subpictus sensu lato (s.l.) and An. sundaicus s.l., collected from the Indian subcontinent, was performed and subjected to genetic distance and molecular phylogenetic analyses.ResultsMolecular characterization of mosquitoes for rDNA revealed the presence of two molecular forms of An. sundaicus s.l. and three molecular forms of An. subpictus s.l. (provisionally designated as Form A, B and C) in the Indian subcontinent. Phylogenetic analyses revealed two distinct clades: (i) subpictus clade, with a single molecular form of An. subpictus (Form A) prevalent in mainland India and Sri Lanka, and (ii) sundaicus clade, comprising of members of Sundaicus Complex, two molecular forms of An. subpictus s.l. (Form B and C), prevalent in coastal areas or islands in Indian subcontinent, and molecular forms of An. subpictus s.l. reported from Thailand and Indonesia. Based on the number of float-ridges on eggs, all An. subpictus molecular Form B were classified as Species B whereas majority (80%) of the molecular Form A were classified as sibling species C. Fixed intragenomic sequence variation in ITS2 with the presence of two haplotypes was found in molecular Form A throughout its distribution.ConclusionA total of three molecular forms of An. subpictus s.l. and two molecular forms of An. sundaicus s.l. were recorded in the Indian subcontinent. Phylogenetically, two forms of An. subpictus s.l. (Form B and C) prevalent in coastal areas or islands in the Indian subcontinent and molecular forms reported from Southeast Asia are members of Sundaicus Complex. Molecular Form A of An. subpictus is distantly related to all other forms and deserve a distinct specific status.

High-throughput sequencing of 5S-IGS in oaks: Exploring intragenomic variation and algorithms to recognize target species in pure and mixed samples.

Measuring biological diversity is a crucial but difficult undertaking, as exemplified in oaks where complex patterns of morphological, ecological, biogeographical and genetic differentiation collide with traditional taxonomy, which measures biodiversity in number of species (or higher taxa). In this pilot study, we generated high-throughput sequencing amplicon data of the intergenic spacer of the 5S nuclear ribosomal DNA cistron (5S-IGS) in oaks, using six mock samples that differ in geographical origin, species composition and pool complexity. The potential of the marker for automated genotaxonomy applications was assessed using a reference data set of 1,770 5S-IGS cloned sequences, covering the entire taxonomic breadth and distribution range of western Eurasian Quercus, and applying similarity (blast) and evolutionary approaches (maximum-likelihood trees and Evolutionary Placement Algorithm). Both methods performed equally well, allowing correct identification of species in sections Ilex and Cerris in the pure and mixed samples, and main lineages shared by species of sect. Quercus. Application of different cut-off thresholds revealed that medium- to high-abundance (>10 or 25) sequences suffice for a net species identification of samples containing one or a few individuals. Lower thresholds identify phylogenetic correspondence with all target species in highly mixed samples (analogous to environmental bulk samples) and include rare variants pointing towards reticulation, incomplete lineage sorting, pseudogenic 5S units and in situ (natural) contamination. Our pipeline is highly promising for future assessments of intraspecific and interpopulation diversity, and of the genetic resources of natural ecosystems, which are fundamental to empower fast and solid biodiversity conservation programmes worldwide.

Genomic Diversity and Hotspot Mutations in 30,983 SARS-CoV-2 Genomes: Moving Toward a Universal Vaccine for the "Confined Virus"?

The COVID-19 pandemic has been ongoing since its onset in late November 2019 in Wuhan, China. Understanding and monitoring the genetic evolution of the virus, its geographical characteristics, and its stability are particularly important for controlling the spread of the disease and especially for the development of a universal vaccine covering all circulating strains. From this perspective, we analyzed 30,983 complete SARS-CoV-2 genomes from 79 countries located in the six continents and collected from 24 December 2019, to 13 May 2020, according to the GISAID database. Our analysis revealed the presence of 3206 variant sites, with a uniform distribution of mutation types in different geographic areas. Remarkably, a low frequency of recurrent mutations has been observed; only 169 mutations (5.27%) had a prevalence greater than 1% of genomes. Nevertheless, fourteen non-synonymous hotspot mutations (>10%) have been identified at different locations along the viral genome; eight in ORF1ab polyprotein (in nsp2, nsp3, transmembrane domain, RdRp, helicase, exonuclease, and endoribonuclease), three in nucleocapsid protein, and one in each of three proteins: Spike, ORF3a, and ORF8. Moreover, 36 non-synonymous mutations were identified in the receptor-binding domain (RBD) of the spike protein with a low prevalence (<1%) across all genomes, of which only four could potentially enhance the binding of the SARS-CoV-2 spike protein to the human ACE2 receptor. These results along with intra-genomic divergence of SARS-CoV-2 could indicate that unlike the influenza virus or HIV viruses, SARS-CoV-2 has a low mutation rate which makes the development of an effective global vaccine very likely.

Haplotype-resolved genome analyses of a heterozygous diploid potato

Potato (Solanum tuberosum L.) is the most important tuber crop worldwide. Efforts are underway to transform the crop from a clonally propagated tetraploid into a seed-propagated, inbred-line-based hybrid, but this process requires a better understanding of potato genome. Here, we report the 1.67-Gb haplotype-resolved assembly of a diploid potato, RH89-039-16, using a combination of multiple sequencing strategies, including circular consensus sequencing. Comparison of the two haplotypes revealed ~2.1% intragenomic diversity, including 22,134 predicted deleterious mutations in 10,642 annotated genes. In 20,583 pairs of allelic genes, 16.6% and 30.8% exhibited differential expression and methylation between alleles, respectively. Deleterious mutations and differentially expressed alleles were dispersed throughout both haplotypes, complicating strategies to eradicate deleterious alleles or stack beneficial alleles via meiotic recombination. This study offers a holistic view of the genome organization of a clonally propagated diploid species and provides insights into technological evolution in resolving complex genomes.

Planktonic foraminifera genomic variations reflect paleoceanographic changes in the Arctic: evidence from sedimentary ancient DNA

Deciphering the evolution of marine plankton is typically based on the study of microfossil groups. Cryptic speciation is common in these groups, and large intragenomic variations occur in ribosomal RNA genes of many morphospecies. In this study, we correlated the distribution of ribosomal amplicon sequence variants (ASVs) with paleoceanographic changes by analyzing the high-throughput sequence data assigned to Neogloboquadrina pachyderma in a 140,000-year-old sediment core from the Arctic Ocean. The sedimentary ancient DNA demonstrated the occurrence of various N. pachyderma ASVs whose occurrence and dominance varied through time. Most remarkable was the striking appearance of ASV18, which was nearly absent in older sediments but became dominant during the last glacial maximum and continues to persist today. Although the molecular ecology of planktonic foraminifera is still poorly known, the analysis of their intragenomic variations through time has the potential to provide new insight into the evolution of marine biodiversity and may lead to the development of new and important paleoceanographic proxies.

Satellitome Analysis in the Ladybird Beetle Hippodamia variegata (Coleoptera, Coccinellidae).

Hippodamia variegata is one of the most commercialized ladybirds used for the biological control of aphid pest species in many economically important crops. This species is the first Coccinellidae whose satellitome has been studied by applying new sequencing technologies and bioinformatics tools. We found that 47% of the H. variegata genome is composed of repeated sequences. We identified 30 satellite DNA (satDNA) families with a median intragenomic divergence of 5.75% and A+T content between 45.6% and 74.7%. This species shows satDNA families with highly variable sizes although the most common size is 100–200 bp. However, we highlight the existence of a satDNA family with a repeat unit of 2 kb, the largest repeat unit described in Coleoptera. PCR amplifications for fluorescence in situ hybridization (FISH) probe generation were performed for the four most abundant satDNA families. FISH with the most abundant satDNA family as a probe shows its pericentromeric location on all chromosomes. This location is coincident with the heterochromatin revealed by C-banding and DAPI staining, also analyzed in this work. Hybridization signals for other satDNA families were located only on certain bivalents and the X chromosome. These satDNAs could be very useful as chromosomal markers due to their reduced location.

Species Boundaries and Molecular Markers for the Classification of 16SrI Phytoplasmas Inferred by Genome Analysis.

Phytoplasmas are plant-pathogenic bacteria that impact agriculture worldwide. The commonly adopted classification system for phytoplasmas is based on the restriction fragment length polymorphism (RFLP) analysis of their 16S rRNA genes. With the increased availability of phytoplasma genome sequences, the classification system can now be refined. This work examined 11 strains in the 16SrI group within the genus ‘Candidatus Phytoplasma’ and investigated the possible species boundaries. We confirmed that the RFLP classification method is problematic due to intragenomic variation of the 16S rRNA genes and uneven weighing of different nucleotide positions. Importantly, our results based on the molecular phylogeny, differentiations in chromosomal segments and gene content, and divergence in homologous sequences, all supported that these strains may be classified into multiple operational taxonomic units (OTUs) equivalent to species. Strains assigned to the same OTU share >97% genome-wide average nucleotide identity (ANI) and >78% of their protein-coding genes. In comparison, strains assigned to different OTUs share < 94% ANI and < 75% of their genes. Reduction in homologous recombination between OTUs is one possible explanation for the discontinuity in genome similarities, and these findings supported the proposal that 95% ANI could serve as a cutoff for distinguishing species in bacteria. Additionally, critical examination of these results and the raw sequencing reads led to the identification of one genome that was presumably mis-assembled by combining two sequencing libraries built from phytoplasmas belonging to different OTUs. This finding provided a cautionary tale for working on uncultivated bacteria. Based on the new understanding of phytoplasma divergence and the current genome availability, we developed five molecular markers that could be used for multilocus sequence analysis (MLSA). By selecting markers that are short yet highly informative, and are distributed evenly across the chromosome, these markers provided a cost-effective system that is robust against recombination. Finally, examination of the effector gene distribution further confirmed the rapid gains and losses of these genes, as well as the involvement of potential mobile units (PMUs) in their molecular evolution. Future improvements on the taxon sampling of phytoplasma genomes will allow further expansions of similar analysis, and thus contribute to phytoplasma taxonomy and diagnostics.

Unambiguous identification of fungi: where do we stand and how accurate and precise is fungal DNA barcoding?

ABSTRACTTrue fungi (Fungi) and fungus-like organisms (e.g. Mycetozoa, Oomycota) constitute the second largest group of organisms based on global richness estimates, with around 3 million predicted species. Compared to plants and animals, fungi have simple body plans with often morphologically and ecologically obscure structures. This poses challenges for accurate and precise identifications. Here we provide a conceptual framework for the identification of fungi, encouraging the approach of integrative (polyphasic) taxonomy for species delimitation, i.e. the combination of genealogy (phylogeny), phenotype (including autecology), and reproductive biology (when feasible). This allows objective evaluation of diagnostic characters, either phenotypic or molecular or both. Verification of identifications is crucial but often neglected. Because of clade-specific evolutionary histories, there is currently no single tool for the identification of fungi, although DNA barcoding using the internal transcribed spacer (ITS) remains a first diagnosis, particularly in metabarcoding studies. Secondary DNA barcodes are increasingly implemented for groups where ITS does not provide sufficient precision. Issues of pairwise sequence similarity-based identifications and OTU clustering are discussed, and multiple sequence alignment-based phylogenetic approaches with subsequent verification are recommended as more accurate alternatives. In metabarcoding approaches, the trade-off between speed and accuracy and precision of molecular identifications must be carefully considered. Intragenomic variation of the ITS and other barcoding markers should be properly documented, as phylotype diversity is not necessarily a proxy of species richness. Important strategies to improve molecular identification of fungi are: (1) broadly document intraspecific and intragenomic variation of barcoding markers; (2) substantially expand sequence repositories, focusing on undersampled clades and missing taxa; (3) improve curation of sequence labels in primary repositories and substantially increase the number of sequences based on verified material; (4) link sequence data to digital information of voucher specimens including imagery. In parallel, technological improvements to genome sequencing offer promising alternatives to DNA barcoding in the future. Despite the prevalence of DNA-based fungal taxonomy, phenotype-based approaches remain an important strategy to catalog the global diversity of fungi and establish initial species hypotheses.

Phylogenomic incongruence in Ceratocystis: a clue to speciation?

BackgroundThe taxonomic history of Ceratocystis, a genus in the Ceratocystidaceae, has been beset with questions and debate. This is due to many of the commonly used species recognition concepts (e.g., morphological and biological species concepts) providing different bases for interpretation of taxonomic boundaries. Species delineation in Ceratocystis primarily relied on genealogical concordance phylogenetic species recognition (GCPSR) using multiple standard molecular markers.ResultsQuestions have arisen regarding the utility of these markers e.g., ITS, BT and TEF1-α due to evidence of intragenomic variation in the ITS, as well as genealogical incongruence, especially for isolates residing in a group referred to as the Latin-American clade (LAC) of the species. This study applied a phylogenomics approach to investigate the extent of phylogenetic incongruence in Ceratocystis. Phylogenomic analyses of a total of 1121 shared BUSCO genes revealed widespread incongruence within Ceratocystis, particularly within the LAC, which was typified by three equally represented topologies. Comparative analyses of the individual gene trees revealed evolutionary patterns indicative of hybridization. The maximum likelihood phylogenetic tree generated from the concatenated dataset comprised of 1069 shared BUSCO genes provided improved phylogenetic resolution suggesting the need for multiple gene markers in the phylogeny of Ceratocystis.ConclusionThe incongruence observed among single gene phylogenies in this study call into question the utility of single or a few molecular markers for species delineation. Although this study provides evidence of interspecific hybridization, the role of hybridization as the source of discordance will require further research because the results could also be explained by high levels of shared ancestral polymorphism in this recently diverged lineage. This study also highlights the utility of BUSCO genes as a set of multiple orthologous genes for phylogenomic studies.

Lost in translation: Population genomics and long-read sequencing reveals relaxation of concerted evolution of the ribosomal DNA cistron

Concerted evolution of the ribosomal DNA array has been studied in numerous eukaryotic taxa, yet is still poorly understood. rDNA genes are repeated dozens to hundreds of times in the eukaryotic genome (Eickbush and Eickbush, 2007) and it is believed that these arrays are homogenized through concerted evolution (Zimmer et al., 1980; Dover, 1993) preventing the accumulation of intragenomic, and intraspecific, variation. However, numerous studies have reported rampant intragenomic and intraspecific variation in the rDNA array (Ganley and Kobayashi, 2011; Naidoo et al., 2013; Hughes and Petersen, 2001; Lindner and Banik, 2011; Li et al., 2013; Lindner et al., 2013; Hughes et al., 2018), contradicting our current understanding of concerted evolution. The internal transcribed spacers (ITS) of the rDNA cistron are the most commonly used DNA barcoding region in Fungi (Schoch et al., 2012), and rely on concerted evolution to homogenize the rDNA array leading to a “barcode gap” (Puillandre et al., 2012). Here we show that in Boletus edulis Bull., ITS intragenomic variation persists at low allele frequencies throughout the rDNA array, this variation does not correlate with genomic relatedness between populations, and rDNA genes may not evolve in a strictly concerted fashion despite the presence of unequal recombination and gene conversion. Under normal assumptions, heterozygous positions found in ITS sequences represent hybridization between populations, yet through allelic mapping of the rDNA array we found numerous heterozygous alleles to be stochastically introgressed throughout, presenting a dishonest signal of gene flow. Moreover, despite the signal of gene flow in ITS, our organisms were highly inbred, indicating a disconnect between true gene flow and barcoding signals. In addition, we show that while the mechanisms of concerted evolution are ongoing in pseudo-heterozygous individuals, they are not fully homogenizing the ITS array. Concerted evolution of the rDNA array may insufficiently homogenize the ITS gene, allowing for misleading signals of gene flow to persist, vastly complicating the use of the ITS locus for DNA barcoding in Fungi.

The impact of intragenomic rRNA variation on metabarcoding‐derived diversity estimates: A case study from marine nematodes

AbstractIntragenomic rRNA variation is a critical concern for eukaryotic metabarcoding studies, due to its potentially confounding effects on species delimitation and biodiversity estimates derived from ‐Omics data. In the present study, we assessed patterns associated with 18S rRNA metabarcoding loci in marine nematodes, including characterization of intragenomic rRNA gene variants (number of variants and abundance profiles) and aspects of datasets that can obscure biological signals (e.g., amplification of nontarget DNA, ambiguous taxonomy assignments). We estimated amplicon sequence variants (ASVs) using DADA2 from an 18S rRNA metabarcoding dataset (Illumina MiSeq) generated from individual marine nematodes. Illumina data were analyzed in conjunction with nematode morphological identifications and nearly full‐length 18S reference sequences (~1,600 bp Sanger barcodes) generated for a subset of the same specimens. Our results indicated that levels of intragenomic rRNA variation appeared to vary widely across nematode taxa (irrespective of phylogenetic clades or ecological feeding groups) and that coamplification of nontarget DNA was common (relic DNA, gut contents, etc.). The DADA2 pipeline appeared to produce a biologically accurate profile of intragenomic rRNA variants in nematodes that was consistent with “Head‐Tail” patterns (of dominant vs. minor rRNA gene variants) identified in previous studies. Although intragenomic rRNA variation appears to be ubiquitous in marine nematodes, nematode identifications were highly congruent across our three methods for species delimitation (traditional morphological taxonomy, Sanger DNA barcoding, and high‐throughput metabarcoding). In spite of pervasive intragenomic variation and high copy number of rRNA genes, the most abundant ASVs in metabarcoding datasets are likely to represent true species barcodes and thus confer an accurate view of extant biodiversity. However, our findings also emphasize the importance of applying bioinformatic filtering techniques and developing well‐curated reference databases in order to better link rRNA molecules with specimen‐level data and alleviate the confounding effects of intragenomic gene variants in studies of microbial eukaryotes.

Characterizing the ribosomal tandem repeat and its utility as a DNA barcode in lichen-forming fungi

BackgroundRegions within the nuclear ribosomal operon are a major tool for inferring evolutionary relationships and investigating diversity in fungi. In spite of the prevalent use of ribosomal markers in fungal research, central features of nuclear ribosomal DNA (nrDNA) evolution are poorly characterized for fungi in general, including lichenized fungi. The internal transcribed spacer (ITS) region of the nrDNA has been adopted as the primary DNA barcode identification marker for fungi. However, little is known about intragenomic variation in the nrDNA in symbiotic fungi. In order to better understand evolution of nrDNA and the utility of the ITS region for barcode identification of lichen-forming fungal species, we generated nearly complete nuclear ribosomal operon sequences from nine species in the Rhizoplaca melanophthalma species complex using short reads from high-throughput sequencing.ResultsWe estimated copy numbers for the nrDNA operon, ranging from nine to 48 copies for members of this complex, and found low levels of intragenomic variation in the standard barcode region (ITS). Monophyly of currently described species in this complex was supported in phylogenetic inferences based on the ITS, 28S, intergenic spacer region, and some intronic regions, independently; however, a phylogenetic inference based on the 18S provided much lower resolution. Phylogenetic analysis of concatenated ITS and intergenic spacer sequence data generated from 496 specimens collected worldwide revealed previously unrecognized lineages in the nrDNA phylogeny.ConclusionsThe results from our study support the general assumption that the ITS region of the nrDNA is an effective barcoding marker for fungi. For the R. melanophthalma group, the limited amount of potential intragenomic variability in the ITS region did not correspond to fixed diagnostic nucleotide position characters separating taxa within this species complex. Previously unrecognized lineages inferred from ITS sequence data may represent undescribed species-level lineages or reflect uncharacterized aspects of nrDNA evolution in the R. melanophthalma species complex.

A comparison between the morphological and molecular approach to identify the benthic diatom community in the St Lucia Estuary, South Africa

The use of diatoms as bioindicators to detect human-induced change is a globally accepted monitoring tool. DNA metabarcoding is an effective complementary tool to standard microscopic methods for species identification, providing reliable and timeous diatom assessments in a range of aquatic ecosystems. This study evaluated the suitability of the DNA metabarcoding approach for diatom monitoring in the St Lucia Estuary, South Africa. Results obtained from the standard microscopic (morphological) and molecular methods were compared and indicated a low similarity in the vegetated (11%), unvegetated (12%) and epiphytic (12%) habitats at species level for the epipelic and epiphytic diatom communities. The difference between the methods was expected given the challenges pertaining to the incomplete molecular reference database, intragenomic variability, physiological differences (number of chloroplasts) and highly variable microphytobenthos biomass. While the molecular method was unable to identify the presence of all the dominant diatom species (routinely used as ecological indicators), the available molecular operational taxonomic units (MOTUs) provided an effective complementary tool to determine the relevant community diversity estimates required for the application of the South African estuarine health index (EHI).

Widespread prevalence but contrasting patterns of intragenomic rRNA polymorphisms in nematodes: Implications for phylogeny, species delimitation and life history inference.

Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free-living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population-level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.

First-tier detection of intragenomic 16S rRNA gene variation in culturable endophytic bacteria from cacao seeds.

BackgroundIntragenomic variability in 16S rDNA is a limiting factor for taxonomic and diversity characterization of Bacteria, and studies on its occurrence in natural/environmental populations are scarce. In this work, direct DNA amplicon sequencing coupled with frequent-cutter restriction analysis allowed detection of intragenomic 16S rDNA variation in culturable endophytic bacteria from cacao seeds in a fast and attractive manner.MethodsTotal genomic DNA from 65 bacterial strains was extracted and the 16S rDNA hyper variable V5–V9 regions were amplified for enzyme digestion and direct Sanger-type sequencing. The resulting electropherograms were visually inspected and compared to the corresponding AluI-restriction profiles, as well as to complete genome sequences in databases. Restriction analysis were employed to substitute the need of amplicon cloning and re-sequencing. A specifically improved polyacrylamide-gradient electrophoresis allowed to resolve 5-bp differences in restriction fragment sizes. Chi-square analysis on 2 × 2 contingency table tested for the independence between the ‘number of AluI bands’ and ‘type of eletropherogram’.ResultsTwo types of electropherograms were obtained: unique template, with single peaks per base (clean chromatograms), and heterogeneous template, with various levels of multiple peaks per base (mixed chromatograms). Statistics revealed significant interaction between number of restriction fragments and type of electropherogram for the same amplicons: clean or mixed ones associated to ≤5 or ≥6 bands, respectively. The mixed-template pattern combined with the AluI-restriction profiles indicated a high proportion of 49% of the culturable endophytes from a tropical environment showing evidence of intragenomic 16S rDNA heterogeneity.ConclusionThe approach presented here was useful for a rapid, first-tier detection of intragenomic variation in culturable isolates, which can be applied in studies of other natural populations; a preliminary view of intragenomic heterogeneity levels can complement culture-dependent and -independent methods. Consequences of these findings in taxonomic and diversity studies in complex bacterial communities are discussed.

Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis

The 16S rRNA gene has been a mainstay of sequence-based bacterial analysis for decades. However, high-throughput sequencing of the full gene has only recently become a realistic prospect. Here, we use in silico and sequence-based experiments to critically re-evaluate the potential of the 16S gene to provide taxonomic resolution at species and strain level. We demonstrate that targeting of 16S variable regions with short-read sequencing platforms cannot achieve the taxonomic resolution afforded by sequencing the entire (~1500 bp) gene. We further demonstrate that full-length sequencing platforms are sufficiently accurate to resolve subtle nucleotide substitutions (but not insertions/deletions) that exist between intragenomic copies of the 16S gene. In consequence, we argue that modern analysis approaches must necessarily account for intragenomic variation between 16S gene copies. In particular, we demonstrate that appropriate treatment of full-length 16S intragenomic copy variants has the potential to provide taxonomic resolution of bacterial communities at species and strain level.

Comparative Analysis of Intra- and Inter-Specific Genomic Variability in the Peach Potato Aphid, Myzus persicae.

The availability of genomic data in the last decade relating to different aphid species has allowed the analysis of the genomic variability occurring among such species, whereas intra-specific variability has hitherto very largely been neglected. In order to analyse the intra-genomic variability in the peach potato aphid, Myzus persicae, comparative analyses were performed revealing several clone-specific gene duplications, together with numerous deletions/rearrangements. Our comparative approach also allowed us to evaluate the synteny existing between the two M. persicae clones tested and between the peach potato aphid and the pea aphid, Acyrthosiphon pisum. Even if part of the observed rearrangements are related to a low quality of some assembled contigs and/or to the high number of contigs present in these aphid genomes, our evidence reveals that aphid clones are genetically more different than expected. These results suggest that the choice of performing genomes sequencing combining different biotypes/populations, as revealed in the case of the soybean aphid, Aphis glycines, is unlikely to be very informative in aphids. Interestingly, it is possible that the holocentric nature of aphid chromosomes favours genome rearrangements that can be successively inherited transgenerationally via the aphid’s apomictic (parthenogenetic) mode of reproduction. Lastly, we evaluated the structure of the cluster of genes coding for the five histones (H1, H2A, H2B, H3 and H4) in order to better understand the quality of the two M. persicae genomes and thereby to improve our knowledge of this functionally important gene family.

Intragenomic variability and extended sequence patterns in the mutational signature of ultraviolet light

Mutational signatures can reveal properties of underlying mutational processes and are important when assessing signals of selection in cancer. Here, we describe the sequence characteristics of mutations induced by ultraviolet (UV) light, a major mutagen in several human cancers, in terms of extended (longer than trinucleotide) patterns as well as variability of the signature across chromatin states. Promoter regions display a distinct UV signature with reduced TCG > TTG transitions, and genome-wide mapping of UVB-induced DNA photoproducts (pyrimidine dimers) showed that this may be explained by decreased damage formation at hypomethylated promoter CpG sites. Further, an extended signature model encompassing additional information from longer contextual patterns improves modeling of UV mutations, which may enhance discrimination between drivers and passenger events. Our study presents a refined picture of the UV signature and underscores that the characteristics of a single mutational process may vary across the genome.