Impact of GC content on gene expression pattern in chicken

Genome-wide mapping of the binding sites of myocyte enhancer factor 2A in chicken primary myoblasts

Myocyte enhancer factor 2A (MEF2A) is a transcription factor that plays a critical role in cell proliferation, differentiation and apoptosis. In contrast to the wide characterization of its regulation mechanism in mammalian skeletal muscle, its role in chickens is limited. Especially, its wide target genes remain to be identified. Therefore, we utilized Cleavage Under Targets and Tagmentation (CUT&Tag) technology to reveal the genome-wide binding profile of MEF2A in chicken primary myoblasts thus gaining insights into its potential role in muscle development. Our results revealed that MEF2A binding sites were primarily distributed in intergenic and intronic regions. Within the promoter region, although only 8.87% of MEF2A binding sites were found, these binding sites were concentrated around the transcription start site (TSS). Following peak annotation, a total of 1903 genes were identified as potential targets of MEF2A. Gene Ontology (GO) enrichment analysis further revealed that MEF2A target genes may be involved in the regulation of embryonic development in multiple organ systems, including muscle development, gland development, and visual system development. Moreover, a comparison of the MEF2A target genes identified in chicken primary myoblasts with those in mouse C2C12 cells revealed 388 target genes are conserved across species, 1515 target genes are chicken specific. Among these conserved genes, ankyrin repeat and SOCS box containing 5 (ASB5), transmembrane protein 182 (TMEM182), myomesin 2 (MYOM2), leucyl and cystinyl aminopeptidase (LNPEP), actinin alpha 2 (ACTN2), sorbin and SH3 domain containing 1 (SORBS1), ankyrin 3 (ANK3), sarcoglycan delta (SGCD), and ORAI calcium release-activated calcium modulator 1 (ORAI1) exhibited consistent expression patterns with MEF2A during embryonic muscle development. Finally, TMEM182, as an important negative regulator of muscle development, has been validated to be regulated by MEF2A by dual-luciferase and quantitative real-time PCR (qPCR) assays. In summary, our study for the first time provides a wide landscape of MEF2A target genes in chicken primary myoblasts, which supports the active role of MEF2A in chicken muscle development.

The Chromosome-level Genome Provides Insights into the Evolution and Adaptation of Extreme Aggression

Extremely aggressive behavior, as the special pattern, is rare in most species and characteristic as contestants severely injured or killed ending the combat. Current studies of extreme aggression are mainly from the perspectives of behavioral ecology and evolution, while lacked the aspects of molecular evolutionary biology. Here, a high-quality chromosome-level genome of the parasitoid Anastatus disparis was provided, in which the males exhibit extreme mate-competition aggression. The integrated multiomics analysis highlighted that neurotransmitter dopamine overexpression, energy metabolism (especially from lipid), and antibacterial activity are likely major aspects of evolutionary formation and adaptation for extreme aggression in A. disparis. Conclusively, our study provided new perspectives for molecular evolutionary studies of extreme aggression as well as a valuable genomic resource in Hymenoptera.

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

The Complete Mitochondrial Genome and Phylogenetic Analyses of To Chicken in Vietnam

Indigenous chicken breeds have both cultural significance and economic value since they possess unique genetic characteristics that enable them to adapt to the local environment and contribute to biodiversity, food security, and sustainable agriculture in Vietnam. To (Tò in Vietnamese) chicken, a Vietnamese indigenous chicken breed, is popularly raised in Thai Binh province; however, little known is about the genetic diversity of this breed. In this study, we sequenced the complete mitochondrial genome of To chicken for a better understanding of the diversity and origin of the breed. The results of sequencing showed that the mitochondrial genome of To chicken spans a total length of 16,784 base pairs and comprises one non-coding control region (known as the displacement-loop (D-loop) region), two ribosomal RNA genes, 13 protein-coding genes, and 22 transfer RNA genes. The phylogenetic tree analyses and estimated genetic distances based on 31 complete mitochondrial genome sequences indicated that To chicken has a close genetic distance with the Laotian native chicken breed, Lv'erwu breed in China, and Nicobari black and Kadaknath breeds in India. The result of the current study might be important for conservation, breeding, and further genetic studies of To chicken.

qPrimerDB: A Powerful and User-Friendly Database for qPCR Primer Design

Real-time quantitative polymerase chain reaction (qPCR) is a powerful tool for analyzing and quantifying gene expression, and its primer design is the first and most important step. In order to improve the efficiency and effectiveness of primer design, we designed a database qPrimerDB, based on thermodynamics gene-specific for multispecies qPCR primers design.In this chapter, we explained the working principle of the database and detailed the step-by-step practical steps with examples. The valuable and time-saving qPrimerDB database is publicly accessible at http://biodb.swu.edu.cn/qprimerdb and will be routinely updated.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.

An interplay between compositional constraint and natural selection dictates the codon usage pattern among select Galliformes

