Levels Of Sequence Divergence Research Articles

Copy number and sequence variation in rDNA of Daphnia pulex from natural populations: insights from whole genome sequencing.

Ribosomal DNA (rDNA) has a vital role in ribosome biogenesis as it contains the genes that encode rRNA separated by intergenic spacers. The rRNA genes occur in hundreds to tens of thousands of copies per haploid genome in eukaryotes and are generally highly conserved with low variation within species. Due to the repetitive nature and large size of rDNA arrays, detecting intraindividual variation can be difficult. In this study, we use whole genome sequences of 169 Daphnia pulex individuals from 10 natural populations to measure the copy number and sequence variation in rDNA. This revealed that variation in rDNA copy number between individuals spans an order of magnitude. We further observed a substantial level of sequence variation within individual genomes. As expected, single nucleotide polymorphisms occurred in regions of lower functional constraint such as the intergenic spacer and expansion segments of the rRNA genes. The presence of strong linkage disequilibrium among variants facilitated identification of haplotypes within each population. Although there was evidence of recombination among haplotypes from different populations, it is insufficient to eliminate linkage disequilibrium within populations. Estimating copy number and haplotype diversity within individuals revealed that the level of intraindividual sequence variation is not strongly correlated with copy number. The observed patterns of variation highlight a complex evolutionary history of rDNA in D. pulex. Future research should explore the functional implications of rDNA copy number and sequence variation on organismal phenotypes.

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

A. 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 non-pathogenic relatives. We hypothesized that A. fumigatus genes exhibit substantial genetic variation in the non-coding 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 non-coding regions showed higher levels of sequence variation compared to 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 non-coding regions and 65 total genes with signatures only in their protein-coding regions. 28 of the 478 non-coding regions and 5 of the 65 protein-coding regions under selection are associated with genes known to modulate A. fumigatus virulence. Non-coding region variation between A. fumigatus strains included single nucleotide polymorphisms and insertions or deletions of at least a few nucleotides. These results show that non-coding 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.

Two new species of the genus Leptomyxa (Goodey, 1915) – Leptomyxa botanica n. sp. and Leptomyxa monrepos n. sp. from the soil habitats of Northwestern Russia

Amoebae of the genus Leptomyxa have variable morphologies and can only be reliably identified using molecular data. However, species distinction based on the 18S rRNA gene sequence is difficult due to the very low level of sequence divergence among morphologically different species. The database for other genes is much smaller, and genomic data on Leptomyxa is almost absent. In this study, we describe two new terrestrial species of the genus Leptomyxa isolated from Northwestern Russia, Leptomyxa botanica n. sp. and Leptomyxa monrepos n. sp. Both species easily adopt an expanded fan-shaped form and have a complex structure of the nucleolar material. Phylogenetic analyses show a derived status of these two species. They form a clade with Leptomyxa valladaresi. Our tree confirms that the 18S rRNA gene sequences of Leptomyxa species are split into two large clades. The morphological synapomorphies of these clades are not obvious. This analysis is complicated by the lack of reliable morphological data on many sequenced strains and probable misidentification of some isolates.

Comparative analysis of Camelina sativa and fungal industrial lipases used for biodiesel production

Aim. To identify the genes of false flax (Camelina sativa) endogenous lipases and to analyze the sequence similarity of their key functional domains with those of commercially available lipases. Methods. A detailed search of the databases was carried out in order to identify the sequences of lipases of various species, as well as their sequences were aligned, conservative sequence motifs were identified, the domain structure of the detected proteins was established, and phylogenetic analysis was carried out. Results. 15 triacylglycerol lipase genes were identified, and corresponding sequences of lipases for commercially available products were identified as well. Their domain structure was analyzed, and the level of sequence divergence of their functional regions was also revealed. Conclusions. The lipases that are most similar in terms of sequence and domain organization to the lipases of false flax have been identified, and therefore could potentially be used for more effective transesterification of oil with a specific fatty acid composition for this species.

Multi-allelic gene editing in an apomictic, tetraploid turf and forage grass (Paspalum notatum Flüggé) using CRISPR/Cas9.

