zPicture: dynamic alignment and visualization tool for analyzing conservation profiles

Regulatory effect of NK homeobox 1 (NKX2.1) on melanocortin 4 receptor (Mc4r) promoter in Mandarin fish

The melanocortin 4 receptor (MC4R) is a G protein-coupled transporter that mediates the regulation of thyroid hormones and leptin on energy balance and food intake. However, the mechanisms of transcriptional regulation of Mc4r by thyroid hormone and leptin in fish have been rarely reported. The messenger RNA expression of Mc4r gene was significantly higher in brain than those in other tissues of mandarin fish. We analyzed the structure and function of a 2029 bp sequence of Mc4r promoter. Meanwhile, overexpression of NKX2.1 and incubation with leptin significantly increased Mc4r promoter activity, but triiodothyronine showed the opposite effect. In addition, mutations in the NKX2.1 binding site abolished not only the activation of Mc4r promoter activity by leptin but also the inhibitory effect of thyroid hormones on Mc4r promoter activity. In summary, these results suggested that thyroid hormones and leptin might regulate the transcriptional expression of Mc4r through NKX2.1.

Chloroplast markers for the Malvaceae and the plastome of Henderson's checkermallow (Sidalcea hendersonii S.Wats.), a rare plant from the Pacific Northwest

OBJECTIVE

Sidalcea is a genus of flowering plants restricted to the west coast of North America, commonly known as checkermallows. Remarkably, of the ~ 30 recognized species, 16 are of conservation concern (vulnerable, imperilled or critically imperilled). To facilitate biological studies in this genus, and in the wider Malvaceae, we have sequenced the whole plastid genome of Sidalcea hendersonii. This will allow us both to check those regions already developed as general Malvaceae markers in a previous study, and to search for new regions.

RESULTS

By comparing the Sidalcea genome to that of Althaea, we have identified a hypervariable circa 1 kb region in the short single copy region. This region shows promise for examining phylogeographic pattern, hybridization and haplotype diversity. Remarkably, considering the conservation of plastome architecture between Sidalcea and Althaea, the former has a 237 bp deletion in the otherwise highly conserved inverted repeat region. Newly designed primers provide a PCR assay to determine presence of this indel across the Malvaceae. Screening of previously designed chloroplast microsatellite markers indicates two markers with variation within S. hendersonii that would be useful in future population conservation genetics.

Mitochondrial genomes do not appear to regulate flowering pattern/reproductive strategy in Cannabis sativa

Currently, the amount of genetic data for Cannabis is lacking due to the illegal nature of the plant. Our study used 73 Cannabis sativa whole-genome shotgun libraries to reveal eight different mtDNA haplotypes. The most common haplotype contained 60 of the 73 samples studied and was composed of only dioecious individuals. However, other haplotypes contained a mix of both mating strategies (i.e. monoecious and dioecious). From these haplotype groupings we further examined the fully annotated mitochondrial genomes of four hemp individuals with different mt haplotypes and recorded gene content, copy number variation and synteny. Our results revealed highly syntenic mitochondrial genomes that contained ~60 identifiable sequences for protein-coding genes, tRNAs and rRNAs and no obvious rearrangements or chimeric genes. We found no clear evidence that modern reproductive patterns are due to simple cytoplasmic male sterility mutations. It is likely the interaction between nuclear genetic components and the X/Y sex chromosomes that determines reproductive strategy. Additionally, we added 50 % more mitochondrial genomes to the publicly available repository.

Gene clusters for biosynthesis of mycosporine-like amino acids in dinoflagellate nuclear genomes: Possible recent horizontal gene transfer between species of Symbiodiniaceae (Dinophyceae)

Global warming increases the temperature of the ocean surface, which can disrupt dinoflagellate-coral symbioses and result in coral bleaching. Photosynthetic dinoflagellates of the family Symbiodiniaceae include bleaching-tolerant and bleaching-sensitive coral symbionts. Therefore, understanding the molecular mechanisms for changing symbiont diversity is potentially useful to assist recovery of coral holobionts (corals and their associated microbes, including multiple species of Symbiodiniaceae), although sexual reproduction has not been observed in the Symbiodiniaceae. Recent molecular phylogenetic analyses estimate that the Symbiodiniaceae appeared 160 million years ago and diversified into 15 groups, five genera of which now have available draft genomes (i.e., Symbiodinium, Durusdinium, Breviolum, Fugacium, and Cladocopium). Comparative genomic analyses have suggested that crown groups have fewer gene families than early-diverging groups, although many genes that were probably acquired via gene duplications and horizontal gene transfers (HGTs) have been found in each decoded genome. Because UV stress is likely a contributor to coral bleaching, and because the highly conserved gene cluster for mycosporine-like amino acid (MAA) biosynthesis has been found in thermal-tolerant symbiont genomes, I reviewed genomic features of the Symbiodiniaceae, focusing on possible acquisition of a biosynthetic gene cluster for MAAs, which absorb UV radiation. On the basis of highly conserved noncoding sequences, I hypothesized that HGTs have occurred among members of the Symbiodiniaceae and have contributed to the diversification of Symbiodiniaceae-host relationships. Finally, I proposed that bleaching tolerance may be strengthened by multiple MAAs from both symbiotic dinoflagellates and corals.

Genomic analyses of human adenoviruses unravel novel recombinant genotypes associated with severe infections in pediatric patients

Human adenoviruses (HAdVs) are highly contagious pathogens of clinical importance, especially among the pediatric population. Studies on comparative viral genomic analysis of cases associated with severe and mild infections due to HAdV are limited. Using whole-genome sequencing (WGS), we investigated whether there were any differences between circulating HAdV strains associated with severe infections (meningitis, sepsis, convulsion, sudden infant death syndrome, death, and hospitalization) and mild clinical presentations in pediatric patients hospitalized between the years 1998 and 2017 in a tertiary care hospital group in Bern, Switzerland covering a population base of approx. 2 million inhabitants. The HAdV species implicated in causing severe infections in this study included HAdV species C genotypes (HAdV1, HAdV2, and HAdV5). Clustering of the HAdV whole-genome sequences of the severe and mild cases did not show any differences except for one sample (isolated from a patient presenting with sepsis, meningitis, and hospitalization) that formed its own cluster with HAdV species C genotypes. This isolate showed intertypic recombination events involving four genotypes, had the highest homology to HAdV89 at complete genome level, but possessed the fiber gene of HAdV1, thereby representing a novel genotype of HAdV species C. The incidence of potential recombination events was higher in severe cases than in mild cases. Our findings confirm that recombination among HAdVs is important for molecular evolution and emergence of new strains. Therefore, further research on HAdVs, particularly among susceptible groups, is needed and continuous surveillance is required for public health preparedness including outbreak investigations.

Enhancer Hijacking Drives Oncogenic BCL11B Expression in Lineage-Ambiguous Stem Cell Leukemia

Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. SIGNIFICANCE: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries.This article is highlighted in the In This Issue feature, p. 2659.

Brassinosteroid-BZR1/2-WAT1 module determines the high level of auxin signalling in vascular cambium during wood formation

The tight regulation of local auxin homeostasis and signalling maxima in xylem precursor cells specifies the organising activity of the vascular cambium and consequently promotes xylem differentiation and wood formation. However, the molecular mechanisms underlying the local auxin signalling maxima in the vascular cambium are largely unknown. Here, we reveal that brassinosteroid (BR)-activated WALLS ARE THIN1 (WAT1) facilitates wood formation by enhancing local auxin signalling in the vascular cambium in Solanum lycopersicum. Growth defects and low auxin signalling readouts in the BR-deficient tomato cultivar, Micro-Tom, were associated with a novel recessive allele, Slwat1-copi, created by the insertion of a retrotransposon in the last exon of the SlWAT1 locus. Molecular and genetic studies by generating the gain-of-function and loss-of-function tomato mutants revealed that SlWAT1 is a critical regulator for fine tuning local auxin homeostasis and signalling outputs in vascular cambium to facilitate secondary growth. Finally, we discovered that BR-regulated SlBZR1/2 directly activated downstream auxin responses by SlWAT1 upregulation in xylem precursor cells to facilitate xylem differentiation and subsequent wood formation. Our data suggest that the BR-SlBZR1/2-WAT1 signalling network contributes to the high level of auxin signalling in the vascular cambium for secondary growth.