Galliformes are believed to be the first avian order that started living in human association and became domesticated. Members of this order ranged from common to rare species. Next-generation sequencing has availed researchers with the whole genome sequences of five Galliformes; chicken, helmeted Guinea fowl, turkey, Japanese quail, and peafowl. Bioinformatic analysis based on codon usage, evolution, and species-specific functional enrichment can provide some crucial information aiding proper understanding of their genomic strategies. In this study, we investigated the genomic features of chicken, helmeted guinea fowl, turkey, and Japanese quail. Their genomes were AT biased although the potentially highly expressed genes contained more GC than AT. Cytosine dominated the third position of frequently used optimal codons. Mutational pressures in the analyzed Galliformes were in the range of 0.2-0.6%. Neutrality plot, translational selection index, and mutational responsive index indicated the dominance of selection pressure over mutational pressure among Galliformes. A pair of di-nucleotides, TpA and CpG, was found to be used less frequently than others in protein-coding genes since both of them are associated with the conversion of euchromatin to heterochromatin. Functional enrichment analysis revealed the dominance of proteins associated with fundamental biological processes. In turkey, chicken and helmeted Guinea fowl proteins with immunity-boosting capacity prevailed along with proteins needed for signal transduction and maintenance of central dogma. Evolutionary analysis indicated a bias towards synonymous substitution than non-synonymous mutation.

Refinement of Draft Genome Assemblies of Pigeonpea (Cajanus cajan)

Genome assembly of short reads from large plant genomes remains a challenge in computational biology despite major developments in next generation sequencing. Of late several draft assemblies have been reported in sequenced plant genomes. The reported draft genome assemblies of Cajanus cajan have different levels of genome completeness, a large number of repeats, gaps, and segmental duplications. Draft assemblies with portions of genome missing are shorter than the referenced original genome. These assemblies come with low map accuracy affecting further functional annotation and the prediction of gene components as desired by crop researchers. Genome coverage, i.e., the number of sequenced raw reads mapped onto a certain location of the genome is an important quality indicator of completeness and assembly quality in draft assemblies. The present work aimed to improve the coverage in reported de novo sequenced draft genomes (GCA_000340665.1 and GCA_000230855.2) of pigeonpea, a legume widely cultivated in India. The two recently sequenced assemblies, A1 and A2 comprised 72% and 75% of the estimated coverage of the genome, respectively. We employed an assembly reconciliation approach to compare the draft assemblies and merge them, filling the gaps by employing an algorithm size sorting mate-pair library to generate a high quality and near complete assembly with enhanced contiguity. The majority of gaps present within scaffolds were filled with right-sized mate-pair reads. The improved assembly reduced the number of gaps than those reported in draft assemblies resulting in an improved genome coverage of 82.4%. Map accuracy of the improved assembly was evaluated using various quality metrics and for the presence of specific trait-related functional genes. Employed pair-end and mate-pair local libraries helped us to reduce gaps, repeats, and other sequence errors resulting in lengthier scaffolds compared to the two draft assemblies. We reported the prediction of putative host resistance genes against Fusarium wilt disease by their performance and evaluated them both in wet laboratory and field phenotypic conditions.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Refinement of Draft Genome Assemblies of Pigeonpea (Cajanus cajan)

Genome assembly of short reads from large plant genomes remains a challenge in computational biology despite major developments in next generation sequencing. Of late several draft assemblies have been reported in sequenced plant genomes. The reported draft genome assemblies of Cajanus cajan have different levels of genome completeness, a large number of repeats, gaps, and segmental duplications. Draft assemblies with portions of genome missing are shorter than the referenced original genome. These assemblies come with low map accuracy affecting further functional annotation and the prediction of gene components as desired by crop researchers. Genome coverage, i.e., the number of sequenced raw reads mapped onto a certain location of the genome is an important quality indicator of completeness and assembly quality in draft assemblies. The present work aimed to improve the coverage in reported de novo sequenced draft genomes (GCA_000340665.1 and GCA_000230855.2) of pigeonpea, a legume widely cultivated in India. The two recently sequenced assemblies, A1 and A2 comprised 72% and 75% of the estimated coverage of the genome, respectively. We employed an assembly reconciliation approach to compare the draft assemblies and merge them, filling the gaps by employing an algorithm size sorting mate-pair library to generate a high quality and near complete assembly with enhanced contiguity. The majority of gaps present within scaffolds were filled with right-sized mate-pair reads. The improved assembly reduced the number of gaps than those reported in draft assemblies resulting in an improved genome coverage of 82.4%. Map accuracy of the improved assembly was evaluated using various quality metrics and for the presence of specific trait-related functional genes. Employed pair-end and mate-pair local libraries helped us to reduce gaps, repeats, and other sequence errors resulting in lengthier scaffolds compared to the two draft assemblies. We reported the prediction of putative host resistance genes against Fusarium wilt disease by their performance and evaluated them both in wet laboratory and field phenotypic conditions.

Genome-wide characterization and analysis of microsatellite sequences in camelid species

Microsatellites or simple sequence repeats (SSRs) are among the genetic markers most widely utilized in research. This includes applications in numerous fields such as genetic conservation, paternity testing, and molecular breeding. Though ordered draft genome assemblies of camels have been announced, including for the Arabian camel, systemic analysis of camel SSRs is still limited. The identification and development of informative and robust molecular SSR markers are essential for marker assisted breeding programs and paternity testing. Here we searched and compared perfect SSRs with 1–6 bp nucleotide motifs to characterize microsatellites for draft genome sequences of the Camelidae. We analyzed and compared the occurrence, relative abundance, relative density, and guanine-cytosine (GC) content in four taxonomically different camelid species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. A total of 546762, 544494, 547974, and 437815 SSRs were mined, respectively. Mononucleotide SSRs were the most frequent in the four genomes, followed in descending order by di-, tetra-, tri-, penta-, and hexanucleotide SSRs. GC content was highest in dinucleotide SSRs and lowest in mononucleotide SSRs. Our results provide further evidence that SSRs are more abundant in noncoding regions than in coding regions. Similar distributions of microsatellites were found in all four species, which indicates that the pattern of microsatellites is conserved in family Camelidae.

qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms

Real-time quantitative polymerase chain reaction (qPCR) is one of the most important methods for analyzing the expression patterns of target genes. However, successful qPCR experiments rely heavily on the use of high-quality primers. Various qPCR primer databases have been developed to address this issue, but these databases target only a few important organisms. Here, we developed the qPrimerDB database, founded on an automatic gene-specific qPCR primer design and thermodynamics-based validation workflow. The qPrimerDB database is the most comprehensive qPCR primer database available to date, with a web front-end providing gene-specific and pre-computed primer pairs across 147 important organisms, including human, mouse, zebrafish, yeast, thale cress, rice and maize. In this database, we provide 3331426 of the best primer pairs for each gene, based on primer pair coverage, as well as 47760359 alternative gene-specific primer pairs, which can be conveniently batch downloaded. The specificity and efficiency was validated for qPCR primer pairs for 66 randomly selected genes, in six different organisms, through qPCR assays and gel electrophoresis. The qPrimerDB database represents a valuable, timesaving resource for gene expression analysis. This resource, which will be routinely updated, is publically accessible at http://biodb.swu.edu.cn/qprimerdb.

Variation of gene expression in plants is influenced by gene architecture and structural properties of promoters

In higher eukaryotes, gene architecture and structural properties of promoters have emerged as significant factors influencing variation in number of transcripts (expression level) and specificity of gene expression in a tissue (expression breadth), which eventually shape the phenotype. In this study, transcriptome data of different tissue types at various developmental stages of A. thaliana, O. sativa, S. bicolor and Z. mays have been used to understand the relationship between properties of gene components and its expression. Our findings indicate that in plants, among all gene architecture and structural properties of promoters, compactness of genes in terms of intron content is significantly linked to gene expression level and breadth, whereas in human an exactly opposite scenario is seen. In plants, for the first time we have carried out a quantitative estimation of effect of a particular trait on expression level and breadth, by using multiple regression analysis and it confirms that intron content of primary transcript (as %) is a powerful determinant of expression breadth. Similarly, further regression analysis revealed that among structural properties of the promoters, stability is negatively linked to expression breadth, while DNase1 sensitivity strongly governs gene expression breadth in monocots and gene expression level in dicots. In addition, promoter regions of tissue specific genes are found to be enriched with TATA box and Y-patch motifs. Finally, multi copy orthologous genes in plants are found to be longer, highly regulated and tissue specific.

Identification of recombination hotspots and quantitative trait loci for recombination rate in layer chickens

Background

The frequency of recombination events varies across the genome and between individuals, which may be related to some genomic features. The objective of this study was to assess the frequency of recombination events and to identify QTL (quantitative trait loci) for recombination rate in two purebred layer chicken lines.

Methods

A total of 1200 white-egg layers (WL) were genotyped with 580 K SNPs and 5108 brown-egg layers (BL) were genotyped with 42 K SNPs (single nucleotide polymorphisms). Recombination events were identified within half-sib families and both the number of recombination events and the recombination rate was calculated within each 0.5 Mb window of the genome. The 10% of windows with the highest recombination rate on each chromosome were considered to be recombination hotspots. A BayesB model was used separately for each line to identify genomic regions associated with the genome-wide number of recombination event per meiosis. Regions that explained more than 0.8% of genetic variance of recombination rate were considered to harbor QTL.

Results

Heritability of recombination rate was estimated at 0.17 in WL and 0.16 in BL. On average, 11.3 and 23.2 recombination events were detected per individual across the genome in 1301 and 9292 meioses in the WL and BL, respectively. The estimated recombination rates differed significantly between the lines, which could be due to differences in inbreeding levels, and haplotype structures. Dams had about 5% to 20% higher recombination rates per meiosis than sires in both lines. Recombination rate per 0.5 Mb window had a strong negative correlation with chromosome size and a strong positive correlation with GC content and with CpG island density across the genome in both lines. Different QTL for recombination rate were identified in the two lines. There were 190 and 199 non-overlapping recombination hotspots detected in WL and BL respectively, 28 of which were common to both lines.

Conclusions

Differences in the recombination rates, hotspot locations, and QTL regions associated with genome-wide recombination were observed between lines, indicating the breed-specific feature of detected recombination events and the control of recombination events is a complex polygenic trait.

Electronic supplementary material

The online version of this article (10.1186/s40104-019-0332-y) contains supplementary material, which is available to authorized users.

Tail Wags the Dog? Functional Gene Classes Driving Genome-Wide GC Content in Plasmodium spp

