Patterns Of Sequence Variation Research Articles

Patterns of occurrence, phenology, and phylogeny of Phloeosinus punctatusLeConte(Coleoptera: Curculionidae, Scolytinae) in giant sequoia.

Here, we describe patterns of reproduction and flight phenology of putative Phloeosinus punctatus in giant sequoia groves and compare morphology and genotypes of beetles from sympatric giant sequoia (Sequoiadendron giganteum) and California incense-cedar (Calocedrus decurrens). Surveys conducted in 2022 revealed that numerous branches fall from giant sequoia crowns (on average ~30 branches/tree), with 20%-50% of trees per site shedding branches, depositing breeding material for beetles on the forest floor that subsequently becomes colonized. When noninfested branches cut from mature giant sequoias were placed at the ground surface, they were colonized by P. punctatus and produced an average of 28 beetles/kg branch. Climbing and examination of sequoia crowns in 2023 showed that 75% of mature trees across 11 groves showed evidence of adult beetle entrance holes in their crowns. In 2021, tests with sticky traps showed that beetles alighted on fallen branches from 20th May to 20th August (peak landing: 2nd July); a logistic model developed from emergence data in 2021 and 2022 predicts the emergence of F1 offspring from branches between 10th July and 1st September (peak emergence: 8th August). Beetles emerging from giant sequoia preferred to settle on giant sequoia, did not reproduce in incense-cedar, and diverged morphologically from beetles emerging from incense-cedar. However, phylogenetic analysis of three genes (28S, CAD, and COI) revealed no clear pattern of sequence divergence, suggesting a single species (P. punctatus) that colonizes both hosts, though cryptic speciation may not be detectable with standard barcoding genes. Ecological and potential management implications are discussed.

Genome-wide patterns of noncoding and protein-coding sequence variation in the major fungal pathogen Aspergillus fumigatus.

Aspergillus fumigatus is a deadly fungal pathogen, responsible for >400,000 infections/year and high mortality rates. A. fumigatus strains exhibit variation in infection-relevant traits, including in their virulence. However, most A. fumigatus protein-coding genes, including those that modulate its virulence, are shared between A. fumigatus strains and closely related nonpathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the noncoding regions immediately upstream to the start codons of genes, which could reflect differences in gene regulation between strains. To begin testing this hypothesis, we identified 5,812 single-copy orthologs across the genomes of 263 A. fumigatus strains. In general, A. fumigatus noncoding regions showed higher levels of sequence variation compared with their corresponding protein-coding regions. Focusing on 2,482 genes whose protein-coding sequence identity scores ranged between 75 and 99%, we identified 478 total genes with signatures of positive selection only in their noncoding regions and 65 total genes with signatures only in their protein-coding regions. Twenty-eight of the 478 noncoding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Noncoding region variation between A. fumigatus strains included single-nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that noncoding regions of A. fumigatus genes harbor greater sequence variation than protein-coding regions, raising the hypothesis that this variation may contribute to A. fumigatus phenotypic heterogeneity.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Conserved patterns of sequence diversification provide insight into the evolution of two-component systems in Enterobacteriaceae.

Two-component regulatory systems (TCSs) are a major mechanism used by bacteria to sense and respond to their environments. Many of the same TCSs are used by biologically diverse organisms with different regulatory needs, suggesting that the functions of TCS must evolve. To explore this topic, we analysed the amino acid sequence divergence patterns of a large set of broadly conserved TCS across different branches of Enterobacteriaceae, a family of Gram-negative bacteria that includes biomedically important genera such as Salmonella, Escherichia, Klebsiella and others. Our analysis revealed trends in how TCS sequences change across different proteins or functional domains of the TCS, and across different lineages. Based on these trends, we identified individual TCS that exhibit atypical evolutionary patterns. We observed that the relative extent to which the sequence of a given TCS varies across different lineages is generally well conserved, unveiling a hierarchy of TCS sequence conservation with EnvZ/OmpR as the most conserved TCS. We provide evidence that, for the most divergent of the TCS analysed, PmrA/PmrB, different alleles were horizontally acquired by different branches of this family, and that different PmrA/PmrB sequence variants have highly divergent signal-sensing domains. Collectively, this study sheds light on how TCS evolve, and serves as a compendium for how the sequences of the TCS in this family have diverged over the course of evolution.