Expansion and Diversification of Fluorescent Protein Genes in Fifteen Acropora Species during the Evolution of Acroporid Corals

In addition to a purple, non-fluorescent chromoprotein (ChrP), fluorescent proteins (FPs) account for the vivid colors of corals, which occur in green (GFP), cyan (CFP), and red (RFP) FPs. To understand the evolution of the coral FP gene family, we examined the genomes of 15 Acropora species and three confamilial taxa. This genome-wide survey identified 219 FP genes. Molecular phylogeny revealed that the 15 Acropora species each have 9–18 FP genes, whereas the other acroporids examined have only two, suggesting a pronounced expansion of the FP genes in the genus Acropora. The data estimates of FP gene duplication suggest that the last common ancestor of the Acropora species that survived in the period of high sea surface temperature (Paleogene period) has already gained 16 FP genes. Different evolutionary histories of lineage-specific duplication and loss were discovered among GFP/CFPs, RFPs, and ChrPs. Synteny analysis revealed core GFP/CFP, RFP, and ChrP gene clusters, in which a tandem duplication of the FP genes was evident. The expansion and diversification of Acropora FPs may have contributed to the present-day richness of this genus.

Conservation of peripheral nervous system formation mechanisms in divergent ascidian embryos

Ascidians with very similar embryos but highly divergent genomes are thought to have undergone extensive developmental system drift. We compared, in four species (Ciona and Phallusia for Phlebobranchia, Molgula and Halocynthia for Stolidobranchia), gene expression and gene regulation for a network of six transcription factors regulating peripheral nervous system (PNS) formation in Ciona. All genes, but one in Molgula, were expressed in the PNS with some differences correlating with phylogenetic distance. Cross-species transgenesis indicated strong levels of conservation, except in Molgula, in gene regulation despite lack of sequence conservation of the enhancers. Developmental system drift in ascidians is thus higher for gene regulation than for gene expression and is impacted not only by phylogenetic distance, but also in a clade-specific manner and unevenly within a network. Finally, considering that Molgula is divergent in our analyses, this suggests deep conservation of developmental mechanisms in ascidians after 390 My of separate evolution.

Bioinformatics analysis reveals four major hexon variants of human adenovirus type-3 (HAdV-3) as the potential strains for development of vaccine and siRNA-based therapeutics against HAdV-3 respiratory infections

Human adenovirus type 3 (HAdV-3) encompasses 15-87% of all adenoviral respiratory infections. The significant morbidity and mortality, especially among the neonates and immunosuppressed patients, demand the need for a vaccine or a targeted antiviral against this type. However, due to the existence of multiple hexon variants (3Hv-1 to 3Hv-25), the selection of vaccine strains of HAdV-3 is challenging. This study was designed to evaluate HAdV-3 hexon variants for the selection of potential vaccine candidates and the use of hexon gene as a target for designing siRNA that can be used as a therapy. Based on the data of worldwide distribution, duration of circulation, co-circulation and their percentage among all the variants, 3Hv-1 to 3Hv-4 were categorized as the major hexon variants. Phylogenetic analysis and the percentage of homology in the hypervariable regions followed by multi-sequence alignment, zPicture analysis and restriction enzyme analysis were carried out. In the phylogram, the variants were arranged in different clusters. The HVR encoding regions of hexon of 3Hv-1 to 3Hv-4 showed 16 point mutations resulting in 12 amino acids substitutions. The homology in HVRs was 81.81-100%. Therefore, the major hexon variants are substantially different from each other which justifies their inclusion as the potential vaccine candidates. Interestingly, despite the significant differences in the DNA sequence, there were many conserved areas in the HVRs, and we have designed functional siRNAs form those locations. We have also designed immunogenic vaccine peptide epitopes from the hexon protein using bioinformatics prediction tool. We hope that our developed siRNAs and immunogenic vaccine peptide epitopes could be used in the future development of siRNA-based therapy and designing a vaccine against HAdV-3.

Complete mitochondrial genomes provide current refined phylogenomic hypotheses for relationships among ten Hirundo species

Hirundo is the most species-rich genus of the passerine swallow family (Hirundinidae) and has a cosmopolitan distribution. Here we report the complete, annotated mitochondrial genomes for 25 individuals from 10 of the 14 extant Hirundo species; these include representatives from four subspecies of the barn swallow, H. rustica. Mitogenomes were conserved in size, ranging from 18,500 to 18,700 base pairs. They all contained 13 protein-coding regions, 22 tRNAs, a control region, and large and small ribosomal subunits. Phylogenetic analysis resolved most of the relationships between the studied species and subspecies which were largely consistent with previously published trees. Several new relationships were observed within the phylogeny that could have only been discovered with the increased amount of genetic material. This study represents the largest Hirundo mitochondrial phylogeny to date, and could serve as a vital tool for other studies focusing on the evolution of the Hirundo genus.

A Zoonotic Adenoviral Human Pathogen Emerged through Genomic Recombination among Human and Nonhuman Simian Hosts

Genomics analysis of a historically intriguing and predicted emergent human adenovirus (HAdV) pathogen, which caused pneumonia and death, provides insight into a novel molecular evolution pathway involving "ping-pong" zoonosis and anthroponosis. The genome of this promiscuous pathogen is embedded with evidence of unprecedented multiple, multidirectional, stable, and reciprocal cross-species infections of hosts from three species (human, chimpanzee, and bonobo). This recombinant genome, typed as HAdV-B76, is identical to two recently reported simian AdV (SAdV) genomes isolated from chimpanzees and bonobos. Additionally, the presence of a critical adenoviral replication element found in HAdV genomes, in addition to genes that are highly similar to counterparts in other HAdVs, reinforces its potential as a human pathogen. Reservoirs in nonhuman hosts may explain periods of apparent absence and then reemergence of human adenoviral pathogens, as well as present pathways for the genesis of those thought to be newly emergent. The nature of the HAdV-D76 genome has implications for the use of SAdVs as gene delivery vectors in human gene therapy and vaccines, selected to avoid preexisting and potentially fatal host immune responses to HAdV.IMPORTANCE An emergent adenoviral human pathogen, HAdV-B76, associated with a fatality in 1965, shows a remarkable degree of genome identity with two recently isolated simian adenoviruses that contain cross-species genome recombination events from three hosts: human, chimpanzee, and bonobo. Zoonosis (nonhuman-to-human transmission) and anthroponosis (human to nonhuman transmission) may play significant roles in the emergence of human adenoviral pathogens.

Identification of neuron-type specific promoters in monkey genome and their functional validation in mice

Viral gene delivery is one of the most versatile techniques for elucidating the mechanisms underlying brain dysfunction, such as neuropsychiatric disorders. Due to the complexity of the brain, expression of genetic tools, such as channelrhodopsin and calcium sensors, often has to be restricted to a specified cell type within a circuit implicated in these disorders. Only a handful of promoters targeting neuronal subtypes are currently used for viral gene delivery. Here, we isolated conserved promoter regions of several subtype-specific genes from the macaque genome and investigated their functionality in the mouse brain when used within lentiviral vectors (LVVs). Immunohistochemical analysis revealed that transgene expression induced by the promoter sequences for somatostatin (SST), cholecystokinin (CCK), parvalbumin (PV), serotonin transporter (SERT), vesicular acetylcholine transporter (vAChT), substance P (SP) and proenkephalin (PENK) was largely colocalized with specific markers for the targeted neuronal populations. Moreover, by combining these results with in silico predictions of transcription factor binding to the isolated sequences, we identified transcription factors possibly underlying cell-type specificity. These findings lay a foundation for the expansion of the current toolbox of promoters suitable for elucidating these neuronal phenotypes.

Whole-Genome Sequence Data Uncover Widespread Heterothallism in the Largest Group of Lichen-Forming Fungi

