Conformational Changes in Three-Dimensional Chromatin Structure in Paulownia fortunei After Phytoplasma Infection.

Higher-order chromatin structures play important roles in regulating multiple biological processes such as growth and development as well as biotic and abiotic stress response. However, little is known about three-dimensional chromatin structures in Paulownia or about whole-genome chromatin conformational changes that occur in response to Paulownia witches' broom (PaWB) disease. We used high-throughput chromosome conformation capture (Hi-C) to obtain genome-wide profiles of chromatin conformation in both healthy and phytoplasma-infected Paulownia fortunei genome. The heat map results indicated that the strongest interactions between chromosomes were in the telomeres. We confirmed that the main structural characteristics of A/B compartments, topologically associated domains, and chromatin loops were prominent in the Paulownia genome and were clearly altered in phytoplasma-infected plants. By combining chromatin immunoprecipitation sequencing, Hi-C signals, and RNA sequencing data, we inferred that the chromatin structure changed and the modification levels of three histones (H3K4me3/K9ac/K36me3) increased in phytoplasma-infected P. fortunei, which was associated with changes of transcriptional activity. We concluded that for epigenetic modifications, transcriptional activity might function in combination to shape chromatin packing in healthy and phytoplasm-infected Paulownia. Finally, 11 genes (e.g., RPN6, Sec61 subunit-α) that were commonly located at specific topologically associated domain boundaries, A/B compartment switching and specific loops, and had been associated with histone marks were identified and considered as closely related to PaWB stress. Our results provide new insights into the nexus between gene regulation and chromatin conformational alterations in nonmodel plants upon phytopathogen infection and plant disease resistance.

Population genomics identifies patterns of genetic diversity and selection in chicken.

BACKGROUND: There are hundreds of phenotypically distinguishable domestic chicken breeds or lines with highly specialized traits worldwide, which provide a unique opportunity to illustrate how selection shapes patterns of genetic variation. There are many local chicken breeds in China. RESULTS: Here, we provide a population genome landscape of genetic variations in 86 domestic chickens representing 10 phenotypically diverse breeds. Genome-wide analysis indicated that sex chromosomes have less genetic diversity and are under stronger selection than autosomes during domestication and local adaptation. We found an evidence of admixture between Tibetan chickens and other domestic population. We further identified strong signatures of selection affecting genomic regions that harbor genes underlying economic traits (typically related to feathers, skin color, growth, reproduction and aggressiveness) and local adaptation (to high altitude). By comparing the genomes of the Tibetan and lowland fowls, we identified genes associated with high-altitude adaptation in Tibetan chickens were mainly involved in energy metabolism, body size maintenance and available food sources. CONCLUSIONS: The work provides crucial insights into the distinct evolutionary scenarios occurring under artificial selection for agricultural production and under natural selection for success at high altitudes in chicken. Several genes were identified as candidates for chicken economic traits and other phenotypic traits.

The temporal expression patterns of brain transcriptome during chicken development and ageing.

BACKGROUND: The transcriptional profiles of mammals during brain development and ageing have been characterized. However the global expression patterns of transcriptome in the chicken brain have not been explored. Here, we systematically investigated the temporal expression profiles of lncRNAs and mRNAs across 8 stages (including 3 embryonic stages, 2 growth stages and 3 adult stages) in the female chicken cerebrum. RESULTS: We identified 39,907 putative lncRNAs and 14,558 mRNAs, investigated the temporal expression patterns by tracking a set of age-dependent genes and predicted potential biological functions of lncRNAs based on co-expression network. The results showed that genes with functions in development, synapses and axons exhibited a progressive decay; genes related to immune response were up-regulated with age. CONCLUSIONS: These results may reflect changes in the regulation of transcriptional networks and provide non-coding RNA gene candidates for further studies and would contribute to a comprehensive understanding of the molecular mechanisms of chicken development and may provide insights or deeper understanding regarding the regulatory mechanisms of age-dependent protein coding and non-protein coding genes in chicken. In addition, as the chicken is an important model organism bridging the evolutionary gap between mammals and other vertebrates, these high resolution data may provide a novel evidence to improve our comprehensive understanding of the brain transcriptome during vertebrate evolution.

Genome-Wide Chromatin Structure Changes During Adipogenesis and Myogenesis.

