BACKGROUND: With the development of third-generation sequencing (TGS) technologies, people are able to obtain DNA sequences with lengths from 10s to 100s of kb. These long reads allow protein domain annotation without assembly, thus can produce important insights into the biological functions of the underlying data. However, the high error rate in TGS data raises a new challenge to established domain analysis pipelines. The state-of-the-art methods are not optimized for noisy reads and have shown unsatisfactory accuracy of domain classification in TGS data. New computational methods are still needed to improve the performance of domain prediction in long noisy reads. RESULTS: In this work, we introduce ProDOMA, a deep learning model that conducts domain classification for TGS reads. It uses deep neural networks with 3-frame translation encoding to learn conserved features from partially correct translations. In addition, we formulate our problem as an open-set problem and thus our model can reject reads not containing the targeted domains. In the experiments on simulated long reads of protein coding sequences and real TGS reads from the human genome, our model outperforms HMMER and DeepFam on protein domain classification. CONCLUSIONS: In summary, ProDOMA is a useful end-to-end protein domain analysis tool for long noisy reads without relying on error correction.

BACKGROUND: Aedes aegypti mosquito, the principal global vector of arboviral diseases, lays eggs and undergoes larval and pupal development to become adult mosquitoes in fresh water (FW). It has recently been observed to develop in coastal brackish water (BW) habitats of up to 50% sea water, and such salinity tolerance shown to be an inheritable trait. Genomics of salinity tolerance in Ae. aegypti has not been previously studied, but it is of fundamental biological interest and important for controlling arboviral diseases in the context of rising sea levels increasing coastal ground water salinity. RESULTS: BW- and FW-Ae. aegypti were compared by RNA-seq analysis on the gut, anal papillae and rest of the carcass in fourth instar larvae (L4), proteomics of cuticles shed when L4 metamorphose into pupae, and transmission electron microscopy of cuticles in L4 and adults. Genes for specific cuticle proteins, signalling proteins, moulting hormone-related proteins, membrane transporters, enzymes involved in cuticle metabolism, and cytochrome P450 showed different mRNA levels in BW and FW L4 tissues. The salinity-tolerant Ae. aegypti were also characterized by altered L4 cuticle proteomics and changes in cuticle ultrastructure of L4 and adults. CONCLUSIONS: The findings provide new information on molecular and ultrastructural changes associated with salinity adaptation in FW mosquitoes. Changes in cuticles of larvae and adults of salinity-tolerant Ae. aegypti are expected to reduce the efficacy of insecticides used for controlling arboviral diseases. Expansion of coastal BW habitats and their neglect for control measures facilitates the spread of salinity-tolerant Ae. aegypti and genes for salinity tolerance. The transmission of arboviral diseases can therefore be amplified in multiple ways by salinity-tolerant Ae. aegypti and requires appropriate mitigating measures. The findings in Ae. aegypti have attendant implications for the development of salinity tolerance in other fresh water mosquito vectors and the diseases they transmit.

BACKGROUND: Plant transitions to land require robust cell walls for regulatory adaptations and to resist changing environments. Cell walls provide essential plasticity for plant cell division and defense, which are often conferred by the expansin superfamily with cell wall-loosening functions. However, the evolutionary mechanisms of expansin during plant terrestrialization are unclear. RESULTS: Here, we identified 323 expansin proteins in 12 genomes from algae to angiosperms. Phylogenetic evolutionary, structural, motif gain and loss and Ka/Ks analyses indicated that highly conserved expansin proteins were already present in algae and expanded and purified after plant terrestrialization. We found that the expansion of the FtEXPA subfamily was caused by duplication events and that the functions of certain duplicated genes may have differentiated. More importantly, we generated space-time expression profiles and finally identified five differentially expressed FtEXPs in both large and small fruit Tartary buckwheat that may regulate fruit size by responding to indoleacetic acid. CONCLUSIONS: A total of 323 expansin proteins from 12 representative plants were identified in our study during terrestrialization, and the expansin family that originated from algae expanded rapidly after the plants landed. The EXPA subfamily has more members and conservative evolution in angiosperms. FtEXPA1, FtEXPA11, FtEXPA12, FtEXPA19 and FtEXPA24 can respond to indole-3-acetic acid (IAA) signals and regulate fruit development. Our study provides a blueprint for improving the agronomic traits of Tartary buckwheat and a reference for defining the evolutionary history of the expansin family during plant transitions to land.

BACKGROUND: Infectious bronchitis virus (IBV) is one of the most relevant viruses affecting the poultry industry, and several studies have investigated the factors involved in its biological cycle and evolution. However, very few of those studies focused on the effect of genome composition and the codon bias of different IBV proteins, despite the remarkable increase in available complete genomes. In the present study, all IBV complete genomes were downloaded (n = 383), and several statistics representative of genome composition and codon bias were calculated for each protein-coding sequence, including but not limited to, the nucleotide odds ratio, relative synonymous codon usage and effective number of codons. Additionally, viral codon usage was compared to host codon usage based on a collection of highly expressed genes in IBV target and nontarget tissues. RESULTS: The results obtained demonstrated a significant difference among structural, non-structural and accessory proteins, especially regarding dinucleotide composition, which appears under strong selective forces. In particular, some dinucleotide pairs, such as CpG, a probable target of the host innate immune response, are underrepresented in genes coding for pp1a, pp1ab, S and N. Although genome composition and dinucleotide bias appear to affect codon usage, additional selective forces may act directly on codon bias. Variability in relative synonymous codon usage and effective number of codons was found for different proteins, with structural proteins and polyproteins being more adapted to the codon bias of host target tissues. In contrast, accessory proteins had a more biased codon usage (i.e., lower number of preferred codons), which might contribute to the regulation of their expression level and timing throughout the cell cycle. CONCLUSIONS: The present study confirms the existence of selective forces acting directly on the genome and not only indirectly through phenotype selection. This evidence might help understanding IBV biology and in developing attenuated strains without affecting the protein phenotype and therefore immunogenicity.

Erratum in BMC Genomics. 2021 Apr 28;22(1):306.

BACKGROUND: Transposable elements are significant components of most organism's genomes, yet the reasons why their abundances vary significantly among species is poorly understood. A recent study has suggested that even in the absence of traditional molecular evolutionary explanations, transposon proliferation may occur through a process known as 'transposon engineering'. However, their model used a fixed beneficial transposon insertion frequency of 20%, which we believe to be unrealistically high. RESULTS: Reducing this beneficial insertion frequency, while keeping all other parameters identical, prevented transposon proliferation. CONCLUSIONS: We conclude that the author's original findings are better explained through the action of positive selection rather than 'transposon engineering', with beneficial insertion effects remaining important during transposon proliferation events.

BACKGROUND: Single nucleotide polymorphism (SNP) array technology has been increasingly used to generate large quantities of SNP data for use in genetic studies. As new arrays are developed to take advantage of new technology and of improved probe design using new genome sequence and panel data, a need to integrate data from different arrays and array platforms has arisen. This study was undertaken in view of our need for an integrated high-quality dataset of Illumina Infinium® 20 K and Affymetrix Axiom® 480 K SNP array data in apple (Malus × domestica). In this study, we qualify and quantify the compatibility of SNP calling, defined as SNP calls that are both accurate and concordant, across both arrays by two approaches. First, the concordance of SNP calls was evaluated using a set of 417 duplicate individuals genotyped on both arrays starting from a set of 10,295 robust SNPs on the Infinium array. Next, the accuracy of the SNP calls was evaluated on additional germplasm (n = 3141) from both arrays using Mendelian inconsistent and consistent errors across thousands of pedigree links. While performing this work, we took the opportunity to evaluate reasons for probe failure and observed discordant SNP calls. RESULTS: Concordance among the duplicate individuals was on average of 97.1% across 10,295 SNPs. Of these SNPs, 35% had discordant call(s) that were further curated, leading to a final set of 8412 (81.7%) SNPs that were deemed compatible. Compatibility was highly influenced by the presence of alternate probe binding locations and secondary polymorphisms. The impact of the latter was highly influenced by their number and proximity to the 3' end of the probe. CONCLUSIONS: The Infinium and Axiom SNP array data were mostly compatible. However, data integration required intense data filtering and curation. This work resulted in a workflow and information that may be of use in other data integration efforts. Such an in-depth analysis of array concordance and accuracy as ours has not been previously described in the literature and will be useful in future work on SNP array data integration and interpretation, and in probe/platform development.

BACKGROUND: The blood feeding poultry red mite (PRM), Dermanyssus gallinae, causes substantial economic damage to the egg laying industry worldwide, and is a serious welfare concern for laying hens and poultry house workers. In this study we have investigated the temporal gene expression across the 6 stages/sexes (egg, larvae, protonymph and deutonymph, adult male and adult female) of this neglected parasite in order to understand the temporal expression associated with development, parasitic lifestyle, reproduction and allergen expression. RESULTS: RNA-seq transcript data for the 6 stages were mapped to the PRM genome creating a publicly available gene expression atlas (on the OrcAE platform in conjunction with the PRM genome). Network analysis and clustering of stage-enriched gene expression in PRM resulted in 17 superclusters with stage-specific or multi-stage expression profiles. The 6 stage specific superclusters were clearly demarked from each other and the adult female supercluster contained the most stage specific transcripts (2725), whilst the protonymph supercluster the fewest (165). Fifteen pairwise comparisons performed between the different stages resulted in a total of 6025 Differentially Expressed Genes (DEGs) (P > 0.99). These data were evaluated alongside a Venn/Euler analysis of the top 100 most abundant genes in each stage. An expanded set of cuticle proteins and enzymes (chitinase and metallocarboxypeptidases) were identified in larvae and underpin cuticle formation and ecdysis to the protonymph stage. Two mucin/peritrophic-A salivary proteins (DEGAL6771g00070, DEGAL6824g00220) were highly expressed in the blood-feeding stages, indicating peritrophic membrane formation during feeding. Reproduction-associated vitellogenins were the most abundant transcripts in adult females whilst, in adult males, an expanded set of serine and cysteine proteinases and an epididymal protein (DEGAL6668g00010) were highly abundant. Assessment of the expression patterns of putative homologues of 32 allergen groups from house dust mites indicated a bias in their expression towards the non-feeding larval stage of PRM. CONCLUSIONS: This study is the first evaluation of temporal gene expression across all stages of PRM and has provided insight into developmental, feeding, reproduction and survival strategies employed by this mite. The publicly available PRM resource on OrcAE offers a valuable tool for researchers investigating the biology and novel interventions of this parasite.

BACKGROUND: Aggressive behavior is an ancient and conserved trait, habitual for most animals in order to eat, protect themselves, compete for mating and defend their territories. Genetic factors have been shown to play an important role in the development of aggression both in animals and humans, displaying moderate to high heritability estimates. Although such types of behaviors have been studied in different animal models, the molecular architecture of aggressiveness remains poorly understood. This study compared gene expression profiles of 16 prefrontal cortex (PFC) samples from aggressive and non-aggressive cattle breeds: Lidia, selected for agonistic responses, and Wagyu, selected for tameness. RESULTS: A total of 918 up-regulated and 278 down-regulated differentially expressed genes (DEG) were identified, representing above-chance overlap with genes previously identified in studies of aggression across species, as well as those implicated in recent human evolution. The functional interpretation of the up-regulated genes in the aggressive cohort revealed enrichment of pathways such as Alzheimer disease-presenilin, integrins and the ERK/MAPK signaling cascade, all implicated in the development of abnormal aggressive behaviors and neurophysiological disorders. Moreover, gonadotropins, are up-regulated as natural mechanisms enhancing aggression. Concomitantly, heterotrimeric G-protein pathways, associated with low reactivity mental states, and the GAD2 gene, a repressor of agonistic reactions associated with PFC activity, are down-regulated, promoting the development of the aggressive responses selected for in Lidia cattle. We also identified six upstream regulators, whose functional activity fits with the etiology of abnormal behavioral responses associated with aggression. CONCLUSIONS: These transcriptional correlates of aggression, resulting, at least in part, from controlled artificial selection, can provide valuable insights into the complex architecture that underlies naturally developed agonistic behaviors. This analysis constitutes a first important step towards the identification of the genes and metabolic pathways that promote aggression in cattle and, providing a novel model species to disentangle the mechanisms underlying variability in aggressive behavior.

BACKGROUND: Antibiotic-producing Streptomyces bacteria are ubiquitous in nature, yet most studies of its diversity have focused on free-living strains inhabiting diverse soil environments and those in symbiotic relationship with invertebrates. RESULTS: We studied the draft genomes of 73 Streptomyces isolates sampled from the skin (wing and tail membranes) and fur surfaces of bats collected in Arizona and New Mexico. We uncovered large genomic variation and biosynthetic potential, even among closely related strains. The isolates, which were initially identified as three distinct species based on sequence variation in the 16S rRNA locus, could be distinguished as 41 different species based on genome-wide average nucleotide identity. Of the 32 biosynthetic gene cluster (BGC) classes detected, non-ribosomal peptide synthetases, siderophores, and terpenes were present in all genomes. On average, Streptomyces genomes carried 14 distinct classes of BGCs (range = 9-20). Results also revealed large inter- and intra-species variation in gene content (single nucleotide polymorphisms, accessory genes and singletons) and BGCs, further contributing to the overall genetic diversity present in bat-associated Streptomyces. Finally, we show that genome-wide recombination has partly contributed to the large genomic variation among strains of the same species. CONCLUSIONS: Our study provides an initial genomic assessment of bat-associated Streptomyces that will be critical to prioritizing those strains with the greatest ability to produce novel antibiotics. It also highlights the need to recognize within-species variation as an important factor in genetic manipulation studies, diversity estimates and drug discovery efforts in Streptomyces.

BACKGROUND: Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far. RESULTS: The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation. CONCLUSIONS: The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

BACKGROUND: Cucumber (Cucumis sativus L.) is cultivated worldwide, and it is essential to produce enough high-quality seeds to meet demand. Pre-harvest sprouting (PHS) in cucumber is a critical problem and causes serious damage to seed production and quality. Nevertheless, the genetic basis and molecular mechanisms underlying cucumber PHS remain unclear. QTL-seq is an efficient approach for rapid quantitative trait loci (QTL) identification that simultaneously takes advantage of bulked-segregant analysis (BSA) and whole-genome resequencing. In the present research, QTL-seq analysis was performed to identify QTLs associated with PHS in cucumber using an F2 segregating population. RESULTS: Two QTLs that spanned 7.3 Mb on Chromosome 4 and 0.15 Mb on Chromosome 5 were identified by QTL-seq and named qPHS4.1 and qPHS5.1, respectively. Subsequently, SNP and InDel markers selected from the candidate regions were used to refine the intervals using the extended F2 populations grown in the 2016 and 2017 seasons. Finally, qPHS4.1 was narrowed to 0.53 Mb on chromosome 4 flanked by the markers SNP-16 and SNP-24 and was found to explain 19-22% of the phenotypic variation in cucumber PHS. These results reveal that qPHS4.1 is a major-effect QTL associated with PHS in cucumber. Based on gene annotations and qRT-PCR expression analyses, Csa4G622760 and Csa4G622800 were proposed as the candidate genes. CONCLUSIONS: These results provide novel insights into the genetic mechanism controlling PHS in cucumber and highlight the potential for marker-assisted selection of PHS resistance breeding.

Erratum in BMC Genomics. 2021 May 19;22(1):365.

