BACKGROUND: Oil flax (linseed, Linum usitatissimum L.) is one of the most important oil crops., However, the increases in drought resulting from climate change have dramatically reduces linseed yield and quality, but very little is known about how linseed coordinates the expression of drought resistance gene in response to different level of drought stress (DS) on the genome-wide level. RESULTS: To explore the linseed transcriptional response of DS and repeated drought (RD) stress, we determined the drought tolerance of different linseed varieties. Then we performed full-length transcriptome sequencing of drought-resistant variety (Z141) and drought-sensitive variety (NY-17) under DS and RD stress at the seedling stage using single-molecule real-time sequencing and RNA-sequencing. Gene Ontology (GO) and reduce and visualize GO (REVIGO) enrichment analysis showed that upregulated genes of Z141 were enriched in more functional pathways related to plant drought tolerance than those of NY-17 were under DS. In addition, 4436 linseed transcription factors were identified, and 1190 were responsive to stress treatments. Moreover, protein-protein interaction (PPI) network analysis showed that the proline biosynthesis pathway interacts with stress response genes through RAD50 (DNA repair protein 50) interacting protein 1 (RIN-1). Finally, proline biosynthesis and DNA repair structural gene expression patterns were verified by RT- PCR. CONCLUSIONS: The drought tolerance of Z141 may be related to its upregulation of drought tolerance genes under DS. Proline may play an important role in linseed drought tolerance by maintaining cell osmotic and protecting DNA from ROS damage. In summary, this study provides a new perspective to understand the drought adaptability of linseed.

BACKGROUND: In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from 'omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments. RESULTS: We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs. CONCLUSION: Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

BACKGROUND: In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported. RESULTS: In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses. CONCLUSIONS: The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

BACKGROUND: The accumulation of intracellular fat depots is a polygenic trait. Therefore, the extent of lipid storage in the individuals of a species covers a broad range and is determined by many genetic factors. Quantitative trait loci analysis can be used to identify those genetic differences between two strains of the same species that are responsible for the differences in a given phenotype. We used this method and complementary approaches to identify genes in the yeast Saccharomyces cerevisiae that are involved in neutral lipid storage. RESULTS: We selected two yeast strains, the laboratory strain BY4741 and the wine yeast AWRI1631, with a more than two-fold difference in neutral lipid content. After crossing, sporulation and germination, we used fluorescence activated cell sorting to isolate a subpopulation of cells with the highest neutral lipid content from the pool of segregants. Whole genome sequencing of this subpopulation and of the unsorted pool of segregants implicated several loci that are involved in lipid accumulation. Three of the identified genes, PIG1, PHO23 and RML2, were investigated in more detail. Deletions of these genes and the exchange of the alleles between the two parental strains confirmed that the encoded proteins contribute to neutral lipid storage in S. cerevisiae and that PIG1, PHO23 and RML2 are the major causative genes. Backcrossing of one of the segregants with the parental strains for seven generations revealed additional regions in the genomes of both strains with potential causative genes for the high lipid accumulation phenotype. CONCLUSIONS: We identified several genes that contribute to the phenotype of lipid accumulation in an allele-specific manner. Surprisingly, no allelic variations of genes with known functions in lipid metabolism were found, indicating that the level of storage lipid accumulation is determined by many cellular processes that are not directly related to lipid metabolism.

BACKGROUND: In evolutionary theory, divergence and speciation can arise from long periods of reproductive isolation, genetic mutation, selection and environmental adaptation. After divergence, alleles can either persist in their initial state (ancestral allele - AA), co-exist or be replaced by a mutated state (derived alleles -DA). In this study, we aligned whole genome sequences of individuals from the Bovinae subfamily to the cattle reference genome (ARS.UCD-1.2) for defining ancestral alleles necessary for selection signatures study. RESULTS: Accommodating independent divergent of each lineage from the initial ancestral state, AA were defined based on fixed alleles on at least two groups of yak, bison and gayal-gaur-banteng resulting in ~ 32.4 million variants. Using non-overlapping scanning windows of 10 Kb, we counted the AA observed within taurine and zebu cattle. We focused on the extreme points, regions with top 0. 1% (high count) and regions without any occurrence of AA (null count). High count regions preserved gene functions from ancestral states that are still beneficial in the current condition, while null counts regions were linked to mutated ones. For both cattle, high count regions were associated with basal lipid metabolism, essential for survival of various environmental pressures. Mutated regions were associated to productive traits in taurine, i.e. higher metabolism, cell development and behaviors and in immune response domain for zebu. CONCLUSIONS: Our findings suggest that retaining and losing AA in some regions are varied and made it species-specific with possibility of overlapping as it depends on the selective pressure they had to experience.

BACKGROUND: Dissolved oxygen (DO) in the water is a vital abiotic factor in aquatic animal farming. A hypoxic environment affects the growth, metabolism, and immune system of fish. Glycolipid metabolism is a vital energy pathway under acute hypoxic stress, and it plays a significant role in the adaptation of fish to stressful environments. In this study, we used multi-omics integrative analyses to explore the mechanisms of hypoxia adaptation in Genetically Improved Farmed Tilapia (GIFT, Oreochromis niloticus). RESULTS: The 96 h median lethal hypoxia (96 h-LH50) for GIFT was determined by linear interpolation. We established control (DO: 5.00 mg/L) groups (CG) and hypoxic stress (96 h-LH50: 0.55 mg/L) groups (HG) and extracted liver tissues for high-throughput transcriptome and metabolome sequencing. A total of 581 differentially expressed (DE) genes and 93 DE metabolites were detected between the CG and the HG. Combined analyses of the transcriptome and metabolome revealed that glycolysis/gluconeogenesis and the insulin signaling pathway were down-regulated, the pentose phosphate pathway was activated, and the biosynthesis of unsaturated fatty acids and fatty acid metabolism were up-regulated in GIFT under hypoxia stress. CONCLUSIONS: The results show that lipid metabolism became the primary pathway in GIFT under acute hypoxia stress. Our findings reveal the changes in metabolites and gene expression that occur under hypoxia stress, and shed light on the regulatory pathways that function under such conditions. Ultimately, this information will be useful to devise strategies to decrease the damage caused by hypoxia stress in farmed fish.

BACKGROUND: Calf mortality generally occurs in calves prior to weaning, which is a serious problem in cattle breeding. Several causative variants of monogenic Mendelian disorders in calf mortality have been identified, whereas genetic factors affecting the susceptibility of calves to death are not well known. To identify variants associated with calf mortality in Japanese Black cattle, we evaluated calf mortality as a categorical trait with a threshold model and performed a genome-wide copy number variation (CNV) association study on calf mortality. RESULTS: We identified a 44-kb deleted-type CNV ranging from 103,317,687 to 103,361,802 bp on chromosome 5, which was associated with the mortality of 1-180-day-old calves. The CNV harbored C1RL, a pseudogene, and an IncRNA localized in the C1R and C1S gene cluster, which is a component of the classical complement activation pathway for immune complexes for infectious pathogens. The average complement activity in CNVR_221 homozygotes at postnatal day 7 was significantly lower than that of wild-type animals and heterozygotes. The frequency of the risk allele in dead calves suffering from diarrhea and pneumonia and in healthy cows was 0.35 and 0.28, respectively (odds ratio = 2.2, P = 0.016), suggesting that CNVR_221 was associated with the mortality of Japanese Black calves suffering from an infectious disease. CONCLUSIONS: This study identified a deleted-type CNV associated with the mortality of 1-180-day-old calves. The complement activity in CNVR_221 homozygotes was significantly lower than that in heterozygotes and wild type animals. The frequency of the risk allele was higher in dead calves suffering from an infectious disease than in healthy cows. These results suggest that the existence of CNVR_221 in calves could be attributed to a reduction in complement activity, which in turn leads to susceptibility to infections. Thus, the risk allele could serve as a useful marker to reduce the mortality of infected Japanese Black calves.