Fungal reproduction is regulated by the mating-type (MAT1) locus, which typically comprises two idiomorphic genes. The presence of one or both allelic variants at the locus determines the reproductive strategy in fungi—homothallism versus heterothallism. It has been hypothesized that self-fertility via homothallism is widespread in lichen-forming fungi. To test this hypothesis, we characterized the MAT1 locus of 41 genomes of lichen-forming fungi representing a wide range of growth forms and reproductive strategies in the class Lecanoromycetes, the largest group of lichen-forming fungi. Our results show the complete lack of genetic homothallism suggesting that lichens evolved from a heterothallic ancestor. We argue that this may be related to the symbiotic lifestyle of these fungi, and may be a key innovation that has contributed to the accelerated diversification rates in this fungal group.

De novo origination of MIRNAs through generation of short inverted repeats in target genes

MIRNA (MIR) gene origin and early evolutionary processes, such as hairpin precursor sequence origination, promoter activity acquirement and the sequence of these two processes, are fundamental and fascinating subjects. Three models, including inverted gene duplication, spontaneous evolution and transposon transposition, have been proposed for de novo origination of hairpin precursor sequence. However, these models still open to discussion. In addition, de novo origination of MIR gene promoters has not been well investigated. Here, I systematically investigated the origin of evolutionarily young polyphenol oxidase gene (PPO)-targeting MIRs, including MIR1444, MIR058 and MIR12112, and a genomic region termed AasPPO-as-hp, which contained a hairpin-forming sequence. I found that MIR058 precursors and the hairpin-forming sequence of AasPPO-as-hp originated in an ancient PPO gene through forming short inverted repeats. Palindromic-like sequences and imperfect inverted repeats in the ancient PPO gene contributed to initiate the generation of short inverted repeats probably by causing errors during DNA duplication. Analysis of MIR058 and AasPPO-as-hp promoters showed that they originated in the 3'-flanking region of the ancient PPO gene. Promoter activities were gained by insertion of a CAAT-box and multiple-copper-response element (CuRE)-containing miniature inverted-repeat transposable element (MITE) in the upstream of AT-rich TATA-box-like sequence. Gain of promoter activities occurred before hairpin-forming sequence origination. Sequence comparison of MIR1444, MIR058 and MIR12112 promoters showed frequent birth and death of CuREs, indicating copper could be vital for the origination and evolution of PPO-targeting MIRs. Based on the evidence obtained, a novel model for plant MIR origination and evolution is proposed.

Halobacterium salinarum virus ChaoS9, a Novel Halovirus Related to PhiH1 and PhiCh1

The unexpected lysis of a large culture of Halobacterium salinarum strain S9 was found to be caused by a novel myovirus, designated ChaoS9. Virus purification from the culture lysate revealed a homogeneous population of caudovirus-like particles. The viral genome is linear, dsDNA that is partially redundant and circularly permuted, has a unit length of 55,145 nt, a G + C% of 65.3, and has 85 predicted coding sequences (CDS) and one tRNA (Arg) gene. The left arm of the genome (0–28 kbp) encodes proteins similar in sequence to those from known caudoviruses and was most similar to myohaloviruses phiCh1 (host: Natrialba magadii) and phiH1 (host: Hbt. salinarum). It carries a tail-fiber gene module similar to the invertible modules present in phiH1 and phiCh1. However, while the tail genes of ChaoS9 were similar to those of phiCh1 and phiH1, the Mcp of ChaoS9 was most similar (36% aa identity) to that of Haloarcula hispanica tailed virus 1 (HHTV-1). Provirus elements related to ChaoS9 showed most similarity to tail/assembly proteins but varied in their similarity with head/assembly proteins. The right arm (29–55 kbp) of ChaoS9 encoded proteins involved in DNA replication (ParA, RepH, and Orc1) but the other proteins showed little similarity to those from phiH1, phiCh1, or provirus elements, and most of them could not be assigned a function. ChaoS9 is probably best classified within the genus Myohalovirus, as it shares many characteristics with phiH1 (and phiCh1), including many similar proteins. However, the head/assembly gene region appears to have undergone a recombination event, and the inferred proteins are different to those of phiH1 and phiCh1, including the major capsid protein. This makes the taxonomic classification of ChaoS9 more ambiguous. We also report a revised genome sequence and annotation of Natrialba virus phiCh1.

Variation in the Promoter Region of the MC4R Gene Elucidates the Association of Body Measurement Traits in Hu Sheep

The melanocortin 4 receptor (MC4R) gene is expressed in the appetite-regulating areas of the brain and is engaged in the leptin signaling pathway. Although previous studies have identified variants in the coding region of the sheep MC4R gene showing significant associations with birth weight, weaning weight, and backfat thickness, no such associations have been reported for the promoter region. Besides, the essential promoter region of the sheep MC4R has not been delineated. In this study, to better understand the transcriptional regulation of MC4R and to elucidate the association between regulatory variants and haplotypes with body measurement traits in sheep, we cloned and characterized the MC4R promoter. We found that the minimal promoter of the gene is located within the region -1207/-880 bp upstream of the first exon. Real-time quantitative PCR (RT-qPCR) data revealed the mRNA expression of the MC4R gene had a significant difference between sex and age. In the association analysis, eight single nucleotide polymorphisms (SNPs) had a significant association with one or more traits (p < 0.05); of these, two SNPs were novel. Notably, individuals with haplotype H1H2 (CT-GA-GT-GA-GT-GA-GA-CG) were heavier in body weight than other haplotypes. Altogether, variations in the MC4R gene promoter, most notably haplotype H1H2, may greatly benefit marker-assisted selection in sheep.

Plant Defenses Against Pests Driven by a Bidirectional Promoter

The plant defense responses to pests results in the synchronized change of a complex network of interconnected genes and signaling pathways. An essential part of this process is mediated by the binding of transcription factors to the specific responsive cis-elements within in the promoters of phytophagous-responsive genes. In this work, it is reported the identification and characterization of a bidirectional promoter that simultaneously co-regulate two divergent genes, At5g10300 and At5g10290, upon arthropod feeding. Computational analysis identified the presence of cis-elements within the intergenic region between two loci, mainly from the DOF but also from the AP2/ERF, Golden 2-like and bHLH families. The function of the bidirectional promoter was analyzed using two enhanced variants of the GFP and CherryFP reporter genes, in both orientations, in transient tobacco and stably transformed Arabidopsis plants. Promoter activity was tested in response to feeding of Tetranychus urticae and Pieris brassicae, as well as wounding, flagellin and chitin treatments. Using RT-qPCR assays and confocal microscopy, it was shown that all treatments resulted in the induction of both reporter genes. Furthermore, our findings revealed the asymmetric character of the promoter with stronger activity in the forward than in the reverse orientation. This study provides an example of a bidirectional promoter with a strong potential to be used in plant biotechnology in pest control that requires stacking of the defense genes.

A draft genome of the striped catfish, Pangasianodon hypophthalmus, for comparative analysis of genes relevant to development and a resource for aquaculture improvement

Background

The striped catfish, Pangasianodon hypophthalmus, is a freshwater and benthopelagic fish common in the Mekong River delta. Catfish constitute a valuable source of dietary protein. Therefore, they are cultured worldwide, and P. hypophthalmus is a food staple in the Mekong area. However, genetic information about the culture stock, is unavailable for breeding improvement, although genetics of the channel catfish, Ictalurus punctatus, has been reported. To acquire genome sequence data as a useful resource for marker-assisted breeding, we decoded a draft genome of P. hypophthalmus and performed comparative analyses.

Results

Using the Illumina platform, we obtained both nuclear and mitochondrial DNA sequences. Molecular phylogeny using the mitochondrial genome confirmed that P. hypophthalmus is a member of the family Pangasiidae and is nested within a clade including the families Cranoglanididae and Ictaluridae. The nuclear genome was estimated at approximately 700 Mb, assembled into 568 scaffolds with an N50 of 14.29 Mbp, and was estimated to contain ~ 28,600 protein-coding genes, comparable to those of channel catfish and zebrafish. Interestingly, zebrafish produce gadusol, but genes for biosynthesis of this sunscreen compound have been lost from catfish genomes. The differences in gene contents between these two catfishes were found in genes for vitamin D-binding protein and cytosolic phospholipase A₂, which have lost only in channel catfish. The Hox cluster in catfish genomes comprised seven paralogous groups, similar to that of zebrafish, and comparative analysis clarified catfish lineage-specific losses of A5a, B10a, and A11a. Genes for insulin-like growth factor (IGF) signaling were conserved between the two catfish genomes. In addition to identification of MHC class I and sex determination-related gene loci, the hypothetical chromosomes by comparison with the channel catfish demonstrated the usefulness of the striped catfish genome as a marker resource.