The recently developed high-throughput chromatin conformation capture (Hi-C) technology enables us to explore the spatial architecture of genomes, which is increasingly considered an important regulator of gene expression. To investigate the changes in three-dimensional (3D) chromatin structure and its mediated gene expression during adipogenesis and myogenesis, we comprehensively mapped 3D chromatin organization for four cell types (3T3-L1 pre-adipocytes, 3T3-L1-D adipocytes, C2C12 myoblasts, and C2C12-D myotubes). We demonstrate that the dynamic spatial genome architecture affected gene expression during cell differentiation. A considerable proportion (~22%) of the mouse genome underwent compartment A/B rearrangement during adipogenic and myogenic differentiation, and most (~80%) upregulated marker genes exhibited an active chromatin state with B to A switch or stable A compartment. More than half (65.4%-73.2%) of the topologically associating domains (TADs) are dynamic. The newly formed TAD and intensified local interactions in the Fabp gene cluster indicated more precise structural regulation of the expression of pro-differentiation genes during adipogenesis. About half (32.39%-59.04%) of the differential chromatin interactions (DCIs) during differentiation are promoter interactions, although these DCIs only account for a small proportion of genome-wide interactions (~9.67% in adipogenesis and ~4.24% in myogenesis). These differential promoter interactions were enriched with promoter-enhancer interactions (PEIs), which were mediated by typical adipogenic and myogenic transcription factors. Differential promoter interactions also included more differentially expressed genes than nonpromoter interactions. Our results provide a global view of dynamic chromatin interactions during adipogenesis and myogenesis and are a resource for studying long-range chromatin interactions mediating the expression of pro-differentiation genes.

Transcriptome differences in frontal cortex between wild boar and domesticated pig.

Animal domestication is a long-term, multistage process that results in modifications of many traits, especially the less aggressive behavior in domesticated animals. In this study, we used the Illumina RNA-seq to compare the transcriptome in brain frontal cortex between wild boar and Rongchang pig, a typical indigenous domestic pig in China, and revealed that 604 genes and 639 genes were specifically detected in wild boar and domesticated pig, respectively, with distinct functional characteristics that may be related to their respective environment. In addition, we identified 60 differentially expressed genes showing an enrichment in immune response-related function. Further comparison of the results against previous studies identified seven genes that are associated with domestication. Our results provide insights for deciphering the mechanism of pig domestication in the future.

Long non-coding RNAs and mRNAs profiling during spleen development in pig.

Genome-wide transcriptomic studies in humans and mice have become extensive and mature. However, a comprehensive and systematic understanding of protein-coding genes and long non-coding RNAs (lncRNAs) expressed during pig spleen development has not been achieved. LncRNAs are known to participate in regulatory networks for an array of biological processes. Here, we constructed 18 RNA libraries from developing fetal pig spleen (55 days before birth), postnatal pig spleens (0, 30, 180 days and 2 years after birth), and the samples from the 2-year-old Wild Boar. A total of 15,040 lncRNA transcripts were identified among these samples. We found that the temporal expression pattern of lncRNAs was more restricted than observed for protein-coding genes. Time-series analysis showed two large modules for protein-coding genes and lncRNAs. The up-regulated module was enriched for genes related to immune and inflammatory function, while the down-regulated module was enriched for cell proliferation processes such as cell division and DNA replication. Co-expression networks indicated the functional relatedness between protein-coding genes and lncRNAs, which were enriched for similar functions over the series of time points examined. We identified numerous differentially expressed protein-coding genes and lncRNAs in all five developmental stages. Notably, ceruloplasmin precursor (CP), a protein-coding gene participating in antioxidant and iron transport processes, was differentially expressed in all stages. This study provides the first catalog of the developing pig spleen, and contributes to a fuller understanding of the molecular mechanisms underpinning mammalian spleen development.

Genome-wide profiling of gene expression and DNA methylation provides insight into low-altitude acclimation in Tibetan pigs.

Efforts have been made to characterize the high-altitude adaption in Tibetan pigs and identified vast of genes or genomic regions undergone natural selection. Nonetheless, information concerning gene expression and DNA methylation changes response to low-altitude acclimation in Tibetan pigs is long overdue. To explore the exceptional mechanisms of gene expression and DNA methylation that are induced by low altitude environments in Tibetan pigs, we performed a comparative transcriptomic and DNA methylation analysis of skeletal muscle between indigenous Tibetan pigs that reside in high altitude regions (~4000m) and their counterparts that migrated to the geographically neighboring low-altitude regions (~500m) for nearly ten generations. Many genes that related to hypoxia response (EGLN3 and FLT1) and energy metabolism (TFB2M) were differentially expressed, but without significant DNA methylation changes. We also found genes embedded in differentially methylated regions were mainly involved in 'Starch and sucrose metabolism', 'glucuronosyltransferase activity' processes. Specifically, our results showed increased SIN3A mRNA expression, with hypomethylation status of its promoter, in longissimus dorsi muscle of low-altitude Tibetan pig. Another gene, CACNG6, showed decreasing expression level with an elevated methylation in its intron 3. These results indicated DNA-methylation-mediated expression alterations in low-altitude acclimation. We envision that this study will serve as a valuable resource for mammal acclimation research and agricultural food industry.