Contrasting patterns of sequence variation in steelhead populations reflect distinct evolutionary processes.

Multiple evolutionary processes influence genome-wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole-genome resequencing of 74 populations of steelhead revealed genome-wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes.

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Transcriptional misexpression in hybrids between species linked by gene flow is associated with patterns of sequence divergence.

The extent to which hybridization disrupts a gene's pattern of expression likely governs its propensity for introgression, while its extent of molecular divergence can itself underlie such disruption. Together, these phenomena shape the landscape of sequence and transcriptional divergence across the genome as species diverge. To understand this process, we characterize gene expression inheritance, regulatory divergence, and molecular divergence in the reproductive transcriptomes of species linked by gene flow: the fruit flies Anastrepha fraterculus and A. obliqua, which show evidence of gene flow despite clear evolutionary divergence. We find that their transcriptional patterns are a mosaic between those typically observed within and between allopatric species. Transcripts showing transgressive expression in hybrids or cis-regulatory divergence between species are associated with greater sequence divergence. This may reflect pleiotropic constraints that make them resistant to gene flow or they may be more likely to experience divergent selection. While these more divergent gene classes are likely to be important contributors to species differences, they are relatively rare. Instead, most differentially regulated transcripts, including those linked to reproduction, show high degrees of dominance in hybrids and trans-regulated divergence between species, suggesting widespread genetic compatibility that potentially allowed for introgression. These findings provide insights into how postzygotic isolating mechanisms might evolve in the presence of gene flow: regions showing cis-regulatory divergence or transgressive expression contribute to reproductive isolation, while regions with dominant expression and trans-regulatory divergence allow for introgression. These patterns create a genomic mosaic of transcriptional regulation that is tied to sequence divergence.

Hepatocytes undergo punctuated expansion dynamics from a periportal stem cell niche in normal human liver

While normal human liver is thought to be generally quiescent, clonal hepatocyte expansions have been observed, though neither their cellular source nor their expansion dynamics have been determined. Knowing the hepatocyte cell of origin, and their subsequent dynamics and trajectory within the human liver will provide an important basis to understand disease-associated dysregulation. Herein, we use invivo lineage tracing and methylation sequence analysis to demonstrate normal human hepatocyte ancestry. We exploit next-generation mitochondrial sequencing to determine hepatocyte clonal expansion dynamics across spatially distinct areas of laser-captured, microdissected, clones, in tandem with computational modelling in morphologically normal human liver. Hepatocyte clones and rare SOX9+ hepatocyte progenitors commonly associate with portal tracts and we present evidence that clones can lineage-trace with cholangiocytes, indicating the presence of a bipotential common ancestor at this niche. Within clones, we demonstrate methylation CpG sequence diversity patterns indicative of periportal not pericentral ancestral origins, indicating a portal to central vein expansion trajectory. Using spatial analysis of mitochondrial DNA variants by next-generation sequencing coupled with mathematical modelling and Bayesian inference across the portal-central axis, we demonstrate that patterns of mitochondrial DNA variants reveal large numbers of spatially restricted mutations in conjunction with limited numbers of clonal mutations. These datasets support the existence of a periportal progenitor niche and indicate that clonal patches exhibit punctuated but slow growth, then quiesce, likely due to acute environmental stimuli. These findings crucially contribute to our understanding of hepatocyte dynamics in the normal human liver. The liver is mainly composed of hepatocytes, but we know little regarding the source of these cells or how they multiply over time within the disease-free human liver. In this study, we determine a source of new hepatocytes by combining many different lab-based methods and computational predictions to show that hepatocytes share a common cell of origin with bile ducts. Both our experimental and computational data also demonstrate hepatocyte clones are likely to expand in slow waves across the liver in a specific trajectory, but often lie dormant for many years. These data show for the first time the expansion dynamics of hepatocytes in normal liver and their cell of origin enabling the accurate measurment of changes to their dynamics that may lead to liver disease. These findings are important for researchers determining cancer risk in human liver.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures.

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

Population Genomics of Plasmodium falciparum and Malaria Control Implications in Abidjan (Cote d’ivoire)

