Genes, Vol. 10, Pages 492: Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.

Genes, Vol. 10, Pages 491: Genome-Wide Homozygosity Patterns and Evidence for Selection in a Set of European and Near Eastern Horse Breeds Intensive artificial and natural selection have shaped substantial variation among European horse breeds. Whereas most equine selection signature studies employ divergent genetic population structures in order to derive specific inter-breed targets of selection, we screened a total of 1476 horses originating from 12 breeds for the loss of genetic diversity by runs of homozygosity (ROH) utilizing a 670,000 single nucleotide polymorphism (SNP) genotyping array. Overlapping homozygous regions (ROH islands) indicating signatures of selection were identified by breed and similarities/dissimilarities between populations were evaluated. In the entire dataset, 180 ROH islands were identified, whilst 100 islands were breed specific, all other overlapped in 36 genomic regions with at least one ROH island of another breed. Furthermore, two ROH hot spots were determined at horse chromosome 3 (ECA3) and ECA11. Besides the confirmation of previously documented target genes involved in selection for coat color (MC1R, STX17, ASIP), body size (LCORL/NCAPG, ZFAT, LASP1, HMGA2), racing ability (PPARGC1A), behavioral traits (GRIN2B, NTM/OPCML) and gait patterns (DMRT3), several putative target genes related to embryonic morphogenesis (HOXB), energy metabolism (IGFBP-1, IGFBP-3), hair follicle morphogenesis (KRT25, KRT27, INTU) and autophagy (RALB) were highlighted. Furthermore, genes were pinpointed which might be involved in environmental adaptation of specific habitats (UVSSA, STXBP4, COX11, HLF, MMD).

Genes, Vol. 10, Pages 490: Novel Molecular Synapomorphies Demarcate Different Main Groups/Subgroups of Plasmodium and Piroplasmida Species Clarifying Their Evolutionary Relationships The class Hematozoa encompasses several clinically important genera, including Plasmodium, whose members cause the major life-threating disease malaria. Hence, a good understanding of the interrelationships of organisms from this class and reliable means for distinguishing them are of much importance. This study reports comprehensive phylogenetic and comparative analyses on protein sequences on the genomes of 28 hematozoa species to understand their interrelationships. In addition to phylogenetic trees based on two large datasets of protein sequences, detailed comparative analyses were carried out on the genomes of hematozoa species to identify novel molecular synapomorphies consisting of conserved signature indels (CSIs) in protein sequences. These studies have identified 79 CSIs that are exclusively present in specific groups of Hematozoa/Plasmodium species, also supported by phylogenetic analysis, providing reliable means for the identification of these species groups and understanding their interrelationships. Of these CSIs, six CSIs are specifically shared by all hematozoa species, two CSIs serve to distinguish members of the order Piroplasmida, five CSIs are uniquely found in all Piroplasmida species except B. microti and two CSIs are specific for the genus Theileria. Additionally, we also describe 23 CSIs that are exclusively present in all genome-sequenced Plasmodium species and two, nine, ten and eight CSIs which are specific for members of the Plasmodium subgenera Haemamoeba, Laverania, Vinckeia and Plasmodium (excluding P. ovale and P. malariae), respectively. Additionally, our work has identified several CSIs that support species relationships which are not evident from phylogenetic analysis. Of these CSIs, one CSI supports the ancestral nature of the avian-Plasmodium species in comparison to the mammalian-infecting groups of Plasmodium species, four CSIs strongly support a specific relationship of species between the subgenera Plasmodium and Vinckeia and three CSIs each that reliably group P. malariae with members of the subgenus Plasmodium and P. ovale within the subgenus Vinckeia, respectively. These results provide a reliable framework for understanding the evolutionary relationships among the Plasmodium/Piroplasmida species. Further, in view of the exclusivity of the described molecular markers for the indicated groups of hematozoa species, particularly large numbers of unique characteristics that are specific for all Plasmodium species, they provide important molecular tools for biochemical/genetic studies and for developing novel diagnostics and therapeutics for these organisms.

