Human-mouse alignments with BLASTZ

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

Appelmans protocol - A directed in vitro evolution enables induction and recombination of prophages with expanded host range

Infections caused by carbapenem-resistant Acinetobacter baumannii (CRAB) present significant healthcare challenges due to limited treatment options. Bacteriophage (phage) therapy offers potential as an alternative treatment. However, the high host specificity of phages poses challenges for their therapeutic application. To broaden the phage spectrum, laboratory-based phage training using the Appelmans protocol was employed in this study. As a result, the protocol successfully expanded the host range of a phage cocktail targeting CRAB. Further analysis revealed that the expanded host range phages isolated from the output cocktail were identified as recombinant derivatives originating from prophages induced from encountered bacterial strains. These findings provide valuable genetic insights into the protocol's mechanism when applied to phages infecting A. baumannii strains that have never been investigated before. However, it is noteworthy that the expanded host range phages obtained from this protocol exhibited limited stability, raising concerns about their suitability for therapeutic purposes.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

Evolution of the Araliaceae family involved rapid diversification of the Asian Palmate group and Hydrocotyle specific mutational pressure

The Araliaceae contain many valuable species in medicinal and industrial aspects. We performed intensive phylogenomics using the plastid genome (plastome) and 45S nuclear ribosomal DNA sequences. A total of 66 plastome sequences were used, 13 of which were newly assembled in this study, 12 from new sequences, and one from existing data. While Araliaceae plastomes showed conserved genome structure, phylogenetic reconstructions based on four different plastome datasets revealed phylogenetic discordance within the Asian Palmate group. The divergence time estimation revealed that splits in two Araliaceae subfamilies and the clades exhibiting phylogenetic discordances in the Asian Palmate group occurred at two climatic optima, suggesting that global warming events triggered species divergence, particularly the rapid diversification of the Asian Palmate group during the Middle Miocene. Nucleotide substitution analyses indicated that the Hydrocotyloideae plastomes have undergone accelerated AT-biased mutations (C-to-T transitions) compared with the Aralioideae plastomes, and the acceleration may occur in their mitochondrial and nuclear genomes as well. This implies that members of the genus Hydrocotyle, the only aquatic plants in the Araliaceae, have experienced a distinct evolutionary history from the other species. We also discussed the intercontinental disjunction in the genus Panax and proposed a hypothesis to complement the previously proposed hypothesis. Our results provide the evolutionary trajectory of Araliaceae and advance our current understanding of the evolution of Araliaceae species.

Enterococcus faecium C171: Modulating the Immune Response to Acute Lethal Viral Challenge

Commensal bacteria modulate acute immune responses to infection in hosts. In this study, Enterococcus faecium C171 was screened and isolated. This strain has similar basic characteristics to the reference probiotic, including strong anti-inflammatory and anti-infective effects. E. faecium C171 inhibits the production of pro-Caspase-1 and significantly reduces the production of interleukin-1β (IL-1β) in vitro. These reactions were confirmed using the Transwell system. Live E. faecium C171 mainly exerted an inhibitory effect on acute inflammation, whereas the anti-infective and immune-activating effects were primarily mediated by the E. faecium C171-produced bacterial extracellular vesicles (Efm-C171-BEVs). Furthermore, in the specific pathogen-free (SPF) chicken model, oral administration of E. faecium C171 increased the relative abundance of beneficial microbiota (Enterococcus and Lactobacillus), particularly Enterococcus, the most important functional bacteria of the gut microbiota. E. faecium C171 significantly inhibited the acute inflammatory response induced by a highly virulent infectious disease, and reduced mortality in SPF chickens by 75%. In addition, E. faecium C171 induced high levels of CD3+, CD4-, and CD8- immunoregulatory cells and CD8+ killer T cells, and significantly improved the proliferative activity of T cells in peripheral blood mononuclear cells, and the secretion of interferon-γ. These findings indicate that E. faecium C171 and Efm-C171-BEVs are promising candidates for adjuvant treatment of acute inflammatory diseases and acute viral infections.

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Learning functional conservation between human and pig to decipher evolutionary mechanisms underlying gene expression and complex traits

Assessment of genomic conservation between humans and pigs at the functional level can improve the potential of pigs as a human biomedical model. To address this, we developed a deep learning-based approach to learn the genomic conservation at the functional level (DeepGCF) between species by integrating 386 and 374 functional profiles from humans and pigs, respectively. DeepGCF demonstrated better prediction performance compared with the previous method. In addition, the resulting DeepGCF score captures the functional conservation between humans and pigs by examining chromatin states, sequence ontologies, and regulatory variants. We identified a core set of genomic regions as functionally conserved that plays key roles in gene regulation and is enriched for the heritability of complex traits and diseases in humans. Our results highlight the importance of cross-species functional comparison in illustrating the genetic and evolutionary basis of complex phenotypes.

A High-Continuity Genome Assembly of Chinese Flowering Cabbage (Brassica rapa var. parachinensis) Provides New Insights into Brassica Genome Structure Evolution

Chinese flowering cabbage (Brassica rapa var. parachinensis) is a popular and widely cultivated leaf vegetable crop in Asia. Here, we performed a high quality de novo assembly of the 384 Mb genome of 10 chromosomes of a typical cultivar of Chinese flowering cabbage with an integrated approach using PacBio, Illumina, and Hi-C technology. We modeled 47,598 protein-coding genes in this analysis and annotated 52% (205.9/384) of its genome as repetitive sequences including 17% in DNA transposons and 22% in long terminal retrotransposons (LTRs). Phylogenetic analysis reveals the genome of the Chinese flowering cabbage has a closer evolutionary relationship with the AA diploid progenitor of the allotetraploid species, Brassica juncea. Comparative genomic analysis of Brassica species with different subgenome types (A, B and C) reveals that the pericentromeric regions on chromosome 5 and 6 of the AA genome have been significantly expanded compared to the orthologous genomic regions in the BB and CC genomes, largely driven by LTR-retrotransposon amplification. Furthermore, we identified a large number of structural variations (SVs) within the B. rapa lines that could impact coding genes, suggesting the functional significance of SVs on Brassica genome evolution. Overall, our high-quality genome assembly of the Chinese flowering cabbage provides a valuable genetic resource for deciphering the genome evolution of Brassica species and it can potentially serve as the reference genome guiding the molecular breeding practice of B. rapa crops.

Diversification and introgression in four chromosomal taxa of the Pearson's horseshoe bat (Rhinolophus pearsoni) group