Introduction: The onset of Plasmodium falciparum (P. falciparum) resistance to antimalarial drugs requires the careful surveillance of African parasite populations. Genomic tools are implemented to detect evolutionary changes that could impact malaria control and elimination strategies. Here, we evaluate the genome-wide pattern of selection and sequence variation in P. falciparum populations in Abidjan, Côte d’Ivoire. Methods: The study was conducted in three localities of Abidjan from 2013 to 2014. We collected 70 blood samples after receiving written informed consent from patients above two years of age. After extracting P. falciparum and human DNA from isolates, we performed whole-genome sequencing and used population genomics approaches to investigate the genetic diversity and complexity of infections and identify loci under positive directional selection. Results: We observed an excess of rare variants in the population, showing a clear mutation process in the isolates. Moderate Fst estimates (0.3) was detected for surfin, an immune invasion gene family. Seven iHS regions that had at least two SNPs with a score > 3.2 were identified. These regions code for genes that have been under strong directional selection. Two of these genes were the chloroquine resistance transporter (crt) on chromosome 7 and the dihydropteroate reductase (dhps) on chromosome 8. Our analyses showed that a recent selective sweep occurred for the erythrocyte membrane protein (Pfemp1). Conclusion: Our analyses identified genes under selective drug pressure and balancing selection on protective immune-specific genes. These findings demonstrate the effectiveness of genomics analyses to follow the evolution of malaria parasites and adopt appropriate strategies to eliminate malaria in Côte d’Ivoire.

Molecular phylogenetics shed light on polyploid speciation in gorses (Ulex, Fabaceae: Genisteae) and on the origin of the invasive Ulex europaeus

Abstract Hybridization and polyploidy are key evolutionary forces in plant diversification, and their co-occurrence in the context of allopolyploid speciation is often associated with increased ability to colonize new environments and invasiveness. In the genus Ulex (Fabaceae), the European gorse (Ulex europaeus subsp. europaeus) is the only invasive and the only polyploid that has recently spread in different eco-geographical regions across the world. Understanding what confers such ecological advantages to this species, compared to its diploid and polyploid congeners, first requires clarification of the ecogeographical and evolutionary context of its formation. To achieve this, the geographical distributions of all Ulex spp. were estimated from species occurrence records, and phylogenetic analyses including all Ulex spp. were performed based on four nuclear (ITS and ETS nrDNA) and plastid (rps12 intron and trnK-matK) regions. The resulting trees were dated using a secondary calibration. Patterns of DNA sequence variation and dated phylogenetic trees were then interpreted in light of previous knowledge of chromosome numbers in Ulex to infer past events of polyploid speciation in the genus. We show that: (1) most current Ulex spp. radiated in the Iberian Peninsula during the past 1–2 Myr; (2) the history of Ulex was punctuated by multiple whole-genome duplication events; and (3) U. europaeus subsp. europaeus is the only gorse taxon that was formed by hybridization of two well-differentiated lineages (which separated c. 5 Mya) with wide climatic ranges (currently represented by Ulex minor and Ulex europaeus subsp. latebracteatus), possibly contributing to the invasive nature and wider climatic range of U. europaeus subsp. europaeus. These findings provide a much-needed evolutionary framework in which to explore the adaptive consequences of genome mergers and duplication in Ulex.

Allelic variation of Escherichia coli outer membrane protein A: Impact on cell surface properties, stress tolerance and allele distribution.

Outer membrane protein A (OmpA) is one of the most abundant outer membrane proteins of Gram-negative bacteria and is known to have patterns of sequence variations at certain amino acids-allelic variation-in Escherichia coli. Here we subjected seven exemplar OmpA alleles expressed in a K-12 (MG1655) ΔompA background to further characterization. These alleles were observed to significantly impact cell surface charge (zeta potential), cell surface hydrophobicity, biofilm formation, sensitivity to killing by neutrophil elastase, and specific growth rate at 42°C and in the presence of acetate, demonstrating that OmpA is an attractive target for engineering cell surface properties and industrial phenotypes. It was also observed that cell surface charge and biofilm formation both significantly correlate with cell surface hydrophobicity, a cell property that is increasingly intriguing for bioproduction. While there was poor alignment between the observed experimental values relative to the known sequence variation, differences in hydrophobicity and biofilm formation did correspond to the identity of residue 203 (N vs T), located within the proposed dimerization domain. The relative abundance of the (I, δ) allele was increased in extraintestinal pathogenic E. coli (ExPEC) isolates relative to environmental isolates, with a corresponding decrease in (I, α) alleles in ExPEC relative to environmental isolates. The (I, α) and (I, δ) alleles differ at positions 203 and 251. Variations in distribution were also observed among ExPEC types and phylotypes. Thus, OmpA allelic variation and its influence on OmpA function warrant further investigation.