BACKGROUND: The polynucleotide phosphorylase (PNPase) is conserved among both Gram-positive and Gram-negative bacteria. As a core part of the Escherichia coli degradosome, PNPase is involved in maintaining proper RNA levels within the bacterial cell. It plays a major role in RNA homeostasis and decay by acting as a 3'-to-5' exoribonuclease. Furthermore, PNPase can catalyze the reverse reaction by elongating RNA molecules in 5'-to-3' end direction which has a destabilizing effect on the prolonged RNA molecule. RNA degradation is often initiated by an endonucleolytic cleavage, followed by exoribonucleolytic decay from the new 3' end. RESULTS: The PNPase mutant from the facultative phototrophic Rhodobacter sphaeroides exhibits several phenotypical characteristics, including diminished adaption to low temperature, reduced resistance to organic peroxide induced stress and altered growth behavior. The transcriptome composition differs in the pnp mutant strain, resulting in a decreased abundance of most tRNAs and rRNAs. In addition, PNPase has a major influence on the half-lives of several regulatory sRNAs and can have both a stabilizing or a destabilizing effect. Moreover, we globally identified and compared differential RNA 3' ends in RNA NGS sequencing data obtained from PNPase, RNase E and RNase III mutants for the first time in a Gram-negative organism. The genome wide RNA 3' end analysis revealed that 885 3' ends are degraded by PNPase. A fair percentage of these RNA 3' ends was also identified at the same genomic position in RNase E or RNase III mutant strains. CONCLUSION: The PNPase has a major influence on RNA processing and maturation and thus modulates the transcriptome of R. sphaeroides. This includes sRNAs, emphasizing the role of PNPase in cellular homeostasis and its importance in regulatory networks. The global 3' end analysis indicates a sequential RNA processing: 5.9% of all RNase E-dependent and 9.7% of all RNase III-dependent RNA 3' ends are subsequently degraded by PNPase. Moreover, we provide a modular pipeline which greatly facilitates the identification of RNA 5'/3' ends. It is publicly available on GitHub and is distributed under ICS license.

Erratum in BMC Genomics. 2021 May 25;22(1):384.

BACKGROUND: The human microbiome forms very complex communities that consist of hundreds to thousands of different microorganisms that not only affect the host, but also participate in disease processes. Several state-of-the-art methods have been proposed for learning the structure of microbial communities and to investigate the relationship between microorganisms and host environmental factors. However, these methods were mainly designed to model and analyze single microbial communities that do not interact with or depend on other communities. Such methods therefore cannot comprehend the properties between interdependent systems in communities that affect host behavior and disease processes. RESULTS: We introduce a novel hierarchical Bayesian framework, called BALSAMICO (BAyesian Latent Semantic Analysis of MIcrobial COmmunities), which uses microbial metagenome data to discover the underlying microbial community structures and the associations between microbiota and their environmental factors. BALSAMICO models mixtures of communities in the framework of nonnegative matrix factorization, taking into account environmental factors. We proposes an efficient procedure for estimating parameters. A simulation then evaluates the accuracy of the estimated parameters. Finally, the method is used to analyze clinical data. In this analysis, we successfully detected bacteria related to colorectal cancer. CONCLUSIONS: These results show that the method not only accurately estimates the parameters needed to analyze the connections between communities of microbiota and their environments, but also allows for the effective detection of these communities in real-world circumstances.

BACKGROUND: Atractylodes DC is the basic original plant of the widely used herbal medicines "Baizhu" and "Cangzhu" and an endemic genus in East Asia. Species within the genus have minor morphological differences, and the universal DNA barcodes cannot clearly distinguish the systemic relationship or identify the species of the genus. In order to solve these question, we sequenced the chloroplast genomes of all species of Atractylodes using high-throughput sequencing. RESULTS: The results indicate that the chloroplast genome of Atractylodes has a typical quadripartite structure and ranges from 152,294 bp (A. carlinoides) to 153,261 bp (A. macrocephala) in size. The genome of all species contains 113 genes, including 79 protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes. Four hotspots, rpl22-rps19-rpl2, psbM-trnD, trnR-trnT(GGU), and trnT(UGU)-trnL, and a total of 42-47 simple sequence repeats (SSR) were identified as the most promising potentially variable makers for species delimitation and population genetic studies. Phylogenetic analyses of the whole chloroplast genomes indicate that Atractylodes is a clade within the tribe Cynareae; Atractylodes species form a monophyly that clearly reflects the relationship within the genus. CONCLUSIONS: Our study included investigations of the sequences and structural genomic variations, phylogenetics and mutation dynamics of Atractylodes chloroplast genomes and will facilitate future studies in population genetics, taxonomy and species identification.

BACKGROUND: With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and FST) to select the most informative SNPs for ancestry inference. RESULTS: Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof. CONCLUSIONS: The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

BACKGROUND: The giant panda (Ailuropoda melanoleuca) is a threatened endemic Chinese species and a flagship species of national and global conservation concern. Life history theory proposes that reproduction and immunity can be mutually constraining and interrelated. Knowledge of immunity changes of male giant pandas during the breeding season is limited. RESULTS: Here, we researched peripheral blood gene expression profiles associated with immunity. Thirteen captive giant pandas, ranging from 9 to 11 years old, were divided into two groups based on their reproductive status. We identified 318 up-regulated DEGs and 43 down-regulated DEGs, which were enriched in 87 GO terms and 6 KEGG pathways. Additionally, we obtained 45 immune-related genes with altered expression, mostly up-regulated, and identified four hub genes HSPA4, SUGT1, SOD1, and IL1B in PPI analysis. These 45 genes were related to pattern recognition receptors, autophagy, peroxisome, proteasome, natural killer cell, antigen processing and presentation. SUGT1 and IL1B were related to pattern recognition receptors. HSP90AA1 was the most up-regulated gene and is a member of heat shock protein 90 family. HSP90 contributes to the translocation of extracellular antigen. KLRD1 encodes CD94, whose complex is an inhibitor of the cytotoxic activity of NK cells, was down-regulated. IGIP, which has the capability of inducing IgA production by B cells, was down-regulated, suggesting low concentration of IgA in male giant pandas. Our results suggest that most immune-related genes were up-regulated and more related to innate immune than adaptive immune. CONCLUSIONS: Our results indicated that breeding male giant pandas presented an immunoenhancement in innate immunity, enhanced antigen presentation and processing in cellular immunity compared to non-breeding males. The humoral immunity of male giant pandas may show a tendency to decrease during the breeding season. This study will provide a foundation for further studies of immunity and reproduction in male giant pandas.

BACKGROUND: Sweetpotato (Ipomoea batatas (L.) Lam.) serves as an important food source for human beings. β-galactosidase (bgal) is a glycosyl hydrolase involved in cell wall modification, which plays essential roles in plant development and environmental stress adaptation. However, the function of bgal genes in sweetpotato remains unclear. RESULTS: In this study, 17 β-galactosidase genes (Ibbgal) were identified in sweetpotato, which were classified into seven subfamilies using interspecific phylogenetic and comparative analysis. The promoter regions of Ibbgals harbored several stress, hormone and light responsive cis-acting elements. Quantitative real-time PCR results displayed that Ibbgal genes had the distinct expression patterns across different tissues and varieties. Moreover, the expression profiles under various hormonal treatments, abiotic and biotic stresses were highly divergent in leaves and root. CONCLUSIONS: Taken together, these findings suggested that Ibbgals might play an important role in plant development and stress responses, which provided evidences for further study of bgal function and sweetpotato breeding.

BACKGROUND: Tartary buckwheat (Fagopyrum tataricum), an important pseudocereal crop, has high economic value due to its nutritional and medicinal properties. However, dehulling of Tartary buckwheat is difficult owing to its thick and tough hull, which has greatly limited the development of the Tartary buckwheat processing industry. The construction of high-resolution genetic maps serves as a basis for identifying quantitative trait loci (QTLs) and qualitative trait genes for agronomic traits. In this study, a recombinant inbred lines (XJ-RILs) population derived from a cross between the easily dehulled Rice-Tartary type and Tartary buckwheat type was genotyped using restriction site-associated DNA (RAD) sequencing to construct a high-density SNP genetic map. Furthermore, QTLs for 1000-grain weight (TGW) and genes controlling hull type were mapped in multiple environments. RESULTS: In total, 4151 bin markers comprising 122,185 SNPs were used to construct the genetic linkage map. The map consisted of 8 linkage groups and covered 1444.15 cM, with an average distance of 0.35 cM between adjacent bin markers. Nine QTLs for TGW were detected and distributed on four loci on chromosome 1 and 4. A major locus detected in all three trials was mapped in 38.2-39.8 cM region on chromosome 1, with an LOD score of 18.1-37.0, and explained for 23.6-47.5% of the phenotypic variation. The genes controlling hull type were mapped to chromosome 1 between marker Block330 and Block331, which was closely followed by the major locus for TGW. The expression levels of the seven candidate genes controlling hull type present in the region between Block330 and Block336 was low during grain development, and no significant difference was observed between the parental lines. Six non-synonymous coding SNPs were found between the two parents in the region. CONCLUSIONS: We constructed a high-density SNP genetic map for the first time in Tartary buckwheat. The mapped major loci controlling TGW and hull type will be valuable for gene cloning and revealing the mechanism underlying grain development and easy dehulling, and marker-assisted selection in Tartary buckwheat.