Chromosomal variation among closely related taxa is common in both plants and animals, and can reduce rates of introgression as well as promote reproductive isolation and speciation. In mammals, studies relating introgression to chromosomal variation have tended to focus on a few model systems and typically characterized levels of introgression using small numbers of loci. Here we took a genome-wide approach to examine how introgression rates vary among four closely related horseshoe bats (Rhinolophus pearsoni group) that possess different diploid chromosome numbers (2n = 42, 44, 46, and 60) resulting from Robertsonian (Rb) changes (fissions/fusions). Using a sequence capture we obtained orthologous loci for thousands of nuclear loci, as well as mitogenomes, and performed phylogenetic and population genetic analyses. We found that the taxon with 2n = 60 was the first to diverge in this group, and that the relationships among the three other taxa (2n = 42, 44 and 46) showed discordance across our different analyses. Our results revealed signatures of multiple ancient introgression events between the four taxa, with evidence of mitonuclar discordance in phylogenetic trees and reticulation events in their evolutionary history. Despite this, we found no evidence of recent and/or ongoing introgression between taxa. Overall, our results indicate that the effects of Rb changes on the reduction of introgression are complicated and that these may contribute to reproductive isolation and speciation in concert with other factors (e.g. phenotypic and genic divergence).

The First Crested Duck Genome Reveals Clues to Genetic Compensation and Crest Cushion Formation

The Chinese crested (CC) duck is a unique indigenous waterfowl breed, which has a crest cushion that affects its survival rate. Therefore, the CC duck is an ideal model to investigate the genetic compensation response to maintain genetic stability. In the present study, we first generated a chromosome-level genome of CC ducks. Comparative genomics revealed that genes related to tissue repair, immune function, and tumors were under strong positive selection, indicating that these adaptive changes might enhance cancer resistance and immune response to maintain the genetic stability of CC ducks. We also assembled a Chinese spot-billed (Csp-b) duck genome, and detected the structural variations (SVs) in the genome assemblies of three ducks (i.e., CC duck, Csp-b duck, and Peking duck). Functional analysis revealed that several SVs were related to the immune system of CC ducks, further strongly suggesting that genetic compensation in the anti-tumor and immune systems supports the survival of CC ducks. Moreover, we confirmed that the CC duck originated from the mallard ducks. Finally, we revealed the physiological and genetic basis of crest traits and identified a causative mutation in TAS2R40 that leads to crest formation. Overall, the findings of this study provide new insights into the role of genetic compensation in adaptive evolution.

Investigation of monacolin K, yellow pigments, and citrinin production capabilities of Monascus purpureus and Monascus ruber (Monascus pilosus)

Red mold rice (RMR) is a traditional Chinese medicine prepared using Monascus fermentation. Monascus ruber ( pilosus) and Monascus purpureus have a long history of use as food and medicine. As an economically important starter culture, the relationship between the taxonomy of Monascus and production capabilities of secondary metabolites is crucial for the Monascus food industry. In this study, monacolin K, monascin, ankaflavin, and citrinin production by M. purpureus and M. ruber were genomically and chemically investigated. Our findings suggest that M. purpureus can produce monascin and ankaflavin in a correlated manner, whereas M. ruber produces monascin with minimum ankaflavin. M. purpureus is capable of producing citrinin; however, it is unlikely able to produce monacolin K. In contrast, M. ruber produces monacolin K, but not citrinin. We suggest that the current monacolin K content-related regulation of Monascus food should be revised, and labeling of Monascus species should be considered.

Characterization of Virulent T4-Like Acinetobacter baumannii Bacteriophages DLP1 and DLP2

The world is currently facing a global health crisis due to the rapid increase in antimicrobial-resistant bacterial infections. One of the most concerning pathogens is Acinetobacter baumannii, which is listed as a Priority 1 pathogen by the World Health Organization. This Gram-negative bacterium has many intrinsic antibiotic resistance mechanisms and the ability to quickly acquire new resistance determinants from its environment. A limited number of effective antibiotics against this pathogen complicates the treatment of A. baumannii infections. A potential treatment option that is rapidly gaining interest is “phage therapy”, or the clinical application of bacteriophages to selectively kill bacteria. The myoviruses DLP1 and DLP2 (vB_AbaM-DLP_1 and vB_AbaM-DLP_2, respectively) were isolated from sewage samples using a capsule minus variant of A. baumannii strain AB5075. Host range analysis of these phages against 107 A. baumannii strains shows a limited host range, infecting 15 and 21 for phages DLP1 and DLP2, respectively. Phage DLP1 has a large burst size of 239 PFU/cell, a latency period of 20 min, and virulence index of 0.93. In contrast, DLP2 has a smaller burst size of 24 PFU/cell, a latency period of 20 min, and virulence index of 0.86. Both phages show potential for use as therapeutics to combat A. baumannii infections.

BLEND: a fast, memory-efficient and accurate mechanism to find fuzzy seed matches in genome analysis

Generating the hash values of short subsequences, called seeds, enables quickly identifying similarities between genomic sequences by matching seeds with a single lookup of their hash values. However, these hash values can be used only for finding exact-matching seeds as the conventional hashing methods assign distinct hash values for different seeds, including highly similar seeds. Finding only exact-matching seeds causes either (i) increasing the use of the costly sequence alignment or (ii) limited sensitivity. We introduce BLEND, the first efficient and accurate mechanism that can identify both exact-matching and highly similar seeds with a single lookup of their hash values, called fuzzy seed matches. BLEND (i) utilizes a technique called SimHash, that can generate the same hash value for similar sets, and (ii) provides the proper mechanisms for using seeds as sets with the SimHash technique to find fuzzy seed matches efficiently. We show the benefits of BLEND when used in read overlapping and read mapping. For read overlapping, BLEND is faster by 2.4×-83.9× (on average 19.3×), has a lower memory footprint by 0.9×-14.1× (on average 3.8×), and finds higher quality overlaps leading to accurate de novo assemblies than the state-of-the-art tool, minimap2. For read mapping, BLEND is faster by 0.8×-4.1× (on average 1.7×) than minimap2. Source code is available at https://github.com/CMU-SAFARI/BLEND.

The chloroplast GATA-motif of Mahonia bealei participates in alkaloid-mediated photosystem inhibition during dark to light transition

