Comparative Sequence Analysis Research Articles

Correlation between the Colony Phenotype and Amino Acid Sequence of the Variable Vaa Antigen in Clinical Isolates of Mycoplasma hominis.

A new Mycoplasma hominis phenotype forming mini-colonies (MC) on agar and distinct from the phenotype forming typical colonies (TC) not only in size, but also in morphology, growth rate, and resistance to adverse factors, has been previously identified. In this study, the phenotype of colonies was determined and a comparative analysis of the amino acid sequence of the main variable antigen Vaa of the laboratory strain N-34 and seven clinical isolates of M. hominis was performed. It is demonstrated that the amino acid sequence of Vaa in clinical isolates forming TC (similar to the laboratory strain N-34) is entirely analogous to that of laboratory strain. Clinical isolates forming MC carry amino acid substitutions in the variable C-terminal region of Vaa, which can contribute to adhesion to eukaryotic cells and immune evasion. The connection between colony phenotype and amino acid sequence of Vaa is established.

Comparative whole genome sequence analysis of mycobacterium tuberculosis isolated from pulmonary tuberculosis and tuberculous lymphadenitis cases in northwest ethiopia: Genomic diversity, transmission dynamics and in-silico molecular drug resistance

Comparative whole genome sequence analysis of mycobacterium tuberculosis isolated from pulmonary tuberculosis and tuberculous lymphadenitis cases in northwest ethiopia: Genomic diversity, transmission dynamics and in-silico molecular drug resistance

Anaerolentibacter hominis gen. nov. sp. nov., Diplocloster hominis sp. nov. and Pilosibacter fragilis gen. nov. sp. nov., isolated from human faeces.

Three Gram-positive, obligately anaerobic bacterial strains, namely CSJ-1T, CSJ-3T, and CSJ-4T, were isolated from faeces of healthy persons. They were characterized through a combination of whole-genome sequencing, phenotypic traits, and metabolomic analysis. The genome sizes of CSJ-1T, CSJ-4T, and CSJ-3T were 3.3, 3.8, and 6.1 Mbp, with DNA G+C contents of 47.2, 48.3, and 48.8 mol%, respectively. Strain CSJ-3T was identified as representing a novel species, Diplocloster hominis (type strain CSJ-3T=CGMCC 1.18033T=JCM 36512T) of the genus Diplocloster. The 16S rRNA gene sequence similarity and whole genome average nucleotide identity (gANI) of CSJ-4T to its closest related species, Diplocloster modestus ASD 4241T, were 98.3 and 91.4 %, respectively. Comparative analysis of 16S rRNA gene sequences showed 91.6 % similarity between CSJ-1T and its closest phylogenetic neighbour, Catenibacillus scindens DSM 106146T, and 93.3 % similarity between CSJ-4T and its closest relative strain, Clostridium fessum SNUG30386T. Based on the polyphasic taxonomic results, we proposed two novel genera and three novel species. Strain CSJ-1T was identified as representing a novel species of novel genus, Anaerolentibacter hominis gen. nov. sp. nov. (type strain CSJ-1T=CGMCC 1.18046T=JCM 36511T) of the family Lachnospiraceae, and strain CSJ-4T was identified as representing a novel species of novel genus Pilosibacter fragilis gen. nov. sp. nov. (type strain CSJ-4T=CGMCC 1.18026T= JCM 36513T) of the family Clostridiaceae.

Competitive Inhibition of Okanin against Plasmodium falciparum Tyrosyl-tRNA Synthetase.