BACKGROUND: The transcriptional changes around the time of embryonic genome activation in pre-implantation embryos indicate that this process is highly dynamic. In vitro produced porcine blastocysts are known to be less competent than in vivo developed blastocysts. To understand the conditions that compromise developmental competence of in vitro embryos, it is crucial to evaluate the transcriptional profile of porcine embryos during pre-implantation stages. In this study, we investigated the transcriptome dynamics in in vivo developed and in vitro produced 4-cell embryos, morulae and hatched blastocysts. RESULTS: In vivo developed and in vitro produced embryos displayed largely similar transcriptome profiles during development. Enriched canonical pathways from the 4-cell to the morula transition that were shared between in vivo developed and in vitro produced embryos included oxidative phosphorylation and EIF2 signaling. The shared canonical pathways from the morula to the hatched blastocyst transition were 14-3-3-mediated signaling, xenobiotic metabolism general signaling pathway, and NRF2-mediated oxidative stress response. The in vivo developed and in vitro produced hatched blastocysts further were compared to identify molecular signaling pathways indicative of lower developmental competence of in vitro produced hatched blastocysts. A higher metabolic rate and expression of the arginine transporter SLC7A1 were found in in vitro produced hatched blastocysts. CONCLUSIONS: Our findings suggest that embryos with compromised developmental potential are arrested at an early stage of development, while embryos developing to the hatched blastocyst stage display largely similar transcriptome profiles, irrespective of the embryo source. The hatched blastocysts derived from the in vitro fertilization-pipeline showed an enrichment in molecular signaling pathways associated with lower developmental competence, compared to the in vivo developed embryos.

BACKGROUND: Osmanthus fragrans is an important economical plant containing multiple secondary metabolites including flavonoids and anthocyanins. During the past years, the roles of miRNAs in regulating the biosynthesis of secondary metabolites in plants have been widely investigated. However, few studies on miRNA expression profiles and the potential roles in regulating flavonoid biosynthesis have been reported in O. fragrans. RESULTS: In this study, we used high-throughput sequencing technology to analyze the expression profiles of miRNAs in leaf and flower tissues of O. fragrans. As a result, 106 conserved miRNAs distributed in 47 families and 88 novel miRNAs were identified. Further analysis showed there were 133 miRNAs differentially expressed in leaves and flowers. Additionally, the potential target genes of miRNAs as well as the related metabolic pathways were predicted. In the end, flavonoid content was measured in flower and leaf tissues and potential role of miR858 in regulating flavonoid synthesis was illustrated in O. fragrans. CONCLUSIONS: This study not only provided the genome-wide miRNA profiles in the flower and leaf tissue of O. fragrans, but also investigated the potential regulatory role of miR858a in flavonoid synthesis in O. fragrans. The results specifically indicated the connection of miRNAs to the regulation of secondary metabolite biosynthesis in non-model economical plant.

BACKGROUND: Isoprenoids are the most ancient and essential class of metabolites produced in all organisms, either via mevalonate (MVA)-and/or methylerythritol phosphate (MEP)-pathways. The MEP-pathway is present in all plastid-bearing organisms and most eubacteria. However, no comprehensive study reveals the origination and evolutionary characteristics of MEP-pathway genes in eukaryotes. RESULTS: Here, detailed bioinformatics analyses of the MEP-pathway provide an in-depth understanding the evolutionary history of this indispensable biochemical route, and offer a basis for the co-existence of the cytosolic MVA- and plastidial MEP-pathway in plants given the established exchange of the end products between the two isoprenoid-biosynthesis pathways. Here, phylogenetic analyses establish the contributions of both cyanobacteria and Chlamydiae sequences to the plant's MEP-pathway genes. Moreover, Phylogenetic and inter-species syntenic block analyses demonstrate that six of the seven MEP-pathway genes have predominantly remained as single-copy in land plants in spite of multiple whole-genome duplication events (WGDs). Substitution rate and domain studies display the evolutionary conservation of these genes, reinforced by their high expression levels. Distinct phenotypic variation among plants with reduced expression levels of individual MEP-pathway genes confirm the indispensable function of each nuclear-encoded plastid-targeted MEP-pathway enzyme in plant growth and development. CONCLUSION: Collectively, these findings reveal the polyphyletic origin and restrict conservation of MEP-pathway genes, and reinforce the potential function of the individual enzymes beyond production of the isoprenoids intermediates.

BACKGROUND: Chloroplast genome resources can provide useful information for the evolution of plant species. Tea plant (Camellia sinensis) is among the most economically valuable member of Camellia. Here, we determined the chloroplast genome of the first natural triploid Chinary type tea ('Wuyi narcissus' cultivar of Camellia sinensis var. sinensis, CWN) and conducted the genome comparison with the diploid Chinary type tea (Camellia sinensis var. sinensis, CSS) and two types of diploid Assamica type teas (Camellia sinensis var. assamica: Chinese Assamica type tea, CSA and Indian Assamica type tea, CIA). Further, the evolutionary mechanism of the chloroplast genome of Camellia sinensis and the relationships of Camellia species based on chloroplast genome were discussed. RESULTS: Comparative analysis showed the evolutionary dynamics of chloroplast genome of Camellia sinensis were the repeats and insertion-deletions (indels), and distribution of the repeats, indels and substitutions were significantly correlated. Chinese tea and Indian tea had significant differences in the structural characteristic and the codon usage of the chloroplast genome. Analysis of sequence characterized amplified region (SCAR) using sequences of the intergenic spacers (trnE/trnT) showed none of 292 different Camellia sinensis cultivars had similar sequence characteristic to triploid CWN, but the other four Camellia species did. Estimations of the divergence time showed that CIA diverged from the common ancestor of two Assamica type teas about 6.2 Mya (CI: 4.4-8.1 Mya). CSS and CSA diverged to each other about 0.8 Mya (CI: 0.4-1.5 Mya). Moreover, phylogenetic clustering was not exactly consistent with the current taxonomy of Camellia. CONCLUSIONS: The repeat-induced and indel-induced mutations were two important dynamics contributed to the diversification of the chloroplast genome in Camellia sinensis, which were not mutually exclusive. Chinese tea and Indian tea might have undergone different selection pressures. Chloroplast transfer occurred during the polyploid evolution in Camellia sinensis. In addition, our results supported the three different domestication origins of Chinary type tea, Chinese Assamica type tea and Indian Assamica type tea. And, the current classification of some Camellia species might need to be further discussed.

BACKGROUND: Mitochondrial transcription termination factor (mTERF) is a large gene family which plays a significant role during plant growth under various environmental stresses. However, knowledge of mTERF genes in grapevine (Vitis L.) is limited. RESULTS: In this research, a comprehensive analysis of grape mTERF (VvmTERF) genes, including chromosome locations, phylogeny, protein motifs, gene structures, gene duplications, synteny analysis and expression profiles, was conducted. As a result, a total of 25 mTERF genes were identified from the grape genome, which are distributed on 13 chromosomes with diverse densities and segmental duplication events. The grape mTERF gene family is classified into nine clades based on phylogenetic analysis and structural characteristics. These VvmTERF genes showed differential expression patterns in response to multiple phytohormone treatments and biotic stresses, including treatments with abscisic acid and methyl jasmonate, and inoculation of Plasmopara viticola and Erysiphe necator. CONCLUSIONS: These research findings, as the first of its kind in grapevine, will provide useful information for future development of new stress tolerant grape cultivars through genetic manipulation of VvmTERF genes.

BACKGROUND: The ladybird beetle Cryptolaemus montrouzieri Mulsant, 1853 (Coleoptera, Coccinellidae) is used worldwide as a biological control agent. It is a predator of various mealybug pests, but it also feeds on alternative prey and can be reared on artificial diets. Relatively little is known about the underlying genetic adaptations of its feeding habits. RESULTS: We report the first high-quality genome sequence for C. montrouzieri. We found that the gene families encoding chemosensors and digestive and detoxifying enzymes among others were significantly expanded or contracted in C. montrouzieri in comparison to published genomes of other beetles. Comparisons of diet-specific larval development, survival and transcriptome profiling demonstrated that differentially expressed genes on unnatural diets as compared to natural prey were enriched in pathways of nutrient metabolism, indicating that the lower performance on the tested diets was caused by nutritional deficiencies. Remarkably, the C. montrouzieri genome also showed a significant expansion in an immune effector gene family. Some of the immune effector genes were dramatically downregulated when larvae were fed unnatural diets. CONCLUSION: We suggest that the evolution of genes related to chemosensing, digestion, and detoxification but also immunity might be associated with diet adaptation of an insect predator. These findings help explain why this predatory ladybird has become a successful biological control agent and will enable the optimization of its mass rearing and use in biological control programs.