Patterns of DNA barcode diversity in butterfly species (Lepidoptera) introduced to the Nearctic

One of the main consequences of globalization is the intensification of biological introductions. Because of their negative impact on environments, the early detection and monitoring of introduced species through molecular approaches is gaining increased uptake. This study assembles 2,278 DNA barcode records to examine contemporary patterns of sequence variation in mitochondrial cytochrome c oxidase I (COI) in five butterfly species introduced to the Nearctic, with a focus on Pieris rapae Linnaeus (Lepidoptera: Pieridae) and Thymelicus lineola Ochsenheimer (Lepidoptera: Hesperiidae). Parameters of genetic diversity were low (i.e., h < 0.606, π < 0.0039) for Nearctic populations of all analyzed species. Those of P. rapae and T. lineola showed marked genetic differentiation from their source populations in the Palearctic. Haplotype distributions in their Nearctic populations exposed a starburst pattern with a few common haplotypes known from Palearctic, and infrequent haplotypes diverging from them at only one or two nucleotide sites. Some uncommon haplotypes were only found in the Nearctic suggesting they originated after invasion, while others also occur in the Palearctic. This study provides an example of genetic paradox of invasion, where species often rapidly expand their distribution and become dominant in the new habitat despite their depleted levels of sequence variation.

A Putative long-range RNA-RNA interaction between ORF8and Spike of SARS-CoV-2.

SARS-CoV-2 has affected people worldwide as the causative agent of COVID-19. The virus is related to the highly lethal SARS-CoV-1 responsible for the 2002-2003 SARS outbreak in Asia. Research is ongoing to understand why both viruses have different spreading capacities and mortality rates. Like other beta coronaviruses, RNA-RNA interactions occur between different parts of the viral genomic RNA, resulting in discontinuous transcription and production of various sub-genomic RNAs. These sub-genomic RNAs are then translated into other viral proteins. In this work, we performed a comparative analysis for novel long-range RNA-RNA interactions that may involve the Spike region. Comparing in-silico fragment-based predictions between reference sequences of SARS-CoV-1 and SARS-CoV-2 revealed several predictions amongst which a thermodynamically stable long-range RNA-RNA interaction between (23660-23703 Spike) and (28025-28060 ORF8) unique to SARS-CoV-2 was observed. The patterns of sequence variation using data gathered worldwide further supported the predicted stability of the sub-interacting region (23679-23690 Spike) and (28031-28042 ORF8). Such RNA-RNA interactions can potentially impact viral life cycle including sub-genomic RNA production rates.

Conservation genomics of urban populations of Streamside Salamander (Ambystoma barbouri).

In Tennessee, populations of the state endangered Streamside Salamander (Ambystoma barbouri) are in decline as their distribution lies mostly within rapidly developing areas in the Nashville Basin. Information regarding the partitioning of genetic variation among populations of A. barbouri and the taxonomic status of these populations relative to northern populations and their congener, the Small-mouthed Salamander (A. texanum), have important implications for management and conservation of this species. Here we combined mitochondrial sequencing and genome-wide single nucleotide polymorphism (SNP) data generated using Genotyping-by-Sequencing (GBS) to investigate patterns of genetic variation within Tennessee populations of A. barbouri, to assess their relationship to populations in Kentucky, and to examine their phylogenetic relationship to the closely related A. texanum. Results from phylogenetic reconstructions reveal a complex history of Tennessee A. barbouri populations with regards to northern populations, unisexual A. barbouri, and A. texanum. Patterns of mitochondrial sequence variation suggest that A. barbouri may have originated within Tennessee and expanded north multiple times into Kentucky, Ohio, Indiana, and West Virginia. Phylogenetic reconstructions based on genome-wide SNP data contradict results based on mitochondrial DNA and correspond to geographic and taxonomic boundaries. Variation in allele frequencies at SNP genotypes, as identified by multivariate analyses and Bayesian assignment tests, identified three evolutionary significant units (ESUs) for A. barbouri within Tennessee. Collectively, these results emphasize the need for prioritizing conservation needs for Tennessee populations of A. barbouri to ensure the long-term persistence of this species.