Conclusions

We developed genomic resources for the striped catfish. Possible conservation of genes for development and marker candidates were confirmed by comparing the assembled genome to that of a model fish, Danio rerio, and to channel catfish. Since the catfish genomic constituent resembles that of zebrafish, it is likely that zebrafish data for gene functions is applicable to striped catfish as well.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5079-x) contains supplementary material, which is available to authorized users.

An Engineered Virus Library as a Resource for the Spectrum-wide Exploration of Virus and Vector Diversity

Adenoviruses (Ads) are large human-pathogenic double-stranded DNA (dsDNA) viruses presenting an enormous natural diversity associated with a broad variety of diseases. However, only a small fraction of adenoviruses has been explored in basic virology and biomedical research, highlighting the need to develop robust and adaptable methodologies and resources. We developed a method for high-throughput direct cloning and engineering of adenoviral genomes from different sources utilizing advanced linear-linear homologous recombination (LLHR) and linear-circular homologous recombination (LCHR). We describe 34 cloned adenoviral genomes originating from clinical samples, which were characterized by next-generation sequencing (NGS). We anticipate that this recombineering strategy and the engineered adenovirus library will provide an approach to study basic and clinical virology. High-throughput screening (HTS) of the reporter-tagged Ad library in a panel of cell lines including osteosarcoma disease-specific cell lines revealed alternative virus types with enhanced transduction and oncolysis efficiencies. This highlights the usefulness of this resource.

BMP3 expression by osteoblast lineage cells is regulated by canonical Wnt signaling

Bone morphogenetic protein (BMP) and canonical Wnt (cWnt) signaling factors are both known to regulate bone mass, fracture risk, fracture repair, and osteoblastogenesis. BMP3 is the most abundant BMP and negatively regulates osteoblastogenesis and bone mass. Thus, identifying the mechanism by which BMP3 acts to depress bone formation may allow for the development of new therapeutics useful in the treatment for osteopenia and osteoporosis. Here, we report that cWnt signaling stimulates BMP3 expression in osteoblast (OB) lineage cells. The expression of BMP3 increases with OB differentiation. Treatment of cells with various cWnt proteins stimulated BMP3 expression. Mice with enhanced cWnt signaling had high expression levels of BMP3. Our data suggest that reduction in BMP3 levels may contribute beneficially to the positive effect of cWnt agonists on bone mass.

The complete mitochondrial genome for Cannabis sativa

The following report details the first annotated mitochondrial genome for the Carmagnola variety of Cannabis sativa, the first reference genome for the Cannabaceae family. The total length is 415,499 bp and contains 54 genes, which sub-divide into 38 protein-coding genes, 15 tRNA genes, and 3 rRNA genes.

An Analysis of Adenovirus Genomes Using Whole Genome Software Tools

The evolution of sequencing technology has lead to an enormous increase in the number of genomes that have been sequenced. This is especially true in the field of virus genomics. In order to extract meaningful biological information from these genomes, whole genome data mining software tools must be utilized. Hundreds of tools have been developed to analyze biological sequence data. However, only some of these tools are user-friendly to biologists. Several of these tools that have been successfully used to analyze adenovirus genomes are described here. These include Artemis, EMBOSS, pDRAW, zPicture, CoreGenes, GeneOrder, and PipMaker. These tools provide functionalities such as visualization, restriction enzyme analysis, alignment, and proteome comparisons that are extremely useful in the bioinformatics analysis of adenovirus genomes.

The peroxisome proliferator-activated receptor α agonist, AZD4619, induces alanine aminotransferase-1 gene and protein expression in human, but not in rat hepatocytes: Correlation with serum ALT levels

Alanine aminotransferase (ALT) in serum is the standard biomarker for liver injury. We have previously described a clinical trial with a novel selective peroxisome proliferator-activated receptor α (PPARα) agonist (AZD4619), which unexpectedly caused increased serum levels of ALT in treated individuals without any other evidence of liver injury. We pinpointed a plausible mechanism through which AZD4619 could increase serum ALT levels; namely through the PPARα-specific activation of the human ALT1 gene at the transcriptional level. In the present study, we present data from the preceding rat toxicity study, demonstrating that AZD4619 had no effect on rat serum ALT activity levels, and further experiments were performed to elucidate the mechanisms responsible for this species-related difference. Our results revealed that AZD4619 increased ALT1 protein expression in a dose-dependent manner in human, but not in rat primary hepatocytes. Cloning of the human and rat ALT1 promoters into luciferase vectors confirmed that AZD4619 induced only the human, but not the rat ALT1 gene promoter in a dose-dependent manner. In PPARα-GAL4 reporter gene assays, AZD4619 was >100-fold more potent on the human vs. rat PPARα levels, explaining the differences in induction of the ALT1 gene between the species at the concentration range tested. These data demonstrate the usefulness of the human and rat ALT1 reporter gene assays for testing future drug candidates at the preclinical stage. In drug discovery projects, these assays elucidate whether elevations in ALT levels observed in vivo or in the clinic are due to metabolic effects rather than a toxic event in the liver.

Muscle regulatory factors regulate T1R3 taste receptor expression

T1R3 is a T1R class of G protein-coupled receptors, composing subunit of the umami taste receptor when complexed with T1R1. T1R3 was originally discovered in gustatory tissue but is now known to be expressed in a wide variety of tissues and cell types such the intestine, pancreatic β-cells, skeletal muscle, and heart. In addition to taste recognition, the T1R1/T1R3 complex functions as an amino acid sensor and has been proposed to be a control mechanism for the secretion of hormones, such as cholecystokinin, insulin, and duodenal HCO3(-) and activates the mammalian rapamycin complex 1 (MTORC1) to inhibit autophagy. T1R3 knockout mice have increased rate of autophagy in the heart, skeletal muscle and liver. Thus, T1R3 has multiple physiological functions and is widely expressed in vivo. However, the exact mechanisms regulating T1R3 expression are largely unknown. Here, we used comparative genomics and functional analyses to characterize the genomic region upstream of the annotated transcriptional start of human T1R3. This revealed that the T1R3 promoter in human and mouse resides in an evolutionary conserved region (ECR). We also identified a repressive element located upstream of the human T1R3 promoter that has relatively high degree of conservation with rhesus macaque. Additionally, the muscle regulatory factors MyoD and Myogenin regulate T1R3 expression and T1R3 expression increases with skeletal muscle differentiation of murine myoblast C2C12 cells. Taken together, our study raises the possibility that MyoD and Myogenin might control skeletal muscle metabolism and homeostasis through the regulation of T1R3 promoter activity.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

Human adenovirus type 8: the major agent of epidemic keratoconjunctivitis (EKC)

Human adenovirus type 8 (HAdV-8) is the most common causative agent of a highly contagious eye disease known as epidemic keratoconjunctivitis (EKC). HAdV-8 strains have been classified into genome types HAdV-8A to 8K and HAdV/D1 to D12 according to restriction endonuclease analysis. This review focuses on the significance of HAdV-8 as an agent of EKC. Molecular analysis of HAdV-8 genome types HAdV-53 and HAdV-54 was performed to reveal potential genetic variation in the hexon and fiber, which might affect the antigenicity and tropism of the virus, respectively. On the basis of the published data, three patterns of HAdV-8 genome type distribution were observed worldwide: (1) genome types restricted to a microenvironment, (2) genome types distributed within a country, and (3) globally dispersed genome types. Simplot and zPicture showed that the HAdV-8 genome types were nearly identical to each other. HAdV-54 is very close to the HAdV-8P, B and E genomes, except in the hexon. In a restriction map, HAdV-8P, B, and E share a very high percentage of restriction sites with each other. Hypervariable regions (HVRs) of the hexon were conserved and were 100% identical among the genome types. The fiber knob of HAdV-8P, A, E, J and HAdV-53 were 100% identical. In phylogeny, HVRs of the hexon and fiber knob of the HAdV-8 genome types segregated into monophyletic clusters. Neutralizing antibodies against one genome type will provide protection against other genome types, and the selection of future vaccine strains would be simple due to the stable HVRs. Molecular analysis of whole genomes, particularly of the capsid proteins of the remaining genome types, would be useful to substantiate our observations.