Mahonia bealei and Mahonia fortunei are important plant resources in Traditional Chinese Medicine that are valued for their high levels of benzylisoquinoline alkaloids (BIAs). Although the phytotoxic activity of BIAs has been recognized, information is limited on the mechanism of action by which these compounds regulate photosynthetic activity. Here, we performed comparative chloroplast genome analysis to examine insertions and deletions in the two species. We found a GATA-motif located in the promoter region of the ndhF gene of only M. bealei. K-mer frequency-based diversity analysis illustrated the close correlation between the GATA-motif and leaf phenotype. We found that the GATA-motif significantly inhibits GUS gene expression in tobacco during the dark-light transition (DLT). The expression of ndhF was downregulated in M. bealei and upregulated in M. fortunei during the DLT. NDH-F activity was remarkably decreased and exhibited a significant negative correlation with BIA levels in M. bealei during the DLT. Furthermore, the NADPH produced through photosynthetic metabolism was found to decrease in M. bealei during the DLT. Taken together, our results indicate that this GATA-motif might act as the functional site by which BIAs inhibit photosynthetic metabolism through downregulating ndhF expression during the DLT.

Complete chloroplast genome sequence of Solanum iopetalum, one of the tuber-bearing wild potato relatives

Solanum iopetalum belongs to the Solanaceae family and is one of the tuber-bearing wild Solanum species. In this study, chloroplast genome sequencing of the species, completed with Illumina sequencing technology, is presented. The length of the chloroplast genome is 155,625 bp with a GC content of 37.86%. It comprises a large single copy (LSC) region of 86,057 bp, a small single copy (SSC) region of 18,382 bp, and two inverted repeat regions (IRa and IRb) of 25,593 bp. Additionally, 158 functional genes in the genome are identified, including 105 protein-coding genes, 8 ribosomal RNAs, and 45 transfer RNAs. Phylogenetic analysis revealed that S. iopetalum is grouped into a large clade with other Solanum species, including cultivated potatoes (S. tuberosum) and is closely related to Mexican Solanum species (S. stoloniferum, S. verrucosum, S. hougasii, S. hjertingii, and S. demissum). This study provides useful genomic information for future breeding and evolutionary studies of S. iopetalum and other Solanum species.

A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

TCEAL1 loss-of-function results in an X-linked dominant neurodevelopmental syndrome and drives the neurological disease trait in Xq22.2 deletions

An Xq22.2 region upstream of PLP1 has been proposed to underly a neurological disease trait when deleted in 46,XX females. Deletion mapping revealed that heterozygous deletions encompassing the smallest region of overlap (SRO) spanning six Xq22.2 genes (BEX3, RAB40A, TCEAL4, TCEAL3, TCEAL1, and MORF4L2) associate with an early-onset neurological disease trait (EONDT) consisting of hypotonia, intellectual disability, neurobehavioral abnormalities, and dysmorphic facial features. None of the genes within the SRO have been associated with monogenic disease in OMIM. Through local and international collaborations facilitated by GeneMatcher and Matchmaker Exchange, we have identified and herein report seven de novo variants involving TCEAL1 in seven unrelated families: three hemizygous truncating alleles; one hemizygous missense allele; one heterozygous TCEAL1 full gene deletion; one heterozygous contiguous deletion of TCEAL1, TCEAL3, and TCEAL4; and one heterozygous frameshift variant allele. Variants were identified through exome or genome sequencing with trio analysis or through chromosomal microarray. Comparison with previously reported Xq22 deletions encompassing TCEAL1 identified a more-defined syndrome consisting of hypotonia, abnormal gait, developmental delay/intellectual disability especially affecting expressive language, autistic-like behavior, and mildly dysmorphic facial features. Additional features include strabismus, refractive errors, variable nystagmus, gastroesophageal reflux, constipation, dysmotility, recurrent infections, seizures, and structural brain anomalies. An additional maternally inherited hemizygous missense allele of uncertain significance was identified in a male with hypertonia and spasticity without syndromic features. These data provide evidence that TCEAL1 loss of function causes a neurological rare disease trait involving significant neurological impairment with features overlapping the EONDT phenotype in females with the Xq22 deletion.

The first complete genome of the simian malaria parasite Plasmodium brasilianum

Naturally occurring human infections by zoonotic Plasmodium species have been documented for P. knowlesi, P. cynomolgi, P. simium, P. simiovale, P. inui, P. inui-like, P. coatneyi, and P. brasilianum. Accurate detection of each species is complicated by their morphological similarities with other Plasmodium species. PCR-based assays offer a solution but require prior knowledge of adequate genomic targets that can distinguish the species. While whole genomes have been published for P. knowlesi, P. cynomolgi, P. simium, and P. inui, no complete genome for P. brasilianum has been available. Previously, we reported a draft genome for P. brasilianum, and here we report the completed genome for P. brasilianum. The genome is 31.4 Mb in size and comprises 14 chromosomes, the mitochondrial genome, the apicoplast genome, and 29 unplaced contigs. The chromosomes consist of 98.4% nucleotide sites that are identical to the P. malariae genome, the closest evolutionarily related species hypothesized to be the same species as P. brasilianum, with 41,125 non-synonymous SNPs (0.0722% of genome) identified between the two genomes. Furthermore, P. brasilianum had 4864 (82.1%) genes that share 80% or higher sequence similarity with 4970 (75.5%) P. malariae genes. This was demonstrated by the nearly identical genomic organization and multiple sequence alignments for the merozoite surface proteins msp3 and msp7. We observed a distinction in the repeat lengths of the circumsporozoite protein (CSP) gene sequences between P. brasilianum and P. malariae. Our results demonstrate a 97.3% pairwise identity between the P. brasilianum and the P. malariae genomes. These findings highlight the phylogenetic proximity of these two species, suggesting that P. malariae and P. brasilianum are strains of the same species, but this could not be fully evaluated with only a single genomic sequence for each species.

A DNA Replication Mechanism Can Explain Structural Variation at the Pigeon Recessive Red Locus

For species to adapt to their environment, evolution must act upon genetic variation that is present in the population. Elucidating the molecular mechanisms that give rise to this variation is thus of crucial importance for understanding how organisms evolve. In addition to variation caused by point mutations, structural variation (deletions, duplications, inversions, translocations) is also an important source of variety. Mechanisms involving recombination, transposition and retrotransposition, and replication have been proposed for generating structural variation, and each are capable of explaining certain rearrangements. In this study, we conduct a detailed analysis of two partially overlapping rearrangements (e1 and e2 allele) in domestic rock pigeon (Columba livia) which are both associated with the recessive red phenotype. We find that a replicative mechanism is best able to explain the complex architecture of the e1 allele, and is also compatible with the simpler architecture of the e2 allele as well.

Analysis of the complete plastomes and nuclear ribosomal DNAs from Euonymus hamiltonianus and its relatives sheds light on their diversity and evolution

