Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents.

BACKGROUND: Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches. RESULTS: We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought. CONCLUSION: The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

Artificial selection footprints in indigenous and commercial chicken genomes.

BACKGROUND: Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China. RESULTS: We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds. CONCLUSION: We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

High quality assemblies of four indigenous chicken genomes and related functional data resources.

Many lines of evidence indicate that red jungle fowl (RJF) is the primary ancestor of domestic chickens. Although multiple versions of RJF (galgal2-galgal5 and GRCg6a) and commercial chickens (GRCg7b/w and Huxu) genomes have been assembled since 2004, no high-quality indigenous chicken genomes have been assembled, hampering the understanding of chicken domestication and evolution. To fill the gap, we sequenced the genomes of four indigenous chickens with distinct morphological traits in southwest China, using a combination of short, long and Hi-C reads. We assembled each genome (~1.0 Gb) into 42 chromosomes with chromosome N50 90.5-90.9 Mb, amongst the highest quality of chicken genome assemblies. To provide resources for gene annotation and functional analysis, we also sequenced transcriptomes of 10 tissues for each of the four chickens. Moreover, we corrected many mis-assemblies and assembled missing micro-chromosomes 29 and 34-39 for GRCg6a. Our assemblies, sequencing data and the correction of GRCg6a can be valuable resources for studying chicken domestication and evolution.

High-quality genome assembly of a C. crossoptilon and related functional and genetics data resources.

There are four species in the Crossoptilon genus inhibiting at from very low to very high altitudes across China, and they are in varying levels of danger of extinction. To better understand the genetic basis of adaptation to high altitudes and genetic changes due to bottleneck, we assembled the genome (~1.02 Gb) of a white eared pheasant (WT) (Crossoptilon crossoptilon) inhibiting at high altitudes (3,000~7,000 m) in northwest of Yunnan province, China, using a combination of Illumina short reads, PacBio long reads and Hi-C reads, with a contig N50 of 19.63 Mb and only six gaps. To further provide resources for gene annotation as well as functional and population genetics analyses, we sequenced transcriptomes of 20 major tissues of the WT individual and re-sequenced another 10 WT individuals and a blue eared pheasant (Crossoptilon auritum) individual inhabiting at intermediate altitudes (1,500~3,000 m). Our assembled WT genome, transcriptome data, and DNA sequencing data can be valuable resources for studying the biology, evolution and developing conservation strategies of these endangered species.

Abnormal saccharides affecting cancer multi-drug resistance (MDR) and the reversal strategies.

Clinically, chemotherapy is the mainstay in the treatment of multiple cancers. However, highly adaptable and activated survival signaling pathways of cancer cells readily emerge after long exposure to chemotherapeutics drugs, resulting in multi-drug resistance (MDR) and treatment failure. Recently, growing evidences indicate that the molecular action mechanisms of cancer MDR are closely associated with abnormalities in saccharides. In this review, saccharides affecting cancer MDR development are elaborated and analyzed in terms of aberrant aerobic glycolysis and its related enzymes, abnormal glycan structures and their associated enzymes, and glycoproteins. The reversal strategies including depletion of ATP, circumventing the original MDR pathway, activation by or inhibition of sugar-related enzymes, combination therapy with traditional cytotoxic agents, and direct modification on the sugar moiety, are ultimately proposed. It follows that abnormal saccharides have a significant effect on cancer MDR development, providing a new perspective for overcoming MDR and improving the outcome of chemotherapy.

Design and synthesis of structurally defined heparan sulfate (HS)-FK506 conjugates as an exogenous approach to investigate biological functions of nucleus HS.

Although heparan sulfate (HS) is widely implicated in numerous physiological and pathological processes, the biological function of nucleus HS remains underexplored, largely due to its complex structure and high hydrophilic property. To supplement these efforts, ideal vehicles are drawing attention as they combine attractive features including lipid solubility for penetrating cell membrane, high affinity binding to its target receptor, metabolic stability, and no cellular actions resulting in toxicity. Herein, we develop a convenient and promising strategy to prepare HS-FK506 conjugates for membrane transport and entry into nucleus, where click chemistry takes easily place between the exocyclic allyl group of a clinic drug FK506 and thiol as a handle incorporated into HS analogues. HS derivatives for constructing the conjugates were synthesized using a cutting-edge chemoenzymatic method. Meantime, [35S] labeled 3'-phosphoadenosine 5'-phosphosulfate (PAP35S) and [14C] glucuronic acid (Glc A) were adopted to label HS-FK506 conjugates, respectively, to evaluate their efficiency of nucleus entry, as a result, 14C Glc A was sensitive, effective and reliable whereas PAP35S gave rise to a mixture of labeled compounds, hampering the understanding of structure-function relationship of nucleus HS. Compared with the corresponding HS, the amount of HS-FK506 conjugates to translocate into nucleus from radioactive assay experiments sharply increased, e.g. tridecasaccharide-FK506 1d increased by approximate 10 folds, offering a simple and robust platform for enabling hydrophilic compounds including carbohydrates to translocate into nucleus and shedding light on their biological functions.

Using MOEA with Redistribution and Consensus Branches to Infer Phylogenies.

In recent years, to infer phylogenies, which are NP-hard problems, more and more research has focused on using metaheuristics. Maximum Parsimony and Maximum Likelihood are two effective ways to conduct inference. Based on these methods, which can also be considered as the optimal criteria for phylogenies, various kinds of multi-objective metaheuristics have been used to reconstruct phylogenies. However, combining these two time-consuming methods results in those multi-objective metaheuristics being slower than a single objective. Therefore, we propose a novel, multi-objective optimization algorithm, MOEA-RC, to accelerate the processes of rebuilding phylogenies using structural information of elites in current populations. We compare MOEA-RC with two representative multi-objective algorithms, MOEA/D and NAGA-II, and a non-consensus version of MOEA-RC on three real-world datasets. The result is, within a given number of iterations, MOEA-RC achieves better solutions than the other algorithms.