BACKGROUND: The behavior of cells in vivo is complex and highly dynamic, as it results from an interplay between intercellular matrix proteins with surface receptors and other microenvironmental cues. Although the effects of the cellular niche have been investigated for a number of cell types using different molecular approaches, comprehensive assessments of how the global transcriptome responds to 3D scaffolds composed of various extracellular matrix (ECM) constituents at different concentrations are still lacking. RESULTS: In this study, we explored the effects of two diverse extracellular matrix (ECM) components, Collagen I and Matrigel, on the transcriptional profile of cells in a cell culture system. Culturing Huh-7 cells on traditional cell culture plates (Control) or on the ECM components at different concentrations to modulate microenvironment properties, we have generated transcriptomics data that may be further explored to understand the differentiation and growth potential of this cell type for the development of 3D cultures. Our analysis infers transcription factors that are most responsible for the transcriptome response to the extracellular cues. CONCLUSION: Our data indicates that the Collagen I substrate induces a robust transcriptional response in the Huh-7 cells, distinct from that induced by Matrigel. Enhanced hepatocyte markers (ALB and miR-122) reveal a potentially robust remodelling towards primary hepatocytes. Our results aid in defining the appropriate culture and transcription pathways while using hepatoma cell lines. As systems mimicking the in vivo structure and function of liver cells are still being developed, our study could potentially circumvent bottlenecks of limited availability of primary hepatocytes for preclinical studies of drug targets.

BACKGROUND: Hypoxic stress plays a critical role in the persistence of Mycobacterium tuberculosis (Mtb) infection, but the mechanisms underlying this adaptive response remain ill defined. MATERIAL AND METHODS: In this study, using M. marinum as a surrogate, we analyzed hypoxic responses at the transcriptional level by Cappable-seq and regular RNA-seq analyses. RESULTS: A total of 6808 transcriptional start sites (TSSs) were identified under normoxic and hypoxic conditions. Among these TSSs, 1112 were upregulated and 1265 were downregulated in response to hypoxic stress. Using SigE-recognized consensus sequence, we identified 59 SigE-dependent promoters and all were upregulated under hypoxic stress, suggesting an important role for SigE in this process. We also compared the performance of Cappable-seq and regular RNA-seq using the same RNA samples collected from normoxic and hypoxic conditions, and confirmed that Cappable-seq is a valuable approach for global transcriptional regulation analyses. CONCLUSIONS: Our results provide insights and information for further characterization of responses to hypoxia in mycobacteria, and prove that Cappable-seq is a valuable approach for global transcriptional studies in mycobacteria.

BACKGROUND: Vernalization is a type of low temperature stress used to promote rapid bolting and flowering in plants. Although rapid bolting and flowering promote the reproduction of Chinese cabbages (Brassica rapa L. ssp. pekinensis), this process causes their commercial value to decline. Clarifying the mechanisms of vernalization is essential for its further application. We performed RNA sequencing of gradient-vernalization in order to explore the reasons for the different bolting process of two Chinese cabbage accessions during vernalization. RESULTS: There was considerable variation in gene expression between different-bolting Chinese cabbage accessions during vernalization. Comparative transcriptome analysis and weighted gene co-expression network analysis (WGCNA) were performed for different-bolting Chinese cabbage during different vernalization periods. The biological function analysis and hub gene annotation of highly relevant modules revealed that shoot system morphogenesis and polysaccharide and sugar metabolism caused early-bolting 'XBJ' to bolt and flower faster; chitin, ABA and ethylene-activated signaling pathways were enriched in late-bolting 'JWW'; and leaf senescence and carbohydrate metabolism enrichment were found in the two Chinese cabbage-related modules, indicating that these pathways may be related to bolting and flowering. The high connectivity of hub genes regulated vernalization, including MTHFR2, CPRD49, AAP8, endoglucanase 10, BXLs, GATLs, and WRKYs. Additionally, five genes related to flower development, BBX32 (binds to the FT promoter), SUS1 (increases FT expression), TSF (the closest homologue of FT), PAO and NAC029 (plays a role in leaf senescence), were expressed in the two Chinese cabbage accessions. CONCLUSION: The present work provides a comprehensive overview of vernalization-related gene networks in two different-bolting Chinese cabbages during vernalization. In addition, the candidate pathways and hub genes related to vernalization identified here will serve as a reference for breeders in the regulation of Chinese cabbage production.

BACKGROUND: For multicellular organisms, much remains unknown about the dynamics of synonymous codon and amino acid use in highly expressed genes, including whether their use varies with expression in different tissue types and sexes. Moreover, specific codons and amino acids may have translational functions in highly transcribed genes, that largely depend on their relationships to tRNA gene copies in the genome. However, these relationships and putative functions are poorly understood, particularly in multicellular systems. RESULTS: Here, we studied codon and amino acid use in highly expressed genes from reproductive and nervous system tissues (male and female gonad, somatic reproductive system, brain and ventral nerve cord, and male accessory glands) in the cricket Gryllus bimaculatus. We report an optimal codon, defined as the codon preferentially used in highly expressed genes, for each of the 18 amino acids with synonymous codons in this organism. The optimal codons were mostly shared among tissue types and both sexes. However, the frequency of optimal codons was highest in gonadal genes. Concordant with translational selection, a majority of the optimal codons had abundant matching tRNA gene copies in the genome, but sometimes obligately required wobble tRNAs. We suggest the latter may comprise a mechanism for slowing translation of abundant transcripts, particularly for cell-cycle genes. Non-optimal codons, defined as those least commonly used in highly transcribed genes, intriguingly often had abundant tRNAs, and had elevated use in a subset of genes with specialized functions (gametic and apoptosis genes), suggesting their use promotes the translational upregulation of particular mRNAs. In terms of amino acids, we found evidence suggesting that amino acid frequency, tRNA gene copy number, and amino acid biosynthetic costs (size/complexity) had all interdependently evolved in this insect model, potentially for translational optimization. CONCLUSIONS: Collectively, the results suggest a model whereby codon use in highly expressed genes, including optimal, wobble, and non-optimal codons, and their tRNA abundances, as well as amino acid use, have been influenced by adaptation for various functional roles in translation within this cricket. The effects of expression in different tissue types and the two sexes are discussed.

BACKGROUND: Spiroplasma is a widely distributed endosymbiont of insects, arthropods, and plants. In insects, Spiroplasma colonizes the gut, hemolymph, and reproductive organs of the host. Previous metagenomic surveys of the domesticated carmine cochineal Dactylopius coccus and the wild cochineal D. opuntiae reported sequences of Spiroplasma associated with these insects. However, there is no analysis of the genomic capabilities and the interaction of this Spiroplasma with Dactylopius. RESULTS: Here we present three Spiroplasma genomes independently recovered from metagenomes of adult males and females of D. coccus, from two different populations, as well as from adult females of D. opuntiae. Single-copy gene analysis showed that these genomes were > 92% complete. Phylogenomic analyses classified these genomes as new members of Spiroplasma ixodetis. Comparative genome analysis indicated that they exhibit fewer genes involved in amino acid and carbon catabolism compared to other spiroplasmas. Moreover, virulence factor-encoding genes (i.e., glpO, spaid and rip2) were found incomplete in these S. ixodetis genomes. We also detected an enrichment of genes encoding the type IV secretion system (T4SS) in S. ixodetis genomes of Dactylopius. A metratranscriptomic analysis of D. coccus showed that some of these T4SS genes (i.e., traG, virB4 and virD4) in addition to the superoxide dismutase sodA of S. ixodetis were overexpressed in the ovaries. CONCLUSION: The symbiont S. ixodetis is a new member of the bacterial community of D. coccus and D. opuntiae. The recovery of incomplete virulence factor-encoding genes in S. ixodetis of Dactylopius suggests that this bacterium is a non-pathogenic symbiont. A high number of genes encoding the T4SS, in the S. ixodetis genomes and the overexpression of these genes in the ovary and hemolymph of the host suggest that S. ixodetis use the T4SS to interact with the Dactylopius cells. Moreover, the transcriptional differences of S. ixodetis among the gut, hemolymph and ovary tissues of D. coccus indicate that this bacterium can respond and adapt to the different conditions (e.g., oxidative stress) present within the host. All this evidence proposes that there is a strong interaction and molecular signaling in the symbiosis between S. ixodetis and the carmine cochineal Dactylopius.

BACKGROUND: Callithrix marmosets are a relatively young primate radiation, whose phylogeny is not yet fully resolved. These primates are naturally para- and allopatric, but three species with highly invasive potential have been introduced into the southeastern Brazilian Atlantic Forest by the pet trade. There, these species hybridize with each other and endangered, native congeners. We aimed here to reconstruct a robust Callithrix phylogeny and divergence time estimates, and identify the biogeographic origins of autochthonous and allochthonous Callithrix mitogenome lineages. We sequenced 49 mitogenomes from four species (C. aurita, C. geoffroyi, C. jacchus, C. penicillata) and anthropogenic hybrids (C. aurita x Callithrix sp., C. penicillata x C. jacchus, Callithrix sp. x Callithrix sp., C. penicillata x C. geoffroyi) via Sanger and whole genome sequencing. We combined these data with previously published Callithrix mitogenomes to analyze five Callithrix species in total. RESULTS: We report the complete sequence and organization of the C. aurita mitogenome. Phylogenetic analyses showed that C. aurita was the first to diverge within Callithrix 3.54 million years ago (Ma), while C. jacchus and C. penicillata lineages diverged most recently 0.5 Ma as sister clades. MtDNA clades of C. aurita, C. geoffroyi, and C. penicillata show intraspecific geographic structure, but C. penicillata clades appear polyphyletic. Hybrids, which were identified by phenotype, possessed mainly C. penicillata or C. jacchus mtDNA haplotypes. The biogeographic origins of mtDNA haplotypes from hybrid and allochthonous Callithrix were broadly distributed across natural Callithrix ranges. Our phylogenetic results also evidence introgression of C. jacchus mtDNA into C. aurita. CONCLUSION: Our robust Callithrix mitogenome phylogeny shows C. aurita lineages as basal and C. jacchus lineages among the most recent within Callithrix. We provide the first evidence that parental mtDNA lineages of anthropogenic hybrid and allochthonous marmosets are broadly distributed inside and outside of the Atlantic Forest. We also show evidence of cryptic hybridization between allochthonous Callithrix and autochthonous C. aurita. Our results encouragingly show that further development of genomic resources will allow to more clearly elucidate Callithrix evolutionary relationships and understand the dynamics of Callithrix anthropogenic introductions into the Brazilian Atlantic Forest.

BACKGROUND: Colossoma macropomum (tambaqui) and Piaractus mesopotamicus (pacu) are good fish species for aquaculture. The tambacu, individuals originating from the induced hybridization of the female tambaqui with the male pacu, present rapid growth and robustness, characteristics which have made the tambacu a good choice for Brazilian fish farms. Here, we used small RNA sequencing to examine global miRNA expression in the genotypes pacu (PC), tambaqui (TQ), and hybrid tambacu (TC), (Juveniles, n = 5 per genotype), to better understand the relationship between tambacu and its parental species, and also to clarify the mechanisms involved in tambacu muscle growth and maintenance based on miRNAs expression. RESULTS: Regarding differentially expressed (DE) miRNAs between the three genotypes, we observed 8 upregulated and 7 downregulated miRNAs considering TC vs. PC; 14 miRNAs were upregulated and 10 were downregulated considering TC vs. TQ, and 15 miRNAs upregulated and 9 were downregulated considering PC vs. TQ. The majority of the miRNAs showed specific regulation for each genotype pair, and no miRNA were shared between the 3 genotype pairs, in both up- and down-regulated miRNAs. Considering only the miRNAs with validated target genes, we observed the miRNAs miR-144-3p, miR-138-5p, miR-206-3p, and miR-499-5p. GO enrichment analysis showed that the main target genes for these miRNAs were grouped in pathways related to oxygen homeostasis, blood vessel modulation, and oxidative metabolism. CONCLUSIONS: Our global miRNA analysis provided interesting DE miRNAs in the skeletal muscle of pacu, tambaqui, and the hybrid tambacu. In addition, in the hybrid tambacu, we identified some miRNAs controlling important molecular muscle markers that could be relevant for the farming maximization.

BACKGROUND: Durum wheat (Triticum turgidum L. ssp. durum Desf. Husn) is the main staple crop used to make pasta products worldwide. Under the current climate change scenarios, genetic variability within a crop plays a crucial role in the successful release of new varieties with high yields and wide crop adaptation. In this study we evaluated a durum wheat collection consisting of 197 genotypes that mainly comprised a historical set of Argentinian germplasm but also included worldwide accessions. RESULTS: We assessed the genetic diversity, population structure and linkage disequilibrium (LD) patterns in this collection using a 35 K SNP array. The level of polymorphism was considered, taking account of the frequent and rare allelic variants. A total of 1547 polymorphic SNPs was located within annotated genes. Genetic diversity in the germplasm collection increased slightly from 1915 to 2010. However, a reduction in genetic diversity using SNPs with rare allelic variants was observed after 1979. However, larger numbers of rare private alleles were observed in the 2000-2009 period, indicating that a high reservoir of rare alleles is still present among the recent germplasm in a very low frequency. The percentage of pairwise loci in LD in the durum genome was low (13.4%) in our collection. Overall LD and the high (r2 > 0.7) or complete (r2 = 1) LD presented different patterns in the chromosomes. The LD increased over three main breeding periods (1915-1979, 1980-1999 and 2000-2020). CONCLUSIONS: Our results suggest that breeding and selection have impacted differently on the A and B genomes, particularly on chromosome 6A and 2A. The collection was structured in five sub-populations and modern Argentinian accessions (cluster Q4) which were clearly differentiated. Our study contributes to the understanding of the complexity of Argentinian durum wheat germplasm and to derive future breeding strategies enhancing the use of genetic diversity in a more efficient and targeted way.

BACKGROUND: The complex life cycle of the coconut crab, Birgus latro, begins when an obligate terrestrial adult female visits the intertidal to hatch zoea larvae into the surf. After drifting for several weeks in the ocean, the post-larval glaucothoes settle in the shallow subtidal zone, undergo metamorphosis, and the early juveniles then subsequently make their way to land where they undergo further physiological changes that prevent them from ever entering the sea again. Here, we sequenced, assembled and analyzed the coconut crab genome to shed light on its adaptation to terrestrial life. For comparison, we also assembled the genomes of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus, and the short-tailed marine-living red king crab, Paralithodes camtschaticus. Our selection of the latter two organisms furthermore allowed us to explore parallel evolution of the crab-like form in anomurans. RESULTS: All three assembled genomes are large, repeat-rich and AT-rich. Functional analysis reveals that the coconut crab has undergone proliferation of genes involved in the visual, respiratory, olfactory and cytoskeletal systems. Given that the coconut crab has atypical mitochondrial DNA compared to other anomurans, we argue that an abundance of kif22 and other significantly proliferated genes annotated with mitochondrial and microtubule functions, point to unique mechanisms involved in providing cellular energy via nuclear protein-coding genes supplementing mitochondrial and microtubule function. We furthermore detected in the coconut crab a significantly proliferated HOX gene, caudal, that has been associated with posterior development in Drosophila, but we could not definitively associate this gene with carcinization in the Anomura since it is also significantly proliferated in the ornate spiny lobster. However, a cuticle-associated coatomer gene, gammacop, that is significantly proliferated in the coconut crab, may play a role in hardening of the adult coconut crab abdomen in order to mitigate desiccation in terrestrial environments. CONCLUSION: The abundance of genomic features in the three assembled genomes serve as a source of hypotheses for future studies of anomuran environmental adaptations such as shell-utilization, perception of visual and olfactory cues in terrestrial environments, and cuticle sclerotization. We hypothesize that the coconut crab exhibits gene proliferation in lieu of alternative splicing as a terrestrial adaptation mechanism and propose life-stage transcriptomic assays to test this hypothesis.