Polyploidy is common among grasses (Poaceae) and poses challenges for conventional breeding. Genome editing technology circumvents crossing and selfing, enabling targeted modifications to multiple gene copies in a single generation while maintaining the heterozygous context of many polyploid genomes. Bahiagrass (Paspalum notatum Flüggé; 2n=4x=40) is an apomictic, tetraploid C4 species that is widely grown in the southeastern United States as forage in beef cattle production and utility turf. The chlorophyll biosynthesis gene magnesium chelatase (MgCh) was selected as a rapid readout target for establishing genome editing in tetraploid bahiagrass. Vectors containing sgRNAs, Cas9 and nptII were delivered to callus cultures by biolistics. Edited plants were characterized through PCR-based assays and DNA sequencing, and mutagenesis frequencies as high as 99% of Illumina reads were observed. Sequencing of wild type (WT) bahiagrass revealed a high level of sequence variation in MgCh likely due to the presence of at least two copies with possibly eight different alleles, including pseudogenes. MgCh mutants exhibited visible chlorophyll depletion with up to 82% reductions in leaf greenness. Two lines displayed progression of editing over time which was linked to somatic editing. Apomictic progeny of a chimeric MgCh editing event were obtained and allowed identification of uniformly edited progeny plants among a range of chlorophyll depletion phenotypes. Sanger sequencing of a highly edited mutant revealed elevated frequency of a WT allele, probably due to frequent homology-directed repair (HDR). To our knowledge these experiments comprise the first report of genome editing applied in perennial, warm-season turf or forage grasses. This technology will accelerate bahiagrass cultivar development.

Testudinimyces gracilis gen. nov, sp. nov. and Astrotestudinimyces divisus gen. nov, sp. nov., two novel, deep-branching anaerobic gut fungal genera from tortoise faeces.

The anaerobic gut fungi (AGF, Neocallimastigomycota) represent a basal zoosporic phylum within the kingdom Fungi. Twenty genera are currently described, all of which were isolated from the digestive tracts of mammalian herbivores. Here, we report on the isolation and characterization of novel AGF taxa from faecal samples of tortoises. Twenty-nine fungal isolates were obtained from seven different tortoise species. Phylogenetic analysis using the D1/D2 region of the LSU rRNA gene, ribosomal internal transcribed spacer 1, and RNA polymerase II large subunit grouped all isolates into two distinct, deep-branching clades (clades T and B), with a high level of sequence divergence to their closest cultured relative (Khoyollomyces ramosus). Average amino acid identity values calculated using predicted peptides from the isolates' transcriptomes ranged between 60.80-66.21 % (clade T), and 61.24-64.83 % (clade B) when compared to all other AGF taxa; values that are significantly below recently recommended thresholds for genus (85%) and family (75%) delineation in the Neocallimastigomycota. Both clades displayed a broader temperature growth range (20-45 °C, optimal 30 °C for clade T, and 30-42 °C, optimal 39 °C for clade B) compared to all other AGF taxa. Microscopic analysis demonstrated that strains from both clades produced filamentous hyphae, polycentric rhizoidal growth patterns, and monoflagellated zoospores. Isolates in clade T were characterized by the production of unbranched, predominantly narrow hyphae, and small zoospores, while isolates in clade B were characterized by the production of multiple sporangiophores and sporangia originating from a single central swelling resulting in large multi-sporangiated structures. Based on the unique phylogenetic positions, AAI values, and phenotypic characteristics, we propose to accommodate these isolates into two novel genera (Testudinimyces and Astrotestudinimyces), and species (T. gracilis and A. divisus) within the order Neocallimastigales. The type species are strains T130AT (T. gracilis) and B1.1T (A. divisus).

Molecular and Morphological Characterization of Tylenchus Zeae n. Sp. (Nematoda: Tylenchida) from Corn (Zea Mays) in South Carolina.