BACKGROUND: Bone marrow mesenchymal stem cells are a potential resource for the clinical therapy of certain diseases. Canine, as a companion animal, living in the same space with human, is an ideal new model for human diseases research. Because of the high prevalence of diabetes, alternative transplantation islets resource (i.e. insulin producing cells) for diabetes treatment will be in urgent need, which makes our research on the transdifferentiation of Bone marrow mesenchymal stem cells into insulin producing cells become more important. RESULT: In this study, we completed the transdifferentiation process and achieved the transcriptome profiling of five samples with two biological duplicates, namely, "BMSCs", "islets", "stage 1", "stage 2" and "stage 3", and the latter three samples were achieved on the second, fifth and eighth day of induction. A total of 11,530 differentially expressed transcripts were revealed in the profiling data. The enrichment analysis of differentially expressed genes revealed several signaling pathways that are essential for regulating proliferation and transdifferentiation, including focal adhesion, ECM-receptor interaction, tight junction, protein digestion and absorption, and the Rap1 signaling pathway. Meanwhile, the obtained protein-protein interaction network and functional identification indicating involvement of three genes, SSTR2, RPS6KA6, and VIP could act as a foundation for further research. CONCLUSION: In conclusion, to the best of our knowledge, this is the first survey of the transdifferentiation of canine BMSCs into insulin-producing cells according with the timeline using next-generation sequencing technology. The three key genes we pick out may regulate decisive genes during the development of transdifferentiation of insulin producing cells.

BACKGROUND: Not all cells in a given individual are identical in their genomic makeup. Mosaicism describes such a phenomenon where a mixture of genotypic states in certain genomic segments exists within the same individual. Mosaicism is a prevalent and impactful class of non-integer state copy number variation (CNV). Mosaicism implies that certain cell types or subset of cells contain a CNV in a segment of the genome while other cells in the same individual do not. Several studies have investigated the impact of mosaicism in single patients or small cohorts but no comprehensive scan of mosaic CNVs has been undertaken to accurately detect such variants and interpret their impact on human health and disease. RESULTS: We developed a tool called Montage to improve the accuracy of detection of mosaic copy number variants in a high throughput fashion. Montage directly interfaces with ParseCNV2 algorithm to establish disease phenotype genome-wide association and determine which genomic ranges had more or less than expected frequency of mosaic events. We screened for mosaic events in over 350,000 samples using 1% allele frequency as the detection limit. Additionally, we uncovered disease associations of multiple phenotypes with mosaic CNVs at several genomic loci. We additionally investigated the allele imbalance observations genome-wide to define non-diploid and non-integer copy number states. CONCLUSIONS: Our novel algorithm presents an efficient tool with fast computational runtime and high levels of accuracy of mosaic CNV detection. A curated mosaic CNV callset of 3716 events in 2269 samples is presented with comparability to previous reports and disease phenotype associations. The new algorithm can be freely accessed via: https://github.com/CAG-CNV/MONTAGE .

BACKGROUND: Saccharomycodes ludwigii belongs to the poorly characterized Saccharomycodeacea family and is known by its ability to spoil wines, a trait mostly attributable to its high tolerance to sulfur dioxide (SO2). To improve knowledge about Saccharomycodeacea our group determined whole-genome sequences of Hanseniaspora guilliermondii (UTAD222) and S. ludwigii (UTAD17), two members of this family. While in the case of H. guilliermondii the genomic information elucidated crucial aspects concerning the physiology of this species in the context of wine fermentation, the draft sequence obtained for S. ludwigii was distributed by more than 1000 contigs complicating extraction of biologically relevant information. In this work we describe the results obtained upon resequencing of S. ludwigii UTAD17 genome using PacBio as well as the insights gathered from the exploration of the annotation performed over the assembled genome. RESULTS: Resequencing of S. ludwigii UTAD17 genome with PacBio resulted in 20 contigs totaling 13 Mb of assembled DNA and corresponding to 95% of the DNA harbored by this strain. Annotation of the assembled UTAD17 genome predicts 4644 protein-encoding genes. Comparative analysis of the predicted S. ludwigii ORFeome with those encoded by other Saccharomycodeacea led to the identification of 213 proteins only found in this species. Among these were six enzymes required for catabolism of N-acetylglucosamine, four cell wall β-mannosyltransferases, several flocculins and three acetoin reductases. Different from its sister Hanseniaspora species, neoglucogenesis, glyoxylate cycle and thiamine biosynthetic pathways are functional in S. ludwigii. Four efflux pumps similar to the Ssu1 sulfite exporter, as well as robust orthologues for 65% of the S. cerevisiae SO2-tolerance genes, were identified in S. ludwigii genome. CONCLUSIONS: This work provides the first genome-wide picture of a S. ludwigii strain representing a step forward for a better understanding of the physiology and genetics of this species and of the Saccharomycodeacea family. The release of this genomic sequence and of the information extracted from it can contribute to guide the design of better wine preservation strategies to counteract spoilage prompted by S. ludwigii. It will also accelerate the exploration of this species as a cell factory, specially in production of fermented beverages where the use of Non-Saccharomyces species (including spoilage species) is booming.

BACKGROUND: Toxoplasma gondii is a globally distributed protozoan parasite that establishes life-long asymptomatic infection in humans, often emerging as a life-threatening opportunistic pathogen during immunodeficiency. As an intracellular microbe, Toxoplasma establishes an intimate relationship with its host cell from the outset of infection. Macrophages are targets of infection and they are important in early innate immunity and possibly parasite dissemination throughout the host. Here, we employ an RNA-sequencing approach to identify host and parasite transcriptional responses during infection of mouse bone marrow-derived macrophages (BMDM). We incorporated into our analysis infection with the high virulence Type I RH strain and the low virulence Type II strain PTG. Because the well-known TLR-MyD88 signaling axis is likely of less importance in humans, we examined transcriptional responses in both MyD88+/+ and MyD88-/- BMDM. Long noncoding (lnc) RNA molecules are emerging as key regulators in infection and immunity, and were, therefore, included in our analysis. RESULTS: We found significantly more host genes were differentially expressed in response to the highly virulent RH strain rather than with the less virulent PTG strain (335 versus 74 protein coding genes for RH and PTG, respectively). Enriched in these protein coding genes were subsets associated with the immune response as well as cell adhesion and migration. We identified 249 and 83 non-coding RNAs as differentially expressed during infection with RH and PTG strains, respectively. Although the majority of these are of unknown function, one conserved lncRNA termed mir17hg encodes the mir17 microRNA gene cluster that has been implicated in down-regulating host cell apoptosis during T. gondii infection. Only a minimal number of transcripts were differentially expressed between MyD88 knockout and wild type cells. However, several immune genes were among the differences. While transcripts for parasite secretory proteins were amongst the most highly expressed T. gondii genes during infection, no differentially expressed parasite genes were identified when comparing infection in MyD88 knockout and wild type host BMDM. CONCLUSIONS: The large dataset presented here lays the groundwork for continued studies on both the MyD88-independent immune response and the function of lncRNAs during Toxoplasma gondii infection.

BACKGROUND: The root system provides nutrient absorption and is closely related to abiotic stress tolerance, but it is difficult to study the roots under field conditions. This study was conducted to identify quantitative trait loci (QTL) associated with primary root length (PRL) during soybean seedling growth in hydroponic conditions. A total of 103 F7 recombinant inbred lines (RILs) derived from a cross between K099 (short primary root) and Fendou 16 (long primary root) were used to identify QTL for PRL in soybean. The RIL population was genotyped with 223 simple sequence repeats markers covering 20 chromosomes. Phenotyping for primary root length was performed for 3-weeks plants grown in hydoponic conditions. The identified QTL was validated in near isogenic lines and in a separate RIL population. RESULTS: QTL analysis using inclusive composite interval mapping method identified a major QTL on Gm16 between SSR markers Sat_165 and Satt621, explaining 30.25 % of the total phenotypic variation. The identified QTL, qRL16.1, was further confirmed in a segregating population derived from a residual heterozygous line (RHLs-98). To validate qRL16.1 in a different genetic background, QTL analysis was performed in another F6 RIL population derived from a cross between Union (medium primary root) and Fendou 16, in which a major QTL was detected again in the same genomic region as qRL16.1, explaining 14 % of the total phenotypic variation for PRL. In addition, the effect of qRL16.1 was confirmed using two pair of near-isogenic lines (NILs). PRL was significantly higher in NILs possessing the qRL16.1 allele from Fendou 16 compared to allele from K099. CONCLUSIONS: The qRL16.1 is a novel QTL for primary root length in soybean which provides important information on the genetic control of root development. Identification of this major QTL will facilitate positional cloning and DNA marker-assisted selection for root traits in soybean.