Euonymus hamiltonianus and its relatives (Celastraceae family) are used for ornamental and medicinal purposes. However, species identification in Euonymus is difficult due to their morphological diversity. Using plastid genome (plastome) data, we attempt to reveal phylogenetic relationship among Euonymus species and develop useful markers for molecular identification. We assembled the plastome and nuclear ribosomal DNA (nrDNA) sequences from five Euonymus lines collected from South Korea: three Euonymus hamiltonianus accessions, E. europaeus, and E. japonicus. We conducted an in-depth comparative analysis using ten plastomes, including other publicly available plastome data for this genus. The genome structures, gene contents, and gene orders were similar in all Euonymus plastomes in this study. Analysis of nucleotide diversity revealed six divergence hotspots in their plastomes. We identified 339 single nucleotide polymorphisms and 293 insertion or deletions among the four E. hamiltonianus plastomes, pointing to abundant diversity even within the same species. Among 77 commonly shared genes, 9 and 33 were identified as conserved genes in the genus Euonymus and E. hamiltonianus, respectively. Phylogenetic analysis based on plastome and nrDNA sequences revealed the overall consensus and relationships between plastomes and nrDNAs. Finally, we developed six barcoding markers and successfully applied them to 31 E. hamiltonianus lines collected from South Korea. Our findings provide the molecular basis for the classification and molecular taxonomic criteria for the genus Euonymus (at least in Korea), which should aid in more objective classification within this genus. Moreover, the newly developed markers will be useful for understanding the species delimitation of E. hamiltonianus and closely related species.

Multiple genome alignment in the telomere-to-telomere assembly era

With the arrival of telomere-to-telomere (T2T) assemblies of the human genome comes the computational challenge of efficiently and accurately constructing multiple genome alignments at an unprecedented scale. By identifying nucleotides across genomes which share a common ancestor, multiple genome alignments commonly serve as the bedrock for comparative genomics studies. In this review, we provide an overview of the algorithmic template that most multiple genome alignment methods follow. We also discuss prospective areas of improvement of multiple genome alignment for keeping up with continuously arriving high-quality T2T assembled genomes and for unlocking clinically-relevant insights.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae)

BACKGROUND

Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species.

AIM

 Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species.

RESULTS

We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample.

CONCLUSIONS

Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

An accurate alignment-free protein sequence comparator based on physicochemical properties of amino acids

Bio-sequence comparators are one of the most basic and significant methods for assessing biological data, and so, due to the importance of proteins, protein sequence comparators are particularly crucial. On the other hand, the complexity of the problem, the growing number of extracted protein sequences, and the growth of studies and data analysis applications addressing protein sequences have necessitated the development of a rapid and accurate approach to account for the complexities in this field. As a result, we propose a protein sequence comparison approach, called PCV, which improves comparison accuracy by producing vectors that encode sequence data as well as physicochemical properties of the amino acids. At the same time, by partitioning the long protein sequences into fix-length blocks and providing encoding vector for each block, this method allows for parallel and fast implementation. To evaluate the performance of PCV, like other alignment-free methods, we used 12 benchmark datasets including classes with homologous sequences which may require a simple preprocessing search tool to select the homologous data. And then, we compared the protein sequence comparison outcomes to those of alternative alignment-based and alignment-free methods, using various evaluation criteria. These results indicate that our method provides significant improvement in sequence classification accuracy, compared to the alternative alignment-free methods and has an average correlation of about 94% with the ClustalW method as our reference method, while considerably reduces the processing time.

Recurrent bacteremia with a hypermucoviscous Escherichia coli isolated from a patient with perihilar cholangiocarcinoma: insights from a comprehensive genome-based analysis

Background

Escherichia coli (E. coli) is a common human pathogen, responsible for a broad spectrum of infections. Sites of infection can vary, but the hepato-biliary system is of particular concern due to the infection-associated formation of gallstones and the spread of pathogens from the bile ducts into the bloodstream.

Case presentation

The presented case is striking, as the detected isolate showed a positive string test. This hypermucoviscous phenotype is atypical for E. coli and a particular feature of hypervirulent Klebsiella pneumoniae (K. pneumoniae) variants.

Objectives

To provide new insights into the genomic background of an E. coli strain with an unusual hypermucoviscous phenotype using hybrid short- and long-read sequencing approaches.

Results

Complete hybrid assemblies of the E. coli genome and plasmids were done and used for genome based typing. Isolate 537–20 was assigned to the multilocus sequence type ST88 and serotype O8:H4. The strain showed a close relationship to avian pathogenic strains. Analysis of the chromosome and plasmids revealed the presence of several virulence factors, such as the Conserved Virulence Plasmidic (CVP) region on plasmid 537-20_1, including several iron acquisition genes (sitABCD, iroABCDEN, iucABCD, hbd) and the iutA gene encoding the receptor of the siderophore aerobactin. The hypermucoviscous phenotype could be caused by encapsulation of putative K. pneumoniae origin.

Conclusions

Hybrid sequencing enabled detailed genomic characterization of the hypermucoviscous E. coli strain, revealing virulence factors that have their putative origin in K. pneumoniae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12941-022-00521-7.

Global, highly specific and fast filtering of alignment seeds

Background

An important initial phase of arguably most homology search and alignment methods such as required for genome alignments is seed finding. The seed finding step is crucial to curb the runtime as potential alignments are restricted to and anchored at the sequence position pairs that constitute the seed. To identify seeds, it is good practice to use sets of spaced seed patterns, a method that locally compares two sequences and requires exact matches at certain positions only.

Results

We introduce a new method for filtering alignment seeds that we call geometric hashing. Geometric hashing achieves a high specificity by combining non-local information from different seeds using a simple hash function that only requires a constant and small amount of additional time per spaced seed. Geometric hashing was tested on the task of finding homologous positions in the coding regions of human and mouse genome sequences. Thereby, the number of false positives was decreased about million-fold over sets of spaced seeds while maintaining a very high sensitivity.

Conclusions

An additional geometric hashing filtering phase could improve the run-time, accuracy or both of programs for various homology-search-and-align tasks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04745-4.

Identification of Two Diamondback Moth Parasitoids, Diadegma fenestrale and Diadegma semiclausum, Using LAMP for Application in Biological Control