Plasmodium parasites are valuable models to understand how nucleotide composition affects mutation, diversification, and adaptation. No other observed eukaryotes have undergone such large changes in genomic Guanine-Cytosine (GC) content as seen in the genus Plasmodium (∼30% within 35-40 Myr). Although mutational biases are known to influence GC content in the human-infective Plasmodium vivax and Plasmodium falciparum; no study has addressed how different gene functional classes contribute to genus-wide compositional changes, or if Plasmodium GC content variation is driven by natural selection. Here, we tested the hypothesis that certain gene processes and functions drive variation in global GC content between Plasmodium species. We performed a large-scale comparative genomic analysis using the genomes and predicted genes of 17 Plasmodium species encompassing a wide genomic GC content range. Genic GC content was sorted and divided into ten equally sized quantiles that were then assessed for functional enrichment classes. In agreement that selection on gene classes may drive genomic GC content, trans-membrane proteins were enriched within extreme GC content quantiles (Q1 and Q10). Specifically, variant surface antigens, which primarily interact with vertebrate immune systems, showed skewed GC content distributions compared with other trans-membrane proteins. Although a definitive causation linking GC content, expression, and positive selection within variant surface antigens from Plasmodium vivax, Plasmodium berghei, and Plasmodium falciparum could not be established, we found that regardless of genomic nucleotide composition, genic GC content and expression were positively correlated during trophozoite stages. Overall, these data suggest that, alongside mutational biases, functional protein classes drive Plasmodium GC content change.

Expression and purification of biologically active recombinant rabbit monocyte chemoattractant protein1 in Escherichia coli

Monocyte chemoattractant protein 1 (MCP1) with recruiting monocytes is an important factor at the beginning of inflammatory disorders such as atherosclerosis which seems its blocking preclude this process and help improvement of related diseases. To perform clinical research in this field, MCP1 protein is required but firstly, animal studies should be done. As the rabbit is a suitable model for many inflammatory disorders, and Escherichia coli BL21(DE3) (BL21) cell is a high-efficiency host for protein expression, we decided to produce recombinant rabbit MCP1 (rRMCP1) in BL21/pET28a system. After codon usage, a construct containing RMCP1 sequence was synthesized, cloned into the pET28a plasmid, and overexpressed in BL21 cells. Followed that, with changing expression condition such as cell concentration before the induction, time period, temperature, shaking rate and inducer concentration (IPTG), rRMCP1 expression was optimized, and purified by Ni-NTA. The biological activity of the expressed protein was verified using monocyte migration assay. Using this expression system, nearly 28 mg/mL rRMCP1 was produced at 26°C/180 rpm for 24 h in LB broth medium with 1 mM IPTG. Therefore, we were succeeded to express the intermediate level of rRMCP1 with this method. This amount of protein is sufficient for biological researches in the laboratory.

Chicken (Gallus gallus) endogenous retrovirus generates genomic variations in the chicken genome

Background

Transposable elements (TEs) comprise ~10% of the chicken (Gallus gallus) genome. The content of TEs is much lower than that of mammalian genomes, where TEs comprise around half of the genome. Endogenous retroviruses are responsible for ~1.3% of the chicken genome. Among them is Gallus gallus endogenous retrovirus 10 (GGERV10), one of the youngest endogenous retrovirus families, which emerged in the chicken genome around 3 million years ago.

Results

We identified a total of 593 GGERV10 elements in the chicken reference genome using UCSC genome database and RepeatMasker. While most of the elements were truncated, 49 GGERV10 elements were full-length retaining 5′ and 3′ LTRs. We examined in detail their structural features, chromosomal distribution, genomic environment, and phylogenetic relationships. We compared LTR sequence among five different GGERV10 subfamilies and found sequence variations among the LTRs. Using a traditional PCR assay, we examined a polymorphism rate of the 49 full-length GGERV10 elements in three different chicken populations of the Korean domestic chicken, Leghorn, and Araucana. The result found a breed-specific GGERV10B insertion locus in the Korean domestic chicken, which could be used as a Korean domestic chicken-specific marker.

Conclusions

GGERV10 family is the youngest ERV family and thus might have contributed to recent genomic variations in different chicken populations. The result of this study showed that one of GGERV10 elements integrated into the chicken genome after the divergence of Korean domestic chicken from other closely related chicken populations, suggesting that GGERV10 could be served as a molecular marker for chicken breed identification.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0085-5) contains supplementary material, which is available to authorized users.

A comparative study on the regulatory region of the PERIOD1 gene among diurnal/nocturnal primates

BACKGROUND

The circadian clock is set up around a 24-h period in humans who are awake in the daytime and sleep in the nighttime, accompanied with physiological and metabolic rhythms. Most haplorhine primates, including humans, are diurnal, while most "primitive" strepsirrhine primates are nocturnal, suggesting primates have evolved from nocturnal to diurnal habits. The mechanisms of physiological changes causing the habits and of genetic changes causing the physiological changes are, however, unknown. To reveal these mechanisms, we focus on the nucleotide sequences of the regulatory region of the PERIOD1 (PER1) gene that is known as one of the key elements of the circadian clock in mammalians.

METHODS

We determined nucleotide sequences of the regulatory region of PER1 concerning the gene expression for six primates and compared those with those of eight primates from the international DNA database. Based on the sequence data, we constructed a phylogenetic tree including both the diurnal/nocturnal species and investigated the guanine and cytosine (GC) content in the regulatory region.

RESULTS

The motif sequences regulating gene expression were evolutionary conservative in the primates examined. The phylogenetic tree simply showed phylogenetic relationship among the species and no branching pattern distinguishable between the diurnal and nocturnal groups. We found two cores showing a statistically significant difference between the diurnal and the nocturnal habits related to the GC contents of the regulatory region of PER1.

CONCLUSION