BACKGROUND: Salvadora persica L. (Toothbrush tree - Miswak; family-Salvadoraceae) grows in the arid-land ecosystem and possesses economic and medicinal importance. The species, genus and the family have no genomic datasets available specifically on chloroplast (cp) genomics and taxonomic evolution. Herein, we have sequenced the complete chloroplast genome of S. persica for the first time and compared it with 11 related specie's cp genomes from the order Brassicales. RESULTS: The S. persica cp genome was 153,379 bp in length containing a sizeable single-copy region (LSC) of 83,818 bp which separated from the small single-copy region (SSC) of 17,683 bp by two inverted repeats (IRs) each 25,939 bp. Among these genomes, the largest cp genome size (160,600 bp) was found in M. oleifera, while in S. persica it was the smallest (153,379 bp). The cp genome of S. persica encoded 131 genes, including 37 tRNA genes, eight rRNA genes and 86 protein-coding genes. Besides, S. persica contains 27 forward, 36 tandem and 19 palindromic repeats. The S. persica cp genome had 154 SSRs with the highest number in the LSC region. Complete cp genome comparisons showed an overall high degree of sequence resemblance between S. persica and related cp genomes. Some divergence was observed in the intergenic spaces of other species. Phylogenomic analyses of 60 shared genes indicated that S. persica formed a single clade with A. tetracantha with high bootstrap values. The family Salvadoraceae is closely related to Capparaceae and Petadiplandraceae rather than to Bataceae and Koberliniacaea. CONCLUSION: The current genomic datasets provide pivotal genetic resources to determine the phylogenetic relationships, genome evolution and future genetic diversity-related studies of S. persica in complex angiosperm families.

BACKGROUND: A body distribution with high intramuscular fat and low abdominal fat is the ideal goal for broiler breeding. Preadipocytes with different origins have differences in terms of metabolism and gene expression. The transcriptome analysis performed in this study of intramuscular preadipocytes (DIMFPs) and adipose tissue-derived preadipocytes (DAFPs) aimed to explore the characteristics of lipid deposition in different chicken preadipocytes by dedifferentiation in vitro. RESULTS: Compared with DAFPs, the total lipid content in DIMFPs was reduced (P < 0.05). Moreover, 72 DEGs related to lipid metabolism were screened, which were involved in adipocyte differentiation, fatty acid transport and fatty acid synthesis, lipid stabilization, and lipolysis. Among the 72 DEGs, 19 DEGs were enriched in the PPAR signaling pathway, indicating its main contribution to the regulation of the difference in lipid deposition between DAFPs and DIMFPs. Among these 19 genes, the representative APOA1, ADIPOQ, FABP3, FABP4, FABP7, HMGCS2, LPL and RXRG genes were downregulated, but the ACSL1, FABP5, PCK2, PDPK1, PPARG, SCD, SCD5, and SLC27A6 genes were upregulated (P < 0.05 or P < 0.01) in the DIMFPs. In addition, the well-known pathways affecting lipid metabolism (MAPK, TGF-beta and calcium) and the pathways related to cell communication were enriched, which may also contribute to the regulation of lipid deposition. Finally, the regulatory network for the difference in lipid deposition between chicken DAFPs and DIMFPs was proposed based on the above information. CONCLUSIONS: Our data suggested a difference in lipid deposition between DIMFPs and DAFPs of chickens in vitro and proposed a molecular regulatory network for the difference in lipid deposition between chicken DAFPs and DIMFPs. The lipid content was significantly increased in DAFPs by the direct mediation of PPAR signaling pathways. These findings provide new insights into the regulation of tissue-specific fat deposition and the optimization of body fat distribution in broilers.

BACKGROUND: The relict gull (Larus relictus), was classified as vulnerable on the IUCN Red List and is a first-class national protected bird in China. Genomic resources for L. relictus are lacking, which limits the study of its evolution and its conservation. RESULTS: In this study, based on the Illumina and PacBio sequencing platforms, we successfully assembled the genome of L. relictus, one of the few known reference genomes in genus Larus. The size of the final assembled genome was 1.21 Gb, with a contig N50 of 8.11 Mb. A total of 18,454 genes were predicted from the assembly results, with 16,967 (91.94%) of these genes annotated. The genome contained 92.52 Mb of repeat sequence, accounting for 7.63% of the assembly. A phylogenetic tree was constructed using 4902 single-copy orthologous genes, which showed L. relictus had closest relative of L. smithsonianus, with divergence time of 14.7 Mya estimated between of them. PSMC analyses indicated that L. relictus had been undergoing a long-term population decline during 0.01-0.1 Mya with a small effective population size fom 8800 to 2200 individuals. CONCLUSIONS: This genome will be a valuable genomic resource for a range of genomic and conservation studies of L. relictus and will help to establish a foundation for further studies investigating whether the breeding population is a complex population. As the species is threatened by habitat loss and fragmentation, actions to protect L. relictus are suggested to alleviate the fragmentation of breeding populations.

BACKGROUND: DNA methyltransferase (DMT) genes contribute to plant stress responses and development by de novo establishment and subsequent maintenance of DNA methylation during replication. The photoperiod and/or temperature-sensitive genic male sterile (P/TGMS) lines play an important role in hybrid seed production of wheat. However, only a few studies have reported on the effect of DMT genes on temperature-sensitive male sterility of wheat. Although DMT genes have been investigated in some plant species, the identification and analysis of DMT genes in wheat (Triticum aestivum L.) based on genome-wide levels have not been reported. RESULTS: In this study, a detailed overview of phylogeny of 52 wheat DMT (TaDMT) genes was presented. Homoeolog retention for TaDMT genes was significantly above the average retention rate for whole-wheat genes, indicating the functional importance of many DMT homoeologs. We found that the strikingly high number of TaDMT genes resulted mainly from the significant expansion of the TaDRM subfamily. Intriguingly, all 5 paralogs belonged to the wheat DRM subfamily, and we speculated that tandem duplications might play a crucial role in the TaDRM subfamily expansion. Through the transcriptional analysis of TaDMT genes in a TGMS line BS366 and its hybrids with the other six fertile lines under sterile and fertile conditions, we concluded that TaCMT-D2, TaMET1-B1, and TaDRM-U6 might be involved in male sterility in BS366. Furthermore, a correlation analysis showed that TaMET1-B1 might negatively regulate the expression of TaRAFTIN1A, an important gene for pollen development, so we speculated regarding an epigenetic regulatory mechanism underlying the male sterility of BS366 via the interaction between TaMET1-B1 and TaRAFTIN1A. CONCLUSIONS: Our findings presented a detailed phylogenic overview of the DMT genes and could provide novel insights into the effects of DMT genes on TGMS wheat.

BACKGROUND: Cotton fibers provide a powerful model for studying cell differentiation and elongation. Each cotton fiber is a singular and elongated cell derived from epidermal-layer cells of a cotton seed. Efforts to understand this dramatic developmental shift have been impeded by the difficulty of separation between fiber and epidermal cells. RESULTS: Here we employed laser-capture microdissection (LCM) to separate these cell types. RNA-seq analysis revealed transitional differences between fiber and epidermal-layer cells at 0 or 2 days post anthesis. Specifically, down-regulation of putative cell cycle genes was coupled with upregulation of ribosome biosynthesis and translation-related genes, which may suggest their respective roles in fiber cell initiation. Indeed, the amount of fibers in cultured ovules was increased by cell cycle progression inhibitor, Roscovitine, and decreased by ribosome biosynthesis inhibitor, Rbin-1. Moreover, subfunctionalization of homoeologs was pervasive in fiber and epidermal cells, with expression bias towards 10% more D than A homoeologs of cell cycle related genes and 40-50% more D than A homoeologs of ribosomal protein subunit genes. Key cell cycle regulators were predicted to be epialleles in allotetraploid cotton. MYB-transcription factor genes displayed expression divergence between fibers and ovules. Notably, many phytohormone-related genes were upregulated in ovules and down-regulated in fibers, suggesting spatial-temporal effects on fiber cell development. CONCLUSIONS: Fiber cell initiation is accompanied by cell cycle arrest coupled with active ribosome biosynthesis, spatial-temporal regulation of phytohormones and MYB transcription factors, and homoeolog expression bias of cell cycle and ribosome biosynthesis genes. These valuable genomic resources and molecular insights will help develop breeding and biotechnological tools to improve cotton fiber production.

BACKGROUND: We previously reported the results of CRISPR/Cas9 knock-out (KO) of type-I and type-III vitellogenins (Vtgs) in zebrafish, which provided the first experimental evidence on essentiality and disparate functioning of Vtgs at different stages during early development. However, the specific contributions of different types of Vtg to major cellular processes remained to be investigated. The present study employed liquid chromatography and tandem mass spectrometry (LC-MS/MS) to meet this deficit. Proteomic profiles of zebrafish eggs lacking three type-I Vtgs simultaneously (vtg1-KO), or lacking only type III Vtg (vtg3-KO) were compared to those of wild type (Wt) eggs. Obtained spectra were searched against a zebrafish proteome database and identified proteins were quantified based on normalized spectral counts. RESULTS: The vtg-KO caused severe changes in the proteome of 1-cell stage zebrafish eggs. These changes were disclosed by molecular signatures that highly resembled the proteomic phenotype of poor quality zebrafish eggs reported in our prior studies. Proteomic profiles of vtg-KO eggs and perturbations in abundances of hundreds of proteins revealed unique, noncompensable contributions of multiple Vtgs to protein and in energy homeostasis. The lack of this contribution appears to have a significant impact on endoplasmic reticulum and mitochondrial functions, and thus embryonic development, even after zygotic genome activation. Increased endoplasmic reticulum stress, Redox/Detox activities, glycolysis/gluconeogenesis, enrichment in cellular proliferation and in human neurodegenerative disease related activities in both vtg1- and vtg3-KO eggs were found to be indicators of the aforementioned conditions. Distinctive increase in apoptosis and Parkinson disease pathways, as well as the decrease in lipid metabolism related activities in vtg3-KO eggs implies compelling roles of Vtg3, the least abundant form of Vtgs in vertebrate eggs, in mitochondrial activities. Several differentially abundant proteins representing the altered molecular mechanisms have been identified as strong candidate markers for studying the details of these mechanisms during early embryonic development in zebrafish and possibly other vertebrates. CONCLUSIONS: These findings indicate that the global egg proteome is subject to extensive modification depending on the presence or absence of specific Vtgs and that these modifications can have a major impact on developmental competence.

Erratum in BMC Genomics. 2021 May 18;22(1):358.

BACKGROUND: Rodent malaria parasites are important models for studying host-malaria parasite interactions such as host immune response, mechanisms of parasite evasion of host killing, and vaccine development. One of the rodent malaria parasites is Plasmodium yoelii, and multiple P. yoelii strains or subspecies that cause different disease phenotypes have been widely employed in various studies. The genomes and transcriptomes of several P. yoelii strains have been analyzed and annotated, including the lethal strains of P. y. yoelii YM (or 17XL) and non-lethal strains of P. y. yoelii 17XNL/17X. Genomic DNA sequences and cDNA reads from another subspecies P. y. nigeriensis N67 have been reported for studies of genetic polymorphisms and parasite response to drugs, but its genome has not been assembled and annotated. RESULTS: We performed genome sequencing of the N67 parasite using the PacBio long-read sequencing technology, de novo assembled its genome and transcriptome, and predicted 5383 genes with high overall annotation quality. Comparison of the annotated genome of the N67 parasite with those of YM and 17X parasites revealed a set of genes with N67-specific orthology, expansion of gene families, particularly the homologs of the Plasmodium chabaudi erythrocyte membrane antigen, large numbers of SNPs and indels, and proteins predicted to interact with host immune responses based on their functional domains. CONCLUSIONS: The genomes of N67 and 17X parasites are highly diverse, having approximately one polymorphic site per 50 base pairs of DNA. The annotated N67 genome and transcriptome provide searchable databases for fast retrieval of genes and proteins, which will greatly facilitate our efforts in studying the parasite biology and gene function and in developing effective control measures against malaria.

BACKGROUND: Phytophthora cinnamomi is an oomycete pathogen of global relevance. It is considered as one of the most invasive species, which has caused irreversible damage to natural ecosystems and horticultural crops. There is currently a lack of a high-quality reference genome for this species despite several attempts that have been made towards sequencing its genome. The lack of a good quality genome sequence has been a setback for various genetic and genomic research to be done on this species. As a consequence, little is known regarding its genome characteristics and how these contribute to its pathogenicity and invasiveness. RESULTS: In this work we generated a high-quality genome sequence and annotation for P. cinnamomi using a combination of Oxford Nanopore and Illumina sequencing technologies. The annotation was done using RNA-Seq data as supporting gene evidence. The final assembly consisted of 133 scaffolds, with an estimated genome size of 109.7 Mb, N50 of 1.18 Mb, and BUSCO completeness score of 97.5%. Genome partitioning analysis revealed that P. cinnamomi has a two-speed genome characteristic, similar to that of other oomycetes and fungal plant pathogens. In planta gene expression analysis revealed up-regulation of pathogenicity-related genes, suggesting their important roles during infection and host degradation. CONCLUSION: This study has provided a high-quality reference genome and annotation for P. cinnamomi. This is among the best assembled genomes for any Phytophthora species assembled to date and thus resulted in improved identification and characterization of pathogenicity-related genes, some of which were undetected in previous versions of genome assemblies. Phytophthora cinnamomi harbours a large number of effector genes which are located in the gene-poor regions of the genome. This unique genomic partitioning provides P. cinnamomi with a high level of adaptability and could contribute to its success as a highly invasive species. Finally, the genome sequence, its annotation and the pathogenicity effectors identified in this study will serve as an important resource that will enable future studies to better understand and mitigate the impact of this important pathogen.

BACKGROUND: Sucrose nonfermenting-1 (SNF1)-related protein kinases (SnRKs) play important roles in regulating metabolism and stress responses in plants, providing a conduit for crosstalk between metabolic and stress signalling, in some cases involving the stress hormone, abscisic acid (ABA). The burgeoning and divergence of the plant gene family has led to the evolution of three subfamilies, SnRK1, SnRK2 and SnRK3, of which SnRK2 and SnRK3 are unique to plants. Therefore, the study of SnRKs in crops may lead to the development of strategies for breeding crop varieties that are more resilient under stress conditions. In the present study, we describe the SnRK gene family of barley (Hordeum vulgare), the widespread cultivation of which can be attributed to its good adaptation to different environments. RESULTS: The barley HvSnRK gene family was elucidated in its entirety from publicly-available genome data and found to comprise 50 genes. Phylogenetic analyses assigned six of the genes to the HvSnRK1 subfamily, 10 to HvSnRK2 and 34 to HvSnRK3. The search was validated by applying it to Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) genome data, identifying 50 SnRK genes in rice (four OsSnRK1, 11 OsSnRK2 and 35 OsSnRK3) and 39 in Arabidopsis (three AtSnRK1, 10 AtSnRK2 and 26 AtSnRK3). Specific motifs were identified in the encoded barley proteins, and multiple putative regulatory elements were found in the gene promoters, with light-regulated elements (LRE), ABA response elements (ABRE) and methyl jasmonate response elements (MeJa) the most common. RNA-seq analysis showed that many of the HvSnRK genes responded to ABA, some positively, some negatively and some with complex time-dependent responses. CONCLUSIONS: The barley HvSnRK gene family is large, comprising 50 members, subdivided into HvSnRK1 (6 members), HvSnRK2 (10 members) and HvSnRK3 (34 members), showing differential positive and negative responses to ABA.

BACKGROUND: Immune responses need to be initiated rapidly, and maintained as needed, to prevent establishment and growth of infections. At the same time, resources need to be balanced with other physiological processes. On the level of transcription, studies have shown that this balancing act is reflected in tight control of the initiation kinetics and shutdown dynamics of specific immune genes. RESULTS: To investigate genome-wide expression dynamics and trade-offs after infection at a high temporal resolution, we performed an RNA-seq time course on D. melanogaster with 20 time points post Imd stimulation. A combination of methods, including spline fitting, cluster analysis, and Granger causality inference, allowed detailed dissection of expression profiles, lead-lag interactions, and functional annotation of genes through guilt-by-association. We identified Imd-responsive genes and co-expressed, less well characterized genes, with an immediate-early response and sustained up-regulation up to 5 days after stimulation. In contrast, stress response and Toll-responsive genes, among which were Bomanins, demonstrated early and transient responses. We further observed a strong trade-off with metabolic genes, which strikingly recovered to pre-infection levels before the immune response was fully resolved. CONCLUSIONS: This high-dimensional dataset enabled the comprehensive study of immune response dynamics through the parallel application of multiple temporal data analysis methods. The well annotated data set should also serve as a useful resource for further investigation of the D. melanogaster innate immune response, and for the development of methods for analysis of a post-stress transcriptional response time-series at whole-genome scale.