Genes, Vol. 10, Pages 489: Further Insights into the Architecture of the PN Promoter That Controls the Expression of the bzd Genes in Azoarcus The anaerobic degradation of benzoate in bacteria involves the benzoyl-CoA central pathway. Azoarcus/Aromatoleum strains are a major group of anaerobic benzoate degraders, and the transcriptional regulation of the bzd genes was extensively studied in Azoarcus sp. CIB. In this work, we show that the bzdR regulatory gene and the PN promoter can also be identified upstream of the catabolic bzd operon in all benzoate-degrader Azoarcus/Aromatoleum strains whose genome sequences are currently available. All the PN promoters from Azoarcus/Aromatoleum strains described here show a conserved architecture including three operator regions (ORs), i.e., OR1 to OR3, for binding to the BzdR transcriptional repressor. Here, we demonstrate that, whereas OR1 is sufficient for the BzdR-mediated repression of the PN promoter, the presence of OR2 and OR3 is required for de-repression promoted by the benzoyl-CoA inducer molecule. Our results reveal that BzdR binds to the PN promoter in the form of four dimers, two of them binding to OR1. The BzdR/PN complex formed induces a DNA loop that wraps around the BzdR dimers and generates a superstructure that was observed by atomic force microscopy. This work provides further insights into the existence of a conserved BzdR-dependent mechanism to control the expression of the bzd genes in Azoarcus strains.

Genes, Vol. 10, Pages 488: Effect of Circadian Clock and Light–Dark Cycles in Onchidium reevesii: Possible Implications for Long-Term Memory The sea slug Onchidium reevesii inhabits the intertidal zone, which is characterized by a changeable environment. Although the circadian modulation of long-term memory (LTM) is well documented, the interaction of the circadian clock with light–dark masking in LTM of intertidal animals is not well understood. We characterized the LTM of Onchidium and tested the expression levels of related genes under a light–dark (LD) cycle and constant darkness (i.e., dark–dark, or DD) cycle. Results indicated that both learning behavior and LTM show differences between circadian time (CT) 10 and zeitgeber time (ZT) 10. In LD, the cry1 gene expressed irregularly, and per2 expression displayed a daily pattern and a peak expression level at ZT 18. OnCREB1 (only in LD conditions) and per2 transcripts cycled in phase with each other. In DD, the cry1 gene had its peak expression at CT 10, and per2 expressed its peak level at CT 18. OnCREB1 had two peak expression levels at ZT 10 or ZT 18 which correspond to the time node of peaks in cry1 and per2, respectively. The obtained results provide an LTM pattern that is different from other model species of the intertidal zone. We conclude that the daily transcriptional oscillations of Onchidium for LTM were affected by circadian rhythms and LD cycle masking.

Genes, Vol. 10, Pages 487: The Complete Mitochondrial Genome of Platysternon megacephalum peguense and Molecular Phylogenetic Analysis Platysternon megacephalum is the only living representative species of Platysternidae and only three subspecies remain: P. m. megalorcephalum, P. m. shiui, and P. m. peguense. However, previous reports implied that P. m. peguense has distinct morphological and molecular features. The characterization of the mitogenome has been accepted as an efficient means of phylogenetic and evolutionary analysis. Hence, this study first determined the complete mitogenome of P. m. peguense with the aim to identify the structure and variability of the P. m. peguense mitogenome through comparative analysis. Furthermore, the phylogenetic relationship of the three subspecies was tested. Based on different tRNA gene loss and degeneration of these three subspecies, their rearrangement pathways have been inferred. Phylogenetic analysis showed that P. m. peguense is a sister group to (P. m. megalorcephalum and P. m. shiui). Furthermore, the divergence time estimation of these three subspecies coincided with the uplift of the Tibetan Plateau. This study shows that the genetic distances between P. m. peguense and the other two subspecies are comparable to interspecific genetic distances, for example within Mauremys. In general, this study provides new and meaningful insights into the evolution of the three Platysternidae subspecies.