Our results suggest the possibility that the two cores in the upstream region of PER1 are related to the regulation of gene expression leading to behavioral differences between diurnal and nocturnal primates.

Codon optimization of antigen coding sequences improves the immune potential of DNA vaccines against avian influenza virus H5N1 in mice and chickens

Background

Highly pathogenic avian influenza viruses are a serious threat to domestic poultry and can be a source of new human pandemic and annual influenza strains. Vaccination is the main strategy of protection against influenza, thus new generation vaccines, including DNA vaccines, are needed. One promising approach for enhancing the immunogenicity of a DNA vaccine is to maximize its expression in the immunized host.

Methods

The immunogenicity of three variants of a DNA vaccine encoding hemagglutinin (HA) from the avian influenza virus A/swan/Poland/305-135V08/2006 (H5N1) was compared in two animal models, mice (BALB/c) and chickens (broilers and layers). One variant encoded the wild type HA while the other two encoded HA without proteolytic site between HA1 and HA2 subunits and differed in usage of synonymous codons. One of them was enriched for codons preferentially used in chicken genes, while in the other modified variant the third position of codons was occupied in almost 100 % by G or C nucleotides.

Results

The variant of the DNA vaccine containing almost 100 % of the GC content in the third position of codons stimulated strongest immune response in two animal models, mice and chickens. These results indicate that such modification can improve not only gene expression but also immunogenicity of DNA vaccine.

Conclusion

Enhancement of the GC content in the third position of the codon might be a good strategy for development of a variant of a DNA vaccine against influenza that could be highly effective in distant hosts, such as birds and mammals, including humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0599-y) contains supplementary material, which is available to authorized users.

Genomics of Entomopathogenic Nematodes and Implications for Pest Control

Entomopathogenic nematodes (EPNs) have been used in biological control but improvement is needed to realize their full potential for broader application in agriculture. Some improvements have been gained through selective breeding and the isolation of additional species and populations. Having genomic sequences for at least six EPNs opens the possibility of genetic improvement, either by facilitating the selection of candidate genes for hypothesis-driven studies of gene-trait relations or by genomics-assisted breeding for desirable traits. However, the genomic data will be of limited use without a more mechanistic understanding of the genes underlying traits that are important for biological control. Additionally, molecular tools are required to fully translate the genomic resources into further functional studies and better biological control.

Characteristics of Antisense Transcript Promoters and the Regulation of Their Activity

Recently, an increasing number of studies on natural antisense transcripts have been reported, especially regarding their classification, temporal and spatial expression patterns, regulatory functions and mechanisms. It is well established that natural antisense transcripts are produced from the strand opposite to the strand encoding a protein. Despite the pivotal roles of natural antisense transcripts in regulating the expression of target genes, the transcriptional mechanisms initiated by antisense promoters (ASPs) remain unknown. To date, nearly all of the studies conducted on this topic have focused on the ASP of a single gene of interest, whereas no study has systematically analyzed the locations of ASPs in the genome, ASP activity, or factors influencing this activity. This review focuses on elaborating on and summarizing the characteristics of ASPs to extend our knowledge about the mechanisms of antisense transcript initiation.

Analysis of 5' gene regions reveals extraordinary conservation of novel non-coding sequences in a wide range of animals

Background

Phylogenetic footprinting is a comparative method based on the principle that functional sequence elements will acquire fewer mutations over time than non-functional sequences. Successful comparisons of distantly related species will thus yield highly important sequence elements likely to serve fundamental biological roles. RNA regulatory elements are less well understood than those in DNA. In this study we use the emerging model organism Nasonia vitripennis, a parasitic wasp, in a comparative analysis against 12 insect genomes to identify deeply conserved non-coding elements (CNEs) conserved in large groups of insects, with a focus on 5’ UTRs and promoter sequences.

Results

We report the identification of 322 CNEs conserved across a broad range of insect orders. The identified regions are associated with regulatory and developmental genes, and contain short footprints revealing aspects of their likely function in translational regulation. The most ancient regions identified in our analysis were all found to overlap transcribed regions of genes, reflecting stronger conservation of translational regulatory elements than transcriptional elements. Further expanding sequence analyses to non-insect species we also report the discovery of, to our knowledge, the two oldest and most ubiquitous CNE’s yet described in the animal kingdom (700 MYA). These ancient conserved non-coding elements are associated with the two ribosomal stalk genes, RPLP1 and RPLP2, and were very likely functional in some of the earliest animals.

Conclusions

We report the identification of the most deeply conserved CNE’s found to date, and several other deeply conserved elements which are without exception, part of 5’ untranslated regions of transcripts, and occur in a number of key translational regulatory genes, highlighting translational regulation of translational regulators as a conserved feature of insect genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0499-6) contains supplementary material, which is available to authorized users.

Adult porcine genome-wide DNA methylation patterns support pigs as a biomedical model

Background

Pigs (Sus scrofa) provide relevant biomedical models to dissect complex diseases due to their anatomical, genetic, and physiological similarities with humans. Aberrant DNA methylation has been linked to many of these diseases and is associated with gene expression; however, the functional similarities and differences between porcine and human DNA methylation patterns are largely unknown.

Methods

