BACKGROUND: Recently, Zika virus (ZIKV) re-emerged in India and was potentially associated with microcephaly. However, the molecular mechanisms underlying ZIKV pathogenesis remain to be explored. RESULTS: Herein, we performed a comprehensive RNA-sequencing analysis on ZIKV-infected JEG-3, U-251 MG, and HK-2 cells versus corresponding uninfected controls. Combined with a series of functional analyses, including gene annotation, pathway enrichment, and protein-protein interaction (PPI) network analysis, we defined the molecular characteristics induced by ZIKV infection in different tissues and invasion time points. Data showed that ZIKV infection and replication in each susceptible organ commonly stimulated interferon production and down-regulated metabolic-related processes. Also, tissue-specific immune responses or biological processes (BPs) were induced after ZIKV infection, including GnRH signaling pathway in JEG-3 cells, MAPK signaling pathway in U-251 MG cells, and PPAR signaling pathway in HK-2 cells. Of note, ZIKV infection induced delayed antiviral interferon responses in the placenta-derived cell lines, which potentially explains the molecular mechanism by which ZIKV replicates rapidly in the placenta and subsequential vertical transmission occurs. CONCLUSIONS: Together, these data may provide a systemic insight into the pathogenesis of ZIKV infection in distinct human tissue-derived cell lines, which is likely to help develop prophylactic and therapeutic strategies against ZIKV infection.

BACKGROUND: The dried stem of Cistanche, is a famous Chinese traditional medicine. The main active pharmacodynamic components are phenylethanol glycosides (PhGs). Cistanche tubulosa produces higher level of PhGs in its stems than that of Cistanche deserticola. However, the key genes in the PhGs biosynthesis pathway is not clear in C. tubulosa. RESULTS: In this study, we performed the full-length transcriptome sequencing and gene expression profiling of C. tubulosa using PacBio combined with BGISEQ-500 RNA-seq technology. Totally, 237,772 unique transcripts were obtained, ranging from 199 bp to 31,857 bp. Among the unique transcripts, 188,135 (79.12%) transcripts were annotated. Interestingly, 1080 transcripts were annotated as 22 enzymes related to PhGs biosynthesis. We measured the content of echinacoside, acteoside and total PhGs at two development stages, and found that the content of PhGs was 46.74% of dry matter in young fleshy stem (YS1) and then decreased to 31.22% at the harvest stage (HS2). To compare with YS1, 13,631 genes were up-regulated, and 15,521 genes were down regulated in HS2. Many differentially expressed genes (DEGs) were identified to be involved in phenylpropanoid biosynthesis pathway, phenylalanine metabolism pathway, and tyrosine metabolism pathway. CONCLUSIONS: This is the first report of transcriptome study of C. tubulosa which provided the foundation for understanding of PhGs biosynthesis. Based on these results, we proposed a potential model for PhGs biosynthesis in C. tubulosa.

BACKGROUND: Previous studies have shown that microtubule actin crosslinking factor 1 (MACF1) can regulate osteoblast proliferation and differentiation through non-coding RNA (ncRNA) in bone-forming osteoblasts. However, the role of MACF1 in targeting the competing endogenous RNA (ceRNA) network to regulate osteoblast differentiation remains poorly understood. Here, we profiled messenger RNA (mRNA), microRNA (miRNA), and long ncRNA (lncRNA) expression in MACF1 knockdown MC3TC‑E1 pre‑osteoblast cells. RESULTS: In total, 547 lncRNAs, 107 miRNAs, and 376 mRNAs were differentially expressed. Significantly altered lncRNAs, miRNAs, and mRNAs were primarily found on chromosome 2. A lncRNA-miRNA-mRNA network was constructed using a bioinformatics computational approach. The network indicated that mir-7063 and mir-7646 were the most potent ncRNA regulators and mef2c was the most potent target gene. Pathway enrichment analysis showed that the fluid shear stress and atherosclerosis, p53 signaling, and focal adhesion pathways were highly enriched and contributed to osteoblast proliferation. Importantly, the fluid shear stress and atherosclerosis pathway was co-regulated by lncRNAs and miRNAs. In this pathway, Dusp1 was regulated by AK079370, while Arhgef2 was regulated by mir-5101. Furthermore, Map3k5 was regulated by AK154638 and mir-466q simultaneously. AK003142 and mir-3082-5p as well as Ak141402 and mir-446 m-3p were identified as interacting pairs that regulate target genes. CONCLUSION: This study revealed the global expression profile of ceRNAs involved in the differentiation of MC3TC‑E1 osteoblasts induced by MACF1 deletion. These results indicate that loss of MACF1 activates a comprehensive ceRNA network to regulate osteoblast proliferation.

BACKGROUND: Coat color is important for registration and maintenance of livestock. Standard coat color of Kumamoto sub-breed of Japanese Brown cattle is solid brown, but individuals with diluted coat color have been observed recently. In this study, we attempted to identify polymorphism(s) responsible for coat color dilution by whole genome analysis. RESULTS: One of the diluted cattle possessed 7302 exonic polymorphisms which could affect genes' function. Among them, 14 polymorphisms in 10 coat color-related genes were assumed to be specific for the diluted cattle. Subsequent genotyping with three diluted cattle and 74 standard cattle elucidated that PMEL p.Leu18del was the causative polymorphism for coat color dilution in this sub-breed. Individuals with del/del type of this polymorphism showed diluted coat color, but coat color of heterozygotes were intermediate with various dilution rates. CONCLUSIONS: Coat color dilution of Kumamoto sub-breed was caused by PMEL p.Leu18del. The causative del allele has been detected in several genetically distant cattle breeds, suggesting that PMEL p.Leu18del can be used as a DNA marker to control cattle coat color.

OBJECTIVES: ST11 KPC-producing Klebsiella pneumoniae (Kp) is highly prevalent in China. We investigated the inter- and intra- host transmission and evolution characteristics of ST11 KPC-producing Kp. METHODS: A retrospective study was conducted in a hospital. The clinical data and antimicrobial resistance (AMR) phenotypes were collected. Whole genome sequencing was performed. The transmission route was reconstructed by combining single nucleotide polymorphisms (SNPs) with the clinical information. Hypervirulent Kp (HvKp) was defined as the presence of some combination of peg-344, iroB, iucA, rmpA, or rmpA2. RESULTS: Fifty-eight Kp strains isolated from thirty-five patients were enrolled. The information of one isolate was missing. The mean age of the patients was 74.3 ± 18.0 years, and 18 (50.0%) were female. Fifteen patients (41.7%, 15/36) presented with poor prognosis. All the strains were identified as ST11, and 57 strains harbored blaKPC-2. Two distinguished clades were identified based on the 1,325 high quality SNPs. In clade 1, carbapenem-resistant (CR)-hvKp accounted for 48.3% of the strains (28/58), which mostly presented as KL64 subclones, whereas CR-classical Klebsiella pneumoniae (cKp) commonly possessing KL47 were clustered in Clade 2. One CR-hvKp strain might have originated from the CR-cKp strain from within-host evolution. Even worse, a prolonged transmission of CR-hvKp has led to its spread into healthcare institutes. CONCLUSION: Two endemic subclones of ST11 KPC-producing Kp, KL64-CR-hvKp and KL47-CR-cKp, were transmitted in parallel within the hospital and/or the healthcare institute, suggesting that the ongoing genomic surveillance should be enhanced.

BACKGROUND: RING (Really Interesting New Gene) zinc finger (RING-zf) proteins belong to an important subclass of zinc fingers superfamily, which play versatile roles during various developmental stages and in abiotic stress responses. Based on the conserved cysteine and histidine residues, the RING-zf domains are classified into RING-HC (C3HC4), RING-H2 (C3H2C3), RING-v, RING-D, RING-S/T, RING-G, and RING-C2. However, little is known about the function of the RING-zfs of wheat. RESULTS: In this study, 129 (93.5%) of 138 members were found in nucleus, indicating TaRING-zf were primarily engaged in the degradation of transcription factors and other nuclear-localized proteins. 138 TaRING-zf domains can be divided into four canonical or modified types (RING-H2, RING-HC, RING-D, and RING-M). The RING-M was newly identified in T. aestivum, and might represent the intermediate other states between RING-zf domain and other modified domains. The consensus sequence of the RING-M domain can be described as M-X2-R-X14-Cys-X1-H-X2-Cys-X2-Cys-X10-Cys-X2-Cys. Further interspecies collinearity analyses showed that TaRING-zfs were more closely related to the genes in Poaceae. According to the public transcriptome data, most of the TaRING-zfs were expressed at different 15 stages of plant growth, development, and some of them exhibited specific responses to drought/heat stress. Moreover, 4 RING-HC (TraesCS2A02G526800.1, TraesCS4A02G290600.1, TraesCS4B02G023600.1 and TraesCS4D02G021200.1) and 2 RING-H2 (TraesCS3A02G288900.1 and TraesCS4A02G174600.1) were significantly expressed at different development stages and under drought stress. These findings provide valuable reference data for further study of their physiological functions in wheat varieties. CONCLUSIONS: Taken together, the characterization and classifications of the TaRING-zf family were extensively studied and some new features about it were revealed. This study could provide some valuable targets for further studies on their functions in growth and development, and abiotic stress responses in wheat.

BACKGROUND: Intramuscular preadipocyte differentiation plays a critical role in bovine intramuscular fat (IMF) deposition. However, the roles of different RNAs, including mRNAs, circRNAs, lncRNAs and miRNAs, in regulating the adipogenic differentiation of intramuscular preadipocytes remain largely unclear. RESULTS: In the present study, a whole transcriptome sequencing and analysis, including the analysis of mRNAs, circRNAs, lncRNAs and miRNAs, during different differentiation stages (0, 3, 6, and 9 d) of intramuscular preadipocytes from Qinchuan cattle was performed. All samples were prepared with 3 biological replicates. Here, a total of 27,153 mRNAs, 14,070 circRNAs, 7035 lncRNAs, and 427 miRNAs were annotated. Among them, we identified 4848 differentially expressed mRNAs (DEMs), 181 DE circRNAs (DECs), 501 DE lncRNAs (DELs) and 77 DE miRNAs (DEmiRs) between 0 d and other differentiation days (3, 6, and 9 d). GO and KEGG functional enrichment analyses showed that these differentially expressed genes were mainly enriched in cell differentiation, fat metabolism and adipogenesis-related pathways. Furthermore, weighted gene coexpression network analysis (WGCNA) and co-expression network analysis screened out multiple important mRNAs, circRNAs and lncRNAs related to intramuscular adipogenesis. Based on the competing endogenous RNA (ceRNA) regulatory mechanism, we finally identified 24 potential ceRNA networks and 31 potential key genes, including FOXO1/miR-330/circRNA2018/MSTRG.20301, GPAM/miR-27b/ciRNA489 and SESN3/miR-433/circRNA2627MSTRG.20342. CONCLUSIONS: This study provides new insights into the differential expression patterns of different transcript types (i.e., mRNAs, circRNAs, lncRNAs and miRNAs) in intramuscular preadipocyte differentiation. Our findings provide data support for studying the molecular mechanism of key mRNAs and noncoding RNAs in IMF deposition, and provide new candidate markers for the molecular breeding of beef cattle.