PsasM2I, a type II restriction-modification system in Pseudomonas savastanoi pv. savastanoi: differential distribution of carrier strains in the environment and the evolutionary history of homologous RM systems in the Pseudomonas syringae complex

A type II restriction-modification system was found in a native plasmid of Pseudomonas savastanoi pv. savastanoi MLLI2. Functional analysis of the methyltransferase showed that the enzyme acts by protecting the DNA sequence CTGCAG from cleavage. Restriction endonuclease expression in recombinant Escherichia coli cells resulted in mutations in the REase sequence or transposition of insertion sequence 1A in the coding sequence, preventing lethal gene expression. Population screening detected homologous RM systems in other P. savastanoi strains and in the Pseudomonas syringae complex. An epidemiological survey carried out by sampling olive and oleander knots in two Italian regions showed an uneven diffusion of carrier strains, whose presence could be related to a selective advantage in maintaining the RM system in particular environments or subpopulations. Moreover, carrier strains can coexist in the same orchards, plants, and knot tissues with non-carriers, revealing unexpected genetic variability on a very small spatial scale. Phylogenetic analysis of the RM system and housekeeping gene sequences in the P. syringae complex demonstrated the ancient acquisition of the RM systems. However, the evolutionary history of the gene complex also showed the involvement of horizontal gene transfer between related strains and recombination events.

Use of single-nucleotide polymorphisms (SNPs) to distinguish gene expression subtypes of chronic fatigue syndrome/myalgic encephalomyelitis (CFS/ME)

AIMS

We have reported gene expression changes in patients with chronic fatigue syndrome/myalgic encephalomyelitis (CFS/ME) and the fact that such gene expression data can be used to identify subtypes of CFS/ME with distinct clinical phenotypes. Due to the difficulties in using a comparative gene expression method as an aid to CFS/ME disease and subtype-specific diagnosis, we have attempted to develop such a method based on single-nucleotide polymorphism (SNP) analysis.

METHODS

To identify SNP allele associations with CFS/ME and CFS/ME subtypes, we tested genomic DNA of patients with CFS/ME (n=108), patients with endogenous depression (n=17) and normal blood donors (n=68) for 504 human SNP alleles located within 88 CFS-associated human genes using the SNP Genotyping GoldenGate Assay (Illumina, San Diego, California, USA). 360 ancestry informative markers (AIM) were also examined.

RESULTS

21 SNPs were significantly associated with CFS/ME compared with depression and normal groups. 148 SNP alleles had a significant association with one or more CFS/ME subtypes. For each subtype, associated SNPs tended to be grouped together within particular genes. AIM SNPs indicated that 4 subjects were of Asian origin while the remainder were Caucasian. Hierarchical clustering of AIM data revealed the relatedness between 2 couples of patients with CFS only and confirmed the overall heterogeneity of all subjects.

CONCLUSIONS

This study provides evidence that human SNPs located within CFS/ME associated genes are associated with particular genomic subtypes of CFS/ME. Further work is required to develop this into a clinically useful subtype-specific diagnostic test.

A piggyBac insertion disrupts Foxl2 expression that mimics BPES syndrome in mice

Blepharophimosis, ptosis, epicanthus inversus syndrome (BPES) is an autosomal dominant genetic disorder characterized by small palpebral fissures and other craniofacial malformations, often with (type I) but could also without (type II) premature ovarian failure. While mutations of the forkhead transcription factor FOXL2 are associated with and likely be responsible for many BPES cases, how FOXL2 affects craniofacial development remain to be understood. Through a large-scale piggyBac (PB) insertion mutagenesis, we have identified a mouse mutant carrying a PB insertion ∼160 kb upstream of the transcription start site (TSS) of Foxl2. The insertion reduces, but not eliminates, the expression of Foxl2. This mutant, but not its revertant, displays BPES-like conditions such as midface hypoplasia, eyelid abnormalities and female subfertility. Further analysis indicates that the mutation does not affect mandible, but causes premature fusion of the premaxilla-maxilla suture, smaller premaxilla and malformed maxilla during midface development. We further identified an evolutionarily conserved fragment near the insertion site and observed enhancer activity of this element in tissue culture cells. Analyses using DNase I hypersensitivity assay and chromosome conformation capture assay in developing maxillary and periocular tissues suggest that the DNA region near the insertion site likely interacts with Foxl2 TSS. Therefore, this mutant presents an excellent animal model for mechanistic study of BPES and regulation of Foxl2.

Evolution of the miR199-214 cluster and vertebrate skeletal development

MicroRNA (miRs) are short non-coding RNAs that fine-tune the regulation of gene expression to coordinate a wide range of biological processes. MicroRNAs are transcribed from miR genes and primary miR transcripts are processed to approximately 22 nucleotide single strand mature forms that function as repressors of transcript translation when bound to the 3'UTR of protein coding transcripts in association with the RISC. Because of their role in the regulation of gene expression, miRs are essential players in development by acting on cell fate determination and progression toward cell differentiation. The miR199 and miR214 genes occupy an intronic cluster located on the opposite strand of the Dynamin3 gene. These miRNAs play major roles in a broad variety of developmental processes and diseases, including skeletal development and several types of cancer. In the work reported here, we first deciphered the origin of the miR199 and miR214 families by following evolution of miR paralogs and their host Dynamin paralogs. We then examined the expression patterns of miR199 and miR214 in developing zebrafish embryos and demonstrated their regulation through a common primary transcript. Results suggest an evolutionarily conserved regulation across vertebrate lineages. Our expression study showed predominant expression patterns for both miR in tissues surrounding developing craniofacial skeletal elements consistent with expression data in mouse and human, thus indicating a conserved role of miR199 and miR214 in vertebrate skeletogenesis.

The vertebrate RCAN gene family: novel insights into evolution, structure and regulation

Recently there has been much interest in the Regulators of Calcineurin (RCAN) proteins which are important endogenous modulators of the calcineurin-NFATc signalling pathway. They have been shown to have a crucial role in cellular programmes such as the immune response, muscle fibre remodelling and memory, but also in pathological processes such as cardiac hypertrophy and neurodegenerative diseases. In vertebrates, the RCAN family form a functional subfamily of three members RCAN1, RCAN2 and RCAN3 whereas only one RCAN is present in the rest of Eukarya. In addition, RCAN genes have been shown to collocate with RUNX and CLIC genes in ACD clusters (ACD21, ACD6 and ACD1). How the RCAN genes and their clustering in ACDs evolved is still unknown. After analysing RCAN gene family evolution using bioinformatic tools, we propose that the three RCAN vertebrate genes within the ACD clusters, which evolved from single copy genes present in invertebrates and lower eukaryotes, are the result of two rounds of whole genome duplication, followed by a segmental duplication. This evolutionary scenario involves the loss or gain of some RCAN genes during evolution. In addition, we have analysed RCAN gene structure and identified the existence of several characteristic features that can be involved in RCAN evolution and gene expression regulation. These included: several transposable elements, CpG islands in the 5′ region of the genes, the existence of antisense transcripts (NAT) associated with the three human genes, and considerable evidence for bidirectional promoters that regulate RCAN gene expression. Furthermore, we show that the CpG island associated with the RCAN3 gene promoter is unmethylated and transcriptionally active. All these results provide timely new insights into the molecular mechanisms underlying RCAN function and a more in depth knowledge of this gene family whose members are obvious candidates for the development of future therapies.

Expression of TLE3 by bone marrow stromal cells is regulated by canonical Wnt signaling