BACKGROUND: In Iran, river buffalo is of great importance. It plays an important role in the economy of the Country, because its adaptation to harsh climate conditions and long productive lifespan permitting its farming across the Country and to convert low-quality feed into valuable milk. The genetic variability in Iranian buffalo breeds have been recently studied using SNPs genotyping data, but a whole genome Copy Number Variants (CNVs) mapping was not available. The aim of this study was to perform a genome wide CNV scan in 361 buffaloes of the three Iranian river breeds (Azeri, Khuzestani and Mazandarani) through the analysis of data obtained using the Axiom® Buffalo Genotyping Array 90 K. RESULTS: CNVs detection resulted in a total of 9550 CNVs and 302 CNVRs identified in at least 5% of samples within breed, covering around 1.97% of the buffalo genome. and A total of 22 CNVRs were identified in all breeds and a different proportion of regions were in common among the three populations. Within the more represented CNVRs (n = 302) mapped a total of 409 buffalo genes, some of which resulted associated with morphological, healthy, milk, meat and reproductive traits, according to Animal Genome Cattle database. CONCLUSIONS: This work provides a step forward in the interpretation of genomic variation within and among the buffalo populations, releasing a first map of CNVs and providing insights about their recent selection and adaptation to environment. The presence of the set of genes and QTL traits harbored in the CNVRs could be possibly linked with the buffalo's natural adaptive history together to a recent selection for milk used as primary food source from this species.

BACKGROUND: The lack of an understanding about the genomic architecture underpinning parental behaviour in subsocial insects displaying simple parental behaviours prevents the development of a full understanding about the evolutionary origin of sociality. Lethrus apterus is one of the few insect species that has biparental care. Division of labour can be observed between parents during the reproductive period in order to provide food and protection for their offspring. RESULTS: Here, we report the draft genome of L. apterus, the first genome in the family Geotrupidae. The final assembly consisted of 286.93 Mbp in 66,933 scaffolds. Completeness analysis found the assembly contained 93.5% of the Endopterygota core BUSCO gene set. Ab initio gene prediction resulted in 25,385 coding genes, whereas homology-based analyses predicted 22,551 protein coding genes. After merging, 20,734 were found during functional annotation. Compared to other publicly available beetle genomes, 23,528 genes among the predicted genes were assigned to orthogroups of which 1664 were in species-specific groups. Additionally, reproduction related genes were found among the predicted genes based on which a reduction in the number of odorant- and pheromone-binding proteins was detected. CONCLUSIONS: These genes can be used in further comparative and functional genomic researches which can advance our understanding of the genetic basis and hence the evolution of parental behaviour.

BACKGROUND: Chloroplasts are important semi-autonomous organelles in plants and algae. Unlike higher plants, the chloroplast genomes of green algal linage have distinct features both in organization and expression. Despite the architecture of chloroplast genome having been extensively studied in higher plants and several model species of algae, little is known about the transcriptional features of green algal chloroplast-encoded genes. RESULTS: Based on full-length cDNA (Iso-Seq) sequencing, we identified widely co-transcribed polycistronic transcriptional units (PTUs) in the green alga Caulerpa lentillifera. In addition to clusters of genes from the same pathway, we identified a series of PTUs of up to nine genes whose function in the plastid is not understood. The RNA data further allowed us to confirm widespread expression of fragmented genes and conserved open reading frames, which are both important features in green algal chloroplast genomes. In addition, a newly fragmented gene specific to C. lentillifera was discovered, which may represent a recent gene fragmentation event in the chloroplast genome. With the newly annotated exon-intron boundary information, gene structural annotation was greatly improved across the siphonous green algae lineages. Our data also revealed a type of non-canonical Group II introns, with a deviant secondary structure and intronic ORFs lacking known splicing or mobility domains. These widespread introns have conserved positions in their genes and are excised precisely despite lacking clear consensus intron boundaries. CONCLUSION: Our study fills important knowledge gaps in chloroplast genome organization and transcription in green algae, and provides new insights into expression of polycistronic transcripts, freestanding ORFs and fragmented genes in algal chloroplast genomes. Moreover, we revealed an unusual type of Group II intron with distinct features and conserved positions in Bryopsidales. Our data represents interesting additions to knowledge of chloroplast intron structure and highlights clusters of uncharacterized genes that probably play important roles in plastids.

BACKGROUND: Tribenuron methyl (TBM) is an herbicide that inhibits sulfonylurea acetolactate synthase (ALS) and is one of the most widely used broad-leaved herbicides for crop production. However, soil residues or drifting of the herbicide spray might affect the germination and growth of rapeseed, Brassica napus, so it is imperative to understand the response mechanism of rape to TBM during germination. The aim of this study was to use transcriptome analysis to reveal the gene responses in herbicide-tolerant rapeseed to TBM stress during seed germination. RESULTS: 2414, 2286, and 1068 differentially expressed genes (DEGs) were identified in TBM-treated resistant vs sensitive lines, treated vs. control sensitive lines, treated vs. control resistant lines, respectively. GO analysis showed that most DEGs were annotated to the oxidation-reduction pathways and catalytic activity. KEGG enrichment was mainly involved in plant-pathogen interactions, α-linolenic acid metabolism, glucosinolate biosynthesis, and phenylpropanoid biosynthesis. Based on GO and KEGG enrichment, a total of 137 target genes were identified, including genes involved in biotransferase activity, response to antioxidant stress and lipid metabolism. Biotransferase genes, CYP450, ABC and GST, detoxify herbicide molecules through physical or biochemical processes. Antioxidant genes, RBOH, WRKY, CDPK, MAPK, CAT, and POD regulate plant tolerance by transmitting ROS signals and triggering antioxidant enzyme expression. Lipid-related genes and hormone-related genes were also found, such as LOX3, ADH1, JAZ6, BIN2 and ERF, and they also played an important role in herbicide resistance. CONCLUSIONS: This study provides insights for selecting TBM-tolerant rapeseed germplasm and exploring the molecular mechanism of TBM tolerance during germination.

BACKGROUND: Paeonia ostii is a potentially important oilseed crop because its seed yield is high, and the seeds are rich in α-linolenic acid (ALA). However, the molecular mechanisms underlying ALA biosynthesis during seed kernel, seed testa, and fruit pericarp development in this plant are unclear. We used transcriptome data to address this knowledge gap. RESULTS: Gas chromatograph-mass spectrometry indicated that ALA content was highest in the kernel, moderate in the testa, and lowest in the pericarp. Therefore, we used RNA-sequencing to compare ALA synthesis among these three tissues. We identified 227,837 unigenes, with an average length of 755 bp. Of these, 1371 unigenes were associated with lipid metabolism. The fatty acid (FA) biosynthesis and metabolism pathways were significantly enriched during the early stages of oil accumulation in the kernel. ALA biosynthesis was significantly enriched in parallel with increasing ALA content in the testa, but these metabolic pathways were not significantly enriched during pericarp development. By comparing unigene transcription profiles with patterns of ALA accumulation, specific unigenes encoding crucial enzymes and transcription factors (TFs) involved in de novo FA biosynthesis and oil accumulation were identified. Specifically, the bell-shaped expression patterns of genes encoding SAD, FAD2, FAD3, PDCT, PDAT, OLE, CLE, and SLE in the kernel were similar to the patterns of ALA accumulation in this tissue. Genes encoding BCCP, BC, KAS I- III, and FATA were also upregulated during the early stages of oil accumulation in the kernel. In the testa, the upregulation of the genes encoding SAD, FAD2, and FAD3 was followed by a sharp increase in the concentrations of ALA. In contrast, these genes were minimally expressed (and ALA content was low) throughout pericarp development. CONCLUSIONS: We used three tissues with high, moderate, and low ALA concentrations as an exemplar system in which to investigate tissue-specific ALA accumulation mechanisms in P. ostii. The genes and TFs identified herein might be useful targets for future studies of ALA accumulation in the tree peony. This study also provides a framework for future studies of FA biosynthesis in other oilseed plants.

BACKGROUND: The newly discovered reversible N6-methyladenosine (m6A) modification plays an important regulatory role in gene expression. Long non-coding RNAs (lncRNAs) participate in Marek's disease virus (MDV) replication but how m6A modifications in lncRNAs are affected during MDV infection is currently unknown. Herein, we profiled the transcriptome-wide m6A modification in lncRNAs in MDV-infected chicken embryo fibroblast (CEF) cells. RESULTS: Methylated RNA immunoprecipitation sequencing results revealed that the lncRNA m6A modification is highly conserved with MDV infection increasing the expression of lncRNA m6A modified sites compared to uninfected cell controls. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that lncRNA m6A modifications were highly associated with signaling pathways associated with MDV infection. CONCLUSIONS: In this study, the alterations seen in transcriptome-wide m6A occurring in lncRNAs following MDV-infection suggest this process plays important regulatory roles during MDV replication. We report for the first time profiling of the alterations in transcriptome-wide m6A modification in lncRNAs of MDV-infected CEF cells.

BACKGROUND: Most Distylium species are endangered. Distylium species mostly display homoplasy in their flowers and fruits, and are classified primarily based on leaf morphology. However, leaf size, shape, and serration vary tremendously making it difficult to use those characters to identify most species and a significant challenge to address the taxonomy of Distylium. To infer robust relationships and develop variable markers to identify Distylium species, we sequenced most of the Distylium species chloroplast genomes. RESULTS: The Distylium chloroplast genome size was 159,041-159,127 bp and encoded 80 protein-coding, 30 transfer RNAs, and 4 ribosomal RNA genes. There was a conserved gene order and a typical quadripartite structure. Phylogenomic analysis based on whole chloroplast genome sequences yielded a highly resolved phylogenetic tree and formed a monophyletic group containing four Distylium clades. A dating analysis suggested that Distylium originated in the Oligocene (34.39 Ma) and diversified within approximately 1 Ma. The evidence shows that Distylium is a rapidly radiating group. Four highly variable markers, matK-trnK, ndhC-trnV, ycf1, and trnT-trnL, and 74 polymorphic simple sequence repeats were discovered in the Distylium plastomes. CONCLUSIONS: The plastome sequences had sufficient polymorphic information to resolve phylogenetic relationships and identify Distylium species accurately.

BACKGROUND: Yersinia ruckeri is a pathogen that can cause enteric redmouth disease in salmonid species, damaging global production of economically important fish including rainbow trout (Oncorhynchus mykiss). Herein, we conducted the transcriptomic profiling of spleen samples from rainbow trout at 24 h post-Y. ruckeri infection via RNA-seq in an effort to more fully understand their immunological responses. RESULTS: We identified 2498 differentially expressed genes (DEGs), of which 2083 and 415 were up- and down-regulated, respectively. We then conducted a more in-depth assessment of 78 DEGs associated with the immune system including CCR9, CXCL11, IL-1β, CARD9, IFN, TNF, CASP8, NF-κB, NOD1, TLR8α2, HSP90, and MAPK11, revealing these genes to be associated with 20 different immunological KEGG pathways including the Cytokine-cytokine receptor interaction, Toll-like receptor signaling, RIG-I-like receptor signaling, NOD-like receptor signaling, and MAPK signaling pathways. Additionally, the differential expression of 8 of these DEGs was validated by a qRT-PCR approach and their immunological importance was then discussed. CONCLUSIONS: Our findings provide preliminary insight on molecular mechanism underlying the immune responses of rainbow trout following Y. ruckeri infection and the base for future studies of host-pathogen interactions in rainbow trout.

BACKGROUND: Sorghum grain mold is the most important disease of the crop. The disease results from simultaneous infection of the grain by multiple fungal species. Host responses to these fungi and the underlying molecular and cellular processes are poorly understood. To understand the genetic, molecular and biochemical components of grain mold resistance, transcriptome profiles of the developing grain of resistant and susceptible sorghum genotypes were studied. RESULTS: The developing kernels of grain mold resistant RTx2911 and susceptible RTx430 sorghum genotypes were inoculated with a mixture of fungal pathogens mimicking the species complexity of the disease under natural infestation. Global transcriptome changes corresponding to multiple molecular and cellular processes, and biological functions including defense, secondary metabolism, and flavonoid biosynthesis were observed with differential regulation in the two genotypes. Genes encoding pattern recognition receptors (PRRs), regulators of growth and defense homeostasis, antimicrobial peptides, pathogenesis-related proteins, zein seed storage proteins, and phytoalexins showed increased expression correlating with resistance. Notably, SbLYK5 gene encoding an orthologue of chitin PRR, defensin genes SbDFN7.1 and SbDFN7.2 exhibited higher expression in the resistant genotype. The SbDFN7.1 and SbDFN7.2 genes are tightly linked and transcribed in opposite orientation with a likely common bidirectional promoter. Interestingly, increased expression of JAZ and other transcriptional repressors were observed that suggested the tight regulation of plant defense and growth. The data suggest a pathogen inducible defense system in the developing grain of sorghum that involves the chitin PRR, MAPKs, key transcription factors, downstream components regulating immune gene expression and accumulation of defense molecules. We propose a model through which the biosynthesis of 3-deoxyanthocynidin phytoalexins, defensins, PR proteins, other antimicrobial peptides, and defense suppressing proteins are regulated by a pathogen inducible defense system in the developing grain. CONCLUSIONS: The transcriptome data from a rarely studied tissue shed light into genetic, molecular, and biochemical components of disease resistance and suggested that the developing grain shares conserved immune response mechanisms but also components uniquely enriched in the grain. Resistance was associated with increased expression of genes encoding regulatory factors, novel grain specific antimicrobial peptides including defensins and storage proteins that are potential targets for crop improvement.

BACKGROUND: The feed conversion ratio (FCR) is an important productive trait that greatly affects profits in the pig industry. Elucidating the genetic mechanisms underpinning FCR may promote more efficient improvement of FCR through artificial selection. In this study, we integrated a genome-wide association study (GWAS) with transcriptome analyses of different tissues in Yorkshire pigs (YY) with the aim of identifying key genes and signalling pathways associated with FCR. RESULTS: A total of 61 significant single nucleotide polymorphisms (SNPs) were detected by GWAS in YY. All of these SNPs were located on porcine chromosome (SSC) 5, and the covered region was considered a quantitative trait locus (QTL) region for FCR. Some genes distributed around these significant SNPs were considered as candidates for regulating FCR, including TPH2, FAR2, IRAK3, YARS2, GRIP1, FRS2, CNOT2 and TRHDE. According to transcriptome analyses in the hypothalamus, TPH2 exhibits the potential to regulate intestinal motility through serotonergic synapse and oxytocin signalling pathways. In addition, GRIP1 may be involved in glutamatergic and GABAergic signalling pathways, which regulate FCR by affecting appetite in pigs. Moreover, GRIP1, FRS2, CNOT2, and TRHDE may regulate metabolism in various tissues through a thyroid hormone signalling pathway. CONCLUSIONS: Based on the results from GWAS and transcriptome analyses, the TPH2, GRIP1, FRS2, TRHDE, and CNOT2 genes were considered candidate genes for regulating FCR in Yorkshire pigs. These findings improve the understanding of the genetic mechanisms of FCR and may help optimize the design of breeding schemes.