Genetic diversity and natural selection in wild fruit flies revealed by whole-genome resequencing.

We characterized 26 wild fruit flies comparative population genomics from six different altitude and latitude locations by whole genome resequencing. Genetic diversity was relatively higher in Ganzi and Chongqing populations. We also found 13 genes showing selection signature between different altitude flies and variants related to hypoxia and temperature stimulus, were preferentially selected during the flies evolution. One of the most striking selective sweeps found in all high altitude flies occurred in the region harboring Hsp70Aa and Hsp70Ab on chromosome 3R. Interestingly, these two genes are involved in GO terms including response to hypoxia, unfolded protein, temperature stimulus, heat, oxygen levels. Mutation in HPH gene, a candidate gene in the hypoxia inducible factor pathway, might contributes to hypoxic high-altitude adaptation. Intriguingly, some of the selected genes, primarily utilized in humans, were involved in the response to hypoxia, which could imply a conserved molecular mechanisms underlying high-altitude adaptation between insects and humans.

Comparative transcriptomics of 5 high-altitude vertebrates and their low-altitude relatives.

Background: Species living at high altitude are subject to strong selective pressures due to inhospitable environments (e.g., hypoxia, low temperature, high solar radiation, and lack of biological production), making these species valuable models for comparative analyses of local adaptation. Studies that have examined high-altitude adaptation have identified a vast array of rapidly evolving genes that characterize the dramatic phenotypic changes in high-altitude animals. However, how high-altitude environment shapes gene expression programs remains largely unknown. Findings: We generated a total of 910 Gb of high-quality RNA-seq data for 180 samples derived from 6 tissues of 5 agriculturally important high-altitude vertebrates (Tibetan chicken, Tibetan pig, Tibetan sheep, Tibetan goat, and yak) and their cross-fertile relatives living in geographically neighboring low-altitude regions. Of these, ∼75% reads could be aligned to their respective reference genomes, and on average ∼60% of annotated protein coding genes in each organism showed FPKM expression values greater than 0.5. We observed a general concordance in topological relationships between the nucleotide alignments and gene expression-based trees. Tissue and species accounted for markedly more variance than altitude based on either the expression or the alternative splicing patterns. Cross-species clustering analyses showed a tissue-dominated pattern of gene expression and a species-dominated pattern for alternative splicing. We also identified numerous differentially expressed genes that could potentially be involved in phenotypic divergence shaped by high-altitude adaptation. Conclusions: These data serve as a valuable resource for examining the convergence and divergence of gene expression changes between species as they adapt or acclimatize to high-altitude environments.

Genomic data for 78 chickens from 14 populations.

Background: Since the domestication of the red jungle fowls ( Gallus gallus ; dating back to ∼10 000 B.P.) in Asia, domestic chickens ( Gallus gallus domesticus ) have been subjected to the combined effects of natural selection and human-driven artificial selection; this has resulted in marked phenotypic diversity in a number of traits, including behavior, body composition, egg production, and skin color. Population genomic variations through diversifying selection have not been fully investigated. The whole genomes of 78 domestic chickens were sequenced to an average of 18-fold coverage for each bird. By combining this data with publicly available genomes of five wild red jungle fowls and eight Xishuangbanna game fowls, we conducted a comprehensive comparative genomics analysis of 91 chickens from 17 populations. After aligning ∼21.30 gigabases (Gb) of high-quality data from each individual to the reference chicken genome, we identified ∼6.44 million (M) single nucleotide polymorphisms (SNPs) for each population. These SNPs included 1.10 M novel SNPs in 17 populations that were absent in the current chicken dbSNP (Build 145) entries. The current data is important for population genetics and further studies in chickens and will serve as a valuable resource for investigating diversifying selection and candidate genes for selective breeding in chickens.

mRNA N6-methyladenosine methylation of postnatal liver development in pig.