DNA and RNA was isolated from eight tissue samples (fat, heart, kidney, liver, lung, lymph node, muscle, and spleen) from the adult female Duroc utilized for the pig genome sequencing project. Reduced representation bisulfite sequencing (RRBS) and RNA-seq were performed on an Illumina HiSeq2000. RRBS reads were aligned using BSseeker2, and only sites with a minimum depth of 10 reads were used for methylation analysis. RNA-seq reads were aligned using Tophat, and expression analysis was performed using Cufflinks. In addition, SNP calling was performed using GATK for targeted control and whole genome sequencing reads for CpG site validation and allelic expression analysis, respectively.

Results

Analysis on the influence of DNA variation in methylation calling revealed a reduced effectiveness of WGS datasets in covering CpG rich regions, as well as the usefulness of a targeted control library for SNP detection. Analysis of over 500,000 CpG sites demonstrated genome wide methylation patterns similar to those observed in humans, including reduced methylation within CpG islands and at transcription start sites (TSS), X chromosome inactivation, and anticorrelation of TSS CpG methylation with gene expression. In addition, a positive correlation between TSS CpG density and expression, and a negative correlation between TSS TpG density and expression were demonstrated. Low but non-random non-CpG methylation (<1%) was also detected in all non-neuronal somatic tissues, with differences in tissue clustering observed based on CpG and non-CpG methylation patterns. Finally, allele specific expression analysis revealed enrichment of genes involved in metabolic and regulatory processes. 

Discussion

These results provide transcriptional and DNA methylation datasets for the biomedical community that are directly relatable to current genomic resources. In addition, the correlation between TSS CpG density and expression suggests increased mutation rates at CpG sites play a significant role in adaptive evolution by reducing CpG density at TSS over time, resulting in higher methylation levels in these regions and more permanent changes to lower gene expression. This is proposed to occur predominantly through deamination of 5-methylcytosine to thymidine, resulting in the replacement of CpG with TpG sites in these regions, as indicated by the increased TSS TpG density observed in non-expressed genes, resulting in a negative correlation between expression and TSS TpG density.

Conclusions

This study provides baseline methylation and gene transcription profiles for a healthy adult pig, reports similar patterns to those observed in humans, and supports future porcine studies related to human disease and development. Additionally, the observed reduced CpG and increased TpG density at TSS of lowly expressed genes suggests DNA methylation plays a significant role in adaptive evolution through more permanent changes to lower gene expression.

Electronic supplementary material

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

Recurrent positive selection and heterogeneous codon usage bias events leading to coexistence of divergent pigeon circoviruses

The capsid genes from 14 pigeon circovirus (PiCV) sequences, collected from Taiwan between 2009 and 2010, were sequenced and compared with 14 PiCV capsid gene sequences from GenBank. Based on pairwise comparison, PiCV strains from Taiwan shared 73.9-100% nucleotide identity and 72-100% amino acid identity with those of the 14 reported PiCV sequences. Phylogenetic analyses revealed that Taiwanese PiCV isolates can be grouped into two clades: clade 1 comprising isolates from Belgium, Australia, USA, Italy and China, and clade 2 showing close relation to isolates from Germany and France. Recurrent positive selection was detected in clade 1 PiCV lineages, which may contribute to the diversification of predominant PiCV sequences in Taiwan. Further observations suggest that synonymous codon usage variations between PiCV clade 1 and clade 2 may reflect the adaptive divergence on translation efficiency of capsid genes in infectious hosts. Variation in selective pressures acting on the evolutionary divergence and codon usage bias of both clades explains the regional coexistence of virus sequences congeners prevented from competitive exclusion within an island such as Taiwan. Our genotyping results also provide insight into the aetiological agents of PiCV outbreak in Taiwan and we present a comparative analysis of the central coding region of PiCV genome. From the sequence comparison results of 28 PiCVs which differs in regard to the geographical origin and columbid species, we identified conserved regions within the capsid gene that are likely to be suitable for primer selection and vaccine development.

Molecular cloning and functional characterisation of chicken Atonal homologue 1: a comparison with human Atoh1

BACKGROUND INFORMATION

The vertebrate basic helix-loop-helix transcription factor Atoh1 is essential for maturation and survival of mechanosensory hair cells of the inner ear, neurogenesis, differentiation of the intestine, homeostasis of the colon and is implicated in cancer progression. Given that mutations in Atoh1 are detected in malignant tumours, study of functionally different Atoh1 alleles and homologues might yield useful avenues for investigation. The predicted sequence of chicken Atoh1 (cAtoh1) has large regions of dissimilarity to that of mammalian Atoh1 homologues. We hypothesise that cAtoh1 might have intrinsic functional differences to mammalian Atoh1.

RESULTS

In this study, we cloned and sequenced the full open reading frame of cAtoh1. In overexpression experiments, we show that this sequence is sufficient to generate a cAtoh1 protein capable of inducing hair cell markers when expressed in nonsensory regions of the developing inner ear, and that morpholino-mediated knock-down using a section of the sequence 5' to the start codon inhibits differentiation of hair cells in the chicken basilar papilla. Furthermore, we compare the behaviour of cAtoh1 and human Atoh1 (hAtoh1) in embryonic mouse cochlear explants, showing that cAtoh1 is a potent inducer of hair cell differentiation and that it can overcome Sox2-mediated repression of hair cell differentiation more effectively than hAtoh1.

CONCLUSIONS

cAtoh1 is both necessary and sufficient for avian mechanosensory hair cell differentiation. The non-conserved regions of the cAtoh1 coding region have functional consequences on its behaviour.