The diamondback moth, Plutella xylostella L., is a lepidopteran pest that damages various vegetable plants belonging to the genus Brassica worldwide. Various biological controls, such as parasitoid wasps, have been used to control this pest. Among these, Diadegma semiclausum and Diadegmafenestrale are widely used globally. In field-based biological control research, the investigation of the population dynamics of parasitoids and the rate of parasitism within the pest population is very important. However, achieving profundity in research is difficult when morphologically similar species coexist in the field. The morphological characteristics of D. semiclausum and D. fenestrale are very similar, and they both parasitize P. xylostella larvae. Therefore, to accurately identify these species, in this study, we developed a molecular diagnostic method by using loop-mediated isothermal amplification (LAMP). The mitochondrial genome of D. fenestrale and partial nucleotide sequences, including the ITS region of D. semiclausum, were analyzed for use as species diagnosis markers. The results showed that the homology of D. fenestrale to D. semiclausum was 94%, due to the excessively low homology of the D loop, but the actual homology was higher than 94%, particularly in the coding region. D. fenestrale species-specific primers for LAMP were designed based on the region encoding COX3, and the optimal diagnostic reaction condition for the four primers (F3, B3, FIP, and BIP) was 63 °C for 35 min. A species-specific primer capable of classifying D. semiclausum was developed based on the ITS2 region, and the optimal reaction condition for diagnosis was 63 °C for 40 min. Under optimal conditions for both species, upon addition of the loop primer LB, the reaction efficiency increased, and the reaction time was shortened by more than 5 min. The diagnostic limit concentration was up to 10 pg under both optimal conditions; therefore, it was possible to detect even very low concentrations. For both species, diagnosis was possible by using LAMP assay with a DNA-releasing technique, without a DNA extraction process, and by incubating a tissue sample or the homogenized whole body at 95 °C for 5 min. In the case of D. fenestrale, it was possible to diagnose the parasitoid in P. xylostella larvae. Therefore, the developed LAMP diagnostic method can be used in a variety of ways to determine whether P. xylostella has been parasitized in the process of field research and mass breeding, and to accurately distinguish the species that are parasitic to P. xylostella larvae. This LAMP-based diagnostic method can be applied to identify various parasitoids that are used for the biological control of P. xylostella.

Complete chloroplast genome sequence of Solanum hjertingii, one of the wild potato relatives

Solanum hjertingii is a wild tuber-bearing species classified in the Solanaceae family. The chloroplast genome of S. hjertingii was completed via de novo assembly using Illumina paired-end sequencing data. Total length of the chloroplast genome of S. hjertingii is 155,545 bp consisting of 85,976 bp in a large single copy, 18,383 bp in a small single copy, and 25,593 bp in a pair of inverted repeat regions. Its structure is circular and typically quadripartite. It contains 158 predicted genes in total, including 105 protein-coding, 45 tRNA, and eight rRNA genes. Maximum likelihood phylogenetic analysis of the sequence along with 33 species in the Solanaceae family revealed that S. hjertingii belongs to a large clade with other Solanum species including S. tuberosum and is most closely grouped in the clade with S. hougasii and S. stoloniferum in the clade.

OUP accepted manuscript