BACKGROUND: As economical traits, food habits domestication can reduce production cost in aquaculture. However, the molecular mechanism underlying food habits domestication has remained elusive. Mandarin fish (Siniperca chuatsi) only feed on live prey fish and refuse artificial diets. In the present study, we domesticated mandarin fish to feed on artificial diets. The two groups were obtained, the fish did not eat artificial diets or ate artificial diets during all of the three domestication processes, named Group W or X, respectively. RESULTS: Using transcriptome and metabolome analysis, we investigated the differentially expressed genes and metabolites between the two groups, and found three common pathways related to food habit domestication, including retinol metabolism, glycerolipid metabolism, and biosynthesis of unsaturated fatty acids pathways. Furthermore, the western blotting and bisulfite sequencing PCR analysis were performed. The gene expression of TFIIF and histone methyltransferase ezh1 were significantly increased and decreased in the fish of Group X, respectively. The total DNA methylation levels of TFIIF gene and tri-methylation of histone H3 at lysine 27 (H3K27me3) were significantly higher and lower in the fish of Group X, respectively. CONCLUSION: It was speculated that mandarin fish which could feed on artificial diets, might be attributed to the lower expression of ezh1, resulting in the decreased level of H3K27me3 and increased level of DNA methylation of TFIIF gene. The high expression of TFIIF gene might up-regulate the expression of genes in retinol metabolism, glycerolipid metabolism and glycerophosphoric metabolism pathways. Our study indicated the relationship between the methylation of DNA and histone and food habits domestication, which might be a novel molecular mechanism of food habits domestication in animals.

BACKGROUND: Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively. RESULTS: To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements. CONCLUSIONS: Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads and the generality of these causes and factors should be tested further in other species.

BACKGROUND: Cataracts are defects of the lens that cause progressive visual impairment and ultimately blindness in many vertebrate species. Most cataracts are age-related, but up to one third have an underlying genetic cause. Cataracts are common in captive zoo animals, but it is often unclear whether these are congenital or acquired (age-related) lesions. RESULTS: Here we used a functional candidate gene screening approach to identify mutations associated with cataracts in a captive giant panda (Ailuropoda melanoleuca). We screened 11 genes often associated with human cataracts and identified a novel missense mutation (c.686G > A) in the MIP gene encoding major intrinsic protein. This is expressed in the lens and normally accumulates in the plasma membrane of lens fiber cells, where it plays an important role in fluid transport and cell adhesion. The mutation causes the replacement of serine with asparagine (p.S229N) in the C-terminal tail of the protein, and modeling predicts that the mutation induces conformational changes that may interfere with lens permeability and cell-cell interactions. CONCLUSION: The c.686G > A mutation was found in a captive giant panda with a unilateral cataract but not in 18 controls from diverse regions in China, suggesting it is most likely a genuine disease-associated mutation rather than a single-nucleotide polymorphism. The mutation could therefore serve as a new genetic marker to predict the risk of congenital cataracts in captive giant pandas.

BACKGROUND: Micrococcus luteus is a group of actinobacteria that is widely used in biotechnology and is being thought as an emerging nosocomial pathogen. With one of the smallest genomes of free-living actinobacteria, it is found in a wide range of environments, but intraspecies genetic diversity and adaptation strategies to various environments remain unclear. Here, comparative genomics, phylogenomics, and genome-wide association studies were used to investigate the genomic diversity, evolutionary history, and the potential ecological differentiation of the species. RESULTS: High-quality genomes of 66 M. luteus strains were downloaded from the NCBI GenBank database and core and pan-genome analysis revealed a considerable intraspecies heterogeneity. Phylogenomic analysis, gene content comparison, and average nucleotide identity calculation consistently indicated that the species has diverged into three well-differentiated clades. Population structure analysis further suggested the existence of an unknown ancestor or the fourth, yet unsampled, clade. Reconstruction of gene gain/loss events along the evolutionary history revealed both early events that contributed to the inter-clade divergence and recent events leading to the intra-clade diversity. We also found convincing evidence that recombination has played a key role in the evolutionary process of the species, with upto two-thirds of the core genes having been affected by recombination. Furthermore, distribution of mammal-associated strains (including pathogens) on the phylogenetic tree suggested that the last common ancestor had a free-living lifestyle, and a few recently diverged lineages have developed a mammal-associated lifestyle separately. Consistently, genome-wide association analysis revealed that mammal-associated strains from different lineages shared genes functionally relevant to the host-associated lifestyle, indicating a recent ecological adaption to the new host-associated habitats. CONCLUSIONS: These results revealed high intraspecies genomic diversity of M. luteus and highlighted that gene gain/loss events and extensive recombination events played key roles in the genome evolution. Our study also indicated that, as a free-living species, some lineages have recently developed or are developing a mammal-associated lifestyle. This study provides insights into the mechanisms that drive the genome evolution and adaption to various environments of a bacterial species.

BACKGROUND: 2-Oxoglutarate and Fe(II)-dependent dioxygenases (2ODDs) belong to the 2-oxoglutarate-dependent dioxygenase (2OGD) superfamily and are involved in various vital metabolic pathways of plants at different developmental stages. These proteins have been extensively investigated in multiple model organisms. However, these enzymes have not been systematically analyzed in tomato. In addition, type I flavone synthase (FNSI) belongs to the 2ODD family and contributes to the biosynthesis of flavones, but this protein has not been characterized in tomato. RESULTS: A total of 131 2ODDs from tomato were identified and divided into seven clades by phylogenetic classification. The Sl2ODDs in the same clade showed similar intron/exon distributions and conserved motifs. The Sl2ODDs were unevenly distributed across the 12 chromosomes, with different expression patterns among major tissues and at different developmental stages of the tomato growth cycle. We characterized several Sl2ODDs and their expression patterns involved in various metabolic pathways, such as gibberellin biosynthesis and catabolism, ethylene biosynthesis, steroidal glycoalkaloid biosynthesis, and flavonoid metabolism. We found that the Sl2ODD expression patterns were consistent with their functions during the tomato growth cycle. These results indicated the significance of Sl2ODDs in tomato growth and metabolism. Based on this genome-wide analysis of Sl2ODDs, we screened six potential FNSI genes using a phylogenetic tree and coexpression analysis. However, none of them exhibited FNSI activity. CONCLUSIONS: Our study provided a comprehensive understanding of the tomato 2ODD family and demonstrated the significant roles of these family members in plant metabolism. We also suggest that no FNSI genes in tomato contribute to the biosynthesis of flavones.

BACKGROUND: The quality and yield of wool determine the economic value of the fine-wool sheep. Therefore, discovering markers or genes relevant to wool traits is the cornerstone for the breeding of fine-wool sheep. In this study, we used the Illumina HiSeq X Ten platform to re-sequence 460 sheep belonging to four different fine-wool sheep breeds, namely, Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Aohan fine-wool sheep (AHS) and Qinghai fine-wool sheep (QHS). Eight wool traits, including fiber diameter (FD), fiber diameter coefficient of variance (FDCV), fiber diameter standard deviation (FDSD), staple length (SL), greasy fleece weight (GFW), clean wool rate (CWR), staple strength (SS) and staple elongation (SE) were examined. A genome-wide association study (GWAS) was performed to detect the candidate genes for the eight wool traits. RESULTS: A total of 8.222 Tb of raw data was generated, with an average of approximately 8.59X sequencing depth. After quality control, 12,561,225 SNPs were available for analysis. And a total of 57 genome-wide significant SNPs and 30 candidate genes were detected for the desired wool traits. Among them, 7 SNPs and 6 genes are related to wool fineness indicators (FD, FDCV and FDSD), 10 SNPs and 7 genes are related to staple length, 13 SNPs and 7 genes are related to wool production indicators (GFW and CWR), 27 SNPs and 10 genes associated with staple elongation. Among these candidate genes, UBE2E3 and RHPN2 associated with fiber diameter, were found to play an important role in keratinocyte differentiation and cell proliferation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results, revealed that multitude significant pathways are related to keratin and cell proliferation and differentiation, such as positive regulation of canonical Wnt signaling pathway (GO:0090263). CONCLUSION: This is the first GWAS on the wool traits by using re-sequencing data in Chinese fine-wool sheep. The newly detected significant SNPs in this study can be used in genome-selective breeding for the fine-wool sheep. And the new candidate genes would provide a good theoretical basis for the fine-wool sheep breeding.

BACKGROUND: Heterosis has been extensively utilized in different crops and made a significant contribution to global food security. Genetic distance (GD) is one of the valuable criteria for selecting parents in hybrid breeding. The objectives of this study were to estimate the GD between parents using both simple sequence repeat (SSR) markers and single nucleotide polymorphism (SNP) markers and to investigate the efficiency of the prediction of hybrid performance based on GD. The experiment comprised of four male parents, 282 female parents and 1128 F1, derived from NCII mating scheme. The hybrids, their parents and two check cultivars were evaluated for two years. Performance of F1, mid-parent heterosis (MPH), and best parent heterosis (BPH) were evaluated for ten agronomic and fiber quality traits, including plant height, boll weight, boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. RESULTS: Heterosis was observed in all hybrids and, the traits like plant height, boll number, boll weight and lint percentage exhibited higher heterosis than the fiber quality traits. Correlations were significant between parental and F1 performances. The F1 performances between three hybrid sets (Elite×Elite, Exotic×Elite, and Historic×Elite) showed significant differences in eight traits, including boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. The correlation of the GD assessed by both SSR and SNP markers was significantly positive. The cluster analysis based on GD results estimated using SNP showed that all the female parents divided into five groups and the F1 performance between these five groups showed significant differences in four traits, including lint percentage, micronaire, fiber strength, and fiber elongation ratio. The correlation between GD and F1 performance, MPH and BPH were significant for lint percentage and micronaire. CONCLUSIONS: Our results suggested that GD between parents could be helpful in heterosis prediction for certain traits. This study reveals that molecular marker analysis can serve as a basis for assigning germplasm into heterotic groups and to provide guidelines for parental selection in hybrid cotton breeding.