BACKGROUND: Ralstonia solanacearum species complex (RSSC) strains are destructive plant pathogenic bacteria and the causative agents of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, their virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape RSSC virulence, the global diversity, distribution, and wider functional gene content of RSSC prophages are unknown. In this study, prophages were identified in a diverse collection of 192 RSSC draft genome assemblies originating from six continents. RESULTS: Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distributions were characterised using metadata on RSSC strain geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 313 intact prophages were identified, forming ten genetically distinct clusters. These included six prophage clusters with similarity to the Inoviridae, Myoviridae, and Siphoviridae phage families, and four uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distributions, being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained many auxiliary genes involved in metabolism and virulence of both phage and bacteria. CONCLUSIONS: Our results show that while RSSC prophages are highly diverse globally, they make lineage-specific contributions to the RSSC accessory genome, which could have resulted from shared coevolutionary history.

BACKGROUND: Hydrocarbons (HCs) are organic compounds composed solely of carbon and hydrogen that are mainly accumulated in oil reservoirs. As the introduction of all classes of hydrocarbons including crude oil and oil products into the environment has increased significantly, oil pollution has become a global ecological problem. However, our perception of pathways for biotic degradation of major HCs and key enzymes in these bioconversion processes has mainly been based on cultured microbes and is biased by uneven taxonomic representation. Here we used Annotree to provide a gene-centric view of the aerobic degradation ability of aliphatic and aromatic HCs in 23,446 genomes from 123 bacterial and 14 archaeal phyla.  RESULTS: Apart from the widespread genetic potential for HC degradation in Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes, genomes from an additional 18 bacterial and 3 archaeal phyla also hosted key HC degrading enzymes. Among these, such degradation potential has not been previously reported for representatives in the phyla UBA8248, Tectomicrobia, SAR324, and Eremiobacterota. Genomes containing whole pathways for complete degradation of HCs were only detected in Proteobacteria and Actinobacteriota. Except for several members of Crenarchaeota, Halobacterota, and Nanoarchaeota that have tmoA, ladA, and alkB/M key genes, respectively, representatives of archaeal genomes made a small contribution to HC degradation. None of the screened archaeal genomes coded for complete HC degradation pathways studied here; however, they contribute significantly to peripheral routes of HC degradation with bacteria. CONCLUSION: Phylogeny reconstruction showed that the reservoir of key aerobic hydrocarbon-degrading enzymes in Bacteria and Archaea undergoes extensive diversification via gene duplication and horizontal gene transfer. This diversification could potentially enable microbes to rapidly adapt to novel and manufactured HCs that reach the environment.

BACKGROUND: Plant-specific BURP domain-containing genes are involved in plant development and stress responses. However, the role of BURP family in jujube (Ziziphus jujuba Mill.) has not been investigated. RESULTS: In this study, 17 BURP genes belonging to four subfamilies were identified in jujube based on homology analysis, gene structures, and conserved motif confirmation. Gene duplication analysis indicated both tandem duplication and segmental duplication had contributed to ZjBURP expansion. The ZjBURPs were extensively expressed in flowers, young fruits, and jujube leaves. Transcriptomic data and qRT-PCR analysis further revealed that ZjBURPs also significantly influence fruit development, and most genes could be induced by low temperature, salinity, and drought stresses. Notably, several BURP genes significantly altered expression in response to low temperature (ZjPG1) and drought stresses (ZjBNM7, ZjBNM8, and ZjBNM9). CONCLUSIONS: These results provided insights into the possible roles of ZjBURPs in jujube development and stress response. These findings would help selecting candidate ZjBURP genes for cold- and drought-tolerant jujube breeding.

BACKGROUND: MicroRNAs (miRNAs) have been confirmed to be inextricably linked to the emergence of human complex diseases. The identification of the disease-related miRNAs has gradually become a routine way to unveil the genetic mechanisms of examined disorders. METHODS: In this study, a method BLNIMDA based on a weighted bi-level network was proposed for predicting hidden associations between miRNAs and diseases. For this purpose, the known associations between miRNAs and diseases as well as integrated similarities between miRNAs and diseases are mapped into a bi-level network. Based on the developed bi-level network, the miRNA-disease associations (MDAs) are defined as strong associations, potential associations and no associations. Then, each miRNA-disease pair (MDP) is assigned two information properties according to the bidirectional information distribution strategy, i.e., associations of miRNA towards disease and vice-versa. Finally, two affinity weights for each MDP obtained from the information properties and the association type are then averaged as the final association score of the MDP. Highlights of the BLNIMDA lie in the definition of MDA types, and the introduction of affinity weights evaluation from the bidirectional information distribution strategy and defined association types, which ensure the comprehensiveness and accuracy of the final prediction score of MDAs. RESULTS: Five-fold cross-validation and leave-one-out cross-validation are used to evaluate the performance of the BLNIMDA. The results of the Area Under Curve show that the BLNIMDA has many advantages over the other seven selected computational methods. Furthermore, the case studies based on four common diseases and miRNAs prove that the BLNIMDA has good predictive performance. CONCLUSIONS: Therefore, the BLNIMDA is an effective method for predicting hidden MDAs.

BACKGROUND: Genetic engineering of crop plants has been successful in transferring traits into elite lines beyond what can be achieved with breeding techniques. Introduction of transgenes originating from other species has conferred resistance to biotic and abiotic stresses, increased efficiency, and modified developmental programs. The next challenge is now to combine multiple transgenes into elite varieties via gene stacking to combine traits. Generating stable homozygous lines with multiple transgenes requires selection of segregating generations which is time consuming and labor intensive, especially if the crop is polyploid. Insertion site effects and transgene copy number are important metrics for commercialization and trait efficiency. RESULTS: We have developed a simple method to identify the sites of transgene insertions using T-DNA-specific primers and high-throughput sequencing that enables identification of multiple insertion sites in the T1 generation of any crop transformed via Agrobacterium. We present an example using the allohexaploid oil-seed plant Camelina sativa to determine insertion site location of two transgenes. CONCLUSION: This new methodology enables the early selection of desirable transgene location and copy number to generate homozygous lines within two generations.

BACKGROUND: MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass (Panicum virgatum), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass. RESULTS: Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession "Alamo", upland accession "Dacotah", and two upland/lowland F1 hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors. CONCLUSION: Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.

BACKGROUND: Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species. RESULTS: From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip. CONCLUSION: This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

BACKGROUND: Recurrent selection is a foundational breeding method for quantitative trait improvement. It typically features rapid breeding cycles that can lead to high rates of genetic gain. Usually, generations are discrete in recurrent selection, which means that breeding candidates are evaluated and considered for selection for only one cycle. Alternately, generations can overlap, with breeding candidates considered for selection as parents for multiple cycles. With recurrent genomic selection but not phenotypic selection, candidates can be re-evaluated by using genomic estimated breeding values without additional phenotyping of the candidates themselves. Therefore, it may be that candidates with true high breeding values discarded in one cycle due to underestimation of breeding value could be identified and selected in subsequent cycles. The consequences of allowing generations to overlap in recurrent selection are unknown. We assessed whether maintaining overlapping and discrete generations led to differences in genetic gain for phenotypic, genomic truncation, and genomic optimum contribution recurrent selection by stochastic simulation. RESULTS: With phenotypic selection, overlapping generations led to decreased genetic gain compared to discrete generations due to increased selection error bias. Selected individuals, which were in the upper tail of the distribution of phenotypic values, tended to also have high absolute error relative to their true breeding value compared to the overall population. Without repeated phenotyping, these individuals erroneously believed to have high value were repeatedly selected across cycles, leading to decreased genetic gain. With genomic truncation selection, overlapping and discrete generations performed similarly as updating breeding values precluded repeatedly selecting individuals with inaccurately high estimates of breeding values in subsequent cycles. Overlapping generations did not outperform discrete generations in the presence of a positive genetic trend with genomic truncation selection, as individuals from previous breeding cycles typically had truly lower breeding values than candidates from the current generation. With genomic optimum contribution selection, overlapping and discrete generations performed similarly, but overlapping generations slightly outperformed discrete generations in the long term if the targeted inbreeding rate was extremely low. CONCLUSION: Maintaining discrete generations in recurrent phenotypic selection leads to increased genetic gain, especially at low heritabilities, by preventing selection error bias. With genomic truncation selection and genomic optimum contribution selection, genetic gain does not differ between discrete and overlapping generations assuming non-genetic effects are not present. Overlapping generations may increase genetic gain in the long term with very low targeted rates of inbreeding in genomic optimum contribution selection.

BACKGROUND: The hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak. RESULTS: We finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development. CONCLUSION: This study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.

BACKGROUND: Soil salinity is a problem in more than 100 countries across all continents. It is one of the abiotic stress that threatens agriculture the most, negatively affecting crops and reducing productivity. Transcriptomics is a technology applied to characterize the transcriptome in a cell, tissue, or organism at a given time via RNA-Seq, also known as full-transcriptome shotgun sequencing. This technology allows the identification of most genes expressed at a particular stage, and different isoforms are separated and transcript expression levels measured. Once determined by this technology, the expression profile of a gene must undergo validation by another, such as quantitative real-time PCR (qRT-PCR). This study aimed to select, annotate, and validate stress-inducible genes-and their promoters-differentially expressed in the leaves of oil palm (Elaeis guineensis) plants under saline stress. RESULTS: The transcriptome analysis led to the selection of 14 genes that underwent structural and functional annotation, besides having their expression validated using the qRT-PCR technique. When compared, the RNA-Seq and qRT-PCR profiles of those genes resulted in some inconsistencies. The structural and functional annotation analysis of proteins coded by the selected genes showed that some of them are orthologs of genes reported as conferring resistance to salinity in other species. There were those coding for proteins related to the transport of salt into and out of cells, transcriptional regulatory activity, and opening and closing of stomata. The annotation analysis performed on the promoter sequence revealed 22 distinct types of cis-acting elements, and 14 of them are known to be involved in abiotic stress. CONCLUSION: This study has helped validate the process of an accurate selection of genes responsive to salt stress with a specific and predefined expression profile and their promoter sequence. Its results also can be used in molecular-genetics-assisted breeding programs. In addition, using the identified genes is a window of opportunity for strategies trying to relieve the damages arising from the salt stress in many glycophyte crops with economic importance.