BACKGROUND: Genes that belong to the same network are frequently co-expressed, but collectively, how the coordination of the whole transcriptome is perturbed during aging remains unclear. To explore this, we calculated the correlation of each gene in the transcriptome with every other, in the brain of young and older outbred deer mice (P. leucopus and P. maniculatus). RESULTS: In about 25 % of the genes, coordination was inversed during aging. Gene Ontology analysis in both species, for the genes that exhibited inverse transcriptomic coordination during aging pointed to alterations in the perception of smell, a known impairment occurring during aging. In P. leucopus, alterations in genes related to cholesterol metabolism were also identified. Among the genes that exhibited the most pronounced inversion in their coordination profiles during aging was THBS4, that encodes for thrombospondin-4, a protein that was recently identified as rejuvenation factor in mice. Relatively to its breadth, abolishment of coordination was more prominent in the long-living P. leucopus than in P. maniculatus but in the latter, the intensity of de-coordination was higher. CONCLUSIONS: There sults suggest that aging is associated with more stringent retention of expression profiles for some genes and more abrupt changes in others, while more subtle but widespread changes in gene expression appear protective. Our findings shed light in the mode of the transcriptional changes occurring in the brain during aging and suggest that strategies aiming to broader but more modest changes in gene expression may be preferrable to correct aging-associated deregulation in gene expression.

BACKGROUND: Preweaned rumen development is vital for animal health and efficient fermentation. In this study, we integrated ruminal transcriptomic and metagenomic data to explore the dynamics of rumen functions, microbial colonization, and their functional interactions during the first 8 weeks of life in goats. RESULTS: The dynamic rumen transcriptomic and microbial profiles both exhibited two distinct phases during early rumen development. The differentially expressed genes of the rumen transcriptome between the two phases showed that the immune-related response was enriched in the first phase and nutrient-related metabolism was enriched in the second phase, whereas the differentially expressed genes of the rumen microbiome were enriched in bacteriocin biosynthesis and glycolysis/gluconeogenesis activities. The developmental shift in the rumen transcriptome (at d 21) was earlier than the feed stimulus (at d 25) and the shift in the rumen microbiome (at d 42). Additionally, 15 temporal dynamic rumen gene modules and 20 microbial modules were revealed by coexpression network analysis. Functional correlations between the rumen and its microbiome were primarily involved in rumen pH homeostasis, nitrogen metabolism and the immune response. Rumen gene modules associated with the microbial alpha diversity index were also enriched in the immune response process. CONCLUSIONS: The present study touched the critical developmental process of rumen functions, microbial colonization and their functional interactions during preweaned development. Taken together, these results demonstrated that rumen development at the first phase is more likely a programmed process rather than stimulation from feed and the microbiome, while the shift of rumen metagenomes was likely regulated by both the diet and host. The intensive functional correlations between rumen genes and the microbiome demonstrated that synergistic processes occurred between them during early rumen development.

BACKGROUND: The key-ancestor approach has been frequently applied to prioritize individuals for whole-genome sequencing based on their marginal genetic contribution to current populations. Using this approach, we selected 70 key ancestors from two lines of the Swiss Large White breed that have been selected divergently for fertility and fattening traits and sequenced their genomes with short paired-end reads. RESULTS: Using pedigree records, we estimated the effective population size of the dam and sire line to 72 and 44, respectively. In order to assess sequence variation in both lines, we sequenced the genomes of 70 boars at an average coverage of 16.69-fold. The boars explained 87.95 and 95.35% of the genetic diversity of the breeding populations of the dam and sire line, respectively. Reference-guided variant discovery using the GATK revealed 26,862,369 polymorphic sites. Principal component, admixture and fixation index (FST) analyses indicated considerable genetic differentiation between the lines. Genomic inbreeding quantified using runs of homozygosity was higher in the sire than dam line (0.28 vs 0.26). Using two complementary approaches, we detected 51 signatures of selection. However, only six signatures of selection overlapped between both lines. We used the sequenced haplotypes of the 70 key ancestors as a reference panel to call 22,618,811 genotypes in 175 pigs that had been sequenced at very low coverage (1.11-fold) using the GLIMPSE software. The genotype concordance, non-reference sensitivity and non-reference discrepancy between thus inferred and Illumina PorcineSNP60 BeadChip-called genotypes was 97.60, 98.73 and 3.24%, respectively. The low-pass sequencing-derived genomic relationship coefficients were highly correlated (r > 0.99) with those obtained from microarray genotyping. CONCLUSIONS: We assessed genetic diversity within and between two lines of the Swiss Large White pig breed. Our analyses revealed considerable differentiation, even though the split into two populations occurred only few generations ago. The sequenced haplotypes of the key ancestor animals enabled us to implement genotyping by low-pass sequencing which offers an intriguing cost-effective approach to increase the variant density over current array-based genotyping by more than 350-fold.

BACKGROUND: Chickpea, pigeonpea, and groundnut are the primary legume crops of semi-arid tropics (SAT) and their global productivity is severely affected by drought stress. The plant-specific NAC (NAM - no apical meristem, ATAF - Arabidopsis transcription activation factor, and CUC - cup-shaped cotyledon) transcription factor family is known to be involved in majority of abiotic stresses, especially in the drought stress tolerance mechanism. Despite the knowledge available regarding NAC function, not much information is available on NAC genes in SAT legume crops. RESULTS: In this study, genome-wide NAC proteins - 72, 96, and 166 have been identified from the genomes of chickpea, pigeonpea, and groundnut, respectively, and later grouped into 10 clusters in chickpea and pigeonpea, while 12 clusters in groundnut. Phylogeny with well-known stress-responsive NACs in Arabidopsis thaliana, Oryza sativa (rice), Medicago truncatula, and Glycine max (soybean) enabled prediction of putative stress-responsive NACs in chickpea (22), pigeonpea (31), and groundnut (33). Transcriptome data revealed putative stress-responsive NACs at various developmental stages that showed differential expression patterns in the different tissues studied. Quantitative real-time PCR (qRT-PCR) was performed to validate the expression patterns of selected stress-responsive, Ca_NAC (Cicer arietinum - 14), Cc_NAC (Cajanus cajan - 15), and Ah_NAC (Arachis hypogaea - 14) genes using drought-stressed and well-watered root tissues from two contrasting drought-responsive genotypes of each of the three legumes. Based on expression analysis, Ca_06899, Ca_18090, Ca_22941, Ca_04337, Ca_04069, Ca_04233, Ca_12660, Ca_16379, Ca_16946, and Ca_21186; Cc_26125, Cc_43030, Cc_43785, Cc_43786, Cc_22429, and Cc_22430; Ah_ann1.G1V3KR.2, Ah_ann1.MI72XM.2, Ah_ann1.V0X4SV.1, Ah_ann1.FU1JML.2, and Ah_ann1.8AKD3R.1 were identified as potential drought stress-responsive candidate genes. CONCLUSION: As NAC genes are known to play role in several physiological and biological activities, a more comprehensive study on genome-wide identification and expression analyses of the NAC proteins have been carried out in chickpea, pigeonpea and groundnut. We have identified a total of 21 potential drought-responsive NAC genes in these legumes. These genes displayed correlation between gene expression, transcriptional regulation, and better tolerance against drought. The identified candidate genes, after validation, may serve as a useful resource for molecular breeding for drought tolerance in the SAT legume crops.

BACKGROUND: Alanine:glyoxylate aminotransferase 2 (AGXT2; EC 2.6.1.44) is the only enzyme that degrades the R-form of 3-aminoisobutyrate, an intermediate metabolite of thymine. AGXT2, as well as diaminoarginine dimethylaminohydrolase 1 (DDAH1; EC 3.5.3.18), works as an enzyme that degrades asymmetric dimethylarginine (ADMA), which competitively inhibits the nitric oxide synthase family. Thus, these two enzyme activities may change vascular vulnerability for a lifetime via the nitric oxide (NO) system. We investigated the association between vascular conditions and diseases such as hypertension and diabetes mellitus and polymorphisms of these two genes in 750 older Japanese subjects (mean age ± standard deviation, 77.0 ± 7.6 years) recruited using the complete enumeration survey method in the Nakayama study. Demographic and biochemical data, such as blood pressure (BP) and casual blood sugar (CBS), were obtained. Four functional single nucleotide polymorphisms (SNPs; rs37370, rs37369, rs180749, and rs16899974) of AGXT2 and one functional insertion/deletion polymorphism in the promotor region with four SNPs (rs307894, rs669173, rs997251, and rs13373844) of DDAH1 were investigated. Plasma ADMA was also analyzed in 163 subjects. RESULTS: The results of multiple regression analysis showed that a loss of the functional haplotype of AGXT2, CAAA, was significantly positively correlated with BP (systolic BP, p = 0.034; diastolic BP, p = 0.025) and CBS (p = 0.021). No correlation was observed between DDAH1 and either BP or CBS. ADMA concentrations were significantly elevated in subjects with two CAAA haplotypes compared with subjects without the CAAA haplotype (p = 0.033). CONCLUSIONS: Missense variants of AGXT2, but not DDAH1, may be related to vulnerability to vascular diseases such as hypertension and DM via the NO system.

BACKGROUND: Analysis of lymphocyte cell lines revealed substantial differences in the expression of mRNA and microRNA (miRNA) among human populations. The extent of such population-associated differences in actual human tissues remains largely unexplored. The placenta is one of the few solid human tissues that can be collected in substantial numbers in a controlled manner, enabling quantitative analysis of transient biomolecules such as RNA transcripts. Here, we analyzed microRNA (miRNA) expression in human placental samples derived from 36 individuals representing four genetically distinct human populations: African Americans, European Americans, South Asians, and East Asians. All samples were collected at the same hospital following a unified protocol, thus minimizing potential biases that might influence the results. RESULTS: Sequence analysis of the miRNA fraction yielded 938 annotated and 70 novel miRNA transcripts expressed in the placenta. Of them, 82 (9%) of annotated and 11 (16%) of novel miRNAs displayed quantitative expression differences among populations, generally reflecting reported genetic and mRNA-expression-based distances. Several co-expressed miRNA clusters stood out from the rest of the population-associated differences in terms of miRNA evolutionary age, tissue-specificity, and disease-association characteristics. Among three non-environmental influenced demographic parameters, the second largest contributor to miRNA expression variation after population was the sex of the newborn, with 32 miRNAs (3% of detected) exhibiting significant expression differences depending on whether the newborn was male or female. Male-associated miRNAs were evolutionarily younger and correlated inversely with the expression of target mRNA involved in neuron-related functions. In contrast, both male and female-associated miRNAs appeared to mediate different types of hormonal responses. Demographic factors further affected reported imprinted expression of 66 placental miRNAs: the imprinting strength correlated with the mother's weight, but not height. CONCLUSIONS: Our results showed that among 12 assessed demographic variables, population affiliation and fetal sex had a substantial influence on miRNA expression variation among human placental samples. The effect of newborn-sex-associated miRNA differences further led to expression inhibition of the target genes clustering in specific functional pathways. By contrast, population-driven miRNA differences might mainly represent neutral changes with minimal functional impacts.

BACKGROUND: Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines. RESULTS: In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2-4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters. CONCLUSIONS: Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

BACKGROUND: Commelinaceae (Commelinales) comprise 41 genera and are widely distributed in both the Old and New Worlds, except in Europe. The relationships among genera in this family have been suggested in several morphological and molecular studies. However, it is difficult to explain their relationships due to high morphological variations and low support values. Currently, many researchers have been using complete chloroplast genome data for inferring the evolution of land plants. In this study, we completed 15 new plastid genome sequences of subfamily Commelinoideae using the Mi-seq platform. We utilized genome data to reveal the structural variations and reconstruct the problematic positions of genera for the first time. RESULTS: All examined species of Commelinoideae have three pseudogenes (accD, rpoA, and ycf15), and the former two might be a synapomorphy within Commelinales. Only four species in tribe Commelineae presented IR expansion, which affected duplication of the rpl22 gene. We identified inversions that range from approximately 3 to 15 kb in four taxa (Amischotolype, Belosynapsis, Murdannia, and Streptolirion). The phylogenetic analysis using 77 chloroplast protein-coding genes with maximum parsimony, maximum likelihood, and Bayesian inference suggests that Palisota is most closely related to tribe Commelineae, supported by high support values. This result differs significantly from the current classification of Commelinaceae. Also, we resolved the unclear position of Streptoliriinae and the monophyly of Dichorisandrinae. Among the ten CDS (ndhH, rpoC2, ndhA, rps3, ndhG, ndhD, ccsA, ndhF, matK, and ycf1), which have high nucleotide diversity values (Pi > 0.045) and over 500 bp length, four CDS (ndhH, rpoC2, matK, and ycf1) show that they are congruent with the topology derived from 77 chloroplast protein-coding genes. CONCLUSIONS: In this study, we provide detailed information on the 15 complete plastid genomes of Commelinoideae taxa. We identified characteristic pseudogenes and nucleotide diversity, which can be used to infer the family evolutionary history. Also, further research is needed to revise the position of Palisota in the current classification of Commelinaceae.

BACKGROUND: Night-break (NB) has been proven to repress flowering of short-day plants (SDPs). Long-noncoding RNAs (lncRNAs) play key roles in plant flowering. However, investigation of the relationship between lncRNAs and NB responses is still limited, especially in Chenopodium quinoa, an important short-day coarse cereal. RESULTS: In this study, we performed strand-specific RNA-seq of leaf samples collected from quinoa seedlings treated by SD and NB. A total of 4914 high-confidence lncRNAs were identified, out of which 91 lncRNAs showed specific responses to SD and NB. Based on the expression profiles, we identified 17 positive- and 7 negative-flowering lncRNAs. Co-expression network analysis indicated that 1653 mRNAs were the common targets of both types of flowering lncRNAs. By mapping these targets to the known flowering pathways in model plants, we found some pivotal flowering homologs, including 2 florigen encoding genes (FT (FLOWERING LOCUS T) and TSF (TWIN SISTER of FT) homologs), 3 circadian clock related genes (EARLY FLOWERING 3 (ELF3), LATE ELONGATED HYPOCOTYL (LHY) and ELONGATED HYPOCOTYL 5 (HY5) homologs), 2 photoreceptor genes (PHYTOCHROME A (PHYA) and CRYPTOCHROME1 (CRY1) homologs), 1 B-BOX type CONSTANS (CO) homolog and 1 RELATED TO ABI3/VP1 (RAV1) homolog, were specifically affected by NB and competed by the positive and negative-flowering lncRNAs. We speculated that these potential flowering lncRNAs may mediate quinoa NB responses by modifying the expression of the floral homologous genes. CONCLUSIONS: Together, the findings in this study will deepen our understanding of the roles of lncRNAs in NB responses, and provide valuable information for functional characterization in future.

BACKGROUND: Cytidine monophospho-n-acetylneuraminic acid hydroxylase (CMAH) gene associated with blood groups in cats encodes CMAH enzyme that converts Neu5Ac to Neu5Gc. Although variations in CMAH gene of pedigree cats have been revealed, the presence/lack of them in non-pedigree stray cats is unknown. Therefore, the present study aimed to investigate the variations in CMAH gene and the quantity of Neu5Ac and Neu5Gc on erythrocytes of non-pedigree stray cats (n:12) living in İzmir, Turkey. Also, the frequency of blood types was determined in 76 stray cats including 12 cats that were used for CMAH and Neu5A/Neu5Gc analysis. RESULTS: In total, 14 SNPs were detected in 5'UTR as well as in exon 2, 4, 9, 10, 11 and 12 of CMAH gene. Among these SNPs, -495 C > T in 5'UTR was detected for the first time as heterozygous in type A and AB cats, and homozygous and heterozygous in type B cats. The remaining 13 that have been detected in previous studies were also found as homozygous or heterozygous. Both Neu5Gc and Neu5Ac were detected in type A and AB cats. In type B cats, only Neu5Ac was detected. Among two type AB cats, the level of Neu5Ac was found higher in cat carrying heterozygous form (T/C) of 1392T > C. The prevalence of type B cats (67.1 %) was higher than others. CONCLUSIONS: The presence of a new SNP as well as previous SNPs indicates that more variations can be found in stray cats with a more comprehensive study in the future. Also, the high prevalence of type B cats demonstrates the possible risk of neonatal isoerythrolysis among stray cats living in İzmir, Turkey.