Transducing-like enhancer of split 3 (TLE3), one of the Groucho/TLE family members, targets Runx2 transcription and suppresses osteoblast differentiation in bone marrow stromal cells (BMSCs). Here, we identify Wnt responsive elements of the TLE3 promoter region through comparative genomic and functional analyses and show that expression of TLE3 is increased by Wnt signaling, which is important for osteoblast differentiation. We also demonstrated that TLE3 is able to suppress canonical Wnt signaling in BMSCs. Taken together, our data suggest that induction of TLE3 by Wnt signaling is part of a negative feedback loop active during osteoblast differentiation.

Mechanisms regulating GLUT4 transcription in skeletal muscle cells are highly conserved across vertebrates

The glucose transporter 4 (GLUT4) plays a key role in glucose uptake in insulin target tissues. This transporter has been extensively studied in many species in terms of its function, expression and cellular traffic and complex mechanisms are involved in its regulation at many different levels. However, studies investigating the transcription of the GLUT4 gene and its regulation are scarce. In this study, we have identified the GLUT4 gene in a teleost fish, the Fugu (Takifugu rubripes), and have cloned and characterized a functional promoter of this gene for the first time in a non-mammalian vertebrate. In silico analysis of the Fugu GLUT4 promoter identified potential binding sites for transcription factors such as SP1, C/EBP, MEF2, KLF, SREBP-1c and GC-boxes, as well as a CpG island, but failed to identify a TATA box. In vitro analysis revealed three transcription start sites, with the main residing 307 bp upstream of the ATG codon. Deletion analysis determined that the core promoter was located between nucleotides -132/+94. By transfecting a variety of 5´deletion constructs into L6 muscle cells we have determined that Fugu GLUT4 promoter transcription is regulated by insulin, PG-J2, a PPARγ agonist, and electrical pulse stimulation. Furthermore, our results suggest the implication of motifs such as PPARγ/RXR and HIF-1α in the regulation of Fugu GLUT4 promoter activity by PPARγ and contractile activity, respectively. These data suggest that the characteristics and regulation of the GLUT4 promoter have been remarkably conserved during the evolution from fish to mammals, further evidencing the important role of GLUT4 in metabolic regulation in vertebrates.

Simian adenovirus type 35 has a recombinant genome comprising human and simian adenovirus sequences, which predicts its potential emergence as a human respiratory pathogen

Emergent human and simian adenoviruses (HAdVs) may arise from genome recombination. Computational analysis of SAdV type 35 reveals a genome comprising a chassis with elements mostly from two simian adenoviruses, SAdV-B21 and -B27, and regions of high sequence similarity shared with HAdV-B21 and HAdV-B16. Although recombination direction cannot be determined, the presence of these regions suggests prior infections of humans by an ancestor of SAdV-B35, and/or vice versa. Absence of this virus in humans may reflect non-optimal conditions for zoonosis or incomplete typing, e.g., limited epitope-based. The presence of both a critical viral replication element found in HAdV genomes and genes that are highly similar to ones in HAdVs suggest the potential to establish in a human host. This allows a prediction that this virus may be a nascent human respiratory pathogen. The recombination potential of human and simian adenovirus genomes should be considered in the use of SAdVs as vectors for gene delivery in humans.

Horizontally acquired divergent O-antigen contributes to escape from cross-immunity in the classical bordetellae

Background

Horizontal gene transfer (HGT) allows for rapid spread of genetic material between species, increasing genetic and phenotypic diversity. Although HGT contributes to adaptation and is widespread in many bacteria, others show little HGT. This study builds on previous work to analyze the evolutionary mechanisms contributing to variation within the locus encoding a prominent antigen of the classical bordetellae.

Results

We observed amongst classical bordetellae discrete regions of the lipopolysaccharide O-antigen locus with higher sequence diversity than the genome average. Regions of this locus had less than 50% sequence similarity, low dN/dS ratios and lower GC content compared to the genome average. Additionally, phylogenetic tree topologies based on genome-wide SNPs were incongruent with those based on genes within these variable regions, suggesting portions of the O-antigen locus may have been horizontally transferred. Furthermore, several predicted recombination breakpoints correspond with the ends of these variable regions. To examine the evolutionary forces that might have selected for this rare example of HGT in bordetellae, we compared in vitro and in vivo phenotypes associated with different O-antigen types. Antibodies against O1- and O2-serotypes were poorly cross-reactive, and did not efficiently kill or mediate clearance of alternative O-type bacteria, while a distinct and poorly immunogenic O-antigen offered no protection against colonization.

Conclusions

This study suggests that O-antigen variation was introduced to the classical bordetellae via HGT through recombination. Additionally, genetic variation may be maintained within the O-antigen locus because it can provide escape from immunity to different O-antigen types, potentially allowing for the circulation of different Bordetella strains within the same host population.

Rapid cohort generation and analysis of disease spectrum of large animal model of cone dystrophy

Large animal models are an important resource for the understanding of human disease and for evaluating the applicability of new therapies to human patients. For many diseases, such as cone dystrophy, research effort is hampered by the lack of such models. Lentiviral transgenesis is a methodology broadly applicable to animals from many different species. When conjugated to the expression of a dominant mutant protein, this technology offers an attractive approach to generate new large animal models in a heterogeneous background. We adopted this strategy to mimic the phenotype diversity encounter in humans and generate a cohort of pigs for cone dystrophy by expressing a dominant mutant allele of the guanylate cyclase 2D (GUCY2D) gene. Sixty percent of the piglets were transgenic, with mutant GUCY2D mRNA detected in the retina of all animals tested. Functional impairment of vision was observed among the transgenic pigs at 3 months of age, with a follow-up at 1 year indicating a subsequent slower progression of phenotype. Abnormal retina morphology, notably among the cone photoreceptor cell population, was observed exclusively amongst the transgenic animals. Of particular note, these transgenic animals were characterized by a range in the severity of the phenotype, reflecting the human clinical situation. We demonstrate that a transgenic approach using lentiviral vectors offers a powerful tool for large animal model development. Not only is the efficiency of transgenesis higher than conventional transgenic methodology but this technique also produces a heterogeneous cohort of transgenic animals that mimics the genetic variation encountered in human patients.

Establishing genomic tools and resources for Guizotia abyssinica (L.f.) Cass.-the development of a library of expressed sequence tags, microsatellite loci, and the sequencing of its chloroplast genome

We present an EST library, chloroplast genome sequence, and nuclear microsatellite markers that were developed for the semi-domesticated oilseed crop noug (Guizotia abyssinica) from Ethiopia. The EST library consists of 25 711 Sanger reads, assembled into 17 538 contigs and singletons, of which 4781 were functionally annotated using the Arabidopsis Information Resource (TAIR). The age distribution of duplicated genes in the EST library shows evidence of two paleopolyploidizations-a pattern that noug shares with several other species in the Heliantheae tribe (Compositae family). From the EST library, we selected 43 microsatellites and then designed and tested primers for their amplification. The number of microsatellite alleles varied between 2 and 10 (average 4.67), and the average observed and expected heterozygosities were 0.49 and 0.54, respectively. The chloroplast genome was sequenced de novo using Illumina's sequencing technology and completed with traditional Sanger sequencing. No large re-arrangements were found between the noug and sunflower chloroplast genomes, but 1.4% of sites have indels and 1.8% show sequence divergence between the two species. We identified 34 tRNAs, 4 rRNA sequences, and 80 coding sequences, including one region (trnH-psbA) with 15% sequence divergence between noug and sunflower that may be particularly useful for phylogeographic studies in noug and its wild relatives.

Computational analysis of four human adenovirus type 4 genomes reveals molecular evolution through two interspecies recombination events

Computational analysis of human adenovirus type 4 (HAdV-E4), a pathogen that is the only HAdV member of species E, provides insights into its zoonotic origin and molecular adaptation. Its genome encodes a domain of the major capsid protein, hexon, from HAdV-B16 recombined into the genome chassis of a simian adenovirus. Genomes of two recent field strains provide a clue to its adaptation to the new host: recombination of a NF-I binding site motif, which is required for efficient viral replication, from another HAdV genome. This motif is absent in the chimpanzee adenoviruses and the HAdV-E4 prototype, but is conserved amongst other HAdVs. This is the first report of an interspecies recombination event for HAdVs, and the first documentation of a lateral partial gene transfer from a chimpanzee AdV. The potential for such recombination events are important when considering chimpanzee adenoviruses as candidate gene delivery vectors for human patients.