BACKGROUND: Polyamines (PAs) are considered promising biostimulants that have diverse key roles during growth and stress responses in plants. Nevertheless, the molecular basis of these roles by PAs has not been completely realized even now, and unfortunately, the transcriptional analyses of the biosynthesis pathway in various wheat tissues have not been investigated under normal or stress conditions. In this research, the findings of genome-wide analyses of genes implicated in the PAs biosynthesis in wheat (ADC, Arginine decarboxylase; ODC, ornithine decarboxylase; AIH, agmatine iminohydrolase; NPL1, Nitrlase like protein 1; SAMDC, S-adenosylmethionine decarboxylase; SPDS, spermidine synthase; SPMS, spermine synthase and ACL5, thermospermine synthase) are shown. RESULTS: In total, thirty PAs biosynthesis genes were identified. Analysis of gene structure, subcellular compartmentation and promoters were discussed. Furthermore, experimental gene expression analyses in roots, shoot axis, leaves, and spike tissues were investigated in adult wheat plants under control and drought conditions. Results revealed structural similarity within each gene family and revealed the identity of two new motifs that were conserved in SPDS, SPMS and ACL5. Analysis of the promoter elements revealed the incidence of conserved elements (STRE, CAAT-box, TATA-box, and MYB TF) in all promoters and highly conserved CREs in >80% of promoters (G-Box, ABRE, TGACG-motif, CGTCA-motif, as1, and MYC). The results of the quantification of PAs revealed higher levels of putrescine (Put) in the leaves and higher spermidine (Spd) in the other tissues. However, no spermine (Spm) was detected in the roots. Drought stress elevated Put level in the roots and the Spm in the leaves, shoots and roots, while decreased Put in spikes and elevated the total PAs levels in all tissues. Interestingly, PA biosynthesis genes showed tissue-specificity and some homoeologs of the same gene family showed differential gene expression during wheat development. Additionally, gene expression analysis showed that ODC is the Put biosynthesis path under drought stress in roots. CONCLUSION: The information gained by this research offers important insights into the transcriptional regulation of PA biosynthesis in wheat that would result in more successful and consistent plant production.

BACKGROUND: Promoters, non-coding DNA sequences located at upstream regions of the transcription start site of genes/gene clusters, are essential regulatory elements for the initiation and regulation of transcriptional processes. Furthermore, identifying promoters in DNA sequences and genomes significantly contributes to discovering entire structures of genes of interest. Therefore, exploration of promoter regions is one of the most imperative topics in molecular genetics and biology. Besides experimental techniques, computational methods have been developed to predict promoters. In this study, we propose iPromoter-Seqvec - an efficient computational model to predict TATA and non-TATA promoters in human and mouse genomes using bidirectional long short-term memory neural networks in combination with sequence-embedded features extracted from input sequences. The promoter and non-promoter sequences were retrieved from the Eukaryotic Promoter database and then were refined to create four benchmark datasets. RESULTS: The area under the receiver operating characteristic curve (AUCROC) and the area under the precision-recall curve (AUCPR) were used as two key metrics to evaluate model performance. Results on independent test sets showed that iPromoter-Seqvec outperformed other state-of-the-art methods with AUCROC values ranging from 0.85 to 0.99 and AUCPR values ranging from 0.86 to 0.99. Models predicting TATA promoters in both species had slightly higher predictive power compared to those predicting non-TATA promoters. With a novel idea of constructing artificial non-promoter sequences based on promoter sequences, our models were able to learn highly specific characteristics discriminating promoters from non-promoters to improve predictive efficiency. CONCLUSIONS: iPromoter-Seqvec is a stable and robust model for predicting both TATA and non-TATA promoters in human and mouse genomes. Our proposed method was also deployed as an online web server with a user-friendly interface to support research communities. Links to our source codes and web server are available at https://github.com/mldlproject/2022-iPromoter-Seqvec .

BACKGROUND: The black-crowned night heron (Nycticorax nycticorax) is an ardeid bird successfully adapted to the nocturnal environment. Previous studies had indicated that the eyes of the night herons have evolved several specialized morphological traits favoring nocturnal vision. However, the molecular mechanisms of the nocturnal vision adaptation of night herons remained inattentions. In this study, the whole genome of N. nycticorax was sequenced and comparative analyses were performed on the vision-related and olfactory receptor (OR) genes to understand the molecular mechanisms of the visual and olfactory adaptation of night herons. RESULTS: The results indicated that a number of vision genes were under positive or relaxed selection in N. nycticorax, whereas a number of other vision genes were under relaxed or intensified selection in the boat-billed heron (Cochlearius cochlearius), which suggested that the two species adapt to nocturnality with different genetic mechanisms. The different selections acting on vision genes are probably associated with the enlargement of eye size and the enhancement of visual sensitivity in night herons. The analyses on olfactory receptor (OR) genes indicated that the total number of OR genes in the genomes of N. nycticorax and C. cochlearius were about half those in the little egret (Egretta garzetta), whereas the diversity of their OR genes was not remarkably different. Additionally, the number of expressed OR genes in the transcriptomes of N. nycticorax was also fewer than that in E. garzetta. These results suggest a reduced olfactory capability in night herons compared with E. garzetta. CONCLUSIONS: Our results provided evidence that several vision genes of the night herons were subjected to different natural selections, which can contribute to a better understanding of the genetic mechanisms of visual adaptions of the night heron. In addition, the finding of the reduced number of total and expressed OR genes in night herons may reflect a trade-off between olfaction and vision.

BACKGROUND: The brown planthopper (Nilaparvata lugens Stål)is a notorious rice pest in many areas of Asia. Study on the molecular mechanisms underlying its development and reproduction will provide scientific basis for effective control. SPARC (Secreted Protein, Acidic and Rich in Cysteine) is one of structural component of the extracellular matrix, which influences a diverse array of biological functions. In this study, the gene for SPARC was identified and functionally analysed from N.lugens. RESULTS: The result showed that the NlSPARC mRNA was highly expressed in fat body, hemolymph and early embryo. The mortality increased significantly when NlSPARC was downregulated after RNA interference (RNAi) in 3 ~ 4th instar nymphs. Downregulation of NlSPARC in adults significantly reduced the number of eggs and offspring, as well as the transcription level of NlSPARC in newly hatched nymphs and survival rate in progeny. The observation with microanatomy on individuals after NlSPARC RNAi showed smaller and less abundant fat body than that in control. No obvious morphological abnormalities in the nymphal development and no differences in development of internal reproductive organ were observed when compared with control. CONCLUSION: NlSPARC is required for oviposition and nymphal development mainly through regulating the tissue of fat body in N.lugens. NlSPARC could be a new candidate target for controlling the rapid propagation of N.lugens population. Our results also demonstrated that the effect of NlSPARC RNAi can transfer to the next generation in N.lugens.

BACKGROUND: Common bean (Phaseolus vulgaris) is considered a staple food in a number of developing countries. Several diseases attack the crop leading to substantial economic losses around the globe. However, the crop has rarely been investigated for multiple disease resistance traits using Meta-analysis approach. RESULTS AND CONCLUSIONS: In this study, in order to identify the most reliable and stable quantitative trait loci (QTL) conveying disease resistance in common bean, we carried out a meta-QTL (MQTL) analysis using 152 QTLs belonging to 44 populations reported in 33 publications within the past 20 years. These QTLs were decreased into nine MQTLs and the average of confidence interval (CI) was reduced by 2.64 folds with an average of 5.12 cM in MQTLs. Uneven distribution of MQTLs across common bean genome was noted where sub-telomeric regions carry most of the corresponding genes and MQTLs. One MQTL was identified to be specifically associated with resistance to halo blight disease caused by the bacterial pathogen Pseudomonas savastanoi pv. phaseolicola, while three and one MQTLs were specifically associated with resistance to white mold and anthracnose caused by the fungal pathogens Sclerotinia sclerotiorum and Colletotrichum lindemuthianum, respectively. Furthermore, two MQTLs were detected governing resistance to halo blight and anthracnose, while two MQTLs were detected for resistance against anthracnose and white mold, suggesting putative genes governing resistance against these diseases at a shared locus. Comparative genomics and synteny analyses provide a valuable strategy to identify a number of well‑known functionally described genes as well as numerous putative novels candidate genes in common bean, Arabidopsis and soybean genomes.