BACKGROUND: The ongoing COVID-19 outbreak has caused devastating mortality and posed a significant threat to public health worldwide. Despite the severity of this illness and 2.3 million worldwide deaths, the disease mechanism is mostly unknown. Previous studies that characterized differential gene expression due to SARS-CoV-2 infection lacked robust validation. Although vaccines are  now available, effective treatment options are still out of reach. RESULTS: To characterize the transcriptional activity of SARS-CoV-2 infection, a gene signature consisting of 25 genes was generated using a publicly available RNA-Sequencing (RNA-Seq) dataset of cultured cells infected with SARS-CoV-2. The signature estimated infection level accurately in bronchoalveolar lavage fluid (BALF) cells and peripheral blood mononuclear cells (PBMCs) from healthy and infected patients (mean 0.001 vs. 0.958; P < 0.0001). These signature genes were investigated in their ability to distinguish the severity of SARS-CoV-2 infection in a single-cell RNA-Sequencing dataset. TNFAIP3, PPP1R15A, NFKBIA, and IFIT2 had shown bimodal gene expression in various immune cells from severely infected patients compared to healthy or moderate infection cases. Finally, this signature was assessed using the publicly available ConnectivityMap database to identify potential disease mechanisms and drug repurposing candidates. Pharmacological classes of tricyclic antidepressants, SRC-inhibitors, HDAC inhibitors, MEK inhibitors, and drugs such as atorvastatin, ibuprofen, and ketoconazole showed strong negative associations (connectivity score < - 90), highlighting the need for further evaluation of these candidates for their efficacy in treating SARS-CoV-2 infection. CONCLUSIONS: Thus, using the 25-gene SARS-CoV-2 infection signature, the SARS-CoV-2 infection status was captured in BALF cells, PBMCs and postmortem lung biopsies. In addition, candidate SARS-CoV-2 therapies with known safety profiles were identified. The signature genes could potentially also be used to characterize the COVID-19 disease severity in patients' expression profiles of BALF cells.

BACKGROUND: Pregnancy-specific glycoprotein (PSG) genes belong to the carcinoembryonic antigen (CEA) gene family, within the immunoglobulin gene superfamily. In humans, 10 PSG genes encode closely related secreted glycoproteins. They are exclusively expressed in fetal syncytiotrophoblast cells and represent the most abundant fetal proteins in the maternal blood. In recent years, a role in modulation of the maternal immune system possibly to avoid rejection of the semiallogeneic fetus and to facilitate access of trophoblast cells to maternal resources via the blood system has been suggested. Alternatively, they could serve as soluble pathogen decoy receptors like other members of the CEA family. Despite their clearly different domain organization, similar functional properties have also been observed for murine and bat PSG. As these species share a hemochorial type of placentation and a seemingly convergent formation of PSG genes during evolution, we hypothesized that hemochorial placentae support the evolution of PSG gene families. RESULTS: To strengthen this hypothesis, we have analyzed PSG genes in 57 primate species which exhibit hemochorial or epitheliochorial placentation. In nearly all analyzed apes some 10 PSG genes each could be retrieved from genomic databases, while 6 to 24 PSG genes were found in Old World monkey genomes. Surprisingly, only 1 to 7 PSG genes could be identified in New World monkeys. Interestingly, no PSG genes were found in more distantly related primates with epitheliochorial placentae like lemurs and lorises. The exons encoding the putative receptor-binding domains exhibit strong selection for diversification in most primate PSG as revealed by rapid loss of orthologous relationship during evolution and high ratios of nonsynonymous and synonymous mutations. CONCLUSION: The distribution of trophoblast-specific PSGs in primates and their pattern of selection supports the hypothesis that PSG are still evolving to optimize fetal-maternal or putative pathogen interactions in mammals with intimate contact of fetal cells with the immune system of the mother like in hemochorial placentation.

BACKGROUND: Insect body coloration often functions as camouflage to survive from predators or mate selection. Transportation of pigment precursors or related metabolites from cytoplasm to subcellular pigment granules is one of the key steps in insect pigmentation and usually executed via such transporter proteins as the ATP-binding cassette (ABC) transmembrane transporters and small G-proteins (e.g. Rab protein). However, little is known about the copy numbers of pigment transporter genes in the butterfly genomes and about the roles of pigment transporters in the development of swallowtail butterflies. RESULTS: Here, we have identified 56 ABC transporters and 58 Rab members in the genome of swallowtail butterfly Papilio xuthus. This is the first case of genome-wide gene copy number identification of ABC transporters in swallowtail butterflies and Rab family in lepidopteran insects. Aiming to investigate the contribution of the five genes which are orthologous to well-studied pigment transporters (ABCG: white, scarlet, brown and ok; Rab: lightoid) of fruit fly or silkworm during the development of swallowtail butterflies, we performed CRISPR/Cas9 gene-editing of these genes using P. xuthus as a model and sequenced the transcriptomes of their morphological mutants. Our results indicate that the disruption of each gene produced mutated phenotypes in the colors of larvae (cuticle, testis) and/or adult eyes in G0 individuals but have no effect on wing color. The transcriptomic data demonstrated that mutations induced by CRISPR/Cas9 can lead to the accumulation of abnormal transcripts and the decrease or dosage compensation of normal transcripts at gene expression level. Comparative transcriptomes revealed 606 ~ 772 differentially expressed genes (DEGs) in the mutants of four ABCG transporters and 1443 DEGs in the mutants of lightoid. GO and KEGG enrichment analysis showed that DEGs in ABCG transporter mutants enriched to the oxidoreductase activity, heme binding, iron ion binding process possibly related to the color display, and DEGs in lightoid mutants are enriched in glycoprotein binding and protein kinases. CONCLUSIONS: Our data indicated these transporter proteins play an important role in body color of P. xuthus. Our study provides new insights into the function of ABC transporters and small G-proteins in the morphological development of butterflies.

BACKGROUND: Heat stress (HS) is a major stress event in the life of an animal, with detrimental upshots in production and health. Long-non-coding RNAs (lncRNAs) play an important role in many biological processes by transcriptional regulation. However, no research has been reported on the characterization and functionality of lncRNAs in heat-stressed rats. RESULTS: We studied expression levels of lncRNAs in rats during HS, using strand-specific RNA sequencing. Six rats, three in each of Control (22 ± 1 °C) and H120 (42 °C for 120 min) experimental groups, were used to screen for lncRNAs in their liver and adrenal glands. Totally, 4498 and 7627 putative lncRNAs were identified in liver and adrenal glands of the Control and H120 groups, respectively. The majority of lncRNAs were relatively shorter and contained fewer exons than protein-coding transcripts. In total, 482 (174 up-regulated and 308 down-regulated) and 271 (126 up-regulated and 145 down-regulated) differentially-expressed lncRNAs (DElncRNAs, P < 0.05) were identified in the liver and adrenal glands of the Control and H120 groups, respectively. Furthermore, 1274, 121, and 73 target differentially-expressed genes (DEGs) in the liver were predicted to interact with DElncRNAs based on trans-/cis- and sequence similarity regulatory modes. Functional annotation analyses indicated that these DEGs were mostly significantly enriched in insulin signalling, myeloid leukaemia, and glucagon signalling pathways. Similarly, 437, 73 and 41 target DEGs in the adrenal glands were mostly significantly enriched in the cell cycle (trans-prediction) and lysosome pathways (cis-prediction). The DElncRNAs interacting with DEGs that encode heat shock proteins (HSPs) may play an important role in HS response, which include Hsf4, Dnaja1, Dnajb4, Hsph1 and Hspb1 in the liver, and Dnajb13 and Hspb8 in the adrenal glands. The strand-specific RNA sequencing findings were also further verified through RT-qPCR. CONCLUSIONS: This study is the first to provide a detailed characterization and functional analysis of expression levels of lncRNAs in liver and adrenal glands of heat-stressed rats, which provides basis for further studies on the biological functions of lncRNAs under heat stress in rats and other mammalian species.