Aryl hydrocarbon receptor negatively regulates expression of the plakoglobin gene (jup)

Plakoglobin is an important component of intercellular junctions, including both desmosomes and adherens junctions, which is known as a tumor suppressor. Although mutations in the plakoglobin gene (Jup) and/or changes in its protein levels have been observed in various disease states, including cancer progression or cardiovascular defects, the information about endogenous or exogenous stimuli orchestrating Jup expression is limited. Here we show that the aryl hydrocarbon receptor (AhR) may regulate Jup expression in a cell-specific manner. We observed a significant suppressive effect of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), a model toxic exogenous activator of the AhR signaling, on Jup expression in a variety of experimental models derived from rodent tissues, including contact-inhibited rat liver progenitor cells (where TCDD induces cell proliferation), rat and mouse hepatoma cell models (where TCDD inhibits cell cycle progression), cardiac cells derived from the mouse embryonic stem cells, or cardiomyocytes isolated from neonatal rat hearts. The small interfering RNA (siRNA)-mediated knockdown of AhR confirmed its role in both basal and TCDD-deregulated Jup expression. The analysis of genomic DNA located ~2.5kb upstream of rat Jup gene revealed a presence of evolutionarily conserved AhR binding motifs, which were confirmed upon their cloning into luciferase reporter construct. The siRNA-mediated knockdown of Jup expression affected both proliferation and attachment of liver progenitor cells. The present data indicate that the AhR may contribute to negative regulation of Jup gene expression in rodent cellular models, which may affect cell adherence and proliferation.

Comparative genomics identifies the mouse Bmp3 promoter and an upstream evolutionary conserved region (ECR) in mammals

The Bone Morphogenetic Protein (BMP) pathway is a multi-member signaling cascade whose basic components are found in all animals. One member, BMP3, which arose more recently in evolution and is found only in deuterostomes, serves a unique role as an antagonist to both the canonical BMP and Activin pathways. However, the mechanisms that control BMP3 expression, and the cis-regulatory regions mediating this regulation, remain poorly defined. With this in mind, we sought to identify the Bmp3 promoter in mouse (M. musculus) through functional and comparative genomic analyses. We found that the minimal promoter required for expression in resides within 0.8 kb upstream of Bmp3 in a region that is highly conserved with rat (R. norvegicus). We also found that an upstream region abutting the minimal promoter acts as a repressor of the minimal promoter in HEK293T cells and osteoblasts. Strikingly, a portion of this region is conserved among all available eutherian mammal genomes (47/47), but not in any non-eutherian animal (0/136). We also identified multiple conserved transcription factor binding sites in the Bmp3 upstream ECR, suggesting that this region may preserve common cis-regulatory elements that govern Bmp3 expression across eutherian mammals. Since dysregulation of BMP signaling appears to play a role in human health and disease, our findings may have application in the development of novel therapeutics aimed at modulating BMP signaling in humans.

Messenger RNA sequence rather than protein sequence determines the level of self-synthesis and antigen presentation of the EBV-encoded antigen, EBNA1

Unique purine-rich mRNA sequences embedded in the coding sequences of a distinct group of gammaherpesvirus maintenance proteins underlie the ability of the latently infected cell to minimize immune recognition. The Epstein-Barr virus nuclear antigen, EBNA1, a well characterized lymphocryptovirus maintenance protein has been shown to inhibit in cis antigen presentation, due in part to a large internal repeat domain encoding glycine and alanine residues (GAr) encoded by a purine-rich mRNA sequence. Recent studies have suggested that it is the purine-rich mRNA sequence of this repeat region rather than the encoded GAr polypeptide that directly inhibits EBNA1 self-synthesis and contributes to immune evasion. To test this hypothesis, we generated a series of EBNA1 internal repeat frameshift constructs and assessed their effects on cis-translation and endogenous antigen presentation. Diverse peptide sequences resulting from alternative repeat reading frames did not alleviate the translational inhibition characteristic of EBNA1 self-synthesis or the ensuing reduced surface presentation of EBNA1-specific peptide-MHC class I complexes. Human cells expressing the EBNA1 frameshift variants were also poorly recognized by antigen-specific T-cells. Furthermore, a comparative analysis of the mRNA sequences of the corresponding repeat regions of different viral maintenance homologues highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. Based on these combined observations, we propose that the cis-translational inhibitory effect of the EBNA1 internal repeat sequence operates mechanistically at the nucleotide level, potentially through RNA secondary structural elements, and is unlikely to be mediated through the GAr polypeptide. The demonstration that the EBNA1 repeat mRNA sequence and not the encoded protein sequence underlies immune evasion in this class of virus suggests a novel approach to therapeutic development through the use of anti-sense strategies or small molecules targeting EBNA1 mRNA structure.

Viruses establishing persistent latent infections have evolved various mechanisms to avoid immune surveillance. The Epstein-Barr virus-encoded nuclear antigen, EBNA1, expressed in all EBV-associated malignancies, modulates its own protein levels at quantities sufficient to maintain viral infection but low enough so as to minimize an immune response by the infected host cell. This evasion mechanism is regulated through an internal purine-rich mRNA repeat sequence encoding glycine and alanine residues. In this study we assess the impact of the repeat's nucleotide versus peptide sequence on inhibiting EBNA1 self-synthesis and antigen presentation. We demonstrate that altered peptide sequences resulting from frameshift mutations within the repeat do not alleviate the immune-evasive function of EBNA1, suggesting that the repetitive purine-rich mRNA sequence itself is responsible for inhibiting EBNA1 synthesis and subsequent poor immunogenicity. Our comparative analysis of the mRNA sequences of the corresponding repeat regions of different gammaherpesvirus maintenance homologues to EBNA1 highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. These studies demonstrate the importance of gammaherpesvirus purine-rich mRNA repeat sequences on antigenic epitope generation and evasion from T-cell mediated immune control, suggesting novel approaches to prevention and treatment of latent infection by this class of virus.

A misplaced lncRNA causes brachydactyly in humans

Translocations are chromosomal rearrangements that are frequently associated with a variety of disease states and developmental disorders. We identified 2 families with brachydactyly type E (BDE) resulting from different translocations affecting chromosome 12p. Both translocations caused downregulation of the parathyroid hormone-like hormone (PTHLH) gene by disrupting the cis-regulatory landscape. Using chromosome conformation capturing, we identified a regulator on chromosome 12q that interacts in cis with PTHLH over a 24.4-megabase distance and in trans with the sex-determining region Y-box 9 (SOX9) gene on chromosome 17q. The element also harbored a long noncoding RNA (lncRNA). Silencing of the lncRNA, PTHLH, or SOX9 revealed a feedback mechanism involving an expression-dependent network in humans. In the BDE patients, the human lncRNA was upregulated by the disrupted chromosomal association. Moreover, the lncRNA occupancy at the PTHLH locus was reduced. Our results document what we believe to be a novel in cis- and in trans-acting DNA and lncRNA regulatory feedback element that is reciprocally regulated by coding genes. Furthermore, our findings provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development.

Comparative genomics of the classical Bordetella subspecies: the evolution and exchange of virulence-associated diversity amongst closely related pathogens

Background

The classical Bordetella subspecies are phylogenetically closely related, yet differ in some of the most interesting and important characteristics of pathogens, such as host range, virulence and persistence. The compelling picture from previous comparisons of the three sequenced genomes was of genome degradation, with substantial loss of genome content (up to 24%) associated with adaptation to humans.

Results

For a more comprehensive picture of lineage evolution, we employed comparative genomic and phylogenomic analyses using seven additional diverse, newly sequenced Bordetella isolates. Genome-wide single nucleotide polymorphism (SNP) analysis supports a reevaluation of the phylogenetic relationships between the classical Bordetella subspecies, and suggests a closer link between ovine and human B. parapertussis lineages than has been previously proposed. Comparative analyses of genome content revealed that only 50% of the pan-genome is conserved in all strains, reflecting substantial diversity of genome content in these closely related pathogens that may relate to their different host ranges, virulence and persistence characteristics. Strikingly, these analyses suggest possible horizontal gene transfer (HGT) events in multiple loci encoding virulence factors, including O-antigen and pertussis toxin (Ptx). Segments of the pertussis toxin locus (ptx) and its secretion system locus (ptl) appear to have been acquired by the classical Bordetella subspecies and are divergent in different lineages, suggesting functional divergence in the classical Bordetellae.