Mouse models have been engineered to reveal the biological mechanisms of human diseases based on an assumption. The assumption is that orthologous genes underlie conserved phenotypes across species. However, genetically modified mouse orthologs of human genes do not often recapitulate human disease phenotypes which might be due to the molecular evolution of phenotypic differences across species from the time of the last common ancestor. Here, we systematically investigated the evolutionary divergence of regulatory relationships between transcription factors (TFs) and target genes in functional modules, and found that the rewiring of gene regulatory networks (GRNs) contributes to the phenotypic discrepancies that occur between humans and mice. We confirmed that the rewired regulatory networks of orthologous genes contain a higher proportion of species-specific regulatory elements. Additionally, we verified that the divergence of target gene expression levels, which was triggered by network rewiring, could lead to phenotypic differences. Taken together, a careful consideration of evolutionary divergence in regulatory networks could be a novel strategy to understand the failure or success of mouse models to mimic human diseases. To help interpret mouse phenotypes in human disease studies, we provide quantitative comparisons of gene expression profiles on our website (http://sbi.postech.ac.kr/w/RN).

The long noncoding RNA H19 regulates tumor plasticity in neuroendocrine prostate cancer

Neuroendocrine (NE) prostate cancer (NEPC) is a lethal subtype of castration-resistant prostate cancer (PCa) arising either de novo or from transdifferentiated prostate adenocarcinoma following androgen deprivation therapy (ADT). Extensive computational analysis has identified a high degree of association between the long noncoding RNA (lncRNA) H19 and NEPC, with the longest isoform highly expressed in NEPC. H19 regulates PCa lineage plasticity by driving a bidirectional cell identity of NE phenotype (H19 overexpression) or luminal phenotype (H19 knockdown). It contributes to treatment resistance, with the knockdown of H19 re-sensitizing PCa to ADT. It is also essential for the proliferation and invasion of NEPC. H19 levels are negatively regulated by androgen signaling via androgen receptor (AR). When androgen is absent SOX2 levels increase, driving H19 transcription and facilitating transdifferentiation. H19 facilitates the PRC2 complex in regulating methylation changes at H3K27me3/H3K4me3 histone sites of AR-driven and NEPC-related genes. Additionally, this lncRNA induces alterations in genome-wide DNA methylation on CpG sites, further regulating genes associated with the NEPC phenotype. Our clinical data identify H19 as a candidate diagnostic marker and predictive marker of NEPC with elevated H19 levels associated with an increased probability of biochemical recurrence and metastatic disease in patients receiving ADT. Here we report H19 as an early upstream regulator of cell fate, plasticity, and treatment resistance in NEPC that can reverse/transform cells to a treatable form of PCa once therapeutically deactivated.

Elevated expression of long noncoding RNA H19 is seen in clinical samples of neuroendocrine prostate cancer (PCa). Here the authors show H19 promotes plasticity from luminal to neuroendocrine by epigenetic reprogramming.

A bioinformatic platform to integrate target capture and whole genome sequences of various read depths for phylogenomics

The increasing availability of short‐read whole genome sequencing (WGS) provides unprecedented opportunities to study ecological and evolutionary processes. Although loci of interest can be extracted from WGS data and combined with target sequence data, this requires suitable bioinformatic workflows. Here, we test different assembly and locus extraction strategies and implement them into secapr, a pipeline that processes short‐read data into multilocus alignments for phylogenetics and molecular ecology analyses. We integrate the processing of data from low‐coverage WGS (<30×) and target sequence capture into a flexible framework, while optimizing de novo contig assembly and loci extraction. Specifically, we test different assembly strategies by contrasting their ability to recover loci from targeted butterfly protein‐coding genes, using four data sets: a WGS data set across different average coverages (10×, 5× and 2×) and a data set for which these loci were enriched prior to sequencing via target sequence capture. Using the resulting de novo contigs, we account for potential errors within contigs and infer phylogenetic trees to evaluate the ability of each assembly strategy to recover species relationships. We demonstrate that choosing multiple sizes of kmer simultaneously for assembly results in the highest yield of extracted loci from de novo assembled contigs, while data sets derived from sequencing read depths as low as 5× recovers the expected species relationships in phylogenetic trees. By making the tested assembly approaches available in the secapr pipeline, we hope to inspire future studies to incorporate complementary data and make an informed choice on the optimal assembly strategy.

Development of a LAMP-Based Molecular Species Diagnosis Method for Four Major Agricultural Pests in the Genus Spodoptera (Lepidoptera: Noctuidae)

Simple Summary

Four major Spodoptera pests, S. exigua, S. frugiperda, S. litura, and S. littoralis, are widely distributed polyphagous pests affecting various crops. Despite different distribution areas, these four species cause serious damage to agriculture worldwide. As these species are morphologically similar at the larval stage, diagnostic methods have been developed and utilized for their identification. Here, we developed a loop-mediated isothermal amplification (LAMP) assay for rapid and effective species diagnosis, along with PCR, to identify Korean field-collected or overseas samples. The optimal conditions for the LAMP assay were 61 °C for 60 min with four LAMP primers. Additional loop primers increased the amplification efficiency in S. exigua, whereas increased non-specific amplification was found in other species. A broad range of DNA concentrations was observed in the LAMP assay, and the minimum detectable DNA concentration was 1 pg. The DNA release method for LAMP involved incubation of larval or adult tissue samples for 5 min at 95 °C, without a DNA extraction step. Considering the gradual diversification invasive pest incidence, this simple and accurate LAMP assay can be used for intensive field monitoring of invasive pests and integrated management of these species.

Abstract

Molecular-based species identification tools are helpful to identify tiny insect and lepidopteran pests that show morphological similarities in the larval stage and are essential for quarantine as well as agricultural research. Here, we focused on four major Spodoptera pests: S. exigua, S. frugiperda, S. litura, and S. littoralis. S. exigua and S. litura mitochondrial genome sequences were newly identified and species-specific sequence regions were identified in the cytochrome c oxidase subunit II and III regions. Species primers were designed and applied in loop-mediated isothermal amplification (LAMP) and PCR to identify Korean field-collected or overseas samples. The optimal incubation conditions for LAMP were 61 °C for 60 min with four LAMP primers. Additional loop primers increased the amplification efficiency for S. exigua, and the nonspecific amplification for other species. The LAMP assay could detect a wide range of DNA concentrations, with the range 1 ng–1 pg in dependence of four LAMP primers. The DNA-releasing technique, without DNA extraction, in the LAMP assay involved larval or adult tissue sample incubation at 95 °C for 5 min. The entire process takes approximately 70 min. This new molecular diagnostic method is simple and accurate, with application in the field and laboratory and for monitoring and ecological studies.

Technology dictates algorithms: recent developments in read alignment

Aligning sequencing reads onto a reference is an essential step of the majority of genomic analysis pipelines. Computational algorithms for read alignment have evolved in accordance with technological advances, leading to today's diverse array of alignment methods. We provide a systematic survey of algorithmic foundations and methodologies across 107 alignment methods, for both short and long reads. We provide a rigorous experimental evaluation of 11 read aligners to demonstrate the effect of these underlying algorithms on speed and efficiency of read alignment. We discuss how general alignment algorithms have been tailored to the specific needs of various domains in biology.

INFIMA leverages multi-omics model organism data to identify effector genes of human GWAS variants

Genome-wide association studies reveal many non-coding variants associated with complex traits. However, model organism studies largely remain as an untapped resource for unveiling the effector genes of non-coding variants. We develop INFIMA, Integrative Fine-Mapping, to pinpoint causal SNPs for diversity outbred (DO) mice eQTL by integrating founder mice multi-omics data including ATAC-seq, RNA-seq, footprinting, and in silico mutation analysis. We demonstrate INFIMA's superior performance compared to alternatives with human and mouse chromatin conformation capture datasets. We apply INFIMA to identify novel effector genes for GWAS variants associated with diabetes. The results of the application are available at http://www.statlab.wisc.edu/shiny/INFIMA/ .

Pseudogenes: Four Decades of Discovery

A pseudogene is defined as a genomic DNA sequence that looks like a mutated or truncated version of a known functional gene. Nearly four decades since their first discovery it has been estimated that between ~12,000 and ~20,000 pseudogenes exist in the human genome. Early efforts to characterize functions for pseudogenes were unsuccessful, thus they were considered functionless relics of evolutionary selection, junk DNA or genetic fossils. Remarkably, an increasing number of pseudogenes have been reported to be expressed as RNA transcripts above and beyond levels considered accidental or spurious transcription. There is emerging evidence that some expressed pseudogene transcripts have biological functions and should be defined as a subclass of functional long noncoding RNAs (lncRNA). In this introductory chapter, I briefly summarize the history and the current knowledge of pseudogenes, and highlight the emerging functions of some pseudogenes in human biology and disease. This second iteration of Pseudogenes in Methods in Molecular Biology highlights new methodological approaches to investigate this intriguing family of lncRNAs and the extent of their biological function.

Numerical Characterization of DNA Sequences for Alignment-free Sequence Comparison - A Review

BACKGROUND

Biological macromolecules namely, DNA, RNA, and protein have their building blocks organized in a particular sequence and the sequential arrangement encodes evolutionary history of the organism (species). Hence, biological sequences have been used for studying evolutionary relationships among the species. This is usually carried out by multiple sequence algorithms (MSA). Due to certain limitations of MSA, alignment-free sequence comparison methods were developed. The present review is on alignment-free sequence comparison methods carried out using numerical characterization of DNA sequences. <P> Discussion: The graphical representation of DNA sequences by chaos game representation and other 2-dimesnional and 3-dimensional methods are discussed. The evolution of numerical characterization from the various graphical representations and the application of the DNA invariants thus computed in phylogenetic analysis is presented. The extension of computing molecular descriptors in chemometrics to the calculation of new set of DNA invariants and their use in alignment-free sequence comparison in a N-dimensional space and construction of phylogenetic tress is also reviewed. <P> Conclusion: The phylogenetic tress constructed by the alignment-free sequence comparison methods using DNA invariants were found to be better than those constructed using alignment-based tools such as PHLYIP and ClustalW. One of the graphical representation methods is now extended to study viral sequences of infectious diseases for the identification of conserved regions to design peptide-based vaccine by combining numerical characterization and graphical representation.

Complete Closed Genome Sequence of the Extremely Heat-Resistant Strain Escherichia coli AW1.7

Escherichia coli isolate AW1.7 is an extremely heat-resistant bacterium and has been widely used as a reference strain in extreme heat resistance studies for almost a decade. Here, we report its complete closed genome sequence.

Bruton's Tyrosine Kinase and Its Isoforms in Cancer

Bruton’s tyrosine kinase (BTK) is a soluble tyrosine kinase with central roles in the development, maturation, and signaling of B cells. BTK has been found to regulate cell proliferation, survival, and migration in various B-cell malignancies. Targeting BTK with recently developed BTK inhibitors has been approved by the Food and Drug Administration (FDA) for the treatment of several hematological malignancies and has transformed the treatment of several B-cell malignancies. The roles that BTK plays in B cells have been appreciated for some time. Recent studies have established that BTK is expressed and plays pro-tumorigenic roles in several epithelial cancers. In this review, we focus on novel isoforms of the BTK protein expressed in epithelial cancers. We review recent work on the expression, function, and signaling of these isoforms and their value as potential therapeutic targets in epithelial tumors.

Molecular mechanisms of mutualistic and antagonistic interactions in a plant-pollinator association

Many insects metamorphose from antagonistic larvae into mutualistic adult pollinators, with reciprocal adaptation leading to specialized insect-plant associations. It remains unknown how such interactions are established at molecular level. Here we assemble high-quality genomes of a fig species, Ficus pumila var. pumila, and its specific pollinating wasp, Wiebesia pumilae. We combine multi-omics with validation experiments to reveal molecular mechanisms underlying this specialized interaction. In the plant, we identify the specific compound attracting pollinators and validate the function of several key genes regulating its biosynthesis. In the pollinator, we find a highly reduced number of odorant-binding protein genes and an odorant-binding protein mainly binding the attractant. During antagonistic interaction, we find similar chemical profiles and turnovers throughout the development of galled ovules and seeds, and a significant contraction of detoxification-related gene families in the pollinator. Our study identifies some key genes bridging coevolved mutualists, establishing expectations for more diffuse insect-pollinator systems.

Development of a simple and accurate molecular tool for Spodoptera frugiperda species identification using LAMP

BACKGROUND

The fall armyworm, Spodoptera frugiperda is a native species of the Americas. First detected in western and central Africa in early 2016, it has become one of the most serious invasive lepidopteran pests in many African and Asian countries. S. frugiperda has spread very quickly; however, there are no molecular-based, simple and accurate diagnostic tools for identification of this species in the field. Methods to identify invasive S. frugiperda are urgently needed because farmers and agricultural managers have no prior experience with this pest.

RESULTS

Based on mitochondrial genome sequence alignment, a S. frugiperda-specific sequence region was identified in the transfer RNA-coding region between NADH dehydrogenase, ND3, and ND5. Using this unique region, species-diagnostic primers were designed and applied in a loop-mediated isothermal amplification (LAMP) assay and a conventional polymerase chain reaction to identify field-collected samples of S. frugiperda. The optimal incubation conditions for the LAMP assay were 61°C for 90 min with four LAMP primers; an additional loop primer increased the amplification efficiency. A response was obtained for a wide range of DNA concentrations in the LAMP assay and the minimum detectable DNA concentration was 10 pg.

CONCLUSIONS

We developed a new LAMP-based molecular diagnostic method that it is easy to use and accurate. The LAMP assay was used with a DNA-releasing technique for larval and adult samples, without a DNA extraction step, by incubating the tissue sample at 95°C for 5 min. This method can be applied in intensive field monitoring of S. frugiperda and its ecological studies. © 2021 Society of Chemical Industry.

Temporal dynamics of teleost populations during the Pleistocene: a report from publicly available genome data

Background

Global climate oscillation, as a selection dynamic, is an ecologically important element resulting in global biodiversity. During the glacial geological periods, most organisms suffered detrimental selection pressures (such as food shortage and habitat loss) and went through population declines. However, during the mild interglacial periods, many species re-flourished. These temporal dynamics of effective population sizes (N_e) provide essential information for understanding and predicting evolutionary outcomes during historical and ongoing global climate changes.

Results

Using high-quality genome assemblies and corresponding sequencing data, we applied the Pairwise Sequentially Markovian Coalescent (PSMC) method to quantify N_e changes of twelve representative teleost species from approximately 10 million years ago (mya) to 10 thousand years ago (kya). These results revealed multiple rounds of population contraction and expansion in most of the examined teleost species during the Neogene and the Quaternary periods. We observed that 83% (10/12) of the examined teleosts had experienced a drastic decline in N_e before the last glacial period (LGP, 110–12 kya), slightly earlier than the reported pattern of N_e changes in 38 avian species. In comparison with the peaks, almost all of the examined teleosts maintained long-term lower N_e values during the last few million years. This is consistent with increasingly dramatic glaciation during this period.

Conclusion

In summary, these findings provide a more comprehensive understanding of the historical N_e changes in teleosts. Results presented here could lead to the development of appropriate strategies to protect species in light of ongoing global climate changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07816-7.

Genomic analysis on Brazilian strains of Anaplasma marginale

Anaplasma marginale is a vector-borne pathogen that causes a disease known as anaplasmosis. No sequenced genomes of Brazilian strains are yet available. The aim of this work was to compare whole genomes of Brazilian strains of A. marginale (Palmeira and Jaboticabal) with genomes of strains from other regions (USA and Australia strains). Genome sequencing of Brazilian strains was performed by means of next-generation sequencing. Reads were mapped using the genome of the Florida strain of A. marginale as a reference sequence. Single nucleotide polymorphisms (SNPs) and insertions/deletions (INDELs) were identified. The data showed that two Brazilian strains grouped together in one particular clade, which grouped in a larger American group together with North American strains. Moreover, some important differences in surface proteins between the two Brazilian isolates can be discerned. These results shed light on the evolutionary history of A. marginale and provide the first genome information on South American isolates. Assessing the genome sequences of strains from different regions is essential for increasing knowledge of the pan-genome of this bacteria.

Learning a genome-wide score of human-mouse conservation at the functional genomics level

Identifying genomic regions with functional genomic properties that are conserved between human and mouse is an important challenge in the context of mouse model studies. To address this, we develop a method to learn a score of evidence of conservation at the functional genomics level by integrating information from a compendium of epigenomic, transcription factor binding, and transcriptomic data from human and mouse. The method, Learning Evidence of Conservation from Integrated Functional genomic annotations (LECIF), trains neural networks to generate this score for the human and mouse genomes. The resulting LECIF score highlights human and mouse regions with shared functional genomic properties and captures correspondence of biologically similar human and mouse annotations. Analysis with independent datasets shows the score also highlights loci associated with similar phenotypes in both species. LECIF will be a resource for mouse model studies by identifying loci whose functional genomic properties are likely conserved.

Vitamin D decreases silencer methylation to downregulate renin gene expression

Renin, encoded by REN, is an essential enzyme in the renin-angiotensin aldosterone system (RAAS) which is responsible for the maintenance of blood pressure homeostasis. Transcriptional regulation of REN has been linked to enhancer-promoter crosstalk, cAMP response element-binding protein (CREB), the active metabolite of vitamin D, 1,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃), and a less well-characterized intronic silencer element. We hypothesized that in addition to these, differential DNA methylation is linked to REN expression and influenced by 1,25(OH)₂D₃. REN expressing cells (HEK293) were used to elucidate the effect of 1,25(OH)₂D₃ on REN methylation and expression as quantified by methylation-sensitive qPCR and RT-qPCR, respectively. In vitro 1,25(OH)₂D₃ supplementation (10 nM) induced significant hypomethylation of the REN silencer (P < 0.050), which was linked to a significant reduction in REN expression (P < 0.010) but had no effect on enhancer methylation. In addition, 1,25(OH)₂D₃ increased VDR (P < 0.05), as well as TET1 (P < 0.05) expression, suggesting an association between 1,25(OH)₂D₃ and DNA methylation. Thus, it appears that the silencer element, which is controlled by DNA methylation and influenced by 1,25(OH)₂D₃, plays an essential role in regulating REN expression.