BACKGROUND: Blueberries (Vaccinium sp.) are native to North America and breeding efforts to improve blueberry fruit quality are focused on improving traits such as increased firmness, enhanced flavor and greater shelf-life. Such efforts require additional genomic resources, especially in southern highbush and rabbiteye blueberries. RESULTS: We generated the first full-length fruit transcriptome for the southern highbush and rabbiteye blueberry using the cultivars, Suziblue and Powderblue, respectively. The transcriptome was generated using the Pacific Biosciences single-molecule long-read isoform sequencing platform with cDNA pooled from seven stages during fruit development and postharvest storage. Raw reads were processed through the Isoseq pipeline and full-length transcripts were mapped to the 'Draper' genome with unmapped reads collapsed using Cogent. Finally, we identified 16,299 and 15,882 non-redundant transcripts in 'Suziblue' and 'Powderblue' respectively by combining the reads mapped to Northern Highbush blueberry 'Draper' genome and Cogent analysis. In both cultivars, > 80% of sequences were longer than 1,000 nt, with the median transcript length around 1,700 nt. Functionally annotated transcripts using Blast2GO were > 92% in both 'Suziblue' and 'Powderblue' with overall equal distribution of gene ontology (GO) terms in the two cultivars. Analyses of alternative splicing events indicated that around 40% non-redundant sequences exhibited more than one isoform. Additionally, long non-coding RNAs were predicted to represent 5.6% and 7% of the transcriptomes in 'Suziblue' and 'Powderblue', respectively. Fruit ripening is regulated by several hormone-related genes and transcription factors. Among transcripts associated with phytohormone metabolism/signaling, the highest number of transcripts were related to abscisic acid (ABA) and auxin metabolism followed by those for brassinosteroid, jasmonic acid and ethylene metabolism. Among transcription factor-associated transcripts, those belonging to ripening-related APETALA2/ethylene-responsive element-binding factor (AP2/ERF), NAC (NAM, ATAF1/2 and CUC2), leucine zipper (HB-zip), basic helix-loop-helix (bHLH), MYB (v-MYB, discovered in avian myeloblastosis virus genome) and MADS-Box gene families, were abundant. Further we measured three fruit ripening quality traits and indicators [ABA, and anthocyanin concentration, and texture] during fruit development and ripening. ABA concentration increased during the initial stages of fruit ripening and then declined at the Ripe stage, whereas anthocyanin content increased during the final stages of fruit ripening in both cultivars. Fruit firmness declined during ripening in 'Powderblue'. Genes associated with the above parameters were identified using the full-length transcriptome. Transcript abundance patterns of these genes were consistent with changes in the fruit ripening and quality-related characteristics. CONCLUSIONS: A full-length, well-annotated fruit transcriptome was generated for two blueberry species commonly cultivated in the southeastern United States. The robustness of the transcriptome was verified by the identification and expression analyses of multiple fruit ripening and quality-regulating genes. The full-length transcriptome is a valuable addition to the blueberry genomic resources and will aid in further improving the annotation. It will also provide a useful resource for the investigation of molecular aspects of ripening and postharvest processes.