Conclusions

Together, these observations, especially in key virulence factors, reveal that multiple mechanisms, such as point mutations, gain or loss of genes, as well as HGTs, contribute to the substantial phenotypic diversity of these versatile subspecies in various hosts.

A machine learning approach for identifying novel cell type-specific transcriptional regulators of myogenesis

Transcriptional enhancers integrate the contributions of multiple classes of transcription factors (TFs) to orchestrate the myriad spatio-temporal gene expression programs that occur during development. A molecular understanding of enhancers with similar activities requires the identification of both their unique and their shared sequence features. To address this problem, we combined phylogenetic profiling with a DNA-based enhancer sequence classifier that analyzes the TF binding sites (TFBSs) governing the transcription of a co-expressed gene set. We first assembled a small number of enhancers that are active in Drosophila melanogaster muscle founder cells (FCs) and other mesodermal cell types. Using phylogenetic profiling, we increased the number of enhancers by incorporating orthologous but divergent sequences from other Drosophila species. Functional assays revealed that the diverged enhancer orthologs were active in largely similar patterns as their D. melanogaster counterparts, although there was extensive evolutionary shuffling of known TFBSs. We then built and trained a classifier using this enhancer set and identified additional related enhancers based on the presence or absence of known and putative TFBSs. Predicted FC enhancers were over-represented in proximity to known FC genes; and many of the TFBSs learned by the classifier were found to be critical for enhancer activity, including POU homeodomain, Myb, Ets, Forkhead, and T-box motifs. Empirical testing also revealed that the T-box TF encoded by org-1 is a previously uncharacterized regulator of muscle cell identity. Finally, we found extensive diversity in the composition of TFBSs within known FC enhancers, suggesting that motif combinatorics plays an essential role in the cellular specificity exhibited by such enhancers. In summary, machine learning combined with evolutionary sequence analysis is useful for recognizing novel TFBSs and for facilitating the identification of cognate TFs that coordinate cell type-specific developmental gene expression patterns.

Ultra-barcoding in cacao (Theobroma spp.; Malvaceae) using whole chloroplast genomes and nuclear ribosomal DNA

PREMISE OF STUDY

To reliably identify lineages below the species level such as subspecies or varieties, we propose an extension to DNA-barcoding using next-generation sequencing to produce whole organellar genomes and substantial nuclear ribosomal sequence. Because this method uses much longer versions of the traditional DNA-barcoding loci in the plastid and ribosomal DNA, we call our approach ultra-barcoding (UBC).

METHODS

We used high-throughput next-generation sequencing to scan the genome and generate reliable sequence of high copy number regions. Using this method, we examined whole plastid genomes as well as nearly 6000 bases of nuclear ribosomal DNA sequences for nine genotypes of Theobroma cacao and an individual of the related species T. grandiflorum, as well as an additional publicly available whole plastid genome of T. cacao.

KEY RESULTS

All individuals of T. cacao examined were uniquely distinguished, and evidence of reticulation and gene flow was observed. Sequence variation was observed in some of the canonical barcoding regions between species, but other regions of the chloroplast were more variable both within species and between species, as were ribosomal spacers. Furthermore, no single region provides the level of data available using the complete plastid genome and rDNA.

CONCLUSIONS

Our data demonstrate that UBC is a viable, increasingly cost-effective approach for reliably distinguishing varieties and even individual genotypes of T. cacao. This approach shows great promise for applications where very closely related or interbreeding taxa must be distinguished.

Regulation of an antisense RNA with the transition of neonatal to IIb myosin heavy chain during postnatal development and hypothyroidism in rat skeletal muscle

Postnatal development of fast skeletal muscle is characterized by a transition in expression of myosin heavy chain (MHC) isoforms, from primarily neonatal MHC at birth to primarily IIb MHC in adults, in a tightly coordinated manner. These isoforms are encoded by distinct genes, which are separated by ∼17 kb on rat chromosome 10. The neonatal-to-IIb MHC transition is inhibited by a hypothyroid state. We examined RNA products [mRNA, pre-mRNA, and natural antisense transcript (NAT)] of developmental and adult-expressed MHC genes (embryonic, neonatal, I, IIa, IIx, and IIb) at 2, 10, 20, and 40 days after birth in normal and thyroid-deficient rat neonates treated with propylthiouracil. We found that a long noncoding antisense-oriented RNA transcript, termed bII NAT, is transcribed from a site within the IIb-Neo intergenic region and across most of the IIb MHC gene. NATs have previously been shown to mediate transcriptional repression of sense-oriented counterparts. The bII NAT is transcriptionally regulated during postnatal development and in response to hypothyroidism. Evidence for a regulatory mechanism is suggested by an inverse relationship between IIb MHC and bII NAT in normal and hypothyroid-treated muscle. Neonatal MHC transcription is coordinately expressed with bII NAT. A comparative phylogenetic analysis also suggests that bII NAT-mediated regulation has been a conserved trait of placental mammals for most of the eutherian evolutionary history. The evidence in support of the regulatory model implicates long noncoding antisense RNA as a mechanism to coordinate the transition between neonatal and IIb MHC during postnatal development.

Genome-wide identification of conserved regulatory function in diverged sequences

Plasticity of gene regulatory encryption can permit DNA sequence divergence without loss of function. Functional information is preserved through conservation of the composition of transcription factor binding sites (TFBS) in a regulatory element. We have developed a method that can accurately identify pairs of functional noncoding orthologs at evolutionarily diverged loci by searching for conserved TFBS arrangements. With an estimated 5% false-positive rate (FPR) in approximately 3000 human and zebrafish syntenic loci, we detected approximately 300 pairs of diverged elements that are likely to share common ancestry and have similar regulatory activity. By analyzing a pool of experimentally validated human enhancers, we demonstrated that 7/8 (88%) of their predicted functional orthologs retained in vivo regulatory control. Moreover, in 5/7 (71%) of assayed enhancer pairs, we observed concordant expression patterns. We argue that TFBS composition is often necessary to retain and sufficient to predict regulatory function in the absence of overt sequence conservation, revealing an entire class of functionally conserved, evolutionarily diverged regulatory elements that we term "covert."

Molecular evolution of the vertebrate TLR1 gene family--a complex history of gene duplication, gene conversion, positive selection and co-evolution

Background

The Toll-like receptors represent a large superfamily of type I transmembrane glycoproteins, some common to a wide range of species and others are more restricted in their distribution. Most members of the Toll-like receptor superfamily have few paralogues; the exception is the TLR1 gene family with four closely related genes in mammals TLR1, TLR2, TLR6 and TLR10, and four in birds TLR1A, TLR1B, TLR2A and TLR2B. These genes were previously thought to have arisen by a series of independent gene duplications. To understand the evolutionary pattern of the TLR1 gene family in vertebrates further, we cloned the sequences of TLR1A, TLR1B, TLR2A and TLR2B in duck and turkey, constructed phylogenetic trees, predicted codons under positive selection and identified co-evolutionary amino acid pairs within the TLR1 gene family using sequences from 4 birds, 28 mammals, an amphibian and a fish.

Results

This detailed phylogenetic analysis not only clarifies the gene gains and losses within the TLR1 gene family of birds and mammals, but also defines orthologues between these vertebrates. In mammals, we predict amino acid sites under positive selection in TLR1, TLR2 and TLR6 but not TLR10. We detect co-evolution between amino acid residues in TLR2 and the other members of this gene family predicted to maintain their ability to form functional heterodimers. In birds, we predict positive selection in the TLR2A and TLR2B genes at functionally significant amino acid residues. We demonstrate that the TLR1 gene family has mostly been subject to purifying selection but has also responded to directional selection at a few sites, possibly in response to pathogen challenge.

Conclusions

Our phylogenetic and structural analyses of the vertebrate TLR1 family have clarified their evolutionary origins and predict amino acid residues likely to be important in the host's defense against invading pathogens.