Genes, Vol. 10, Pages 486: The Genome of the Steller Sea Lion (Eumetopias jubatus) The Steller sea lion is the largest member of the Otariidae family and is found in the coastal waters of the northern Pacific Rim. Here, we present the Steller sea lion genome, determined through DNA sequencing approaches that utilized microfluidic partitioning library construction, as well as nanopore technologies. These methods constructed a highly contiguous assembly with a scaffold N50 length of over 14 megabases, a contig N50 length of over 242 kilobases and a total length of 2.404 gigabases. As a measure of completeness, 95.1% of 4104 highly conserved mammalian genes were found to be complete within the assembly. Further annotation identified 19,668 protein coding genes. The assembled genome sequence and underlying sequence data can be found at the National Center for Biotechnology Information (NCBI) under the BioProject accession number PRJNA475770.

Genes, Vol. 10, Pages 485: The Challenge of the Sponge Suberites domuncula (Olivi, 1792) in the Presence of a Symbiotic Bacterium and a Pathogen Bacterium Sponges, which are in close contact with numerous bacteria in prey/predator, symbiotic and pathogenic relationships, must provide an appropriate response in such situations. This starts with a discriminating recognition of the partner either by a physical contact or through secreted molecules or both. We investigated the expression of the Toll-like receptor, Caspase 3/7, Tumor Necrosis Factor receptor-associated factor 6, Bcl-2 homology protein-2 and macrophage expressed genes of axenic sponge cells in the presence of a symbiotic bacterium (Endozoicomonas sp. Hex311), a pathogen bacterium (Pseudoalteromonas sp. 1A1), their exoproducts and lipopolysaccharides. The vast majority of answers are in line with what could be observed with the symbiotic bacterium. The pathogenic bacterium seems to profit from the eukaryotic cell: suppression of the production of the antibacterial compound, inhibition of the apoptosis caspase-dependent pathway, deregulation of bacterial recognition. This work contributes new scientific knowledge in the field of immunology and apoptosis in early branching metazoan harboring within its tissue and cells a large number of symbiotic bacteria.

Genes, Vol. 10, Pages 484: Population Genetic Divergence and Environment Influence the Gut Microbiome in Oregon Threespine Stickleback Much of animal-associated microbiome research has been conducted in species for which little is known of their natural ecology and evolution. Microbiome studies that combine population genetic, environment, and geographic data for wild organisms can be very informative, especially in situations where host genetic variation and the environment both influence microbiome variation. The few studies that have related population genetic and microbiome variation in wild populations have been constrained by observation-based kinship data or incomplete genomic information. Here we integrate population genomic and microbiome analyses in wild threespine stickleback fish distributed throughout western Oregon, USA. We found that gut microbiome diversity and composition partitioned more among than within wild host populations and was better explained by host population genetic divergence than by environment and geography. We also identified gut microbial taxa that were most differentially abundant across environments and across genetically divergent populations. Our findings highlight the benefits of studies that investigate host-associated microbiomes in wild organisms.

Genes, Vol. 10, Pages 483: Host-Associated Bacterial Succession during the Early Embryonic Stages and First Feeding in Farmed Gilthead Sea Bream (Sparus aurata) One of the most widely reared fish in the Mediterranean Sea is Sparus aurata. The succession of S. aurata whole-body microbiota in fertilized eggs, five, 15, 21 and 71 days post hatch (dph) larvae and the contribution of the rearing water and the provided feed (rotifers, Artemia sp. and commercial diet) to the host’s microbiota was investigated by 454 pyrosequencing of the 16S rRNA gene diversity. In total, 1917 bacterial operational taxonomic units (OTUs) were found in all samples. On average, between 93 ± 2.1 and 366 ± 9.2 bacterial OTUs per sample were found, with most of them belonging to Proteobacteria and Bacteroidetes. Ten OTUs were shared between all S. aurata stages and were also detected in the rearing water or diet. The highest OTU richness occurred at the egg stage and the lowest at the yolk sac stage (5 dph). The rearing water and diet microbial communities contributed in S. aurata microbiota without overlaps in their microbial composition and structure. The commercial diet showed higher contribution to the S. aurata microbiota than the rearing water. After stage D71 the observed microbiota showed similarities with that of adult S. aurata as indicated by the increased number of OTUs associated with γ-Proteobacteria and Firmicutes.