BACKGROUND: Histone lysine methylation plays an important role in plant development and stress responses by activating or repressing gene expression. Histone lysine methylation is catalyzed by a class of SET-domain group proteins (SDGs). Although an increasing number of studies have shown that SDGs play important regulatory roles in development and stress responses, the functions of SDGs in apple remain unclear. RESULTS: A total of 67 SDG members were identified in the Malus×domestica genome. Syntenic analysis revealed that most of the MdSDG duplicated gene pairs were associated with a recent genome-wide duplication event of the apple genome. These 67 MdSDG members were grouped into six classes based on sequence similarity and the findings of previous studies. The domain organization of each MdSDG class was characterized by specific patterns, which was consistent with the classification results. The tissue-specific expression patterns of MdSDGs among the 72 apple tissues in the different apple developmental stages were characterized to provide insight into their potential functions in development. The expression profiles of MdSDGs were also investigated in fruit development, the breaking of bud dormancy, and responses to abiotic and biotic stress; the results indicated that MdSDGs might play a regulatory role in development and stress responses. The subcellular localization and putative interaction network of MdSDG proteins were also analyzed. CONCLUSIONS: This work presents a fundamental comprehensive analysis of SDG histone methyltransferases in apple and provides a basis for future studies of MdSDGs involved in apple development and stress responses.

BACKGROUND: Methionine (Met) supply during late-pregnancy enhances fetal development in utero and leads to greater rates of growth during the neonatal period. Due to its central role in coordinating nutrient and one-carbon metabolism along with immune responses of the newborn, the liver could be a key target of the programming effects induced by dietary methyl donors such as Met. To address this hypothesis, liver biopsies from 4-day old calves (n = 6/group) born to Holstein cows fed a control or the control plus ethyl-cellulose rumen-protected Met for the last 28 days prepartum were used for DNA methylation, transcriptome, metabolome, proteome, and one-carbon metabolism enzyme activities. RESULTS: Although greater withers and hip height at birth in Met calves indicated better development in utero, there were no differences in plasma systemic physiological indicators. RNA-seq along with bioinformatics and transcription factor regulator analyses revealed broad alterations in 'Glucose metabolism', 'Lipid metabolism, 'Glutathione', and 'Immune System' metabolism due to enhanced maternal Met supply. Greater insulin sensitivity assessed via proteomics, and efficiency of transsulfuration pathway activity suggested beneficial effects on nutrient metabolism and metabolic-related stress. Maternal Met supply contributed to greater phosphatidylcholine synthesis in calf liver, with a role in very low density lipoprotein secretion as a mechanism to balance metabolic fates of fatty acids arising from the diet or adipose-depot lipolysis. Despite a lack of effect on hepatic amino acid (AA) transport, a reduction in metabolism of essential AA within the liver indicated an AA 'sparing effect' induced by maternal Met. CONCLUSIONS: Despite greater global DNA methylation, maternal Met supply resulted in distinct alterations of hepatic transcriptome, proteome, and metabolome profiles after birth. Data underscored an effect on maintenance of calf hepatic Met homeostasis, glutathione, phosphatidylcholine and taurine synthesis along with greater efficiency of nutrient metabolism and immune responses. Transcription regulators such as FOXO1, PPARG, E2F1, and CREB1 appeared central in the coordination of effects induced by maternal Met. Overall, maternal Met supply induced better immunometabolic status of the newborn liver, conferring the calf a physiologic advantage during a period of metabolic stress and suboptimal immunocompetence.

BACKGROUND: Heterosis has been exploited for decades in different crops due to resulting in dramatic increases in yield, but relatively little molecular evidence on this topic was reported in cotton. RESULTS: The elite cotton hybrid variety 'Huaza Mian H318' (H318) and its parental lines were used to explore the source of its yield heterosis. A four-year investigation of yield-related traits showed that the boll number of H318 showed higher stability than that of its two parents, both in suitable and unsuitable climate years. In addition, the hybrid H318 grew faster and showed higher fresh and dry weights than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differentially expressed genes (DEGs) between H318 and its parental lines, and 3490 extremely changed DEGs were screened out for later analysis. Most DEGs (3472/3490) were gathered between H318 and its paternal line (4-5), and only 64 DEGs were found between H318 and its maternal line (B0011), which implied that H318 displays more similar transcriptional patterns to its maternal parent at the seedling stage. GO and KEGG analyses showed that these DEGs were highly enriched in photosynthesis, lipid metabolic, carbohydrate metabolic and oxidation-reduction processes, and the expression level of these DEGs was significantly higher in H318 relative to its parental lines, which implied that photosynthesis, metabolism and stress resistances were enhanced in H318. CONCLUSION: The enhanced photosynthesis, lipid and carbohydrate metabolic capabilities contribute to the heterosis of H318 at the seedling stage, and establishes a material foundation for subsequent higher boll-setting rates in complex field environments.

BACKGROUND: Mammalian centromeres are satellite-rich chromatin domains that execute conserved roles in kinetochore assembly and chromosome segregation. Centromere satellites evolve rapidly between species, but little is known about population-level diversity across these loci. RESULTS: We developed a k-mer based method to quantify centromere copy number and sequence variation from whole genome sequencing data. We applied this method to diverse inbred and wild house mouse (Mus musculus) genomes to profile diversity across the core centromere (minor) satellite and the pericentromeric (major) satellite repeat. We show that minor satellite copy number varies more than 10-fold among inbred mouse strains, whereas major satellite copy numbers span a 3-fold range. In contrast to widely held assumptions about the homogeneity of mouse centromere repeats, we uncover marked satellite sequence heterogeneity within single genomes, with diversity levels across the minor satellite exceeding those at the major satellite. Analyses in wild-caught mice implicate subspecies and population origin as significant determinants of variation in satellite copy number and satellite heterogeneity. Intriguingly, we also find that wild-caught mice harbor dramatically reduced minor satellite copy number and elevated satellite sequence heterogeneity compared to inbred strains, suggesting that inbreeding may reshape centromere architecture in pronounced ways. CONCLUSION: Taken together, our results highlight the power of k-mer based approaches for probing variation across repetitive regions, provide an initial portrait of centromere variation across Mus musculus, and lay the groundwork for future functional studies on the consequences of natural genetic variation at these essential chromatin domains.

BACKGROUND: Rice, which serves as a staple food for more than half of the world's population, is grown worldwide. The hybridization of wild and cultivated rice has enabled the incorporation of resistance to varying environmental conditions. Endophytic microbiota are known to be transferred with their host plants. Although some studies have reported on the endophytic microbiota of wild and cultivated rice, the inheritance from wild and cultivated rice accessions in next generations, in terms of endophytic microbiota, has not been examined. RESULTS: In the present study, the endophytic microbial community structures of Asian and African wild and cultivated rice species were compared with those of their F1 offspring. High-throughput sequencing data of bacterial 16S rDNA and fungal internal transcribed spacer regions were used to classify the endophytic microbiota of collected samples of rice. Results indicated that when either African or Asian wild rice species were crossed with cultivated rice accessions, the first generation harbored a greater number of root endophytic fungi than the cultivated parent used to make the crosses. Network analysis of the bacterial and fungal operational taxonomic units revealed that Asian and African wild rice species clustered together and exhibited a greater number of significant correlations between fungal taxa than cultivated rice. The core bacterial genus Acidovorax and the core fungal order Pleosporales, and genera Myrothecium and Bullera connected African and Asian wild rice accessions together, and both the wild rice accessions with their F1 offspring. On the other hand, the core bacterial genus Bradyrhizobium and the core fungal genera Dendroclathra linked the African and Asian cultivated rice accessions together. CONCLUSIONS: This study has theoretical significance for understanding the effect of breeding on the inheritance of endophytic microbiota of rice and identifying beneficial endophytic bacteria and fungi among wild and cultivated rice species, and their F1 offspring.

BACKGROUND: Sepsis is the major cause of death in Intensive Care Unit (ICU) globally. Molecular detection enables rapid diagnosis that allows early intervention to minimize the death rate. Recent studies showed that long non-coding RNAs (lncRNAs) regulate proinflammatory genes and are related to the dysfunction of organs in sepsis. Identifying lncRNA signature with absolute abundance is challenging because of the technical variation and the systematic experimental bias. RESULTS: Cohorts (n = 768) containing whole blood lncRNA profiling of sepsis patients in the Gene Expression Omnibus (GEO) database were included. We proposed a novel diagnostic strategy that made use of the relative expressions of lncRNA pairs, which are reversed between sepsis patients and normal controls (eg. lncRNAi > lncRNAj in sepsis patients and lncRNAi < lncRNAj in normal controls), to identify 14 lncRNA pairs as a sepsis diagnostic signature. The signature was then applied to independent cohorts (n = 644) to evaluate its predictive performance across different ages and normalization methods. Comparing to common machine learning models and existing signatures, SepSigLnc consistently attains better performance on the validation cohorts from the same age group (AUC = 0.990 & 0.995 in two cohorts) and across different groups (AUC = 0.878 on average), as well as cohorts processed by an alternative normalization method (AUC = 0.953 on average). Functional analysis demonstrates that the lncRNA pairs in SepsigLnc are functionally similar and tend to implicate in the same biological processes including cell fate commitment and cellular response to steroid hormone stimulus. CONCLUSION: Our study identified 14 lncRNA pairs as signature that can facilitate the diagnosis of septic patients at an intervenable point when clinical manifestations are not dramatic. Also, the computational procedure can be generalized to a standard procedure for discovering diagnostic molecule signatures.

BACKGROUND: Stem rot caused by Sclerotium rolfsii is a very important soil-borne disease of peanut. S. rolfsii is a necrotrophic plant pathogenic fungus with an extensive host range and worldwide distribution. It can infect peanut stems, roots, pegs and pods, leading to varied yield losses. S. rolfsii strains GP3 and ZY collected from peanut in different provinces of China exhibited a significant difference in aggressiveness on peanut plants by artificial inoculation test. In this study, de-novo genome sequencing of these two distinct strains was performed aiming to reveal the genomic basis of difference in aggressiveness. RESULTS: Scleotium rolfsii strains GP3 and ZY, with weak and high aggressiveness on peanut plants, exhibited similar growth rate and oxalic acid production in laboratory. The genomes of S. rolfsii strains GP3 and ZY were sequenced by Pacbio long read technology and exhibited 70.51 Mb and 70.61 Mb, with contigs of 27 and 23, and encoded 17,097 and 16,743 gene models, respectively. Comparative genomic analysis revealed that the pathogenicity-related gene repertoires, which might be associated with aggressiveness, differed between GP3 and ZY. There were 58 and 45 unique pathogen-host interaction (PHI) genes in GP3 and ZY, respectively. The ZY strain had more carbohydrate-active enzymes (CAZymes) in its secretome than GP3, especially in the glycoside hydrolase family (GH), the carbohydrate esterase family (CBM), and the polysaccharide lyase family (PL). GP3 and ZY also had different effector candidates and putative secondary metabolite synthetic gene clusters. These results indicated that differences in PHI, secreted CAZymes, effectors and secondary metabolites may play important roles in aggressive difference between these two strains. CONCLUSIONS: The data provided a further understanding of the S. rolfsii genome. Genomic comparison provided clues to the difference in aggressiveness of S. rolfsii strains.

BACKGROUND: Ectoparasites from the family Diplozoidae (Platyhelminthes, Monogenea) belong to obligate haematophagous helminths of cyprinid fish. Current knowledge of these worms is for the most part limited to their morphological, phylogenetic, and population features. Information concerning the biochemical and molecular nature of physiological processes involved in host-parasite interaction, such as evasion of the immune system and its regulation, digestion of macromolecules, suppression of blood coagulation and inflammation, and effect on host tissue and physiology, is lacking. In this study, we report for the first time a comprehensive transcriptomic/secretome description of expressed genes and proteins secreted by the adult stage of Eudiplozoon nipponicum (Goto, 1891) Khotenovsky, 1985, an obligate sanguivorous monogenean which parasitises the gills of the common carp (Cyprinus carpio). RESULTS: RNA-seq raw reads (324,941 Roche 454 and 149,697,864 Illumina) were generated, de novo assembled, and filtered into 37,062 protein-coding transcripts. For 19,644 (53.0%) of them, we determined their sequential homologues. In silico functional analysis of E. nipponicum RNA-seq data revealed numerous transcripts, pathways, and GO terms responsible for immunomodulation (inhibitors of proteolytic enzymes, CD59-like proteins, fatty acid binding proteins), feeding (proteolytic enzymes cathepsins B, D, L1, and L3), and development (fructose 1,6-bisphosphatase, ferritin, and annexin). LC-MS/MS spectrometry analysis identified 721 proteins secreted by E. nipponicum with predominantly immunomodulatory and anti-inflammatory functions (peptidyl-prolyl cis-trans isomerase, homolog to SmKK7, tetraspanin) and ability to digest host macromolecules (cathepsins B, D, L1). CONCLUSIONS: In this study, we integrated two high-throughput sequencing techniques, mass spectrometry analysis, and comprehensive bioinformatics approach in order to arrive at the first comprehensive description of monogenean transcriptome and secretome. Exploration of E. nipponicum transcriptome-related nucleotide sequences and translated and secreted proteins offer a better understanding of molecular biology and biochemistry of these, often neglected, organisms. It enabled us to report the essential physiological pathways and protein molecules involved in their interactions with the fish hosts.

BACKGROUND: Messenger RNA sequencing is becoming more common in studies of non-model species and is most often used for gene expression-based investigations. However, the method holds potential for numerous other applications as well-including analyses of alternative splicing, population structure, and signatures of selection. To maximize the utility of mRNA data sets, distinct analyses may be combined such as by exploring dynamics between gene expression with signatures of selection in the context of population structure. Here, we compare two published data sets describing two populations of a minnow species endemic to the San Francisco Estuary (Sacramento splittail, Pogonichthys macrolepidotus): a microsatellite data set showing population structure, and an mRNA whole transcriptome data set obtained after the two populations were exposed to a salinity challenge. We compared measures of population structure and genetic variation using single nucleotide polymorphisms (SNPs) called from mRNA from the whole transcriptome sequencing study with those patterns determined from microsatellites. For investigating plasticity and evolution, intra- and inter-population transcriptome plasticity was investigated with differential gene expression, differential exon usage, and gene expression variation. Outlier SNP analysis was also performed on the mRNA data set and signatures of selection and phenotypic plasticity were investigated on an individual-gene basis. RESULTS: We found that mRNA sequencing revealed patterns of population structure consistent with those found with microsatellites, but with lower magnitudes of genetic variation and population differentiation consistent with widespread purifying selection expected when using mRNA. In addition, within individual genes, phenotypic plasticity or signatures of selection were found in almost mutual exclusion (except heatr6, nfu1, slc22a6, sya, and mmp13). CONCLUSIONS: These results show that an mRNA sequencing data set may have multiple uses, including describing population structure and for investigating the mechanistic interplay of evolution and plasticity in adaptation. MRNA sequencing thus complements traditional sequencing methods used for population genetics, in addition to its utility for describing phenotypic plasticity.

BACKGROUND: Infestation by tea green leafhoppers (Empoasca (Matsumurasca) onukii) can cause a series of biochemical changes in tea leaves. As a typical cell-rupture feeder, E. onukii secretes proteases while using its stylet to probe the tender shoots of tea plants (Camellia sinensis). This study identified and analyzed proteases expressed specifically in the salivary gland (SG) and gut of E. onukii through enzymatic activity assays complemented with an integrated analysis of transcriptomic and proteomic data. RESULTS: In total, 129 contigs representing seven types of putative proteases were identified. Transcript abundance of digestive proteases and enzymatic activity assays showed that cathepsin B-like protease, cathepsin L-like protease, and serine proteases (trypsin- and chymotrypsin-like protease) were highly abundant in the gut but moderately abundant in the SG. The abundance pattern of digestive proteases in the SG and gut of E. onukii differed from that of other hemipterans, including Nilaparvata lugens, Laodelphax striatellus, Acyrthosiphum pisum, Halyomorpha halys and Nephotettix cincticeps. Phylogenetic analysis showed that aminopeptidase N-like proteins and serine proteases abundant in the SG or gut of hemipterans formed two distinct clusters. CONCLUSIONS: Altogether, this study provides insightful information on the digestive system of E. onukii. Compared to five other hemipteran species, we observed different patterns of proteases abundant in the SG and gut of E. onukii. These results will be beneficial in understanding the interaction between tea plants and E. onukii.