Malaria is a severe disease that presents a significant threat to human health. As resistance to current drugs continues to increase, there is an urgent need for new antimalarial medications. Aminoacyl-tRNA synthetases (aaRSs) represent promising targets for drug development. In this study, we identified Plasmodium falciparum tyrosyl-tRNA synthetase (PfTyrRS) as a potential target for antimalarial drug development through a comparative analysis of the amino acid sequences and three-dimensional structures of human and plasmodium TyrRS, with particular emphasis on differences in key amino acids at the aminoacylation site. A total of 2141 bioactive compounds were screened using a high-throughput thermal shift assay (TSA). Okanin, known as an inhibitor of LPS-induced TLR4 expression, exhibited potent inhibitory activity against PfTyrRS, while showing limited inhibition of human TyrRS. Furthermore, bio-layer interferometry (BLI) confirmed the high affinity of okanin for PfTyrRS. Molecular dynamics (MD) simulations highlighted the stable conformation of okanin within PfTyrRS and its sustained binding to the enzyme. A molecular docking analysis revealed that okanin binds to both the tyrosine and partial ATP binding sites of the enzyme, preventing substrate binding. In addition, the compound inhibited the production of Plasmodium falciparum in the blood stage and had little cytotoxicity. Thus, okanin is a promising lead compound for the treatment of malaria caused by P. falciparum.

Symbiont Diversity in Imidacloprid-Resistant and Imidacloprid-Susceptible Populations of Nilaparvata lugens (Hemiptera: Delphacidae)1

Abstract Imidacloprid, a neonicotinoid insecticide, is efficacious against hemipterans, including the brown planthopper, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae). However, frequent use of this insecticide has resulted in the development of high levels of resistance among brown planthopper populations. Endosymbionts of insects have contributed to host physiology and evolution and play a role in resistance to chemical toxins. In this study, polymerase chain reaction–denaturing gradient gel electrophoresis was used to analyze the bacterial and yeast-like symbiont communities of imidacloprid-resistant and -susceptible brown planthopper populations. The Shannon-Weaver diversity index and the evenness index indicated no differences in the richness and the expression of overall species distribution of the symbiotic communities of resistant and susceptible populations. The similarity coefficients of 0.53 and 0.56 for bacterial and yeast-like symbionts, respectively, indicated the main types of differences among microorganisms in resistant and susceptible populations. Sequence comparison analysis indicated the bacterial species in the susceptible population were members of the Enterobacteriaceae and Moraxellaceae, and those in the resistant population were members of the Enterobacteriaceae, Oxalobacteriaceae, Rhodobacteriaceae, and Sphingomonadaceae. Differences also were found in the composition of yeast-like symbionts of the two populations; Cryptococcus luteolus, Pseudozyma aphidis, and Pseudozyma antarctica were detected in the susceptible population, and Cladosporium perangustum was detected in the resistant population.

Staphylococcus hsinchuensis sp. nov., Isolated from Soymilk.

A novel coagulase-negative Staphylococcus strain (H164T) was isolated from soymilk in Taiwan. Comparative sequence analysis of the 16S rRNA gene revealed that the H164T strain is a member of the genus Staphylococcus. We used multilocus sequence analysis (MLSA) and phylogenomic analyses to demonstrate that the novel strain was closely related to Staphylococcus gallinarum, Staphylococcus nepalensis, Staphylococcus cohnii, and Staphylococcus urealyuticus. The average nucleotide identity and digital DNA-DNA hybridization values between H164T and its closest relatives were <95% and <70%, respectively. The H164T strain could also be distinguished from its closest relatives by the fermentation of d-fructose, d-maltose, d-trehalose, and d-mannitol, as well as by the activities of α-glucosidase and alkaline phosphatase. The major cellular fatty acids were C15:0 iso and C15:0 anteiso, and the predominant menaquinones were MK-7 and MK-8, respectively. The major cellular fatty acids and predominant menaquinones were C15:0 iso and C15:0 anteiso and MK-7 and MK-8, respectively. In conclusion, this strain represents a novel species, named Staphylococcus hsinchuensis sp. nov., with the type strain H164T (=BCRC 81404T = NBRC 116174T).

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development.

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Molecular characterization of lumpy skin disease virus in North Central Vietnam during 2021 and early 2022.

