Search

Harmful algal blooms (HABs) of the toxic haptophytePrymnesium parvumare a recurrent problem in many inland and estuarine waters around the world. Strains ofP. parvumvary in the toxins they produce and in other physiological traits associated with HABs, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for fifteen phylogenetically and geographically diverse strains ofP. parvumincluding Hi-C guided, near-chromosome level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 Mbp to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across strains identified several functional categories associated with metabolism, including candidates for the biosynthesis of toxic metabolites, as well as genome size variation, including recent proliferations of transposable elements. Together, our results indicate thatP. parvumis comprised of multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present inP. parvumand demonstrate the need for similar resources for other HAB-forming morphospecies.SIGNIFICANCE STATEMENTHarmful algal blooms (HABs) are a global concern. Efforts to understand the genetic basis of traits associated with the success of HAB-forming species are limited by a dearth of genomic resources. In this paper we present genomes for fifteen strains ofPrymnesium parvum, a toxic alga that causes ecosystem and societally disruptive HABs around the world. We uncover an unprecedented amount of sequence-level, gene family, and genome architecture evolution inP. parvumand provide evidence for both cryptic speciation and hybridization. These results illustrate how both inter- and intraspecific genetic variation can be dramatically underestimated in a protist morphospecies. More work is needed to understand the eco-physiological consequences of hidden genetic diversity inP. parvumand HAB-forming species more generally.

Significance Advancements in long-read DNA sequencing technologies provide more comprehensive views of genomes. We used long-read sequences to assemble a Great Dane dog genome that provides several improvements over the existing reference derived from a Boxer. Assembly comparisons revealed that gaps in the Boxer assembly often occur at the beginning of protein-coding genes and have a high-GC content, which likely reflects limitations of previous technologies in resolving GC-rich sequences. Dimorphic LINE-1 and SINEC retrotransposons represent the predominant differences between the Great Dane and Boxer assemblies. Proof-of-principle experiments demonstrated that expression of a canine LINE-1 could promote the retrotransposition of itself and a SINEC_Cf consensus sequence in cultured human cells. Thus, ongoing retrotransposon activity is a major contributor to canine genetic diversity., Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3′ end of LINE-1_Cfs (i.e., LINE-1_Cf 3′-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

Lithospermum erythrorhizon (red gromwell; zicao) is a medicinal and economically valuable plant belonging to the Boraginaceae family. Roots from L. erythrorhizon have been used for centuries based on the antiviral and wound-healing properties produced from the bioactive compound shikonin and its derivatives. More recently, shikonin, its enantiomer alkannin, and several other shikonin/alkannin derivatives have collectively emerged as valuable natural colorants and as novel drug scaffolds. Despite several transcriptomes and proteomes having been generated from L. erythrorhizon, a reference genome is still unavailable. This has limited investigations into elucidating the shikonin/alkannin pathway and understanding its evolutionary and ecological significance. In this study, we obtained a de novo genome assembly for L. erythrorhizon using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies. The resulting genome is ∼367.41 Mb long, with a contig N50 size of 314.31 kb and 27,720 predicted protein-coding genes. Using the L. erythrorhizon genome, we identified several additional p-hydroxybenzoate:geranyltransferase (PGT) homologs and provide insight into their evolutionary history. Phylogenetic analysis of prenyltransferases suggests that PGTs originated in a common ancestor of modern shikonin/alkannin-producing Boraginaceous species, likely from a retrotransposition-derived duplication event of an ancestral prenyltransferase gene. Furthermore, knocking down expression of LePGT1 in L. erythrorhizon hairy root lines revealed that LePGT1 is predominantly responsible for shikonin production early in culture establishment. Taken together, the reference genome reported in this study and the provided analysis on the evolutionary origin of shikonin/alkannin biosynthesis will guide elucidation of the remainder of the pathway.

Cdc14 protein phosphatases play an important role in plant infection by several fungal pathogens. This and other properties of Cdc14 enzymes make them an intriguing target for development of new antifungal crop treatments. Active site architecture and substrate specificity of Cdc14 from the model fungus Saccharomyces cerevisiae (ScCdc14) are well-defined and unique among characterized phosphatases. Cdc14 appears absent from some model plants. However, the extent of conservation of Cdc14 sequence, structure, and specificity in fungal plant pathogens is unknown. We addressed this by performing a comprehensive phylogenetic analysis of the Cdc14 family and comparing the conservation of active site structure and specificity among a sampling of plant pathogen Cdc14 homologs. We show that Cdc14 was lost in the common ancestor of angiosperm plants but is ubiquitous in ascomycete and basidiomycete fungi. The unique substrate specificity of ScCdc14 was invariant in homologs from eight diverse species of dikarya, suggesting it is conserved across the lineage. A synthetic substrate mimetic inhibited diverse fungal Cdc14 homologs with similar low µM Ki values, but had little effect on related phosphatases. Our results justify future exploration of Cdc14 as a broad spectrum antifungal target for plant protection.

Mammalian genomes contain a fossilized record of ancient retroviral infections in the form of endogenous retroviruses (ERVs). We used whole genome sequence data to assess the origin and evolution of the recently active ERV-Fc gammaretroviral lineage based on the record of past infections retained in the genome of the domestic dog,Canis lupus familiaris.We identified 165 loci, including 58 insertions absent from the dog reference assembly, and characterized element polymorphism across 332 canids from nine species. Insertions were found throughout the dog genome including within and near gene models. Analysis of 19 proviral sequences identified shared disruptive mutations indicating defective proviruses were spread via complementation. The patterns of ERV polymorphism and sequence variation indicate multiple circulating viruses infected canid ancestors within the last 20 million to within 1.6 million years with a recent bust of germline invasion in the lineage leading to wolves and dogs.

BackgroundDogs (Canis lupus familiaris) were domesticated from gray wolves between 10-40 kya in Eurasia, yet details surrounding the process of domestication remain unclear. The vast array of phenotypes exhibited by dogs mirror other domesticated animal species, a phenomenon known as the Domestication Syndrome. Here, we use signatures persisting in the dog genome to identify genes and pathways altered by the intensive selective pressures of domestication.ResultsWe identified 246 candidate domestication regions containing 10.8Mb of genome sequence and 178 genes through whole-genome SNP analysis of 43 globally distributed village dogs and 10 wolves. Comparisons with ancient dog genomes suggest that these regions reflect signatures of domestication rather than breed formation. The strongest hit is located in theRetinoic Acid-Induced 1(RAI1) gene, mutations of which cause Smith-Magenis syndrome. The identified regions contain a significant enrichment of genes linked to neural crest cell migration, differentiation and development. Read depth analysis suggests that copy number variation played a minor role in dog domestication.ConclusionOur results indicate that phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development, are determined by genes which act early in embryogenesis. These differences are all phenotypes of the Domestication Syndrome, which can be explained by decreases in neural crest cells at these sites. We propose that initial selection during early dog domestication was for behavior, a trait also influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.