N6-methyladenosine (m6A) is a ubiquitous reversible epigenetic RNA modification that plays an important role in the regulation of post-transcriptional protein coding gene expression. Liver is a vital organ and plays a major role in metabolism with numerous functions. Information concerning the dynamic patterns of mRNA m6A methylation during postnatal development of liver has been long overdue and elucidation of this information will benefit for further deciphering a multitude of functional outcomes of mRNA m6A methylation. Here, we profile transcriptome-wide m6A in porcine liver at three developmental stages: newborn (0 day), suckling (21 days) and adult (2 years). About 33% of transcribed genes were modified by m6A, with 1.33 to 1.42 m6A peaks per modified gene. m6A was distributed predominantly around stop codons. The consensus motif sequence RRm6ACH was observed in 78.90% of m6A peaks. A negative correlation (average Pearson's r = -0.45, P < 10-16) was found between levels of m6A methylation and gene expression. Functional enrichment analysis of genes consistently modified by m6A methylation at all three stages showed genes relevant to important functions, including regulation of growth and development, regulation of metabolic processes and protein catabolic processes. Genes with higher m6A methylation and lower expression levels at any particular stage were associated with the biological processes required for or unique to that stage. We suggest that differential m6A methylation may be important for the regulation of nutrient metabolism in porcine liver.

Comprehensive variation discovery and recovery of missing sequence in the pig genome using multiple de novo assemblies.

Uncovering genetic variation through resequencing is limited by the fact that only sequences with similarity to the reference genome are examined. Reference genomes are often incomplete and cannot represent the full range of genetic diversity as a result of geographical divergence and independent demographic events. To more comprehensively characterize genetic variation of pigs (Sus scrofa), we generated de novo assemblies of nine geographically and phenotypically representative pigs from Eurasia. By comparing them to the reference pig assembly, we uncovered a substantial number of novel SNPs and structural variants, as well as 137.02-Mb sequences harboring 1737 protein-coding genes that were absent in the reference assembly, revealing variants left by selection. Our results illustrate the power of whole-genome de novo sequencing relative to resequencing and provide valuable genetic resources that enable effective use of pigs in both agricultural production and biomedical research.

A de novo silencer causes elimination of MITF-M expression and profound hearing loss in pigs.

BACKGROUND: Genesis of novel gene regulatory modules is largely responsible for morphological and functional evolution. De novo generation of novel cis-regulatory elements (CREs) is much rarer than genomic events that alter existing CREs such as transposition, promoter switching or co-option. Only one case of de novo generation has been reported to date, in fish and without involvement of phenotype alteration. Yet, this event likely occurs in other animals and helps drive genetic/phenotypic variation. RESULTS: Using a porcine model of spontaneous hearing loss not previously characterized we performed gene mapping and mutation screening to determine the genetic foundation of the phenotype. We identified a mutation in the non-regulatory region of the melanocyte-specific promoter of microphthalmia-associated transcription factor (MITF) gene that generated a novel silencer. The consequent elimination of expression of the MITF-M isoform led to early degeneration of the intermediate cells of the cochlear stria vascularis and profound hearing loss, as well as depigmentation, all of which resemble the typical phenotype of Waardenburg syndrome in humans. The mutation exclusively affected MITF-M and no other isoforms. The essential function of Mitf-m in hearing development was further validated using a knock-out mouse model. CONCLUSIONS: Elimination of the MITF-M isoform alone is sufficient to cause deafness and depigmentation. To our knowledge, this study provides the first evidence of a de novo CRE in mammals that produces a systemic functional effect.

Transcriptomic analysis provides insight into high-altitude acclimation in domestic goats.

Domestic goats are distributed in a wide range of habitats and have acclimated to their local environmental conditions. To investigate the gene expression changes of goats that are induced by high altitude stress, we performed RNA-seq on 27 samples from the three hypoxia-sensitive tissues (heart, lung, and skeletal muscle) in three indigenous populations from distinct altitudes (600 m, 2000 m, and 3000 m). We generated 129Gb of high-quality sequencing data (~4Gb per sample) and catalogued the expression profiles of 12,421 annotated hircine genes in each sample. The analysis showed global similarities and differences of high-altitude transcriptomes among populations and tissues as well as revealed that the heart underwent the most high-altitude induced expression changes. We identified numerous differentially expressed genes that exhibited distinct expression patterns, and nonsynonymous single nucleotide variant-containing genes that were highly differentiated between the high- and low-altitude populations. These genes have known or potential roles in hypoxia response and were enriched in functional gene categories potentially responsible for high-altitude stress. Therefore, they are appealing candidates for further investigation of the gene expression and associated regulatory mechanisms related to high-altitude acclimation.