In October 2020, the first outbreaks of lumpy skin disease (LSD) in Lang Son Province, Vietnam were reported by our laboratory. The disease had rapidly spread to the South, and it was reported in 55 of 63 provinces and cities of Vietnam by the end of 2021. The most economic loss caused by this disease occurred in the north-central region in 2021 where approximately 46,788 LSD virus (LSDV) infected cattle and buffaloes have been reported and 8,976 animals have been culled. However, the information on this pathogen circulating in this region is missing. Here, we describe the molecular characterization of LSDV circulating in north-central Vietnam in 2021 and early 2022. In total, 155 LSDV samples were collected during this period and three of these samples from each province were further characterized by Sanger sequencing analysis based on three key maker genes (GPCR, RPO30, and p32). Sequence comparison and phylogenetic analysis based on GPCR, RPO30, and p32 genes indicated that LSDV strains circulating in north-central Vietnam are closely related to previously reported strains in Vietnam regions which bordered China and all LSDV strains were 100% identical. These results show the importance of continuous monitoring and characterization of circulating LSDV strains and are important for vaccine development for the control and eradication of LSD in Vietnam.

Analysis of genetic variations of a recombinant strain derived from the three lineages of porcine reproductive and respiratory syndrome virus-2

Since the rise of the porcine reproductive and respiratory syndrome virus (PRRSV) in China, gene mutations have frequently occurred. To understand the current prevalence and evolution of PRRSV in Shandong Province, 1,528 samples suspected of PRRSV were collected from local pig farms of different sizes. The complete genome sequence of the PRRSV strain SDLY-27 was determined by next-generation sequencing (NGS) technology. The genomic sequence of SDLY-27 was 15,363 nucleotides (nt) in length, comparative analysis of the whole genome sequence suggested that the homology between SDLY 27 and 81 PRRSV strains from China and other countries in genbank was 61.9 ~ 96.4%. This study is the first to detect recombinants from multiple recombination events among the Lineage 8 (JXA1-like strains), Lineage 5 (RespPRRSV-MLV and VR2332 strains) and Sublineage 1.5 (NADC34-like strains) in Shandong, China, and provides new data for the epidemiological study of PRRSV in China. This study enriches the epidemiological data on PRRSV in Shandong Province, China. It provides an important reference for the development of new vaccines and for the prevention and control of PRRSV in China.

Electrophilic descriptor from third‐order Taylor expansion: The role of hyperhardness

AbstractIn this study, we introduce a novel definition for the electrophilic descriptor. Following the proposal of other researchers, the foundation of our approach lies in expanding the change in the energy of an electronic system using a third‐order Taylor series, leading to the integration of hyperhardness into the formulation of the electrophilic descriptor. To validate our results, we conducted a comparative analysis of reactivity sequences among 35 electrophilic systems with known reactivities, as reported in external experimental studies. The electrophilic descriptor demonstrated a remarkable ability to accurately replicate the relative reactivities of the compounds under consideration, surpassing predictions based solely on the electrophilicity index. These findings underscore the significance of incorporating hyperhardness as a crucial reactivity descriptor. Additionally, our study provides interesting insights into the interpretation of this concept.

Comparative analysis and classification of highly divergent mouse rDNA units based on their intergenic spacer (IGS) variability.

Ribosomal DNA (rDNA) repeat units are organized into tandem clusters in eukaryotic cells. In mice, these clusters are located on at least eight chromosomes and show extensive variation in the number of repeats between mouse genomes. To analyze intra- and inter-genomic variation of mouse rDNA repeats, we selectively isolated 25 individual rDNA units using Transformation-Associated Recombination (TAR) cloning. Long-read sequencing and subsequent comparative sequence analysis revealed that each full-length unit comprises an intergenic spacer (IGS) and a ∼13.4 kb long transcribed region encoding the three rRNAs, but with substantial variability in rDNA unit size, ranging from ∼35to ∼46 kb. Within the transcribed regions of rDNA units, we found 209 variants, 70 of which are in external transcribed spacers (ETSs); but the rDNA size differences are driven primarily by IGS size heterogeneity, due to indels containing repetitive elements and some functional signals such as enhancers. Further evolutionary analysis categorized rDNA units into distinct clusters with characteristic IGS lengths; numbers of enhancers; and presence/absence of two common SNPs in promoter regions, one of which is located within promoter (p)RNA and may influence pRNA folding stability. These characteristic features of IGSs also correlated significantly with 5'ETS variant patterns described previously and associated with differential expression of rDNA units. Our results suggest that variant rDNA units are differentially regulated and open a route to investigate the role of rDNA variation on nucleolar formation and possible associations with pathology.