Evolution of genome structure in the Drosophila simulans species complex

The rapid evolution of repetitive DNA sequences, including satellite DNA, tandem duplications, and transposable elements, underlies phenotypic evolution and contributes to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented and misassembled in most contemporary genome assemblies. We generated highly contiguous de novo reference genomes for the Drosophila simulans species complex (D. simulans, D. mauritiana, and D. sechellia), which speciated ∼250,000 yr ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive sequences between these four species. We find that at least 15% of the D. simulans complex species genomes fail to align uniquely to D. melanogaster owing to structural divergence-twice the number of single-nucleotide substitutions. We also find rapid turnover of satellite DNA and extensive structural divergence in heterochromatic regions, whereas the euchromatic gene content is mostly conserved. Despite the overall preservation of gene synteny, euchromatin in each species has been shaped by clade- and species-specific inversions, transposable elements, expansions and contractions of satellite and tRNA tandem arrays, and gene duplications. We also find rapid divergence among Y-linked genes, including copy number variation and recent gene duplications from autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.

Topologically associating domains and their role in the evolution of genome structure and function in Drosophila

Topologically associating domains (TADs) were recently identified as fundamental units of three-dimensional eukaryotic genomic organization, although our knowledge of the influence of TADs on genome evolution remains preliminary. To study the molecular evolution of TADs in Drosophila species, we constructed a new reference-grade genome assembly and accompanying high-resolution TAD map for D. pseudoobscura Comparison of D. pseudoobscura and D. melanogaster, which are separated by ∼49 million years of divergence, showed that ∼30%-40% of their genomes retain conserved TADs. Comparative genomic analysis of 17 Drosophila species revealed that chromosomal rearrangement breakpoints are enriched at TAD boundaries but depleted within TADs. Additionally, genes within conserved TADs show lower expression divergence than those located in nonconserved TADs. Furthermore, we found that a substantial proportion of long genes (>50 kbp) in D. melanogaster (42%) and D. pseudoobscura (26%) constitute their own TADs, implying transcript structure may be one of the deterministic factors for TAD formation. By using structural variants (SVs) identified from 14 D. melanogaster strains, its three closest sibling species from the D. simulans species complex, and two obscura clade species, we uncovered evidence of selection acting on SVs at TAD boundaries, but with the nature of selection differing between SV types. Deletions are depleted at TAD boundaries in both divergent and polymorphic SVs, suggesting purifying selection, whereas divergent tandem duplications are enriched at TAD boundaries relative to polymorphism, suggesting they are adaptive. Our findings highlight how important TADs are in shaping the acquisition and retention of structural mutations that fundamentally alter genome organization.