BACKGROUND: Structural variants (SVs) play important roles in adaptation evolution and species diversification. Especially, in plants, many phenotypes of response to the environment were found to be associated with SVs. Despite the prevalence and significance of SVs, long insertions remain poorly detected and studied in all but model species. RESULTS: We used whole-genome resequencing of paired reads from 80 Asian butternuts to detect long insertions and further analyse their characteristics and potential functional effects. By combining of mapping-based and de novo assembly-based methods, we obtained a multiple related species pangenome representing higher taxonomic groups. We obtained 89,312 distinct contigs totaling 147,773,999 base pair (bp) of new sequences, of which 347 were putative long insertions placed in the reference genome. Most of the putative long insertions appeared in multiple species; in contrast, only 62 putative long insertions appeared in one species, which may be involved in the response to the environment. 65 putative long insertions fell into 61 distinct protein-coding genes involved in plant development, and 105 putative long insertions fell into upstream of 106 distinct protein-coding genes involved in cellular respiration. 3,367 genes were annotated in 2,606 contigs. We propose PLAINS ( https://github.com/CMB-BNU/PLAINS.git ), a streamlined, comprehensive pipeline for the prediction and analysis of long insertions using whole-genome resequencing. CONCLUSIONS: Our study lays down an important foundation for further whole-genome long insertion studies, allowing the investigation of their effects by experiments.

BACKGROUND: The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F1 populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes. RESULTS: Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285-4674 markers per map, with sizes ranging from c. 1643-2292 cM, and mean marker intervals of 0.7-2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1-3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F1 population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits. CONCLUSIONS: Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits.

BACKGROUND: The rapid growth of annual shoots is detrimental to peach production. While gibberellin (GA) promotes the rapid growth of peach shoots, there is limited information on the identity and expression profiles of GA-metabolism genes for this species. RESULTS: All six GA biosynthetic gene families were identified in the peach genome, and the expression profiles of these family members were determined in peach shoots. The upstream biosynthetic gene families have only one or two members (1 CPS, 2 KSs, and 1 KO), while the downstream gene families have multiple members (7 KAOs, 6 GA20oxs, and 5 GA3oxs). Between the two KS genes, PpKS1 showed a relatively high transcript level in shoots, while PpKS2 was undetectable. Among the seven KAO genes, PpKAO2 was highly expressed in shoots, while PpKAO1 and - 6 were weakly expressed. For the six GA20ox genes, both PpGA20ox1 and - 2 were expressed in shoots, but PpGA20ox1 levels were higher than PpGA20ox2. For the five GA3ox genes, only PpGA3ox1 was highly expressed in shoots. Among these biosynthesis genes, PpGA20ox1 and PpGA3ox1 showed a gradual decrease in transcript level along shoots from top to bottom, and a similar trend was observed in bioactive GA1 and GA4 distribution. Among the GA-deactivation genes, PpGA2ox6 was highly expressed in peach shoots. PpGA2ox1 and - 5 transcripts were relatively lower and showed a similar pattern to PpGA20ox1 and PpGA3ox1 in peach shoots. Overexpression of PpGA20ox1, - 2, or PpGA2ox6 in Arabidopsis or tobacco promoted or depressed the plant growth, respectively, while PpGA3ox1 did not affect plant height. Transient expression of PpGA20ox1 in peach leaves significantly increased bioactive GA1 content. CONCLUSIONS: Our results suggest that PpGA20ox and PpGA2ox expression are closely associated with the distribution of active GA1 and GA4 in peach annual shoots. Our research lays a foundation for future studies into ways to effectively repress the rapid growth of peach shoot.

BACKGROUND: Crustose coralline algae (CCA) are calcifying red macroalgae that play important ecological roles including stabilisation of reef frameworks and provision of settlement cues for a range of marine invertebrates. Previous research into the responses of CCA to ocean warming (OW) and ocean acidification (OA) have found magnitude of effect to be species-specific. Response to OW and OA could be linked to divergent underlying molecular processes across species. RESULTS: Here we show Sporolithon durum, a species that exhibits low sensitivity to climate stressors, had little change in metabolic performance and did not significantly alter the expression of any genes when exposed to temperature and pH perturbations. In contrast, Porolithon onkodes, a major coral reef builder, reduced photosynthetic rates and had a labile transcriptomic response with over 400 significantly differentially expressed genes, with differential regulation of genes relating to physiological processes such as carbon acquisition and metabolism. The differential gene expression detected in P. onkodes implicates possible key metabolic pathways, including the pentose phosphate pathway, in the stress response of this species. CONCLUSIONS: We suggest S. durum is more resistant to OW and OA than P. onkodes, which demonstrated a high sensitivity to climate stressors and may have limited ability for acclimatisation. Understanding changes in gene expression in relation to physiological processes of CCA could help us understand and predict how different species will respond to, and persist in, future ocean conditions predicted for 2100.

BACKGROUND: Isolates of Serratia entomophila and S. proteamaculans (Yersiniaceae) cause disease specific to the endemic New Zealand pasture pest, Costelytra giveni (Coleoptera: Scarabaeidae). Previous genomic profiling has shown that S. entomophila isolates appear to have conserved genomes and, where present, conserved plasmids. In the absence of C. giveni larvae, S. entomophila prevalence reduces in the soil over time, suggesting that S. entomophila has formed a host-specific relationship with C. giveni. To help define potential genetic mechanisms driving retention of the chronic disease of S. entomophila, the genome of the isolate 626 was sequenced, enabling the identification of unique chromosomal properties, and defining the gain/loss of accessory virulence factors relevant to pathogenicity to C. giveni larvae. RESULTS: We report the complete sequence of S. entomophila isolate 626, a causal agent of amber disease in C. giveni larvae. The genome of S. entomophila 626 is 5,046,461 bp, with 59.1% G + C content and encoding 4,695 predicted CDS. Comparative analysis with five previously sequenced Serratia species, S. proteamaculans 336X, S. marcescens Db11, S. nematodiphila DH-S01, S. grimesii BXF1, and S. ficaria NBRC 102596, revealed a core of 1,165 genes shared. Further comparisons between S. entomophila 626 and S. proteamaculans 336X revealed fewer predicted phage-like regions and genomic islands in 626, suggesting less horizontally acquired genetic material. Genomic analyses revealed the presence of a four-gene itaconate operon, sharing a similar gene order as the Yersinia pestis ripABC complex. Assessment of a constructed 626::RipC mutant revealed that the operon confer a possible metabolic advantage to S. entomophila in the initial stages of C. giveni infection. CONCLUSIONS: Evidence is presented where, relative to S. proteamaculans 336X, S. entomophila 626 encodes fewer genomic islands and phages, alluding to limited horizontal gene transfer in S. entomophila. Bioassay assessments of a S. entomophila-mutant with a targeted mutation of the itaconate degradation region unique to this species, found the mutant to have a reduced capacity to replicate post challenge of the C. giveni larval host, implicating the itaconate operon in establishment within the host.

BACKGROUND: As one of the microelements, nitrogen play essential roles in cereal production. Although the use of chemical fertilizers has significantly improved the yield of wheat, it has also caused increasingly adverse environmental pollution. Revealing the molecular mechanism manipulating wheat nitrogen use efficiency (NUE), and cultivating wheat germplasms with high nitrogen use efficiency has become important goals for wheat researchers. In this study, we investigated the physiological and transcriptional differences of three wheat cultivars with different NUE under low nitrogen stress. RESULTS: The results showed that, under low nitrogen conditions, the activities of nitrogen metabolism-related enzymes (GS, NR, GDH), antioxidant enzymes (SOD, POD, CAT) and soluble protein contents of ZM366 (high NUE cultivar) were higher than those of JD8 (low NUE cultivar). The hybrid cultivar of ZM366 and JD8 showed mid-parent or over-parent heterosis. Transcriptome analysis revealed that 'alanine, aspartate and glutamate metabolism', 'terpenoid backbone biosynthesis' and 'vitamin B6 metabolism' pathways play key roles in nitrogen use efficiency in wheat. The significant enhancement of the 'Calvin cycle' and 'photorespiration' in ZM366 contributed to its higher level of carbon metabolism under low nitrogen stress, which is an important attribute differs from the other two varieties. In addition, the activation of ABA signal transduction and biosynthesis pathways also helps to maintain NUE under low- nitrogen conditions. Moreover, bHLH transcription factors were also found to play a positive role in wheat NUE. CONCLUSIONS: In conclusion, these results enriched our knowledge of the mechanism of wheat NUE, and provided a theoretical basis for improving wheat NUE and breeding new cultivars.

BACKGROUND: Eukaryotic genome is compartmentalized into structural and functional domains. One of the concepts of higher order organization of chromatin posits that the DNA is organized in constrained loops that behave as independent functional domains. Nuclear Matrix (NuMat), a ribo-proteinaceous nucleoskeleton, provides the structural basis for this organization. DNA sequences located at base of the loops are known as the Matrix Attachment Regions (MARs). NuMat relates to multiple nuclear processes and is partly cell type specific in composition. It is a biochemically defined structure and several protocols have been used to isolate the NuMat where some of the steps have been critically evaluated. These sequences play an important role in genomic organization it is imperative to know their dynamics during development and differentiation. RESULTS: Here we look into the dynamics of MARs when the preparation process is varied and during embryonic development of D. melanogaster. A subset of MARs termed as "Core-MARs" present abundantly in pericentromeric heterochromatin, are constant unalterable anchor points as they associate with NuMat through embryonic development and are independent of the isolation procedure. Euchromatic MARs are dynamic and reflect the transcriptomic profile of the cell. New MARs are generated by nuclear stabilization, and during development, mostly at paused RNA polymerase II promoters. Paused Pol II MARs depend on RNA transcripts for NuMat association. CONCLUSIONS: Our data reveals the role of MARs in functionally dynamic nucleus and contributes to the current understanding of nuclear architecture in genomic context.

BACKGROUND: Insect females undergo a huge transition in energy homeostasis after mating to compensate for nutrient investment during reproduction. To manage with this shift in metabolism, mated females experience extensive morphological, behavioral and physiological changes, including increased food intake and altered digestive processes. However, the mechanisms by which the digestive system responds to mating in females remain barely characterized. Here we performed transcriptomic analysis of the main digestive organ, the midgut, to investigate how gene expression varies with female mating status in Drosophila suzukii, a destructive and invasive soft fruit pest. RESULTS: We sequenced 15,275 unique genes with an average length of 1,467 bp. In total, 652 differentially expressed genes (DEGs) were detected between virgin and mated D. suzukii female midgut libraries. The DEGs were functionally annotated utilizing the GO and KEGG pathway annotation methods. Our results showed that the major GO terms associated with the DEGs from the virgin versus mated female midgut were largely appointed to the metabolic process, response to stimulus and immune system process. We obtained a mass of protein and lipid metabolism genes which were up-regulated and carbohydrate metabolism and immune-related genes which were down-regulated at different time points after mating in female midgut by qRT-PCR. These changes in metabolism and immunity may help supply the female with the nutrients and energy required to sustain egg production. CONCLUSION: Our study characterizes the transcriptional mechanisms driven by mating in the D. suzukii female midgut. Identification and characterization of the DEGs between virgin and mated females midgut will not only be crucial to better understand molecular research related to intestine plasticity during reproduction, but may also provide abundant target genes for the development of effective and ecofriendly pest control strategies against this economically important species.

BACKGROUND: Salt stress is one of the most destructive environmental factors limiting crop growth and development. MicroRNAs (miRNAs) are a class of conserved endogenous small non-coding RNAs, playing the crucial role in regulating salt response and tolerance in plants. However, the miRNAs in wild emmer wheat, especially the key and specific salt-responsive miRNAs are not well studied. RESULTS: Here, we performed small RNA, transcriptome, and degradome sequencing of both of salt-tolerance (ST) and salt-sensitive (SS) wild emmer genotypes to identify the miRNA-mRNA modules associating with salt tolerance. Totally, 775 miRNAs, including 361 conserved known miRNAs and 414 novel miRNAs were detected. Differential expression analysis identified 93 salt-responsive miRNAs under salt stress. Combined with RNA-seq and degradome sequencing analysis, 224 miRNA-mRNA modules displayed the complete opposite expression trends between ST and SS genotypes, most of which functionally enriched into ROS homeostasis maintaining, osmotic pressure modulating, and root growth and development. Finally, the qRT-PCR and a large-scale yeast functional screening were also performed to initially validate the expression pattern and function of candidate genes. CONCLUSIONS: This study reported the key and specific miRNA-mRNA modules associated with salt tolerance in wild emmer, which lay the foundation for improving the salt tolerance in cultivated emmer and bread wheat through miRNA engineering approach.

BACKGROUND: During embryogenesis, the developmental potential of initially pluripotent cells becomes progressively restricted as they transit to lineage restricted states. The pluripotent cells of Xenopus blastula-stage embryos are an ideal system in which to study cell state transitions during developmental decision-making, as gene expression dynamics can be followed at high temporal resolution. RESULTS: Here we use transcriptomics to interrogate the process by which pluripotent cells transit to four different lineage-restricted states: neural progenitors, epidermis, endoderm and ventral mesoderm, providing quantitative insights into the dynamics of Waddington's landscape. Our findings provide novel insights into why the neural progenitor state is the default lineage state for pluripotent cells and uncover novel components of lineage-specific gene regulation. These data reveal an unexpected overlap in the transcriptional responses to BMP4/7 and Activin signaling and provide mechanistic insight into how the timing of signaling inputs such as BMP are temporally controlled to ensure correct lineage decisions. CONCLUSIONS: Together these analyses provide quantitative insights into the logic and dynamics of developmental decision making in early embryos. They also provide valuable lineage-specific time series data following the acquisition of specific lineage states during development.

BACKGROUND: Homalodisca vitripennis Germar, the glassy-winged sharpshooter, is an invasive insect in California and a critical threat to agriculture through its transmission of the plant pathogen, Xylella fastidiosa. Quarantine, broad-spectrum insecticides, and biological control have been used for population management of H. vitripennis since its invasion and subsequent proliferation throughout California. Recently wide-spread neonicotinoid resistance has been detected in populations of H. vitripennis in the southern portions of California's Central Valley. In order to better understand potential mechanisms of H. vitripennis neonicotinoid resistance, we performed RNA sequencing on wild-caught insecticide-resistant and relatively susceptible sharpshooters to profile their transcriptome and population structure. RESULTS: We identified 81 differentially expressed genes with higher expression in resistant individuals. The significant largest differentially expressed candidate gene linked to resistance status was a cytochrome P450 gene with similarity to CYP6A9. Furthermore, we observed an over-enrichment of GO terms representing functions supportive of roles in resistance mechanisms (cytochrome P450s, M13 peptidases, and cuticle structural proteins). Finally, we saw no evidence of broad-scale population structure, perhaps due to H. vitripennis' relatively recent introduction to California or due to the relatively small geographic scale investigated here. CONCLUSIONS: In this work, we characterized the transcriptome of insecticide-resistant and susceptible H. vitripennis and identified candidate genes that may be involved in resistance mechanisms for this species. Future work should seek to build on the transcriptome profiling performed here to confirm the role of the identified genes, particularly the cytochrome P450, in resistance in H. vitripennis. We hope this work helps aid future population management strategies for this and other species with growing insecticide resistance.

BACKGROUND: Merino sheep exhibit high wool production and excellent wool quality. The fleece of Merino sheep is predominantly composed of wool fibers grown from hair follicles (HFs). The HF is a complex biological system involved in a dynamic process governed by gene regulation, and gene expression is regulated by microRNAs (miRNAs). miRNA inhibits posttranscriptional gene expression by specifically binding to target messenger RNA (mRNA) and plays an important role in regulating gene expression, the cell cycle and biological development sequences. The purpose of this study was to examine mRNA and miRNA binding to identify key miRNAs and target genes related to HF development. This will provide new and important insights into fundamental mechanisms that regulate cellular activity and cell fate decisions within and outside of the skin. RESULTS: We analyzed miRNA data in skin tissues collected from 18 Merino sheep on four embryonic days (E65, E85, E105 and E135) and two postnatal days (D7 and D30) and identified 87 differentially expressed miRNAs (DE-miRNAs). These six stages were further divided into two longer developmental stages based on heatmap cluster analysis, and the results showed that DE-mRNAs in Stage A were closely related to HF morphogenesis. A coanalysis of Stage A DE-mRNAs and DE-miRNAs revealed that 9 DE-miRNAs and 17 DE-mRNAs presented targeting relationships in Stage A. We found that miR-23b and miR-133 could target and regulate ACVR1B and WNT10A. In dermal fibroblasts, the overexpression of miR-133 significantly reduced the mRNA and protein expression levels of ACVR1B. The overexpression of miR-23b significantly reduced the mRNA and protein expression levels of WNT10A. CONCLUSION: This study provides a new reference for understanding the molecular basis of HF development and lays a foundation for further improving sheep HF breeding. miRNAs and target genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine-wool sheep.

BACKGROUND: The mechanistic target of rapamycin (mTOR) signal pathway plays a critical regulating role in the occurrence and development of cataract. However, the role of mTORC1 downstream proteins, including ribosomal protein S6K (RP-S6K), eukaryotic initiation factor 4E-binding protein (EIF4EBP), eukaryotic initiation factor 4G (EIF-4G), eukaryotic initiation factor 4E (EIF-4E), and eukaryotic initiation factor 4A (EIF-4A), in regulating cataract development is still unknown. Herein, we conducted a mendelian randomization (MR) study to understand the function of mTORC1 signaling in the process of cataract development. RESULTS: The causal estimate was evaluated with inverse-variance weighted (IVW) estimate, weighted median estimator, MR-Egger and MR robust adjusted profile score (MR. RAPS). The single-nucleotide polymorphisms (SNPs), P<5 × 10- 6 and r2<0.05, were selected to genetically predict the RP-S6K, EIF4EBP, EIF-4E, EIF-4A, and EIF-4G. We included a total of 26,758 cases and 189,604 controls in this MR study. The study revealed causal association between circulating EIF4EBP (OR 1.09, 95% confidence interval 1.03,1.16, P = 0.004), RP-S6K (OR 1.04, 95% confidence interval 1.01, 1.08, P = 0.02) and cataract formation with IVW estimate. Whereas after correcting outliers, MR robust adjusted profile score (MR. RAPS) shows consistent result with IVW for EIF4EBP (OR = 1.08, 95%CI:1.05-1.11, P = 0.007). The observation strengthened the confidence in the true causal associations. However, no association was found for circulating EIF-4E (OR 1.03, 95% confidence interval 0.97, 1.09, P = 0.31), EIF-4A (OR 1.02, 95% confidence interval 0.98, 1.07, P = 0.34), and EIF-4G (OR 1.02, 95% confidence interval 0.94, 1.01, P = 0.64) levels with cataract formation. No evidence of heterogeneity and unbalanced horizontal pleiotropy was detected. CONCLUSION: The MR study suggests that EIF4EBP is a high-risk factor for cataract development. There may be a potential causal association between the mTORC1/EIF4EBP axis and cataract. This research highlights the potential mechanism for cataract development and a genetic target to prevent as well as treat cataracts.

BACKGROUND: Semiothisa cinerearia belongs to Geometridae, which is one of the most species-rich families of lepidopteran insects. It is also one of the most economically significant pests of the Chinese scholar tree (Sophora japonica L.), which is an important urban greenbelt trees in China due to its high ornamental value. A genome assembly of S. cinerearia would facilitate study of the control and evolution of this species. RESULTS: We present a reference genome for S. cinerearia; the size of the genome was ~ 580.89 Mb, and it contained 31 chromosomes. Approximately 43.52% of the sequences in the genome were repeat sequences, and 21,377 protein-coding genes were predicted. Some important gene families involved in the detoxification of pesticides (P450) have expanded in S. cinerearia. Cytochrome P450 gene family members play key roles in mediating relationships between plants and insects, and they are important in plant secondary metabolite detoxification and host-plant selection. Using comparative analysis methods, we find positively selected gene, Sox15 and TipE, which may play important roles during the larval-pupal metamorphosis development of S. cinerearia. CONCLUSION: This assembly provides a new genomic resource that will aid future comparative genomic studies of Geometridae species and facilitate future evolutionary studies on the S. cinerearia.

Mantidis Ootheca (Sangpiaoxiao, mantis egg case) is a typical multi-origin Chinese medicinal material. The Chinese Pharmacopoeia stipulates that the Mantidis Ootheca originates from three species of Mantis: Tenodera sinensis, Statilia maculate, and Hierodula patellifera. However, Mantidis Ootheca mainly relies on field collection, which leads to confusion of its actual origin in the market. As the clinical use of Mantidis Ootheca with unknown original mantis species will pose potential risks to drug safety, it is necessary to survey the commercially available Mantidis Ootheca origin species. However, as the egg case of Mantis, the morphological characters of Mantidis Ootheca are limited and usually cannot serve as accurate identification tool. DNA barcoding, which is widely used in taxonomic studies of animals, is severely affected by the impact of storage pests and DNA degradation. Thus, this study collected a total of 4580 Mantidis Ootheca and pooled separately Mantidis Ootheca samples according to 18 different sources as DNA samples to analyze the origin diversity of Mantidis Ootheca individuals contaminated by common store pests collected in in the market using DNA metabarcoding, and to provide a basis for quality control of Mantidis Ootheca. 37 Mantis ASVs and 9 Mantis MOTUs were identified through species delimitation, and the high-level intraspecific diversity was depicted as haplotype network plot. Besides Tenodera sinensis and Hierodula patellifera as genuine original mantis species defined in the Chinese Pharmacopoeia, Tenodera angustipennis was also the origin species of these Mantidis Ootheca samples.

BACKGROUND: It is important to understand the functional impact of somatic mutation and methylation aberration at an individual level to implement precision medicine. Recent studies have demonstrated that the perturbation of gene interaction networks can provide a fundamental link between genotype (or epigenotype) and phenotype. However, it is unclear how individual mutations affect the function of biological networks, especially for individual methylation aberration. To solve this, we provided a sample-specific driver module construction method using the 2-order network theory and hub-gene theory to identify individual perturbation networks driven by mutations or methylation aberrations. RESULTS: Our method integrated multi-omics of breast cancer, including genomics, transcriptomics, epigenomics and interactomics, and provided new insight into the synergistic collaboration between methylation and mutation at an individual level. A common driver pattern of breast cancer was identified from a novel perspective of a driver module, which is correlated to the occurrence and development of breast cancer. The constructed driver module reflects the survival prognosis and degree of malignancy among different subtypes of breast cancer. Additionally, subtype-specific driver modules were identified. CONCLUSIONS: This study explores the driver module of individual cancer, and contributes to a better understanding of the mechanism of breast cancer driven by the mutations and methylation variations from the point of view of the driver network. This work will help identify new therapeutic combinations of gene mutations and drugs in humans.

BACKGROUND: Mycobacteroides abscessus complex (MABC), an emerging pathogen, causes human infections resistant to multiple antibiotics. In this study, the genome data of 1,581 MABC strains were downloaded from NCBI database for phylogenetic relatedness inference, resistance profile identification and the estimation of evolutionary pressure on resistance genes in silico. RESULTS: From genes associated with resistance to 28 antibiotic classes, 395 putative proteins (ARPs) were identified, based on the information in two antibiotic resistance databases (CARD and ARG-ANNOT). The ARPs most frequently identified in MABC were those associated with resistance to multiple antibiotic classes, beta-lactams and aminoglycosides. After excluding ARPs that had undergone recombination, two ARPs were predicted to be under diversifying selection and 202 under purifying selection. This wide occurrence of purifying selection suggested that the diversity of commonly shared ARPs in MABC have been reduced to achieve stability. The unequal distribution of ARPs in members of the MABC could be due to horizontal gene transfer or ARPs pseudogenization events. Most (81.5%) of the ARPs were observed in the accessory genome and 72.2% ARPs were highly homologous to proteins associated with mobile genetic elements such as plasmids, prophages and viruses. On the other hand, with TBLASTN search, only 18 of the ARPs were identified as pseudogenes. CONCLUSION: Altogether, our results suggested an important role of horizontal gene transfer in shaping the resistome of MABC.

BACKGROUND: The subgenus Gynopodium belonging to genus Magnolia have high ornamental, economic, and ecological value. Subgenus Gynopodium contains eight species, but six of these species are threatened. No studies to date have characterized the characteristics of the chloroplast genomes (CPGs) within subgenus Gynopodium species. In this study, we compared the structure of CPGs, identified the mutational hotspots and resolved the phylogenetic relationship of subgenus Gynopodium. RESULTS: The CPGs of six subgenus Gynopodium species ranged in size from 160,027 bp to 160,114 bp. A total of 131 genes were identified, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. We detected neither major expansions or contractions in the inverted repeat region, nor rearrangements or insertions in the CPGs of six subgenus Gynopodium species. A total of 300 large repeat sequences (forward, reverse, and palindrome repeats), 847 simple sequence repeats, and five highly variable regions were identified. One gene (ycf1) and four intergenic regions (psbA-trnH-GUG, petA-psbJ, rpl32-trnL-UAG, and ccsA-ndhD) were identified as mutational hotspots by their high nucleotide diversity (Pi) values (≥ 0.004), which were useful for species discrimination. Maximum likelihood and Bayesian inference trees were concordant and indicated that Magnoliaceae consisted of two genera Liriodendron and Magnolia. Six species of subgenus Gynopodium clustered as a monophyletic clade, forming a sister clade with subgenus Yulania (BS = 100%, PP = 1.00). Due to the non-monophyly of subgenus Magnolia, subgenus Gynopodium should be treated as a section of Magnolia. Within section Gynopodium, M. sinica diverged first (posterior probability = 1, bootstrap = 100), followed by M. nitida, M. kachirachirai and M. lotungensis. M. omeiensis was sister to M. yunnanensis (posterior probability = 0.97, bootstrap = 50). CONCLUSION: The CPGs and characteristics information provided by our study could be useful in species identification, conservation genetics and resolving phylogenetic relationships of Magnoliaceae species.

BACKGROUND: The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada. RESULTS: We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data. CONCLUSIONS: Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination.

BACKGROUND: Preeclampsia is a pregnancy-related condition that causes high blood pressure and proteinuria after 20 weeks of pregnancy. It is linked to increased maternal mortality, organ malfunction, and foetal development limitation. In this view, there is a need critical to identify biomarkers for the early detection of preeclampsia. The objective of this study is to discover critical genes and explore medications for preeclampsia treatment that may influence these genes. METHODS: Four datasets, including GSE10588, GSE25906, GSE48424 and GSE60438 were retrieved from the Gene Expression Omnibus database. The GSE10588, GSE25906, and GSE48424 datasets were then removed the batch effect using the "sva" R package and merged into a complete dataset. The differentially expressed genes (DEGs) were identified using the "limma" R package. The potential small-molecule agents for the treatment of PE was further screened using the Connective Map (CMAP) drug database based on the DEGs. Further, Weight gene Co-expression network (WGNCA) analysis was performed to identified gene module associated with preeclampsia, hub genes were then identified using the logistic regression analysis. Finally, the immune cell infiltration level of genes was evaluated through the single sample gene set enrichment analysis (ssGSEA). RESULTS: A total of 681 DEGs (376 down-regulated and 305 up-regulated genes) were identified between normal and preeclampsia samples. Then, Dexamethasone, Prednisone, Rimexolone, Piretanide, Trazodone, Buflomedil, Scoulerin, Irinotecan, and Camptothecin drugs were screened based on these DEGs through the CMAP database. Two modules including yellow and brown modules were the most associated with disease through the WGCNA analysis. KEGG analysis revealed that the chemokine signaling pathway, Th1 and Th2 cell differentiation, B cell receptor signalling pathway and oxytocin signalling pathway were significantly enriched in these modules. Moreover, two key genes, PLEK and LEP were evaluated using the univariate and multivariate logistic regression analysis from the hub modules. These two genes were further validated in the external validation cohort GSE60438 and qRT-PCR experiment. Finally, we evaluated the relationship between immune cell and two genes. CONCLUSION: In conclusion, the present study investigated key genes associated with PE pathogenesis that may contribute to identifying potential biomarkers, therapeutic agents and developing personalized treatment for PE.

BACKGROUND: Bacteria are an essential component of the earth`s biota and affect circulation of matters through their metabolic activity. They also play an important role in the carbon and nitrogen cycle in the deep-sea environment. In this paper, two strains from deep-sea sediments were investigated in order to understand nitrogen cycling involved in the deep-sea environment. RESULTS: In this paper, the basic genomic information of two strains was obtained by whole genome sequencing. The Cobetia amphilecti N-80 and Halomonas profundus 13 genome sizes are 4,160,095 bp with a GC content of 62.5% and 5,251,450 bp with a GC content of 54.84%. Through a comparison of functional analyses, we predicted the possible C and N metabolic pathways of the two strains and determined that Halomonas profundus 13 could use more carbon sources than Cobetia amphilecti N-80. The main genes associated with N metabolism in Halomonas profundus 13 are narG, narY, narI, nirS, norB, norC, nosZ, and nirD. On the contrast, nirD, using NH4+ for energy, plays a main role in Cobetia amphilecti N-80. Both of them have the same genes for fixing inorganic carbon: icd, ppc, fdhA, accC, accB, accD, and accA. CONCLUSION: In this study, the whole genomes of two strains were sequenced to clarify the basic characteristics of their genomes, laying the foundation for further studying nitrogen-metabolizing bacteria. Halomonas profundus 13 can utilize more carbon sources than Cobetia amphilecti N-80, as indicated by API as well as COG and KEGG prediction results. Finally, through the analysis of the nitrification and denitrification abilities as well as the inorganic carbon fixation ability of the two strains, the related genes were identified, and the possible metabolic pathways were predicted. Together, these results provide molecular markers and theoretical support for the mechanisms of inorganic carbon fixation by deep-sea microorganisms.

BACKGROUND: Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome. RESULTS: Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively. CONCLUSIONS: Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

BACKGROUND: Single primer enrichment technology (SPET) is an emerging and increasingly popular solution for high-throughput targeted genotyping in plants. Although SPET requires a priori identification of polymorphisms for probe design, this technology has potentially higher reproducibility and transferability compared to other reduced representation sequencing (RRS) approaches, also enabling the discovery of closely linked polymorphisms surrounding the target one. RESULTS: The potential for SPET application in fruit trees was evaluated by developing a 25K target SNPs assay to genotype a panel of apricot accessions and progenies. A total of 32,492 polymorphic sites were genotyped in 128 accessions (including 8,188 accessory non-target SNPs) with extremely low levels of missing data and a significant correlation of allelic frequencies compared to whole-genome sequencing data used for array design. Assay performance was further validated by estimating genotyping errors in two biparental progenies, resulting in an overall 1.8% rate. SPET genotyping data were used to infer population structure and to dissect the architecture of fruit maturity date (MD), a quantitative reproductive phenological trait of great agronomical interest in apricot species. Depending on the year, GWAS revealed loci associated to MD on several chromosomes. The QTLs on chromosomes 1 and 4 (the latter explaining most of the phenotypic variability in the panel) were the most consistent over years and were further confirmed by linkage mapping in two segregating progenies. CONCLUSIONS: Besides the utility for marker assisted selection and for paving the way to in-depth studies to clarify the molecular bases of MD trait variation in apricot, the results provide an overview of the performance and reliability of SPET for fruit tree genetics.

BACKGROUND: The sucrose non-fermenting 1 (SNF1)-related protein kinases (SnRKs) play a vivid role in regulating plant metabolism and stress response, providing a pathway for regulation between metabolism and stress signals. Conducting identification and stress response studies on SnRKs in plants contributes to the development of strategies for tree species that are more tolerant to stress conditions. RESULTS: In the present study, a total of 30 LcSnRKs were identified in Liriodendron chinense (L. chinense) genome, which was distributed across 15 chromosomes and 4 scaffolds. It could be divided into three subfamilies: SnRK1, SnRK2, and SnRK3 based on phylogenetic analysis and domain types. The LcSnRK of the three subfamilies shared the same Ser/Thr kinase structure in gene structure and motif composition, while the functional domains, except for the kinase domain, showed significant differences. A total of 13 collinear gene pairs were detected in L. chinense and Arabidopsis thaliana (A. thaliana), and 18 pairs were detected in L. chinense and rice, suggesting that the LcSnRK family genes may be evolutionarily more closely related to rice. Cis-regulation element analysis showed that LcSnRKs were LTR and TC-rich, which could respond to different environmental stresses. Furthermore, the expression patterns of LcSnRKs are different at different times under low-temperature stress. LcSnRK1s expression tended to be down-regulated under low-temperature stress. The expression of LcSnRK2s tended to be up-regulated under low-temperature stress. The expression trend of LcSnRK3s under low-temperature stress was mainly up-or down-regulated. CONCLUSION: The results of this study will provide valuable information for the functional identification of the LcSnRK gene in the future.

Erratum in BMC Genomics. 2023 Jan 11;24(1):15.

BACKGROUND: Livestock animals is of great significance in agricultural production. However, the role of specific gene expression, especially alternative splicing in determining phenotype, is not well understood. The livestock research community needs a gene expression and alternative splicing database contributing to livestock genetic improvement. DESCRIPTION: We report the construction of LivestockExp ( https://bioinfo.njau.edu.cn/livestockExp ), a web-based database server for the exploration of gene expression and alternative splicing using 43,710 uniformly processed RNA-seq samples from livestock animals and several relative species across six orders. The database is equipped with basic querying functions and multiple online analysis modules including differential/specific expression analysis, co-expression network analysis, and cross-species gene expression conservation analysis. In addition to the re-analysis of public datasets, users can upload personal datasets to perform co-analysis with public datasets. The database also offers a wide range of visualization tools and diverse links to external databases enabling users to efficiently explore the results and to gain additional insights. CONCLUSION: LivestockExp covers by far the largest number of livestock animal species and RNA-seq samples and provides a valuable data resource and analysis platform for the convenient utilization of public RNA-seq datasets.

BACKGROUND: Elongases of very long chain fatty acids (ELOVLs), a family of first rate-limiting enzymes in the synthesis of long-chain fatty acids, play an essential role in the biosynthesis of complex lipids. Disrupting any of ELOVLs affects normal growth and development in mammals. Genetic variations in ELOVLs are associated with backfat or intramuscular fatty acid composition in livestock. However, the effects of ELOVL gene family on breeding selection and lipid deposition in different tissues are still unknown in chickens. RESULTS: Genetic variation patterns and genetic associations analysis showed that the genetic variations of ELOVL genes were contributed to breeding selection of commercial varieties in chicken, and 14 SNPs in ELOVL2-6 were associated with body weight, carcass or fat deposition traits. Especially, one SNP rs17631638T > C in the promoter of ELOVL3 was associated with intramuscular fat content (IMF), and its allele frequency was significantly higher in native and layer breeds compared to that in commercial broiler breeds. Quantitative real-time PCR (qRT-PCR) determined that the ELOVL3 expressions in pectoralis were affected by the genotypes of rs17631638T > C. In addition, the transcription levels of ELOVL genes except ELOVL5 were regulated by estrogen in chicken liver and hypothalamus with different regulatory pathways. The expression levels of ELOVL1-6 in hypothalamus, liver, abdominal fat and pectoralis were correlated with abdominal fat weight, abdominal fat percentage, liver lipid content and IMF. Noteworthily, expression of ELOVL3 in pectoralis was highly positively correlated with IMF and glycerophospholipid molecules, including phosphatidyl choline, phosphatidyl ethanolamine, phosphatidyl glycerol and phospholipids inositol, rich in ω-3 and ω-6 long-chain unsaturated fatty acids, suggesting ELOVL3 could contribute to intramuscular fat deposition by increasing the proportion of long-chain unsaturated glycerophospholipid molecules in pectoralis. CONCLUSIONS: In summary, we demonstrated the genetic contribution of ELOVL gene family to breeding selection for specialized varieties, and revealed the expression regulation of ELOVL genes and their potential roles in regulating lipid deposition in different tissues. This study provides new insights into understanding the functions of ELOVL family on avian growth and lipid deposition in different tissues and the genetic variation in ELOVL3 may aid the marker-assisted selection of meat quality in chicken.

BACKGROUND: Biotin carboxyl carrier protein (BCCP) is a subunit of Acetyl CoA-carboxylase (ACCase) which catalyzes the conversion of acetyl-CoA to malonyl-CoA in a committed step during the de novo biosynthesis of fatty acids. Lipids, lipid metabolites, lipid-metabolizing and -modifying enzymes are known to play a role in biotic and abiotic stress tolerance in plants. In this regard, an understanding of the Brassica napus BCCP genes will aid in the improvement of biotic and abiotic stress tolerance in canola. RESULTS: In this study, we identified 43 BCCP genes in five Brassica species based on published genome data. Among them, Brassica rapa, Brassica oleracea, Brassica nigra, Brassica napus and Brassica juncea had six, seven, seven, 10 and 13 BCCP homologs, respectively. Phylogenetic analysis categorized them into five classes, each with unique conserved domains. The promoter regions of all BCCP genes contained stress-related cis-acting elements as determined by cis-element analysis. We identified four and three duplicated gene pairs (segmental) in B. napus and B. juncea respectively, indicating the role of segmental duplication in the expansion of this gene family. The Ka/Ks ratios of orthologous gene pairs between Arabidopsis thaliana and five Brassica species were mostly less than 1.0, implying that purifying selection, i.e., selective removal of deleterious alleles, played a role during the evolution of Brassica genomes. Analysis of 10 BnaBCCP genes using qRT-PCR showed a different pattern of expression because of exposure of the plants to biotic stresses, such as clubroot and sclerotinia diseases, and abiotic stresses such as drought, low temperature and salinity stresses. CONCLUSIONS: The identification and functional analysis of the Brassica BCCPs demonstrated that some of these genes might play important roles in biotic and abiotic stress responses. Results from this study could lay the foundation for a better understanding of these genes for the improvement of Brassica crops for stress tolerance.

BACKGROUND: The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines. RESULTS: We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria. CONCLUSIONS: Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

BACKGROUND: Influenza is a contagious disease that affects people of all ages and is linked to considerable mortality during epidemics and occasional outbreaks. Moreover, effective immunological biomarkers are needed for elucidating aetiology and preventing and treating severe influenza. Herein, we aimed to evaluate the key genes linked with the disease severity in influenza patients needing invasive mechanical ventilation (IMV). Three gene microarray data sets (GSE101702, GSE21802, and GSE111368) from blood samples of influenza patients were made available by the Gene Expression Omnibus (GEO) database. The GSE101702 and GSE21802 data sets were combined to create the training set. Hub indicators for IMV patients with severe influenza were determined using differential expression analysis and Weighted correlation network analysis (WGCNA) from the training set. The receiver operating characteristic curve (ROC) was also used to evaluate the hub genes from the test set's diagnostic accuracy. Different immune cells' infiltration levels in the expression profile and their correlation with hub gene markers were examined using single-sample gene set enrichment analysis (ssGSEA). RESULTS: In the present study, we evaluated a total of 447 differential genes. WGCNA identified eight co-expression modules, with the red module having the strongest correlation with IMV patients. Differential genes were combined to obtain 3 hub genes (HLA-DPA1, HLA-DRB3, and CECR1). The identified genes were investigated as potential indicators for patients with severe influenza who required IMV using the least absolute shrinkage and selection operator (LASSO) approach. The ROC showed the diagnostic value of the three hub genes in determining the severity of influenza. Using ssGSEA, it has been revealed that the expression of key genes was negatively correlated with neutrophil activation and positively associated with adaptive cellular immune response. CONCLUSION: We evaluated three novel hub genes that could be linked to the immunopathological mechanism of severe influenza patients who require IMV treatment and could be used as potential biomarkers for severe influenza prevention and treatment.

BACKGROUND: Cellular events during meiosis can differ between inbred lines in maize. Substantial differences in the average numbers of chiasmata and double-strand breaks (DSBs) per meiotic cell have been documented among diverse inbred lines of maize: CML228, a tropical maize inbred line, B73 and Mo17, temperate maize lines. To determine if gene expression might explain these observed differences, an RNA-Seq experiment was performed on CML228 male meiocytes which was compared to B73 and Mo17 male meiocytes, where plants were grown in the same controlled environment. RESULTS: We found that a few DSB-repair/meiotic genes which promote class I crossovers (COs) and the Zyp1 gene which limits newly formed class I COs were up-regulated, whereas Mus81 homolog 2 which promotes class II COs was down-regulated in CML228. Although we did not find enriched gene ontology (GO) categories directly related to meiosis, we found that GO categories in membrane, localization, proteolysis, energy processes were up-regulated in CML228, while chromatin remodeling, epigenetic regulation, and cell cycle related processes including meiosis related cell cycle processes were down-regulated in CML228. The degree of similarity in expression patterns between the three maize lines reflect their genetic relatedness: B73 and Mo17 had similar meiotic expressions and CML228 had a more distinct expression profile. CONCLUSIONS: We found that meiotic related genes were mostly conserved among the three maize inbreds except for a few DSB-repair/meiotic genes. The findings that the molecular players in limiting class I CO formation (once CO assurance is achieved) were up-regulated and those involved in promoting class II CO formation were down-regulated in CML228 agree with the lower chiasmata number observed in CML228 previously. In addition, epigenetics such as chromatin remodeling and histone modification might play a role. Transport and energy-related processes was up-regulated and Cyclin13 was down-regulated in CML228. The direction of gene expression of these processes agree with that previously found in meiotic tissues compared with vegetative tissues. In summary, we used different natural maize inbred lines from different climatic conditions and have shown their differences in expression landscape in male meiocytes.

Long non-coding RNAs (lncRNAs) play a vital role in regulating adipogenesis. However, the associated regulatory mechanisms have yet to be described in detail in pig. In this study, we demonstrate a critical role for lncMYOZ2 in adipogenesis from porcine preadipocytes. Specifically, lncMYOZ2 was more abundant in the adipose tissue of Mashen (fat-type) pigs than for Large White (lean-type) pigs, and knockdown of this lncRNA significantly inhibited the differentiation of porcine preadipocytes into adipocytes. Mechanistically, we used RNA pull-down and RIP assays to establish that lncMYOZ2 interacts with adenosylhomocysteinase (AHCY). Moreover, lncMYOZ2 knockdown increased promoter methylation of the target gene MYOZ2 and lowered its expression. Finally, we describe a positive regulatory role for MYOZ2 in adipogenesis. Collectively, these findings establish lncMYOZ2 as an important epigenetic regulator of adipogenesis via the aforementioned AHCY/MYOZ2 pathway, and provide insights into the role of lncRNAs in porcine adipose development.

BACKGROUND: Forest musk deer is an endangered species globally. The death of captive forest musk deer can be caused by certain respiratory system diseases. Acute respiratory distress syndrome (ARDS) is a huge threat to the life of forest muck deer that breed in our department. METHODS: Lung histopathologic analysis was conducted by hematoxylin and eosin (HE) staining. The lung gene changes triggered by ARDS were examined by RNA sequencing and related bioinformatics analysis in forest musk deer. The potential functions of unigenes were investigated by NR, SwissProt KOG, GO, and KEGG annotation analyses. Vital biological processes or pathways in ARDS were examined by GO and KEGG enrichment analyses. RESULTS: A total of 3265 unigenes were differentially expressed (|log2fold-change|> 2 and adjusted P value < 0.01) in lung tissues of 3 forest musk deer with ARDS compared with normal lung tissues of the non-ARDS group. These differentially expressed unigenes (DEGs) played crucial roles in immunity and defense responses to pathogens. Moreover, we identified the DEGs related to one or more of the following biological processes: lung development, immunity, and bacterial/viral/fungal infection. And six DEGs that might be involved in lung injury caused by immune dysregulation or viral/fungal infection were identified. CONCLUSION: ARDS-mediated lung gene alterations were identified in forest musk deer. Moreover, multiple genes involved in lung development and lung defense responses to bacteria/viruses/fungi in ARDS were filtered out in forest musk deer.

BACKGROUND: One of the most complex prokaryotic organelles are magnetosomes, which are formed by magnetotactic bacteria as sensors for navigation in the Earth's magnetic field. In the alphaproteobacterium Magnetospirillum gryphiswaldense magnetosomes consist of chains of magnetite crystals (Fe3O4) that under microoxic to anoxic conditions are biomineralized within membrane vesicles. To form such an intricate structure, the transcription of > 30 specific structural genes clustered within the genomic magnetosome island (MAI) has to be coordinated with the expression of an as-yet unknown number of auxiliary genes encoding several generic metabolic functions. However, their global regulation and transcriptional organization in response to anoxic conditions most favorable for magnetite biomineralization are still unclear. RESULTS: Here, we compared transcriptional profiles of anaerobically grown magnetosome forming cells with those in which magnetosome biosynthesis has been suppressed by aerobic condition. Using whole transcriptome shotgun sequencing, we found that transcription of about 300 of the > 4300 genes was significantly enhanced during magnetosome formation. About 40 of the top upregulated genes are directly or indirectly linked to aerobic and anaerobic respiration (denitrification) or unknown functions. The mam and mms gene clusters, specifically controlling magnetosome biosynthesis, were highly transcribed, but constitutively expressed irrespective of the growth condition. By Cappable-sequencing, we show that the transcriptional complexity of both the MAI and the entire genome decreased under anaerobic conditions optimal for magnetosome formation. In addition, predominant promoter structures were highly similar to sigma factor σ70 dependent promoters in other Alphaproteobacteria. CONCLUSIONS: Our transcriptome-wide analysis revealed that magnetite biomineralization relies on a complex interplay between generic metabolic processes such as aerobic and anaerobic respiration, cellular redox control, and the biosynthesis of specific magnetosome structures. In addition, we provide insights into global regulatory features that have remained uncharacterized in the widely studied model organism M. gryphiswaldense, including a comprehensive dataset of newly annotated transcription start sites and genome-wide operon detection as a community resource (GEO Series accession number GSE197098).

BACKGROUND: The importance of fathers' engagement in care and its critical role in the offspring's cognitive and emotional development is now well established. Yet, little is known on the underlying neurobiology due to the lack of appropriate animal models. In the socially monogamous and bi-parental prairie vole (Microtus ochrogaster), while 60-80% of virgin males show spontaneous paternal behaviors (Paternal), others display pup-directed aggression (Attackers). Here we took advantage of this phenotypic dichotomy and used RNA-sequencing in three important brain areas to characterize gene expression associated with paternal behaviors of Paternal males and compare it to experienced Fathers and Mothers. RESULTS: While Paternal males displayed the same range and extent of paternal behaviors as experienced Fathers, we observed structure-specific transcriptomic differences between parental behaviors phenotypes. Using differential expression, gene set expression, as well as co-expression network analyses, we found that phenotypic differences between Paternal males and Attackers were mainly reflected by the lateral septum (LS), and to a lower extent, the nucleus accumbens (NAc), transcriptomes. In the medial preoptic area (MPOA), the profiles of gene expression mainly reflected differences between females and males regardless of their parental behaviors phenotype. Functional enrichment analyses of those gene sets associated with Paternal males or Attackers in the LS and the NAc revealed the involvement of processes related to the mitochondria, RNA translation, protein degradation processes, as well as epigenetic regulation of gene expression. CONCLUSIONS: By leveraging the natural phenotypic differences in parental behaviors in virgin male prairie voles alongside fathers and mothers, we identified a marked structure- and phenotype-specific pattern of gene expression associated with spontaneous paternal behaviors independently from fatherhood and pair-bonding. The LS transcriptome related to the mitochondria, RNA translation, and protein degradation processes was thus highlighted as a primary candidate associated with the spontaneous display of paternal behaviors. Altogether, our observations further characterize the behavioral and transcriptomic signature of parental behaviors in the socially monogamous prairie vole and lay the groundwork to further our understanding of the molecular underpinnings of paternal behavior.

BACKGROUND: Considering relatives' health history in logistic regression for case-control genome-wide association studies (CC-GWAS) may provide new information that increases accuracy and power to detect disease associated genetic variants. We conducted simulations and analyzed type 2 diabetes (T2D) data from the Framingham Heart Study (FHS) to compare two methods, liability threshold model conditional on both case-control status and family history (LT-FH) and Fam-meta, which incorporate family history into CC-GWAS. RESULTS: In our simulation scenario of trait with modest T2D heritability (h2 = 0.28), variant minor allele frequency ranging from 1% to 50%, and 1% of phenotype variance explained by the genetic variants, Fam-meta had the highest overall power, while both methods incorporating family history were more powerful than CC-GWAS. All three methods had controlled type I error rates, while LT-FH was the most conservative with a lower-than-expected error rate. In addition, we observed a substantial increase in power of the two familial history methods compared to CC-GWAS when the prevalence of the phenotype increased with age. Furthermore, we showed that, when only the phenotypes of more distant relatives were available, Fam-meta still remained more powerful than CC-GWAS, confirming that leveraging disease history of both close and distant relatives can increase power of association analyses. Using FHS data, we confirmed the well-known association of TCF7L2 region with T2D at the genome-wide threshold of P-value < 5 × 10-8, and both familial history methods increased the significance of the region compared to CC-GWAS. We identified two loci at 5q35 (ADAMTS2) and 5q23 (PRR16), not previously reported for T2D using CC-GWAS and Fam-meta; both genes play a role in cardiovascular diseases. Additionally, CC-GWAS detected one more significant locus at 13q31 (GPC6) reported associated with T2D-related traits. CONCLUSIONS: Overall, LT-FH and Fam-meta had higher power than CC-GWAS in simulations, especially using phenotypes that were more prevalent in older age groups, and both methods detected known genetic variants with lower P-values in real data application, highlighting the benefits of including family history in genetic association studies.