Computational antigenic insights into the novel NADC-34-like Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) isolate YC-2020

In this computational study, we advanced the understanding of the antigenic properties of the NADC-34-like isolate of the Porcine Reproductive and Respiratory Syndrome Virus (PRRSV), named YC-2020, relevant in veterinary pathology. We utilized sequence comparison analyses of the M and N proteins, comparing them with those of NADC34, identifying substantial amino acid homology that allowed us to highlight conserved epitopes and crucial variants. Through the application of Clustal Omega for multiple sequence alignment and platforms like Vaxijen and AllerTOP for predicting antigenic and allergenic potential, our analyses revealed important insights into the conservation and variation of epitopes essential for the development of effective diagnostic tools and vaccines. Our findings, aligned with initial experimental studies, underscore the importance of these epitopes in the development of targeted immunodiagnostic platforms and significantly contribute to the management and control of PRRSV. However, further studies are required to validate the computational predictions of antigenicity for this new viral isolate. This approach underscores the potential of computational models to enable ongoing monitoring and control of PRRSV evolution in swine. While this study provides valuable insights into the antigenic properties of the novel PRRSV isolate YC-2020 through computational analysis, it is important to acknowledge the limitations inherent to in silico predictions, specifically, the absence of laboratory validation.

Porphyromonas levii – A Boon for Periodontists?

ABSTRACT The majority of species previously categorized as Bacteroides have been reassigned into new genera. Bacteroides levii (Holdeman, Cato, and Mooretaxonomic)’s status has remained uncertain. This species shares a high degree of similarity with members of the genus Porphyromonas based on biochemical, chemical, and comparative 16s rRNA sequence analysis. As a result, Bacteroides levii (Holdeman, Cato, and Moore) was reclassified as Porphyromonas levii comb. now under the genus Porphyromonas.

Comparative single-cell RNA sequencing analysis of immune response to inactivated vaccine and natural SARS-CoV-2 infection.

Uncovering the immune response to an inactivated SARS-CoV-2 vaccine (In-Vac) and natural infection is crucial for comprehending COVID-19 immunology. Here we conducted an integrated analysis of single-cell RNA sequencing (scRNA-seq) data from serial peripheral blood mononuclear cell (PBMC) samples derived from 12 individuals receiving In-Vac compared with those from COVID-19 patients. Our study reveals that In-Vac induces subtle immunological changes in PBMC, including cell proportions and transcriptomes, compared with profound changes for natural infection. In-Vac modestly upregulates IFN-α but downregulates NF-κB pathways, while natural infection triggers hyperactive IFN-α and NF-κB pathways. Both In-Vac and natural infection alter T/B cell receptor repertoires, but COVID-19 has more significant change in preferential VJ gene, indicating a vigorous immune response. Our study reveals distinct patterns of cellular communications, including a selective activation of IL-15RA/IL-15 receptor pathway after In-Vac boost, suggesting its potential role in enhancing In-Vac-induced immunity. Collectively, our study illuminates multifaceted immune responses to In-Vac and natural infection, providing insights for optimizing SARS-CoV-2 vaccine efficacy.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5’ regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

Peptostreptococcus equinus sp. nov., isolated from horse faeces.