Specimens of a tylenchid nematode were recovered in 2019 from soil samples collected from a corn field, located in Pickens County, South Carolina, USA. A moderate number of Tylenchus sp. adults (females and males) were recovered. Extracted nematodes were examined morphologically and molecularly for species identification, which indicated that the specimens of the tylenchid adults were a new species, described herein as Tylenchus zeae n. sp. Morphological examination and the morphometric details of the specimens were very close to the original descriptions of Tylenchus sherianus and T. rex. However, females of the new species can be differentiated from these species by body shape and length, shape of excretory duct, distance between anterior end and esophageal intestinal valve, and a few other characteristics given in the diagnosis. Males of the new species can be differentiated from the two closely related species by tail, spicules, and gubernaculum length. Cryo-scanning electron microscopy confirmed head bearing five or six annules; four to six cephalic sensilla represented by small pits at the rounded corners of the labial plate; a small, round oral plate; and a large, pit-like amphidial opening confined to the labial plate and extending three to four annules beyond it. Phylogenetic analysis of 18S rRNA gene sequences placed Tylenchus zeae n. sp. in a clade with Tylenchus arcuatus and several Filenchus spp., and the mitochondrial cytochrome oxidase c subunit 1 (COI) gene region separated the new species from T. arcuatus and other tylenchid species. In the 28S tree, T. zeae n. sp. showed a high level of sequence divergence and was positioned outside of the main Tylenchus-Filenchus clade.

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.

Identifying host-specific amino acid signatures for influenza A viruses using an adjusted entropy measure

BackgroundInfluenza A viruses (IAV) exhibit vast genetic mutability and have great zoonotic potential to infect avian and mammalian hosts and are known to be responsible for a number of pandemics. A key computational issue in influenza prevention and control is the identification of molecular signatures with cross-species transmission potential. We propose an adjusted entropy-based host-specific signature identification method that uses a similarity coefficient to incorporate the amino acid substitution information and improve the identification performance. Mutations in the polymerase genes (e.g., PB2) are known to play a major role in avian influenza virus adaptation to mammalian hosts. We thus focus on the analysis of PB2 protein sequences and identify host specific PB2 amino acid signatures.ResultsValidation with a set of H5N1 PB2 sequences from 1996 to 2006 results in adjusted entropy having a 40% false negative discovery rate compared to a 60% false negative rate using unadjusted entropy. Simulations across different levels of sequence divergence show a false negative rate of no higher than 10% while unadjusted entropy ranged from 9 to 100%. In addition, under all levels of divergence adjusted entropy never had a false positive rate higher than 9%. Adjusted entropy also identifies important mutations in H1N1pdm PB2 previously identified in the literature that explain changes in divergence between 2008 and 2009 which unadjusted entropy could not identify.ConclusionsBased on these results, adjusted entropy provides a reliable and widely applicable host signature identification approach useful for IAV monitoring and vaccine development.

Medicago truncatula (model legume), Medicago sativa (alfalfa), Medicago polymorpha (bur clover), and Medicago ruthenica

Alfalfa (Medicago sativa) is one of the most widely grown forage crops, providing a major protein source for livestock. In the tetraploid alfalfa genome, the four allelic chromosomes have a low level of sequence divergence, with similar sizes, numbers of genes, and repeat element contents, proving that autotetraploid alfalfa is a stable, random-pairing species.

Complete mitochondrial genomes of endemic charr of the genus Salvelinus from Lake Nachikinskoe (Kamchatka, Russia)