BACKGROUND: Autophagy, meaning 'self-eating', is required for the degradation and recycling of cytoplasmic constituents under stressful and non-stressful conditions, which helps to maintain cellular homeostasis and delay aging and longevity in eukaryotes. To date, the functions of autophagy have been heavily studied in yeast, mammals and model plants, but few studies have focused on economically important crops, especially tea plants (Camellia sinensis). The roles played by autophagy in coping with various environmental stimuli have not been fully elucidated to date. Therefore, investigating the functions of autophagy-related genes in tea plants may help to elucidate the mechanism governing autophagy in response to stresses in woody plants. RESULTS: In this study, we identified 35 C. sinensis autophagy-related genes (CsARGs). Each CsARG is highly conserved with its homologues from other plant species, except for CsATG14. Tissue-specific expression analysis demonstrated that the abundances of CsARGs varied across different tissues, but CsATG8c/i showed a degree of tissue specificity. Under hormone and abiotic stress conditions, most CsARGs were upregulated at different time points during the treatment. In addition, the expression levels of 10 CsARGs were higher in the cold-resistant cultivar 'Longjing43' than in the cold-susceptible cultivar 'Damianbai' during the CA period; however, the expression of CsATG101 showed the opposite tendency. CONCLUSIONS: We performed a comprehensive bioinformatic and physiological analysis of CsARGs in tea plants, and these results may help to establish a foundation for further research investigating the molecular mechanisms governing autophagy in tea plant growth, development and response to stress. Meanwhile, some CsARGs could serve as putative molecular markers for the breeding of cold-resistant tea plants in future research.

BACKGROUND: Genetic connectedness is a critical component of genetic evaluation as it assesses the comparability of predicted genetic values across units. Genetic connectedness also plays an essential role in quantifying the linkage between reference and validation sets in whole-genome prediction. Despite its importance, there is no user-friendly software tool available to calculate connectedness statistics. RESULTS: We developed the GCA R package to perform genetic connectedness analysis for pedigree and genomic data. The software implements a large collection of various connectedness statistics as a function of prediction error variance or variance of unit effect estimates. The GCA R package is available at GitHub and the source code is provided as open source. CONCLUSIONS: The GCA R package allows users to easily assess the connectedness of their data. It is also useful to determine the potential risk of comparing predicted genetic values of individuals across units or measure the connectedness level between training and testing sets in genomic prediction.

BACKGROUND: The homologous recombination (HR) pathway is largely inactive in early embryos prior to the first cell division, making it difficult to achieve targeted gene knock-ins. The homology-mediated end joining (HMEJ)-based strategy has been shown to increase knock-in efficiency relative to HR, non-homologous end joining (NHEJ), and microhomology-mediated end joining (MMEJ) strategies in non-dividing cells. RESULTS: By introducing gRNA/Cas9 ribonucleoprotein complex and a HMEJ-based donor template with 1 kb homology arms flanked by the H11 safe harbor locus gRNA target site, knock-in rates of 40% of a 5.1 kb bovine sex-determining region Y (SRY)-green fluorescent protein (GFP) template were achieved in Bos taurus zygotes. Embryos that developed to the blastocyst stage were screened for GFP, and nine were transferred to recipient cows resulting in a live phenotypically normal bull calf. Genomic analyses revealed no wildtype sequence at the H11 target site, but rather a 26 bp insertion allele, and a complex 38 kb knock-in allele with seven copies of the SRY-GFP template and a single copy of the donor plasmid backbone. An additional minor 18 kb allele was detected that looks to be a derivative of the 38 kb allele resulting from the deletion of an inverted repeat of four copies of the SRY-GFP template. CONCLUSION: The allelic heterogeneity in this biallelic knock-in calf appears to have resulted from a combination of homology directed repair, homology independent targeted insertion by blunt-end ligation, NHEJ, and rearrangement following editing of the gRNA target site in the donor template. This study illustrates the potential to produce targeted gene knock-in animals by direct cytoplasmic injection of bovine embryos with gRNA/Cas9, although further optimization is required to ensure a precise single-copy gene integration event.

BACKGROUND: High-pressure processing (HPP) is a commonly used technique in the food industry to inactivate pathogens, including L. monocytogenes. It has been shown that L. monocytogenes is able to recover from HPP injuries and can start to grow again during long-term cold storage. To date, the gene expression profiling of L. monocytogenes during HPP damage recovery at cooling temperature has not been studied. In order identify key genes that play a role in recovery of the damage caused by HPP treatment, we performed RNA-sequencing (RNA-seq) for two L. monocytogenes strains (barotolerant RO15 and barosensitive ScottA) at nine selected time points (up to 48 h) after treatment with two pressure levels (200 and 400 MPa). RESULTS: The results showed that a general stress response was activated by SigB after HPP treatment. In addition, the phosphotransferase system (PTS; mostly fructose-, mannose-, galactitol-, cellobiose-, and ascorbate-specific PTS systems), protein folding, and cobalamin biosynthesis were the most upregulated genes during HPP damage recovery. We observed that cell-division-related genes (divIC, dicIVA, ftsE, and ftsX) were downregulated. By contrast, peptidoglycan-synthesis genes (murG, murC, and pbp2A) were upregulated. This indicates that cell-wall repair occurs as a part of HPP damage recovery. We also observed that prophage genes, including anti-CRISPR genes, were induced by HPP. Interestingly, a large amount of RNA-seq data (up to 85%) was mapped to Rli47, which is a non-coding RNA that is upregulated after HPP. Thus, we predicted that Rli47 plays a role in HPP damage recovery in L. monocytogenes. Moreover, gene-deletion experiments showed that amongst peptidoglycan biosynthesis genes, pbp2A mutants are more sensitive to HPP. CONCLUSIONS: We identified several genes and mechanisms that may play a role in recovery from HPP damage of L. monocytogenes. Our study contributes to new information on pathogen inactivation by HPP.