A Gram-stain-positive, non-spore-forming, and obligate anaerobic bacteria designated strain CBA3647T was isolated from a horse faecal sample in Jeju, Republic of Korea. Phylogenetic analysis based on 16S rRNA gene sequences showed that strain CBA3647T formed a distinct phyletic lineage from closely related species within the genus Peptostreptococcus. Based on comparative analysis of 16S rRNA gene sequences, Peptostreptococcus anaerobius ATCC 27337T is most closely related to strain CBA3647T with a 16S rRNA gene similarity of 98.31 %, while similarity to other type strains is below 98.0 %. The genomic DNA G+C content of strain CBA3647T was 30.0 mol%. The digital DNA-DNA hybridization values between strain CBA3647T and the six Peptostreptococcus species were equal to or less than 24 %. Cells were non-motile and oval-shaped cocci with catalase-positive and oxidase-negative activities. Growth occurred at 20-40 °C (optimum, 35 °C), pH 6-8 (optimum, pH 7), and in the presence of 0-2 % (w/v) NaCl (optimum, 1 %). Strain CBA3647T contained C14 : 0 iso and C16 : 0 as major fatty acids. Phenotypic, chemotaxonomic, and molecular properties of strain CBA3647T suggest that it represents a novel species in the genus Peptostreptococcus, which has been named Peptostreptococcus equinus sp. nov. The type strain is CBA3647T (=KACC 22891T= JCM 35846T).

Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis.

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

COSAP: Comparative Sequencing Analysis Platform

BackgroundRecent improvements in sequencing technologies enabled detailed profiling of genomic features. These technologies mostly rely on short reads which are merged and compared to reference genome for variant identification. These operations should be done with computers due to the size and complexity of the data. The need for analysis software resulted in many programs for mapping, variant calling and annotation steps. Currently, most programs are either expensive enterprise software with proprietary code which makes access and verification very difficult or open-access programs that are mostly based on command-line operations without user interfaces and extensive documentation. Moreover, a high level of disagreement is observed among popular mapping and variant calling algorithms in multiple studies, which makes relying on a single algorithm unreliable. User-friendly open-source software tools that offer comparative analysis are an important need considering the growth of sequencing technologies.ResultsHere, we propose Comparative Sequencing Analysis Platform (COSAP), an open-source platform that provides popular sequencing algorithms for SNV, indel, structural variant calling, copy number variation, microsatellite instability and fusion analysis and their annotations. COSAP is packed with a fully functional user-friendly web interface and a backend server which allows full independent deployment for both individual and institutional scales. COSAP is developed as a workflow management system and designed to enhance cooperation among scientists with different backgrounds. It is publicly available at https://cosap.bio and https://github.com/MBaysanLab/cosap/. The source code of the frontend and backend services can be found at https://github.com/MBaysanLab/cosap-webapi/ and https://github.com/MBaysanLab/cosap_frontend/ respectively. All services are packed as Docker containers as well. Pipelines that combine algorithms can be customized and new algorithms can be added with minimal coding through modular structure.ConclusionsCOSAP simplifies and speeds up the process of DNA sequencing analyses providing commonly used algorithms for SNV, indel, structural variant calling, copy number variation, microsatellite instability and fusion analysis as well as their annotations. COSAP is packed with a fully functional user-friendly web interface and a backend server which allows full independent deployment for both individual and institutional scales. Standardized implementations of popular algorithms in a modular platform make comparisons much easier to assess the impact of alternative pipelines which is crucial in establishing reproducibility of sequencing analyses.

CSF1R inhibitor PLX3397 depletes microglia in Mongolian gerbil Meriones unguiculatus, but not in syrian hamster Mesocricetus auratus

Microglia are the residential immune cells in the central nervous system. Their roles as innate immune cells and regulators of synaptic remodeling are critical to the development and the maintenance of the brain. Numerous studies have depleted microglia to elucidate their involvement in healthy and pathological conditions. PLX3397, a blocker of colony stimulating factor 1 receptor (CSF1R), is widely used to deplete mouse microglia due to its non-invasiveness and convenience. Recently, other small rodents, including Syrian hamsters (Mesocricetus auratus) and Mongolian gerbils (Meriones unguiculatus), have been recognized as valuable animal models for studying brain functions and diseases. However, whether microglia depletion via PLX3397 is feasible in these species remains unclear. Here, we administered PLX3397 orally via food pellets to hamsters and gerbils. PLX3397 successfully depleted gerbil microglia but had no effect on microglial density in hamsters. Comparative analysis of the CSF1R amino acid sequence in different species hints that amino acid substitutions in the juxtamembrane domain may potentially contribute to the inefficacy of PLX3397 in hamsters.

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5′ to 3′). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.