BACKGROUND: Human cancer cell line profiling and drug sensitivity studies provide valuable information about the therapeutic potential of drugs and their possible mechanisms of action. The goal of those studies is to translate the findings from in vitro studies of cancer cell lines into in vivo therapeutic relevance and, eventually, patients' care. Tremendous progress has been made. RESULTS: In this work, we built predictive models for 453 drugs using data on gene expression and drug sensitivity (IC50) from cancer cell lines. We identified many known drug-gene interactions and uncovered several potentially novel drug-gene associations. Importantly, we further applied these predictive models to ~ 17,000 bulk RNA-seq samples from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) database to predict drug sensitivity for both normal and tumor tissues. We created a web site for users to visualize and download our predicted data ( https://manticore.niehs.nih.gov/cancerRxTissue ). Using trametinib as an example, we showed that our approach can faithfully recapitulate the known tumor specificity of the drug. CONCLUSIONS: We demonstrated that our approach can predict drugs that 1) are tumor-type specific; 2) elicit higher sensitivity from tumor compared to corresponding normal tissue; 3) elicit differential sensitivity across breast cancer subtypes. If validated, our prediction could have relevance for preclinical drug testing and in phase I clinical design.

BACKGROUND: Bacterial plant pathogens of the Pectobacterium genus are responsible for a wide spectrum of diseases in plants, including important crops such as potato, tomato, lettuce, and banana. Investigation of the genetic diversity underlying virulence and host specificity can be performed at genome level by using a comprehensive comparative approach called pangenomics. A pangenomic approach, using newly developed functionalities in PanTools, was applied to analyze the complex phylogeny of the Pectobacterium genus. We specifically used the pangenome to investigate genetic differences between virulent and avirulent strains of P. brasiliense, a potato blackleg causing species dominantly present in Western Europe. RESULTS: Here we generated a multilevel pangenome for Pectobacterium, comprising 197 strains across 19 species, including type strains, with a focus on P. brasiliense. The extensive phylogenetic analysis of the Pectobacterium genus showed robust distinct clades, with most detail provided by 452,388 parsimony-informative single-nucleotide polymorphisms identified in single-copy orthologs. The average Pectobacterium genome consists of 47% core genes, 1% unique genes, and 52% accessory genes. Using the pangenome, we zoomed in on differences between virulent and avirulent P. brasiliense strains and identified 86 genes associated to virulent strains. We found that the organization of genes is highly structured and linked with gene conservation, function, and transcriptional orientation. CONCLUSION: The pangenome analysis demonstrates that evolution in Pectobacteria is a highly dynamic process, including gene acquisitions partly in clusters, genome rearrangements, and loss of genes. Pectobacterium species are typically not characterized by a set of species-specific genes, but instead present themselves using new gene combinations from the shared gene pool. A multilevel pangenomic approach, fusing DNA, protein, biological function, taxonomic group, and phenotypes, facilitates studies in a flexible taxonomic context.

BACKGROUND: The number of myofiber is determined during the embryonic stage and does not increase during the postnatal period for birds, including goose. Thus, muscle production of adult goose is pre-determined during embryogenesis. Previous studies show N6-methyladenosine (m6A) is an important regulator for skeletal muscle development of birds and miRNAs play as a co-regulator for the skeletal muscle development in birds. Herein, we sequenced m6A and miRNA transcriptomes to investigate the profiles of m6A and their potential mechanism of regulating breast muscle development in Dingan Goose. RESULTS: We selected embryonic 21th day (E21) and embryonic 30th day (E30) to investigate the roles of transcriptome-wide m6A modification combining with mRNAs and miRNAs in goose breast muscle development. In this study, m6A peaks were mainly enriched in coding sequence (CDS) and start codon and397 genes were identified as differentially methylated genes (DMGs). GO and KEGG analysis showed that DMGs were highly related to cellular and metabolic process and that most DMGs were enriched in muscle-related pathways including Wnt signaling pathway, mTOR signaling and FoxO signaling pathway. Interestingly, a negative correlation between m6A methylation level and mRNA abundance was found through the analysis of m6A-RNA and RNA-seq data. Besides, we found 26 muscle-related genes in 397 DMGs. We also detected 228 differentially expressed miRNAs (DEMs), and further found 329 genes shared by the target genes of DEMs and DMGs (m6A-miRNA-genes), suggesting a tightly relationship between DEMs and DMGs. Among the m6A-miRNA-genes, we found 10 genes are related to breast muscle development. We further picked out an m6A-miRNA-gene, PDK3, from the 10 genes to visualize it and the result showed differentially methylated peaks on the mRNA transcript consistent with our m6A-seq results. CONCLUSION: GO and KEGG of DMGs between E21 and E30 showed most DMGs were muscle-related. In total, 228 DEMs were found, and the majority of DMGs were overlapped with the targets of DEGs. The differentially methylated peaks along with an m6A-miRNA-gene, PDK3, showed the similar results with m6A-seq results. Taken together, the results presented here provide a reference for further investigation of embryonic skeletal muscle development mechanism in goose.

BACKGROUND: The back plays a vital role in horse locomotion, where the spine functions as a spring during the stride cycle. A complex interaction between the spine and the muscles of the back contribute to locomotion soundness, gait ability, and performance of riding and racehorses. Conformation is commonly used to select horses for breeding and performance in multiple horse breeds, where the back and croup conformation plays a significant role. The conformation of back and croup plays an important role on riding ability in Icelandic horses. However, the genes behind this trait are still unknown. Therefore, the aim of this study was to identify genomic regions associated with conformation of back and croup in Icelandic horses and to investigate their effects on riding ability. One hundred seventy-seven assessed Icelandic horses were included in the study. A genome-wide association analysis was performed using the 670 K+ Axiom Equine Genotyping Array, and the effects of different haplotypes in the top associated region were estimated for riding ability and additional conformation traits assessed during breeding field tests. RESULTS: A suggestive quantitative trait loci (QTL) for the score of back and croup was detected on Equus caballus (ECA) 22 (p-value = 2.67 × 10- 7). Haplotype analysis revealed two opposite haplotypes, which resulted in higher and lower scores of the back and croup, respectively (p-value < 0.001). Horses with the favorable haplotype were more inclined to have a well-balanced backline with an uphill conformation and had, on average, higher scores for the lateral gaits tölt (p-value = 0.02) and pace (p-value = 0.004). This genomic region harbors three genes: C20orf85, ANKRD60 and LOC100056167. ANKRD60 is associated with body height in humans. C20orf85 and ANKRD60 are potentially linked to adolescent idiopathic scoliosis in humans. CONCLUSIONS: Our results show that the detected QTL for conformation of back and croup is of importance for quality of lateral gaits in Icelandic horses. These findings could result in a genetic test to aid in the selection of breeding horses, thus they are of major interest for horse breeders. The results may also offer a gateway to comparative functional genomics by potentially linking both motor laterality and back inclination in horses with scoliosis in humans.

BACKGROUND: Fertilization in flowering plants depends on the early contact and acceptance of pollen grains by the receptive papilla cells of the stigma. Deciphering the specific transcriptomic response of both pollen and stigmatic cells during their interaction constitutes an important challenge to better our understanding of this cell recognition event. RESULTS: Here we describe a transcriptomic analysis based on single nucleotide polymorphisms (SNPs) present in two Arabidopsis thaliana accessions, one used as female and the other as male. This strategy allowed us to distinguish 80% of transcripts according to their parental origins. We also developed a tool which predicts male/female specific expression for genes without SNP. We report an unanticipated transcriptional activity triggered in stigma upon incompatible pollination and show that following compatible interaction, components of the pattern-triggered immunity (PTI) pathway are induced on the female side. CONCLUSIONS: Our work unveils the molecular signatures of compatible and incompatible pollinations both at the male and female side. We provide invaluable resource and tools to identify potential new molecular players involved in pollen-stigma interaction.

BACKGROUND: The pathogen Listeria (L.) monocytogenes is known to survive heat, cold, high pressure, and other extreme conditions. Although the response of this pathogen to pH, osmotic, temperature, and oxidative stress has been studied extensively, its reaction to the stress produced by high pressure processing HPP (which is a preservation method in the food industry), and the activated gene regulatory network (GRN) in response to this stress is still largely unknown. RESULTS: We used RNA sequencing transcriptome data of L. monocytogenes (ScottA) treated at 400 MPa and 8∘C, for 8 min and combined it with current information in the literature to create a transcriptional regulation database, depicting the relationship between transcription factors (TFs) and their target genes (TGs) in L. monocytogenes. We then applied network component analysis (NCA), a matrix decomposition method, to reconstruct the activities of the TFs over time. According to our findings, L. monocytogenes responded to the stress applied during HPP by three statistically different gene regulation modes: survival mode during the first 10 min post-treatment, repair mode during 1 h post-treatment, and re-growth mode beyond 6 h after HPP. We identified the TFs and their TGs that were responsible for each of the modes. We developed a plausible model that could explain the regulatory mechanism that L. monocytogenes activated through the well-studied CIRCE operon via the regulator HrcA during the survival mode. CONCLUSIONS: Our findings suggest that the timely activation of TFs associated with an immediate stress response, followed by the expression of genes for repair purposes, and then re-growth and metabolism, could be a strategy of L. monocytogenes to survive and recover extreme HPP conditions. We believe that our results give a better understanding of L. monocytogenes behavior after exposure to high pressure that may lead to the design of a specific knock-out process to target the genes or mechanisms. The results can help the food industry select appropriate HPP conditions to prevent L. monocytogenes recovery during food storage.

BACKGROUND: Uncovering the mechanisms underlying rapid genetic adaptation can provide insight into adaptive evolution and shed light on conservation, invasive species control, and natural resource management. However, it can be difficult to experimentally explore rapid adaptation due to the challenges associated with propagating and maintaining species in captive environments for long periods of time. By contrast, many introduced species have experienced strong selection when colonizing environments that differ substantially from their native range and thus provide a "natural experiment" for studying rapid genetic adaptation. One such example occurred when sea lamprey (Petromyzon marinus), native to the northern Atlantic, naturally migrated into Lake Champlain and expanded their range into the Great Lakes via man-made shipping canals. RESULTS: Utilizing 368,886 genome-wide single nucleotide polymorphisms (SNPs), we calculated genome-wide levels of genetic diversity (i.e., heterozygosity and π) for sea lamprey collected from native (Connecticut River), native but recently colonized (Lake Champlain), and invasive (Lake Michigan) populations, assessed genetic differentiation between all populations, and identified candidate genes that responded to selection imposed by the novel environments. We observed a 14 and 24% reduction in genetic diversity in Lake Michigan and Lake Champlain populations, respectively, compared to individuals from the Connecticut River, suggesting that sea lamprey populations underwent a genetic bottleneck during colonization. Additionally, we identified 121 and 43 outlier genes in comparisons between Lake Michigan and Connecticut River and between Lake Champlain and Connecticut River, respectively. Six outlier genes that contained synonymous SNPs in their coding regions and two genes that contained nonsynonymous SNPs may underlie the rapid evolution of growth (i.e., GHR), reproduction (i.e., PGR, TTC25, STARD10), and bioenergetics (i.e., OXCT1, PYGL, DIN4, SLC25A15). CONCLUSIONS: By identifying the genomic basis of rapid adaptation to novel environments, we demonstrate that populations of invasive species can be a useful study system for understanding adaptive evolution. Furthermore, the reduction in genome-wide levels of genetic diversity associated with colonization coupled with the identification of outlier genes underlying key life history traits known to have changed in invasive sea lamprey populations (e.g., growth, reproduction) illustrate the utility in applying genomic approaches for the successful management of introduced species.

BACKGROUND: Protein content determines the state of cells. The variation in protein abundance is crucial when organisms are in the early stages of heat stress, but the reasons affecting their changes are largely unknown. RESULTS: We quantified 47,535 mRNAs and 3742 proteins in the filling grains of wheat in two different thermal environments. The impact of mRNA abundance and sequence features involved in protein translation and degradation on protein expression was evaluated by regression analysis. Transcription, codon usage and amino acid frequency were the main drivers of changes in protein expression under heat stress, and their combined contribution explains 58.2 and 66.4% of the protein variation at 30 and 40 °C (20 °C as control), respectively. Transcription contributes more to alterations in protein content at 40 °C (31%) than at 30 °C (6%). Furthermore, the usage of codon AAG may be closely related to the rapid alteration of proteins under heat stress. The contributions of AAG were 24 and 13% at 30 and 40 °C, respectively. CONCLUSION: In this study, we analyzed the factors affecting the changes in protein expression in the early stage of heat stress and evaluated their influence.

Erratum in BMC Genomics. 2021 May 20;22(1):368.

BACKGROUND: Polyploid species often originate recurrently. While this is well known, there is little information on the extent to which distinct allotetraploid species formed from the same parent species differ in gene expression. The tetraploid yarrow species Achillea alpina and A. wilsoniana arose independently from allopolyploidization between diploid A. acuminata and A. asiatica. The genetics and geography of these origins are clear from previous studies, providing a solid basis for comparing gene expression patterns of sibling allopolyploid species that arose independently. RESULTS: We conducted comparative RNA-sequencing analyses on the two Achillea tetraploid species and their diploid progenitors to evaluate: 1) species-specific gene expression and coexpression across the four species; 2) patterns of inheritance of parental gene expression; 3) parental contributions to gene expression in the allotetraploid species, and homeolog expression bias. Diploid A. asiatica showed a higher contribution than diploid A. acuminata to the transcriptomes of both tetraploids and also greater homeolog bias in these transcriptomes, possibly reflecting a maternal effect. Comparing expressed genes in the two allotetraploids, we found expression of ca. 30% genes were species-specific in each, which were most enriched for GO terms pertaining to "defense response". Despite species-specific and differentially expressed genes between the two allotetraploids, they display similar transcriptome changes in comparison to their diploid progenitors. CONCLUSION: Two independently originated Achillea allotetraploid species exhibited difference in gene expression, some of which must be related to differential adaptation during their post-speciation evolution. On the other hand, they showed similar expression profiles when compared to their progenitors. This similarity might be expected when pairs of merged diploid genomes in tetraploids are similar, as is the case in these two particular allotetraploids.

This correspondence responds to the critique by Butler et al. (BMC Genomics 22:241, 2021) of our recent paper on transposable element (TE) persistence. We address the three main objections raised by Butler et al. After running a series of additional simulations that were inspired by the authors' criticisms, we are able to present a more nuanced understanding of the conditions that generate long-term persistence.

BACKGROUND: Leaf color mutants are the ideal materials to explore the pathways of chlorophyll (Chl) metabolism, chloroplast development, and photosynthesis system. In this study, a spontaneous yellow-green leaf wucai (Brassica campestris L.) mutant "WY16-13" was identified, which exhibited yellow-green leaf color during its entire growth period. However, current understanding of the molecular mechanism underlying Chl metabolism and chloroplast development of "WY16-13" is limited. RESULTS: Total Chl and carotenoid content in WY16-13 was reduced by 60.92 and 58.82%, respectively, as compared with its wild type parental line W16-13. Electron microscopic investigation revealed fewer chloroplasts per cell and looser stroma lamellae in WY16-13 than in W16-13. A comparative transcriptome profiling was performed using leaves from the yellow-green leaf type (WY16-13) and normal green-leaf type (W16-13). A total of 54.12 million (M) (WY16-13) and 56.17 M (W16-13) reads were generated. A total of 40,578 genes were identified from the mapped libraries. We identified 3882 differentially expressed genes (DEGs) in WY16-13 compared with W16-13 (i.e., 1603 upregulated genes and 2279 downregulated genes). According to the Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, these DEGs are involved in porphyrin and Chl metabolism [i.e., chlorophyllase (CLH), heme oxygenase (HO), chlorophyll (ide) b reductase (NYC), and protochlorophyllide oxidoreductase (POR) genes], carbohydrate metabolism, photosynthesis, and carbon fixation in photosynthetic organisms. Moreover, deficiency in Chl biosynthetic intermediates in WY16-13 revealed that the formation of the yellow-green phenotype was related to the disorder of heme metabolism. CONCLUSIONS: Our results provide valuable insights into Chl deficiency in the yellow-green leaf mutant and a bioinformatics resource for further functional identification of key allelic genes responsible for differences in Chl content.