The complete mitochondrial genomes were sequenced from three individuals of the endemic charr Salvelinus sp. 4 [ from Lake Nachikinskoe (Kamchatka, Russia). The complete mitochondrial genomes were 16,654 bp in size; the genome organization and GC content (45.6%) was consistent with charr mitochondrial genomes published previously. The low level of sequence divergence detected between the Salvelinus sp. 4 and the GenBank Salvelinus taranetzi genomes indicated recent divergence and their origin from a common ancestor. Our results could play an essential role in resolving the conflict over current taxonomic status of endemic charr in genus Salvelinus.

Implications of a highly divergent dengue virus strain for cross-neutralization, protection, and vaccine immunity

Although divergent dengue viruses (DENVs) have been isolated in insects, nonhuman primates, and humans, their relationships to the four canonical serotypes (DENV 1-4) are poorly understood. One virus isolated from a dengue patient, DKE-121, falls between genotype and serotype levels of sequence divergence to DENV-4. To examine its antigenic relationship to DENV-4, we assessed serum neutralizing and protective activity. Whereas DENV-4-immune mouse sera neutralize DKE-121 infection, DKE-121-immune sera inhibit DENV-4 less efficiently. Passive transfer of DENV-4 or DKE-121-immune sera protects mice against homologous, but not heterologous, DENV-4 or DKE-121 challenge. Antigenic cartography suggests that DENV-4 and DKE-121 are related but antigenically distinct. However, DENV-4 vaccination confers protection against DKE-121 in nonhuman primates, and serum from humans immunized with a tetravalent vaccine neutralize DENV-4 and DKE-121 infection equivalently. As divergent DENV strains, such as DKE-121, may meet criteria for serotype distinction, monitoring their capacity to impact dengue disease and vaccine efficacy appears warranted.

Integration of transcriptomic and proteomic approaches for snake venom profiling

ABSTRACT Introduction Snake venoms contain many protein and peptide isoforms with high levels of sequence variation, even within a single species. Areas covered In this review, we highlight several examples, from both published and unpublished work in our lab, demonstrating how a combined venom gland transcriptome and proteome methodology allows for comprehensive characterization of venoms, including those from understudied rear-fanged snake species, and we provide recommendations for using these approaches. Expert Opinion When characterizing venoms, peptide mass fingerprinting using databases built predominately from protein sequences originating from model organisms can be disadvantageous, especially when the intention is to document protein diversity. Therefore, the use of species-specific venom gland transcriptomes corrects for the absence of these unique peptide sequences in databases. The integration of transcriptomics and proteomics improves the accuracy of either approach alone for venom profiling.

The complete chloroplast genome sequence of Korean Neolitsea sericea (Lauraceae)

The complete chloroplast (cp) genome sequence of Neolitsea sericea was determined by Illumina sequencing. The complete cp genome was 152,446bp in length, containing a large single-copy region of 93,796 bp and a small single-copy region of 18,506bp, which were separated by a pair of 20,072bp inverted repeats. A total of 112 unique genes were annotated, including 78 protein-coding genes (PCGs), 30 transfer RNAs, and four ribosomal RNAs. Among the PCGs, 18 genes contained one or two introns. A very low level of sequence variation between two cp genomes of N. sericea was found with seven insertions or deletions and only one single nucleotide polymorphism. An analysis using the maximum likelihood method showed that N. sericea was closely related to Actinodaphne trichocarpa.

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes.

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

Androgenic-Induced Transposable Elements Dependent Sequence Variation in Barley.

A plant genome usually encompasses different families of transposable elements (TEs) that may constitute up to 85% of nuclear DNA. Under stressful conditions, some of them may activate, leading to sequence variation. In vitro plant regeneration may induce either phenotypic or genetic and epigenetic changes. While DNA methylation alternations might be related, i.e., to the Yang cycle problems, DNA pattern changes, especially DNA demethylation, may activate TEs that could result in point mutations in DNA sequence changes. Thus, TEs have the highest input into sequence variation (SV). A set of barley regenerants were derived via in vitro anther culture. High Performance Liquid Chromatography (RP-HPLC), used to study the global DNA methylation of donor plants and their regenerants, showed that the level of DNA methylation increased in regenerants by 1.45% compared to the donors. The Methyl-Sensitive Transposon Display (MSTD) based on methylation-sensitive Amplified Fragment Length Polymorphism (metAFLP) approach demonstrated that, depending on the selected elements belonging to the TEs family analyzed, varying levels of sequence variation were evaluated. DNA sequence contexts may have a different impact on SV generated by distinct mobile elements belonged to various TE families. Based on the presented study, some of the selected mobile elements contribute differently to TE-related SV. The surrounding context of the TEs DNA sequence is possibly important here, and the study explained some part of SV related to those contexts.

Evolution and function of nitrogen fixation gene clusters in sugarcane associated Bradyrhizobium strains.

Bradyrhizobium spp. are well known to mediate biological nitrogen fixation (BNF) as microsymbionts inhabiting nodules on leguminous plants. However, they may also contribute to plant growth via free-living N2 fixation (FLNF) in association with non-legumes. Notably, several Bradyrhizobium strains from sugarcane roots display FLNF activity. Among them, Bradyrhizobium sacchari is a legume symbiotic species, whereas strains AG48 and M12 are non-symbiotic. In the present study, a phylogenomic approach was applied to study peculiarities of these and other Bradyrhizobium strains with respect to N fixation (nif) gene content in order to reveal genetic features that enable FNLF in Bradyrhizobium spp. All FLNF strains carry an ancestral 'non-symbiotic' nif-gene cluster (NSC). B. sacchari also contains a second 'symbiotic' nif-gene cluster (SC), a characteristic observed in only three of 156 evaluated genomes. B. sacchari stood out and presented a high level of sequence divergence between individual nif-gene homologues and we discuss scenarios for the evolutionary origin of these clusters. The transcript level of NSC nifH gene increased during FLNF, when compared to symbiotic conditions. The data suggest that sugarcane roots harbor diverse Bradyrhizobium spp. that are genetically adapted to a dynamic environment where leguminous and non-leguminous host plants are alternately available.

Taxonomic assessment of the genus Procamallanus (Nematoda) in Middle American cichlids (Osteichthyes) with molecular data, and the description of a new species from Nicaragua and Costa Rica.

Procamallanus is a species-rich genus of parasitic nematodes of marine, brackish, and freshwater fishes, occurring also occasionally in amphibians and reptiles. In the Neotropical region, this genus is highly diverse, with species described from a wide range of fish families. In this study, we reassess the taxonomic status of Procamallanus rebecae with molecular and morphological data and describe a new species endemic to Nicaragua and Costa Rica. We analyzed all Procamallanus isolated from fish from the Nicaraguan lakes and some rivers in Costa Rica after an exhaustive analysis of their freshwater fish endoparasite fauna. Procamallanus rebecae is a host-specific parasite of Middle American cichlids, previously reported in southern Mexico, Nicaragua, and Costa Rica. We therefore compared these Central American specimens with individuals of P. rebecae collected in cichlids from southeastern Mexico using two genomic regions (28S rDNA and mitochondrial cytochrome oxidase subunit 1, COI). We found high levels of sequence divergence between Procamallanus from the two geographical regions, with up to 9.8 and 10.5% for both genetic markers, respectively. We also analyzed their morphology and found conspicuous differences in the shape of the mouth and the structure of the female cauda. We therefore describe the specimens of Procamallanus from Central American cichlids as a new species. Both Procamallanus species occur in different cichlid species and are allopatrically distributed. The host specificity and ancient association patterns between cichlids and Procamallanus and the jointly colonization of both hosts and parasites during their northern dispersal from South America are briefly discussed.

A Survey of Transposon Landscapes in the Putative Ancient Asexual Ostracod Darwinula stevensoni

How asexual reproduction shapes transposable element (TE) content and diversity in eukaryotic genomes remains debated. We performed an initial survey of TE load and diversity in the putative ancient asexual ostracod Darwinula stevensoni. We examined long contiguous stretches of DNA in clones from a genomic fosmid library, totaling about 2.5 Mb, and supplemented these data with results on TE abundance and diversity from an Illumina draft genome. In contrast to other TE studies in putatively ancient asexuals, which revealed relatively low TE content, we found that at least 19% of the fosmid dataset and 26% of the genome assembly corresponded to known transposons. We observed a high diversity of transposon families, including LINE, gypsy, PLE, mariner/Tc, hAT, CMC, Sola2, Ginger, Merlin, Harbinger, MITEs and helitrons, with the prevalence of DNA transposons. The predominantly low levels of sequence diversity indicate that many TEs are or have recently been active. In the fosmid data, no correlation was found between telomeric repeats and non-LTR retrotransposons, which are present near telomeres in other taxa. Most TEs in the fosmid data were located outside of introns and almost none were found in exons. We also report an N-terminal Myb/SANT-like DNA-binding domain in site-specific R4/Dong non-LTR retrotransposons. Although initial results on transposable loads need to be verified with high quality draft genomes, this study provides important first insights into TE dynamics in putative ancient asexual ostracods.

Genetic divergence of the Sunda ratsnake (Coelognathus subradiatus) across the Lesser Sunda Islands (Squamata: Colubridae)

Abstract The Lesser Sunda Archipelago consists of hundreds of oceanic islands located in southern Wallacea. The Sunda ratsnake, Coelognathus subradiatus, is endemic to the Lesser Sundas and is found on most of the major islands. Mitochondrial DNA was sequenced from snakes representing five of the major islands revealing that levels of sequence divergence between islands range from 2-7%. Phylogenetic analyses recover what can be interpreted as a three-lineage polytomy consisting of lineages from 1) Alor, 2) Sumbawa + Flores, and 3) Timor + Wetar. The archipelago was colonized from the Sunda Shelf ∼7 Myr with subsequent population divergence occurring ∼4.5 Myr, likely resulting in insular species formation.