Methods to Identify and Study the Evolution of Pseudogenes Using a Phylogenetic Approach

The discovery that pseudogenes are involved in important biological processes has excited enthusiasm and increased the research interest on them. An accurate detection and analysis of pseudogenes can be achieved using comparative methods, but only the use of phylogenetic tools can provide accurate information about their birth, their evolution and their death, hence about the impact that they have on genes and genomes. Here, phylogenetic methods that allow for studying pseudogene history are described.

The complete chloroplast genome sequence of the Alstroemeria hybrid variety 'Hanhera'

Alstroemeria, a member of the Alstroemeriaceae family, is a species from South America. The chloroplast genome of Alstroemeria spp. was completed by de novo assembly using a small amount of whole genome sequencing data. The chloroplast genome of Alstroemeria spp. was 155,672 bp in length consisting of 84,379 bp of large single copy, 17,815 bp of small single copy, and 26,739 bp of a pair of inverted repeat regions. A total of 157 genes were annotated including 103 protein-coding genes (PCGs), 46 tRNA genes, and eight rRNA genes. Maximum likelihood phylogenetic analysis with seven species belonging to the Alstroemeriaceae or Liliaceae family revealed that Alstroemeria spp. is grouped with the species in the Alstroemeriaceae family.

Genomic signatures of recombination in a natural population of the bdelloid rotifer Adineta vaga

Sexual reproduction is almost ubiquitous among extant eukaryotes. As most asexual lineages are short-lived, abandoning sex is commonly regarded as an evolutionary dead end. Still, putative anciently asexual lineages challenge this view. One of the most striking examples are bdelloid rotifers, microscopic freshwater invertebrates believed to have completely abandoned sexual reproduction tens of Myr ago. Here, we compare whole genomes of 11 wild-caught individuals of the bdelloid rotifer Adineta vaga and present evidence that some patterns in its genetic variation are incompatible with strict clonality and lack of genetic exchange. These patterns include genotype proportions close to Hardy-Weinberg expectations within loci, lack of linkage disequilibrium between distant loci, incongruent haplotype phylogenies across the genome, and evidence for hybridization between divergent lineages. Analysis of triallelic sites independently corroborates these findings. Our results provide evidence for interindividual genetic exchange and recombination in A. vaga, a species previously thought to be anciently asexual.