Stepped Geomorphology Shaped the Phylogeographic Structure of a Widespread Tree Species (Toxicodendron vernicifluum, Anacardiaceae) in East Asia

Species’ phylogeographic patterns reflect the interplay between landscape features, climatic forces, and evolutionary processes. Here, we used two chloroplast DNA (cpDNA) markers (trnL and trnL-F) to explore the role of stepped geomorphology in shaping the phylogeographic structure of Toxicodendron vernicifluum, an economically important tree species widely distributed in East Asia. The range-wide pattern of sequence variation was analyzed based on a dataset including 357 individuals from China, together with published sequences of 92 individuals mainly from Japan and South Korea. We identified five chloroplast haplotypes based on seven substitutions across the 717-bp alignment. A clear east-west phylogeographic break was recovered according to the stepped landforms of mainland China. The wild trees of the western clade were found to be geographically restricted to the “middle step”, which is characterized by high mountains and plateaus, while those of the eastern clade were confined to the “low step”, which is mainly made up of hills and plains. The two major clades were estimated to have diverged during the Early Pleistocene, suggesting that the cool glacial climate may have caused the ancestral population to retreat to at least two glacial refugia, leading to allopatric divergence in response to long-term geographic isolation. Migration vector analyses based on the outputs of ecological niche models (ENMs) supported a gradual range expansion since the Last Interglacial. Mountain ranges in western China and the East China Sea land bridge were inferred to be dispersal corridors in the western and eastern distributions of T. vernicifluum, respectively. Overall, our study provides solid evidence for the role of stepped geomorphology in shaping the phylogeographic patterns of T. vernicifluum. The resulting east-west genetic discontinuities could persist for a long time, and could occur at a much larger scale than previously reported, extending from subtropical (e.g., the Xuefeng Mountain) to warm-temperate China (e.g., the Taihang Mountain).

Development of anthocyanin markers: gene mapping, genomic analysis and genetic diversity studies in Ipomoea species.

The anthocyanins are pigments responsible for a wide range of colours in plants, from blue, red and purple, play essential biological roles as well as their genes are evolutionarily conserved. Purple sweet potatoes have anthocyanins as the predominant colour, even though they are present in orange roots masked by carotenoids. Several studies have focused on molecular aspects of anthocyanin genes, mainly in wild Ipomoea species, although the structure and segregation analysis of those genes in sweet potato hexaploid species are still unknown. Based on an "exon-primed intron-crossing" (EPIC) approach, fourteen pairs of primers were designed, on five structural anthocyanin genes as candidates. The strategy exploits the Intron Length Polymorphism (ILP) from Candidate Genes (CG), resulting in 93% of successful markers giving scorable and reproducible alleles. The results allowed to define partial structure and sequence of the introns and exons from the selected CG, and to determine patterns of sequence variation. The evaluation of marker dosage and allelic segregations in an Ipomoea batatas (L.) Lam mapping population identified several alleles for linkage analysis. The study validated the utility of ILP-CG markers for genetic diversity and conservation applicability and a successful amplification gradient across wild Ipomoea species validated their transferability.

A new minisatellite VNTR marker, Pscp1, discovered for the identification of opium poppy