Hydrophobicity and aromaticity are primary factors shaping variation in amino acid usage of chicken proteome

Amino acids are utilized with different frequencies both among species and among genes within the same genome. Up to date, no study on the amino acid usage pattern of chicken has been performed. In the present study, we carried out a systematic examination of the amino acid usage in the chicken proteome. Our data indicated that the relative amino acid usage is positively correlated with the tRNA gene copy number. GC contents, including GC1, GC2, GC3, GC content of CDS and GC content of the introns, were correlated with the most of the amino acid usage, especially for GC rich and GC poor amino acids, however, multiple linear regression analyses indicated that only approximately 10–40% variation of amino acid usage can be explained by GC content for GC rich and GC poor amino acids. For other intermediate GC content amino acids, only approximately 10% variation can be explained. Correspondence analyses demonstrated that the main factors responsible for the variation of amino acid usage in chicken are hydrophobicity, aromaticity and genomic GC content. Gene expression level also influenced the amino acid usage significantly. We argued that the amino acid usage of chicken proteome likely reflects a balance or near balance between the action of selection, mutation, and genetic drift.

Abundance, arrangement, and function of sequence motifs in the chicken promoters

Background

Eukaryotic promoters are regions containing various sequence motifs necessary to control gene transcription. Much evidence has emerged showing that structural and/or contextual changes in regulatory elements can critically affect cis-regulatory activity. As sequence motifs can be key factors in maintaining complex promoter architectures, one effective approach to further understand the evolution of promoter regions in vertebrates is to compare the abundance and distribution patterns of sequence motifs in these regions between divergent species. When compared with mammals, the chicken (Gallus gallus) has a very different genome composition and sufficient genomic information to make it a good model for the exploration of promoter structure and evolution.

Results

More than 10% of chicken genes contained short tandem repeat (STR) in the region 2 kb upstream of promoters, but the total number of STRs observed in chicken is approximately half of that detected in human promoters. In terms of the STR motif frequencies, chicken promoter regions were more similar to other avian and mammalian promoters than these were to the entire chicken genome. Unlike other STRs, nearly half of the trinucleotide repeats found in promoters partly or entirely overlapped with CpG islands, indicating potential association with nucleosome positions. Moreover, the chicken promoters are abundant with sequence motifs such as poly-A, poly-G and G-quadruplexes, especially in the core region, that are otherwise rare in the genome. Most of sequence motifs showed strong functional enrichment for particular gene ontology (GO) categories, indicating roles in regulation of transcription and gene expression, as well as immune response and cognition.

Conclusions

Chicken promoter regions share some, but not all, of the structural features observed in mammalian promoters. The findings presented here provide empirical evidence suggesting that the frequencies and locations of STR motifs have been conserved through promoter evolution in a lineage-specific manner. Correlation analysis between GO categories and sequence motifs suggests motif-specific constraints acting on gene function.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-900) contains supplementary material, which is available to authorized users.

Epigenetic structure and the role of polymorphism in the shaping of DNA methylation patterns of equine OAS1 locus

DNA methylation patterns and their relation with genetic polymorphisms were determined in the equine OAS1 locus. Genetic variants of OAS1 were previously found to be associated with susceptibility to West Nile virus infections in horses. The subject of the study were white blood cells of 13 juvenile and 13 old horses from AA and HC breed and a set of solid tissues from a single adult horse. The aim was to determine the degree of variation of CpG methylation profiles with concern for tissue type, horse breed and age. Results of direct BSPCR and cloned BSPCR sequencing revealed that all of determined CpG islands (CGIs) were hypermethylated in exception to CGI covering OAS1 promoter and exon 1. One of intragenic CGIs displayed variability of methylation patterns across eight tissue types. The variability of particular sub-types of white blood cells between AA and HC horses were considered as the possible cause of interbreed differences of methylation levels. Comparison of sequence variability between converted and unconverted DNAs of both horse breeds showed polymorphisms of CpG sites to be the source of monoallelic methylation in exception to the polymorphic CpGs located in the OAS1 promoter. Two of them are new polymorphic variants in the OAS1 promoter region. Application of methylation data in conjunction with genetic variation detected at the OAS1 locus might be useful to deepen the knowledge about mechanisms underlying immunity to viral infections in the horse.

Length and GC content variability of introns among teleostean genomes in the light of the metabolic rate hypothesis