BACKGROUND: Processing and analyzing whole genome sequencing (WGS) is computationally intense: a single Illumina MiSeq WGS run produces ~ 1 million 250-base-pair reads for each of 24 samples. This poses significant obstacles for smaller laboratories, or laboratories not affiliated with larger projects, which may not have dedicated bioinformatics staff or computing power to effectively use genomic data to protect public health. Building on the success of the cloud-based Galaxy bioinformatics platform ( http://galaxyproject.org ), already known for its user-friendliness and powerful WGS analytical tools, the Center for Food Safety and Applied Nutrition (CFSAN) at the U.S. Food and Drug Administration (FDA) created a customized 'instance' of the Galaxy environment, called GalaxyTrakr ( https://www.galaxytrakr.org ), for use by laboratory scientists performing food-safety regulatory research. The goal was to enable laboratories outside of the FDA internal network to (1) perform quality assessments of sequence data, (2) identify links between clinical isolates and positive food/environmental samples, including those at the National Center for Biotechnology Information sequence read archive ( https://www.ncbi.nlm.nih.gov/sra/ ), and (3) explore new methodologies such as metagenomics. GalaxyTrakr hosts a variety of free and adaptable tools and provides the data storage and computing power to run the tools. These tools support coordinated analytic methods and consistent interpretation of results across laboratories. Users can create and share tools for their specific needs and use sequence data generated locally and elsewhere. RESULTS: In its first full year (2018), GalaxyTrakr processed over 85,000 jobs and went from 25 to 250 users, representing 53 different public and state health laboratories, academic institutions, international health laboratories, and federal organizations. By mid-2020, it has grown to 600 registered users and processed over 450,000 analytical jobs. To illustrate how laboratories are making use of this resource, we describe how six institutions use GalaxyTrakr to quickly analyze and review their data. Instructions for participating in GalaxyTrakr are provided. CONCLUSIONS: GalaxyTrakr advances food safety by providing reliable and harmonized WGS analyses for public health laboratories and promoting collaboration across laboratories with differing resources. Anticipated enhancements to this resource will include workflows for additional foodborne pathogens, viruses, and parasites, as well as new tools and services.

BACKGROUND: Despite several RNA-Seq and microarray studies on differentially expressed genes (DEGs) between high- and low-abdominal fat deposition in different broiler lines, to our knowledge, gene coexpression analysis across multiple broiler lines has rarely been reported. Here, we constructed a consensus gene coexpression network focused on identifying consensus gene coexpression modules associated with abdominal fat deposition across multiple broiler lines using two public RNA-Seq datasets (GSE42980 and GSE49121). RESULTS: In the consensus gene coexpression network, we identified eight consensus modules significantly correlated with abdominal fat deposition across four broiler lines using the consensus module analysis function in the weighted gene coexpression network analysis (WGCNA) package. The eight consensus modules were moderately to strongly preserved in the abdominal fat RNA-Seq dataset of another broiler line (SRP058295). Furthermore, we identified 5462 DEGs between high- and low-abdominal fat lines (FL and LL) (GSE42980) and 6904 DEGs between high- and low-growth (HG and LG) (GSE49121), including 1828 overlapping DEGs with similar expression profiles in both datasets, which were clustered into eight consensus modules. Pyruvate metabolism, fatty acid metabolism, and steroid biosynthesis were significantly enriched in the green, yellow, and medium purple 3 consensus modules. The PPAR signaling pathway and adipocytokine signaling pathway were significantly enriched in the green and purple consensus modules. Autophagy, mitophagy, and lysosome were significantly enriched in the medium purple 3 and yellow consensus modules. CONCLUSION: Based on lipid metabolism pathways enriched in eight consensus modules and the overexpression of numerous lipogenic genes in both FL vs. LL and HG vs. LG, we hypothesize that more fatty acids, triacylglycerols (TAGs), and cholesterol might be synthesized in broilers with high abdominal fat than in broilers with low abdominal fat. According to autophagy, mitophagy, and lysosome enrichment in eight consensus modules, we inferred that autophagy might participate in broiler abdominal fat deposition. Altogether, these studies suggest eight consensus modules associated with abdominal fat deposition in broilers. Our study also provides an idea for investigating the molecular mechanism of abdominal fat deposition across multiple broiler lines.

BACKGROUND: Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses. RESULTS: Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max. CONCLUSIONS: The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

BACKGROUND: Intramuscular fat (IMF) content has become one of the most important indicators for measuring meat quality, and levels of IMF are affected by various genes. Long non-coding RNAs (lncRNAs) are widely expressed non-coding RNAs that play an important regulatory role in a variety of biological processes; however, research on the lncRNAs involved in sheep IMF deposition is still in its infancy. Aohan fine-wool sheep (AFWS), one of China's most important meat-hair, dual-purpose sheep breed, provides a great model for studying the role of lncRNAs in the regulation of IMF deposition. We identified lncRNAs by RNA sequencing in Longissimus thoracis et lumborum (LTL) samples of sheep at two ages: 2 months (Mth-2) and 12 months (Mth-12). RESULTS: We identified a total of 26,247 genes and 6935 novel lncRNAs in LTL samples of sheep. Among these, 199 mRNAs and 61 lncRNAs were differentially expressed. We then compared the structural characteristics of lncRNAs and mRNAs. We obtained target genes of differentially expressed lncRNAs (DELs) and performed enrichment analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). We found that target mRNAs were enriched in metabolic processes and developmental pathways. One pathway was significantly enriched, namely tight junction. Based on the analysis of critical target genes, we obtained seven candidate lncRNAs that potentially regulated lipid deposition and constructed a lncRNA-mRNA co-expression network that included MSTRG.4051.3-FZD4, MSTRG.16157.3-ULK1, MSTRG.21053.3-PAQR3, MSTRG.19941.2-TPI1, MSTRG.12864.1-FHL1, MSTRG.2469.2-EXOC6 and MSTRG.21381.1-NCOA1. We speculated that these candidate lncRNAs might play a role by regulating the expression of target genes. We randomly selected five mRNAs and five lncRNAs to verify the accuracy of the sequencing data by qRT-PCR. CONCLUSIONS: Our study identified the differentially expressed mRNAs and lncRNAs during intramuscular lipid deposition in Aohan fine-wool sheep. The work may widen the knowledge about the annotation of the sheep genome and provide a working basis for investigating intramuscular fat deposition in sheep.

BACKGROUND: Paeonia lactiflora 'Hangshao' is widely cultivated in China as a traditional Chinese medicine 'Radix Paeoniae Alba'. Due to the abundant unsaturated fatty acids in its seed, it can also be regarded as a new oilseed plant. However, the process of the biosynthesis of unsaturated fatty acids in it has remained unknown. Therefore, transcriptome analysis is helpful to better understand the underlying molecular mechanisms. RESULTS: Five main fatty acids were detected, including stearic acid, palmitic acid, oleic acid, linoleic acid and α-linolenic acid, and their absolute contents first increased and then decreased during seed development. A total of 150,156 unigenes were obtained by transcriptome sequencing. There were 15,005 unigenes annotated in the seven functional databases, including NR, NT, GO, KOG, KEGG, Swiss-Prot and InterPro. Based on the KEGG database, 1766 unigenes were annotated in the lipid metabolism. There were 4635, 12,304, and 18,291 DEGs in Group I (60 vs 30 DAF), Group II (90 vs 60 DAF) and Group III (90 vs 30 DAF), respectively. A total of 1480 DEGs were detected in the intersection of the three groups. In 14 KEGG pathways of lipid metabolism, 503 DEGs were found, belonging to 111 enzymes. We screened out 123 DEGs involved in fatty acid biosynthesis (39 DEGs), fatty acid elongation (33 DEGs), biosynthesis of unsaturated fatty acid (24 DEGs), TAG assembly (17 DEGs) and lipid storage (10 DEGs). Furthermore, qRT-PCR was used to analyze the expression patterns of 16 genes, including BBCP, BC, MCAT, KASIII, KASII, FATA, FATB, KCR, SAD, FAD2, FAD3, FAD7, GPAT, DGAT, OLE and CLO, most of which showed the highest expression at 45 DAF, except for DGAT, OLE and CLO, which showed the highest expression at 75 DAF. CONCLUSIONS: We predicted that MCAT, KASIII, FATA, SAD, FAD2, FAD3, DGAT and OLE were the key genes in the unsaturated fatty acid biosynthesis and oil accumulation in herbaceous peony seed. This study provides the first comprehensive genomic resources characterizing herbaceous peony seed gene expression at the transcriptional level. These data lay the foundation for elucidating the molecular mechanisms of fatty acid biosynthesis and oil accumulation for herbaceous peony.

BACKGROUND: Despite the increasing use of RNAseq for transcriptome analysis, microarrays remain a widely-used methodology for genomic studies. The latest generation of Affymetrix/Thermo-Fisher microarrays, the ClariomD/XTA and ClariomS array, provide a sensitive and facile method for complex transcriptome expression analysis. However, existing methods of analysis for these high-density arrays do not leverage the statistical power contained in having multiple oligonucleotides representing each gene/exon, but rather summarize probes into a single expression value. We previously developed a methodology, the Sscore algorithm, for probe-level identification of differentially expressed genes (DEGs) between treatment and control samples with oligonucleotide microarrays. The Sscore algorithm was validated for sensitive detection of DEGs by comparison with existing methods. However, the prior version of the Sscore algorithm and a R-based implementation software, sscore, do not function with the latest generations of Affymetrix/Fisher microarrays due to changes in microarray design that eliminated probes previously used for estimation of non-specific binding. RESULTS: Here we describe the GCSscore algorithm, which utilizes the GC-content of a given oligonucleotide probe to estimate non-specific binding using antigenomic background probes found on new generations of arrays. We implemented this algorithm in an improved GCSscore R package for analysis of modern oligonucleotide microarrays. GCSscore has multiple methods for grouping individual probes on the ClariomD/XTA chips, providing the user with differential expression analysis at the gene-level and the exon-level. By utilizing the direct probe-level intensities, the GCSscore algorithm was able to detect DEGs under stringent statistical criteria for all Clariom-based arrays. We demonstrate that for older 3'-IVT arrays, GCSscore produced very similar differential gene expression analysis results compared to the original Sscore method. However, GCSscore functioned well for both the ClariomS and ClariomD/XTA newer microarrays and outperformed existing analysis approaches insofar as the number of DEGs and cognate biological functions identified. This was particularly striking for analysis of the highly complex ClariomD/XTA based arrays. CONCLUSIONS: The GCSscore package represents a powerful new application for analysis of the newest generation of oligonucleotide microarrays such as the ClariomS and ClariomD/XTA arrays produced by Affymetrix/Fisher.

BACKGROUND: Intramuscular fat (IMF) is associated with meat quality and insulin resistance in animals. Research on genetic mechanism of IMF decomposition has positive meaning to pork quality and diseases such as obesity and type 2 diabetes treatment. In this study, an IMF trait segregation population was used to perform RNA sequencing and to analyze the joint or independent effects of genes and long intergenic non-coding RNAs (lincRNAs) on IMF. RESULTS: A total of 26 genes including six lincRNA genes show significantly different expression between high- and low-IMF pigs. Interesting, one lincRNA gene, named IMF related lincRNA (IRLnc) not only has a 292-bp conserved region in 100 vertebrates but also has conserved up and down stream genes (< 10 kb) in pig and humans. Real-time quantitative polymerase chain reaction (RT-qPCR) validation study indicated that nuclear receptor subfamily 4 group A member 3 (NR4A3) which located at the downstream of IRLnc has similar expression pattern with IRLnc. RNAi-mediated loss of function screens identified that IRLnc silencing could inhibit both of the RNA and protein expression of NR4A3. And the in-situ hybridization co-expression experiment indicates that IRLnc may directly binding to NR4A3. As the NR4A3 could regulate the catecholamine catabolism, which could affect insulin sensitivity, we inferred that IRLnc influence IMF decomposition by regulating the expression of NR4A3. CONCLUSIONS: In conclusion, a novel functional noncoding variation named IRLnc has been found contribute to IMF by regulating the expression of NR4A3. These findings suggest novel mechanistic approach for treatment of insulin resistance in human beings and meat quality improvement in animal.