BACKGROUND: Prior work demonstrated that female rats (but not their male littermates) exposed to methamphetamine become hypersensitive to myocardial ischemic injury. Importantly, this sex-dependent effect persists following 30 days of subsequent abstinence from the drug, suggesting that it may be mediated by long term changes in gene expression that are not rapidly reversed following discontinuation of methamphetamine use. The goal of the present study was to determine whether methamphetamine induces sex-dependent changes in myocardial gene expression and whether these changes persist following subsequent abstinence from methamphetamine. RESULTS: Methamphetamine induced changes in the myocardial transcriptome were significantly greater in female hearts than male hearts both in terms of the number of genes affected and the magnitude of the changes. The largest changes in female hearts involved genes that regulate the circadian clock (Dbp, Per3, Per2, BMal1, and Npas2) which are known to impact myocardial ischemic injury. These genes were unaffected by methamphetamine in male hearts. All changes in gene expression identified at day 11 returned to baseline by day 30. CONCLUSIONS: These data demonstrate that female rats are more sensitive than males to methamphetamine-induced changes in the myocardial transcriptome and that methamphetamine does not induce changes in myocardial transcription that persist long term after exposure to the drug has been discontinued.

BACKGROUND: Increases in ocean temperatures and in the frequency and severity of hypoxic events are expected with climate change, and may become a challenge for cultured Atlantic salmon and negatively affect their growth, immunology and welfare. Thus, we examined how an incremental temperature increase alone (Warm & Normoxic-WN: 12 → 20 °C; 1 °C week- 1), and in combination with moderate hypoxia (Warm & Hypoxic-WH: ~ 70% air saturation), impacted the salmon's hepatic transcriptome expr\ession compared to control fish (CT: 12 °C, normoxic) using 44 K microarrays and qPCR. RESULTS: Overall, we identified 2894 differentially expressed probes (DEPs, FDR < 5%), that included 1111 shared DEPs, while 789 and 994 DEPs were specific to WN and WH fish, respectively. Pathway analysis indicated that the cellular mechanisms affected by the two experimental conditions were quite similar, with up-regulated genes functionally associated with the heat shock response, ER-stress, apoptosis and immune defence, while genes connected with general metabolic processes, proteolysis and oxidation-reduction were largely suppressed. The qPCR assessment of 41 microarray-identified genes validated that the heat shock response (hsp90aa1, serpinh1), apoptosis (casp8, jund, jak2) and immune responses (apod, c1ql2, epx) were up-regulated in WN and WH fish, while oxidative stress and hypoxia sensitive genes were down-regulated (cirbp, cyp1a1, egln2, gstt1, hif1α, prdx6, rraga, ucp2). However, the additional challenge of hypoxia resulted in more pronounced effects on heat shock and immune-related processes, including a stronger influence on the expression of 14 immune-related genes. Finally, robust correlations between the transcription of 19 genes and several phenotypic traits in WH fish suggest that changes in gene expression were related to impaired physiological and growth performance. CONCLUSION: Increasing temperature to 20 °C alone, and in combination with hypoxia, resulted in the differential expression of genes involved in similar pathways in Atlantic salmon. However, the expression responses of heat shock and immune-relevant genes in fish exposed to 20 °C and hypoxia were more affected, and strongly related to phenotypic characteristics (e.g., growth). This study provides valuable information on how these two environmental challenges affect the expression of stress-, metabolic- and immune-related genes and pathways, and identifies potential biomarker genes for improving our understanding of fish health and welfare.

BACKGROUND: Jasmonates (JAs) are important for plants to coordinate growth, reproduction, and defense responses. In JA signaling, jasmonate ZIM-domain (JAZ) proteins serve as master regulators at the initial stage of herbivores attacks. Although discovered in many plant species, little in-depth characterization of JAZ gene expression has been reported in the agronomically important crop, maize (Zea mays L.). RESULTS: In this study 16 JAZ genes from the maize genome were identified and classified. Phylogenetic analyses were performed from maize, rice, sorghum, Brachypodium, and Arabidopsis using deduced protein sequences, total six clades were proposed and conservation was observed in each group, such as similar gene exon/intron structures. Synteny analysis across four monocots indicated these JAZ gene families had a common ancestor, and duplication events in maize genome may drive the expansion of JAZ gene family, including genome-wide duplication (GWD), transposon, and/or tandem duplication. Strong purifying selection acted on all JAZ genes except those in group 4, which were under neutral selection. Further, we cloned three paralogous JAZ gene pairs from two maize inbreds differing in JA levels and insect resistance, and gene polymorphisms were observed between two inbreds. CONCLUSIONS: Here we analyzed the composition and evolution of JAZ genes in maize with three other monocot plants. Extensive phylogenetic and synteny analysis revealed the expansion and selection fate of maize JAZ. This is the first study comparing the difference between two inbreds, and we propose genotype-specific JAZ gene expression might be present in maize plants. Since genetic redundancy in JAZ gene family hampers our understanding of their role in response to specific elicitors, we hope this research could be pertinent to elucidating the defensive responses in plants.

BACKGROUND: The differentiation of bone marrow mesenchymal stem cells is a complex and dynamic process. The gene expression pattern and mechanism of different periods of adipogenic and osteogenic differentiation remain unclear. Additionally, the interaction between these two lineage determination requires further exploration. RESULTS: Five modules that were most significantly associated with osteogenic or adipogenic differentiation of BMSCs were selected for further investigation. Biological terms (e.g. ribosome biogenesis, TNF-α signalling pathway, glucose import and fatty acid metabolism) along with hub transcription factors (e.g. PPARG and YY1) and hub miRNAs (e.g. hsa-mir-26b-5p) were enriched in different modules. The expression pattern of 6 hub genes, ADIPOQ, FABP4, SLC7A5, SELPLG, BIRC3, and KLHL30 was validated by RT-qPCR. Finally, cell staining experiments extended the findings of bioinformatics analysis. CONCLUSION: This study identified the key genes, biological functions, and regulators of each time point of adipogenic and osteogenic differentiation of BMSCs and provided novel evidence and ideas for further research on the differentiation of BMSCs.

BACKGROUND: Lysine succinylation is a naturally occurring post-translational modification (PTM) that is ubiquitous in organisms. Lysine succinylation plays important roles in regulating protein structure and function as well as cellular metabolism. Global lysine succinylation at the proteomic level has been identified in a variety of species; however, limited information on lysine succinylation in plant species, especially paper mulberry, is available. Paper mulberry is not only an important plant in traditional Chinese medicine, but it is also a tree species with significant economic value. Paper mulberry is found in the temperate and tropical zones of China. The present study analyzed the effects of lysine succinylation on the growth, development, and physiology of paper mulberry. RESULTS: A total of 2097 lysine succinylation sites were identified in 935 proteins associated with the citric acid cycle (TCA cycle), glyoxylic acid and dicarboxylic acid metabolism, ribosomes and oxidative phosphorylation; these pathways play a role in carbon fixation in photosynthetic organisms and may be regulated by lysine succinylation. The modified proteins were distributed in multiple subcellular compartments and were involved in a wide variety of biological processes, such as photosynthesis and the Calvin-Benson cycle. CONCLUSION: Lysine-succinylated proteins may play key regulatory roles in metabolism, primarily in photosynthesis and oxidative phosphorylation, as well as in many other cellular processes. In addition to the large number of succinylated proteins associated with photosynthesis and oxidative phosphorylation, some proteins associated with the TCA cycle are succinylated. Our study can serve as a reference for further proteomics studies of the downstream effects of succinylation on the physiology and biochemistry of paper mulberry.

BACKGROUND: Non-Drosophila insects provide diverse aging types and important complementary systems for studies of aging biology. However, little attention has been paid to the special roles of non-Drosophila insects in aging research. Here, the aging-related features of the migratory locust, Locusta migratoria, were determined at the physiological, cellular, and transcriptional levels. RESULTS: In physiological assessments, the flight performance and sperm state of locusts displayed clear aging-related decline in male adults. Transcriptional analyses demonstrated locusts have similar aging-related genes with model species. However, different from those of Drosophila and mammals, the organ-specific aging transcriptional features of locusts were characterized by intensive expression changes in flight muscle and fat body and little transcriptional changes in brain. The predominant transcriptional characteristics of flight muscle and fat body aging were changes in expression of mitochondrion-related genes and detoxification and phagocytosis genes, respectively. Cellular assessments revealed the incidence of mitochondrial abnormalities significantly increased in aged flight muscle, and apoptotic signals and nuclear abnormalities were enhanced in aged fat body but not in brain. In addition, some well-known aging genes and locust aging-related genes (i.e., IAP1, PGRP-SA, and LIPT1), whose roles in aging regulation were rarely reported, were demonstrated to affect lifespan, metabolism, and flight ability of locusts after RNAi. CONCLUSION: This study revealed multi-level aging signatures of locust, thus laying a foundation for further investigation of aging mechanisms in this famous insect in the future.

BACKGROUND: Olive orchards are threatened by a wide range of pathogens. Of these, Verticillium dahliae has been in the spotlight for its high incidence, the difficulty to control it and the few cultivars that has increased tolerance to the pathogen. Disease resistance not only depends on detection of pathogen invasion and induction of responses by the plant, but also on barriers to avoid the invasion and active resistance mechanisms constitutively expressed in the absence of the pathogen. In a previous work we found that two healthy non-infected plants from cultivars that differ in V. dahliae resistance such as 'Frantoio' (resistant) and 'Picual' (susceptible) had a different root morphology and gene expression pattern. In this work, we have addressed the issue of basal differences in the roots between Resistant and Susceptible cultivars. RESULTS: The gene expression pattern of roots from 29 olive cultivars with different degree of resistance/susceptibility to V. dahliae was analyzed by RNA-Seq. However, only the Highly Resistant and Extremely Susceptible cultivars showed significant differences in gene expression among various groups of cultivars. A set of 421 genes showing an inverse differential expression level between the Highly Resistant to Extremely Susceptible cultivars was found and analyzed. The main differences involved higher expression of a series of transcription factors and genes involved in processes of molecules importation to nucleus, plant defense genes and lower expression of root growth and development genes in Highly Resistant cultivars, while a reverse pattern in Moderately Susceptible and more pronounced in Extremely Susceptible cultivars were observed. CONCLUSION: According to the different gene expression patterns, it seems that the roots of the Extremely Susceptible cultivars focus more on growth and development, while some other functions, such as defense against pathogens, have a higher expression level in roots of Highly Resistant cultivars. Therefore, it seems that there are constitutive differences in the roots between Resistant and Susceptible cultivars, and that susceptible roots seem to provide a more suitable environment for the pathogen than the resistant ones.

BACKGROUND: Barley is one of the founder crops of Neolithic agriculture and is among the most-grown cereals today. The only trait that universally differentiates the cultivated and wild subspecies is 'non-brittleness' of the rachis (the stem of the inflorescence), which facilitates harvesting of the crop. Other phenotypic differences appear to result from facultative or regional selective pressures. The population structure resulting from these regional events has been interpreted as evidence for multiple domestications or a mosaic ancestry involving genetic interaction between multiple wild or proto-domesticated lineages. However, each of the three mutations that confer non-brittleness originated in the western Fertile Crescent, arguing against multiregional origins for the crop. RESULTS: We examined exome data for 310 wild, cultivated and hybrid/feral barley accessions and showed that cultivated barley is structured into six genetically-defined groups that display admixture, resulting at least in part from two or more significant passages of gene flow with distinct wild populations. The six groups are descended from a single founding population that emerged in the western Fertile Crescent. Only a few loci were universally targeted by selection, the identity of these suggesting that changes in seedling emergence and pathogen resistance could represent crucial domestication switches. Subsequent selection operated on a regional basis and strongly contributed to differentiation of the genetic groups. CONCLUSIONS: Identification of genetically-defined groups provides clarity to our understanding of the population history of cultivated barley. Inference of population splits and mixtures together with analysis of selection sweeps indicate descent from a single founding population, which emerged in the western Fertile Crescent. This founding population underwent relatively little genetic selection, those changes that did occur affecting traits involved in seedling emergence and pathogen resistance, indicating that these phenotypes should be considered as 'domestication traits'. During its expansion out of the western Fertile Crescent, the crop underwent regional episodes of gene flow and selection, giving rise to a modern genetic signature that has been interpreted as evidence for multiple domestications, but which we show can be rationalized with a single origin.

BACKGROUND: The Chinese giant salamander Andrias davidianus is an important amphibian species in China because of its increasing economic value, protection status and special evolutionary position from aquatic to terrestrial animal. Its large genome presents challenges to genetic research. Genetic linkage mapping is an important tool for genome assembly and determination of phenotype-related loci. RESULTS: In this study, we constructed a high-density genetic linkage map using ddRAD sequencing technology to obtain SNP genotyping data of members from an full-sib family which sex had been determined. A total of 10,896 markers were grouped and oriented into 30 linkage groups, representing 30 chromosomes of A. davidianus. The genetic length of LGs ranged from 17.61 cM (LG30) to 280.81 cM (LG1), with a mean inter-locus distance ranging from 0.11(LG3) to 0.48 cM (LG26). The total genetic map length was 2643.10 cM with an average inter-locus distance of 0.24 cM. Three sex-related loci and four sex-related markers were found on LG6 and LG23, respectively. CONCLUSION: We constructed the first High-density genetic linkage map and identified three sex-related loci in the Chinese giant salamander. Current results are expected to be a useful tool for future genomic studies aiming at the marker-assisted breeding of the species.

BACKGROUND: The red panda (Ailurus fulgens) is a riddle of morphology, making it hard to tell whether it is an ursid, a procyonid, a mustelid, or a member of its own family. Previous genetic studies have given quite contradictory results as to its phylogenetic placement. RESULTS: A recently developed whole genome-based algorithm, the Whole Genome K-mer Signature algorithm was used to analyze the genomes of 28 species of Carnivora, including A. fulgens and several felid, ursid, mustelid, one mephitid species. This algorithm has the advantage of holistically using all the information in the genomes of these species. Being a genomics-based algorithm, it also reduces stochastic error to a minimum. Besides the whole genome, the mitochondrial DNA from 52 mustelids, mephitids, ursids, procyonids and A. fulgens were aligned to draw further phylogenetic inferences. The results from the whole genome study suggested that A. fulgens is a member of the mustelid clade (p = 9·10- 97). A. fulgens also separates from the mephitid Spilogala gracilis. The giant panda, Ailuropoda melanoleuca also clusters away from A. fulgens, together with other ursids (p = 1.2·10- 62). This could be due to the geographic isolation of A. fulgens from other mustelid species. However, results from the mitochondrial study as well as neighbor-joining methods based on the sequence identity matrix suggests that A. fulgens forms a monophyletic group. A Maximum Likelihood tree suggests that A. fulgens and Ursidae form a monophyletic group, although the bootstrap value is weak. CONCLUSIONS: The main conclusion that we can draw from this study is that on a whole genome level A. fulgens possibly belongs to the mustelid clade, and not an ursid or a mephitid. This despite the fact that previously some researchers classified A. fulgens and A. melanoleuca as relatives. Since the genotype determines the phenotype, molecular-based classification takes precedence over morphological classifications. This affirms the results of some previous studies, which studied smaller portions of the genome. However, mitochondrial analyses based on neighbor-joining and maximum likelihood methods suggest otherwise.