A comparative analysis of five teleostean genomes, namely zebrafish, medaka, three-spine stickleback, fugu and pufferfish was performed with the aim to highlight the nature of the forces driving both length and base composition of introns (i.e., bpi and GCi). An inter-genome approach using orthologous intronic sequences was carried out, analyzing independently both variables in pairwise comparisons. An average length shortening of introns was observed at increasing average GCi values. The result was not affected by masking transposable and repetitive elements harbored in the intronic sequences. The routine metabolic rate (mass specific temperature-corrected using the Boltzmann's factor) was measured for each species. A significant correlation held between average differences of metabolic rate, length and GC content, while environmental temperature of fish habitat was not correlated with bpi and GCi. Analyzing the concomitant effect of both variables, i.e., bpi and GCi, at increasing genomic GC content, a decrease of bpi and an increase of GCi was observed for the significant majority of the intronic sequences (from ∼40% to ∼90%, in each pairwise comparison). The opposite event, concomitant increase of bpi and decrease of GCi, was counter selected (from <1% to ∼10%, in each pairwise comparison). The results further support the hypothesis that the metabolic rate plays a key role in shaping genome architecture and evolution of vertebrate genomes.

Identification, characterization, and utilization of single copy genes in 29 angiosperm genomes

Background

Single copy genes are common across angiosperm genomes. With the sufficiently high quality sequenced genomes, the identification of large-scale single copy genes among multiple species is possible. Although some characteristics have been reported, our study provides novel insights into single copy genes.

Results

We identified single copy genes across 29 angiosperm genomes. A significant negative correlation was found between the number of duplicate blocks and the number of single copy genes. We found that a considerable number of single copy genes are located in organelles, showing a preference for binding and catalytic activity. The analysis of effective number of codons (Nc) illustrates that single copy genes have a stronger codon bias than non-single copy genes in eudicots. The relative high expression level of single copy genes was partially confirmed by the RNA-seq data, rather than the Codon Adaptation Index (CAI). Unlike in most other species, a strongly negatively correlation occurs between Nc and GC3 among single copy genes in grass genomes. When compared to all non-single copy genes, single copy genes indicate more conservation (as indicated by Ka and Ks values). But our alternative splicing (AS) results reveal that selective constraints are weaker in single copy genes than in low copy family genes (1–10 in-paralogs) and stronger than high copy family genes (>10 in-paralogs). Using concatenated shared single copy genes, we obtained a well-resolved phylogenetic tree. With the addition of intron sequences, the branch support is improved, but striking incongruences are also evident. Therefore, it is noteworthy that inclusion of intron sequences seems more appropriate for the phylogenetic reconstruction at lower taxonomic levels.

Conclusions

Our analysis provides insight into the evolutionary characteristics of single copy genes across 29 angiosperm genomes. The results suggest that there are key differences in evolutionary constraints between single copy genes and non-single copy genes. And to some extent, these evolutionary constraints show some species-specific differences, especially between eudicots and monocots. Our preliminary evidence also suggests that the concatenated shared single copy genes are well suited for use in resolving phylogenetic relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-504) contains supplementary material, which is available to authorized users.

Molecular evolution of NASP and conserved histone H3/H4 transport pathway

BACKGROUND

NASP is an essential protein in mammals that functions in histone transport pathways and maintenance of a soluble reservoir of histones H3/H4. NASP has been studied exclusively in Opisthokonta lineages where some functional diversity has been reported. In humans, growing evidence implicates NASP miss-regulation in the development of a variety of cancers. Although a comprehensive phylogenetic analysis is lacking, NASP-family proteins that possess four TPR motifs are thought to be widely distributed across eukaryotes.

RESULTS

We characterize the molecular evolution of NASP by systematically identifying putative NASP orthologs across diverse eukaryotic lineages ranging from excavata to those of the crown group. We detect extensive silent divergence at the nucleotide level suggesting the presence of strong purifying selection acting at the protein level. We also observe a selection bias for high frequencies of acidic residues which we hypothesize is a consequence of their critical function(s), further indicating the role of functional constraints operating on NASP evolution. Our data indicate that TPR1 and TPR4 constitute the most rapidly evolving functional units of NASP and may account for the functional diversity observed among well characterized family members. We also show that NASP paralogs in ray-finned fish have different genomic environments with clear differences in their GC content and have undergone significant changes at the protein level suggesting functional diversification.

CONCLUSION

We draw four main conclusions from this study. First, wide distribution of NASP throughout eukaryotes suggests that it was likely present in the last eukaryotic common ancestor (LECA) possibly as an important innovation in the transport of H3/H4. Second, strong purifying selection operating at the protein level has influenced the nucleotide composition of NASP genes. Further, we show that selection has acted to maintain a high frequency of functionally relevant acidic amino acids in the region that interrupts TPR2. Third, functional diversity reported among several well characterized NASP family members can be explained in terms of quickly evolving TPR1 and TPR4 motifs. Fourth, NASP fish specific paralogs have significantly diverged at the protein level with NASP2 acquiring a NNR domain.

Genome-enabled prediction of quantitative traits in chickens using genomic annotation

BACKGROUND

Genome-wide association studies have been deemed successful for identifying statistically associated genetic variants of large effects on complex traits. Past studies have found enrichment of trait-associated SNPs in functionally annotated regions, while depletion was reported for intergenic regions (IGR). However, no systematic examination of connections between genomic regions and predictive ability of complex phenotypes has been carried out.

RESULTS

In this study, we partitioned SNPs based on their annotation to characterize genomic regions that deliver low and high predictive power for three broiler traits in chickens using a whole-genome approach. Additive genomic relationship kernels were constructed for each of the genic regions considered, and a kernel-based Bayesian ridge regression was employed as prediction machine. We found that the predictive performance for ultrasound area of breast meat from using genic regions marked by SNPs was consistently better than that from SNPs in IGR, while IGR tagged by SNPs were better than the genic regions for body weight and hen house egg production. We also noted that predictive ability delivered by the whole battery of markers was close to the best prediction achieved by one of the genomic regions.

CONCLUSIONS

Whole-genome regression methods use all available quality filtered SNPs into a model, contrary to accommodating only validated SNPs from exonic or coding regions. Our results suggest that, while differences among genomic regions in terms of predictive ability were observed, the whole-genome approach remains as a promising tool if interest is on prediction of complex traits.