Opium poppy, a member of the Papaveraceae family, is an ancient herbaceous plant and well-known medical resource in the pharmaceutical industry. However, opium poppies are grown worldwide for producing illicit drugs, significantly increasing the incidence of narcotic drug abuse. Since the narcotic poppy has not yet been genetically investigated, we characterized a novel variable number tandem repeat (VNTR) marker of forensically important poppy species based on the genetic analysis of 164 samples collected from two locations spanning the Jeolla province and Jeju island of South Korea. Comparing analysis of the chloroplast (cp) genome sequences for four representative species of Papaver (Papaver somniferum, Papaver somniferum subs. setigerum, Papaver orientale, and Papaver rhoeas) revealed a unique region with 1–3 repeats for 16 nucleotide motifs in the genome inverted repeat A (IRA, positions 128,651 to 128,698) region. For 16 nucleotide motifs, 3 repeats were found in P. somniferum, and 2 repeats were found in P. somniferum subs. setigerum. Therefore, 10 known and the 133 unknown, seized Papaver species were compared to determine whether the species could be identified via variations in the repeat units. The sizes of a novel VNTR ranged from 181 to 252 bp between the species. Phylogenetic analysis confirmed that a novel VNTR, which we named Pscp1, could clearly distinguish between the narcotic and non-narcotic types of Papaver species based on the patterns of sequence variation. Interestingly, we found that Pscp1 could also distinguish between P. somniferum and P. somniferum subs. setigerum. The regions of eight non-narcotic species displayed similar patterns and also differences were found due to the nucleotide substitution and deletion events. The structural differences of Pscp1 were observed within the two narcotic species or between the narcotic and non-narcotic species, suggesting that these variations may act as a genetic marker. We, therefore, developed a new Pscp1 PCR-capillary electrophoresis (CE) method that can reliably identify the narcotic type of Papaver species. Taken together, our findings suggest that the newly developed Pscp1 can be used as an identification marker of opium poppy, and establish that the Pscp1 genotyping method by PCR-CE is an effective primary screening tool that can also contribute to species discrimination in the field of forensic diagnosis and applications.

A DNA barcode library for the butterflies of North America

Although the butterflies of North America have received considerable taxonomic attention, overlooked species and instances of hybridization continue to be revealed. The present study assembles a DNA barcode reference library for this fauna to identify groups whose patterns of sequence variation suggest the need for further taxonomic study. Based on 14,626 records from 814 species, DNA barcodes were obtained for 96% of the fauna. The maximum intraspecific distance averaged 1/4 the minimum distance to the nearest neighbor, producing a barcode gap in 76% of the species. Most species (80%) were monophyletic, the others were para- or polyphyletic. Although 15% of currently recognized species shared barcodes, the incidence of such taxa was far higher in regions exposed to Pleistocene glaciations than in those that were ice-free. Nearly 10% of species displayed high intraspecific variation (>2.5%), suggesting the need for further investigation to assess potential cryptic diversity. Aside from aiding the identification of all life stages of North American butterflies, the reference library has provided new perspectives on the incidence of both cryptic and potentially over-split species, setting the stage for future studies that can further explore the evolutionary dynamics of this group.

Biogeography of ITS Variation in the Helenium autumnale (Asteraceae) Species Complex

Abstract—We used ITS sequences as species barcodes to evaluate 127 samples of 12 Helenium species from the eastern and central USA, focusing on a species complex of H. autumnale (76 sequences from 11 U.S. states and two Canadian provinces), the federally “threatened” H. virginicum, now recommended for delisting (30 sequences from three U.S. states), and H. flexuosum (11 sequences from four U.S. states). ITS sequences confirmed most species identifications and supported the presence of the first population of the “threatened” endemic Missouri-Virginia disjunct, Helenium virginicum, in Indiana. Because the Indiana plants grow in a restored wetland, have a cpDNA haplotype previously known only from Missouri, with a morphology similar to Missouri variants, and an herbarium search for additional populations in Indiana found none, it is unclear whether the Indiana population is natural or planted. The presence of a putative sister lineage to H. virginicum thought to exist on the Bruce Peninsula, Ontario, Canada, was not supported after sequencing 36 plants with the sister morphology that grew there along 18 km of beach fens and finding they had H. autumnale sequences. Fine-scale biogeographic patterns of intraspecific sequence variation were found mostly in H. autumnale, with centers of different base site polymorphisms found in northern North America and the Missouri Ozarks. As in a previous study, we found DNA evidence of hybridization between Helenium species in Missouri. We offer hypotheses to explain the biogeography of North American Helenium, focusing on the three species that compose the H. autumnale species complex and suggesting that H. autumnale may be a compilospecies showing incomplete lineage sorting. We encourage exploration of more Helenium species and their conspecific populations in search of fine-scale ITS base site polymorphisms to reveal emerging lineages and resolve the origins and evolutionary implications of these biogeographic patterns.