Ensembl Genomes 2016: more genomes, more complexity

Allelic Diversity and Development of Breeder-Friendly Marker Specific to floury2 Gene Regulating the Accumulation of α-Zeins and Essential Amino Acids in Maize Kernel

Maize zeins lack essential amino acids, such as methionine, lysine, and tryptophan. The floury2 (fl2) mutation reduces zein synthesis and increases methionine and lysine content in kernels. In this study, fl2 gene (1612 bp) was sequenced in eight wild-type and two mutant inbreds and detected 218 SNPs and 18 InDels. Transversion of C to T at 343 bp position caused the substitution of alanine by valine in the fl2 mutant. A PCR-based marker (FL-SNP-CT) was developed, which distinguished the favorable mutant fl2 allele (T) from the wild-type (C) Fl2 allele. Gene-based diversity analysis using seven gene-based InDel markers grouped 48 inbred lines into three major clusters, with an average genetic dissimilarity coefficient of 0.534. The average major allele frequency, gene diversity, heterozygosity, and polymorphism information content of the InDel markers were 0.701, 0.392, 0.039, and 0.318, respectively. Haplotype analysis revealed 29 haplotypes of fl2 gene among these 48 inbreds. Amino acid substitution (Ala-Val) at the signal peptide cleavage site produced unprocessed 24-kDa mutant protein instead of 22-kDa zein found in normal genotype. Eight paralogues of fl2 detected in the study showed variation in exon lengths (616-1170 bp) and translation lengths (135-267 amino acids). Orthologue analysis among 15 accessions of Sorghum bicolor and two accessions of Saccharum spontaneum revealed a single exon in fl2 gene, ranging from 267 to 810 bp. The study elucidated the molecular basis of fl2 mutation and reported a breeder-friendly marker for molecular breeding programs. This is the first study to characterize fl2 gene in a set of subtropically adapted inbreds.

Genome-wide characterization and expression profiling of FARL (FHY3/FAR1) family genes in Zea mays

A significant role of the plant is played by the transcription factor FARL, which is light signal transduction as well as plant growth and development. Despite being transposases, FARL has developed a variety of dominant biological actions in evolution and speciation. On the other hand, little is known about the Zea mays FARL protein family. This study identifies and characterizes fifteen ZmFARL genes genome-wide, and RNA sequencing data was used to profile their expression. 105 FARL proteins from five plant species were classified into five groups based on sequence alignment and phylogeny. The ZmFARL genes' exon-intron and motif distribution were conserved based on their evolutionary group. The fifteen ZmFARL genes were distributed over seven of the ten Z. mays chromosomes, although no duplication was discovered. Cis-element analysis reveals that ZmFARL genes play a variety of activities, including tissue-specific, stress- and hormone-responsive expressions. Furthermore, the results of the RNA sequencing used to profile expression showed that the genes ZmFARL2 and ZmFARL5 were much more expressed than other genes in various tissues, particularly in leaf characteristics. The identification of likely genes involved in cellular activity in Z. mays and related species will be aided by the characterization of the FARL genes.

Molecular Diagnostics for Invasive Fungal Diseases: Current and Future Approaches

Invasive fungal diseases (IFDs) comprise a growing healthcare burden, especially given the expanding population of immunocompromised hosts. Early diagnosis of IFDs is required to optimise therapy with antifungals, especially in the setting of rising rates of antifungal resistance. Molecular techniques including nucleic acid amplification tests and whole genome sequencing have potential to offer utility in overcoming limitations with traditional phenotypic testing. However, standardisation of methodology and interpretations of these assays is an ongoing undertaking. The utility of targeted Aspergillus detection has been well-defined, with progress in investigations into the role of targeted assays for Candida, Pneumocystis, Cryptococcus, the Mucorales and endemic mycoses. Likewise, whilst broad-range polymerase chain reaction assays have been in use for some time, pathology stewardship and optimising diagnostic yield is a continuing exercise. As costs decrease, there is also now increased access and experience with whole genome sequencing, including metagenomic sequencing, which offers unparalleled resolution especially in the investigations of potential outbreaks. However, their role in routine diagnostic use remains uncommon and standardisation of techniques and workflow are required for wider implementation.

MarR family proteins sense sulfane sulfur in bacteria

Members of the multiple antibiotic resistance regulator (MarR) protein family are ubiquitous in bacteria and play critical roles in regulating cellular metabolism and antibiotic resistance. MarR family proteins function as repressors, and their interactions with modulators induce the expression of controlled genes. The previously characterized modulators are insufficient to explain the activities of certain MarR family proteins. However, recently, several MarR family proteins have been reported to sense sulfane sulfur, including zero-valent sulfur, persulfide (R-SSH), and polysulfide (R-SnH, n ≥ 2). Sulfane sulfur is a common cellular component in bacteria whose levels vary during bacterial growth. The changing levels of sulfane sulfur affect the expression of many MarR-controlled genes. Sulfane sulfur reacts with the cysteine thiols of MarR family proteins, causing the formation of protein thiol persulfide, disulfide bonds, and other modifications. Several MarR family proteins that respond to reactive oxygen species (ROS) also sense sulfane sulfur, as both sulfane sulfur and ROS induce the formation of disulfide bonds. This review focused on MarR family proteins that sense sulfane sulfur. However, the sensing mechanisms reviewed here may also apply to other proteins that detect sulfane sulfur, which is emerging as a modulator of gene regulation.

Real-time emulation of future global warming reveals realistic impacts on the phenological response and quality deterioration in rice

Decreased production of crops due to climate change has been predicted scientifically. While climate-resilient crops are necessary to ensure food security and support sustainable agriculture, predicting crop growth under future global warming is challenging. Therefore, we aimed to assess the impact of realistic global warming conditions on rice cultivation. We developed a crop evaluation platform, the agro-environment (AE) emulator, which generates diverse environments by implementing the complexity of natural environmental fluctuations in customized, fully artificial lighting growth chambers. We confirmed that the environmental responsiveness of rice obtained in the fluctuation of artificial environments is similar to those exhibited in natural environments by validating our AE emulator using publicly available meteorological data from multiple years at the same location and multiple locations in the same year. Based on the representative concentration pathway, real-time emulation of severe global warming unveiled dramatic advances in the rice life cycle, accompanied by a 35% decrease in grain yield and an 85% increase in quality deterioration, which is higher than the recently reported projections. The transcriptome dynamism showed that increasing temperature and CO2 concentrations synergistically changed the expression of various genes and strengthened the induction of flowering, heat stress adaptation, and CO2 response genes. The predicted severe global warming greatly alters rice environmental adaptability and negatively impacts rice production. Our findings offer innovative applications of artificial environments and insights for enhancing varietal potential and cultivation methods in the future.

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

BACKGROUND

Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors.

METHODS

We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers.

RESULTS

We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors.

CONCLUSIONS

This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Assessing sequence variation, haplotype analysis and molecular characterisation of aspartate kinase2 (ask2) gene regulating methionine biosynthesis in diverse maize inbreds

Traditional maize grain is deficient in methionine, an essential amino acid required for proper growth and development in humans and poultry birds. Thus, development of high methionine maize (HMM) assumes great significance in alleviating malnutrition through sustainable and cost-effective approach. Of various genetic loci, aspartate kinase2 (ask2) gene plays a pivotal role in regulating methionine accumulation in maize. Here, we sequenced the entire ask2 gene of 5394 bp with 13 exons in five wild and five mutant maize inbreds to understand variation at nucleotide level. Sequence analysis revealed that an SNP in exon-13 caused thymine to adenine transversion giving rise to a favourable mutant allele associated with leucine to glutamine substitution in mutant ASK2 protein. Gene-based diversity analysis with 11 InDel markers grouped 48 diverse inbreds into three major clusters with an average genetic dissimilarity of 0.570 (range, 0.0-0.9). The average major allele frequency, gene diversity and PIC are 0.693, 0.408 and 0.341, respectively. A total of 45 haplotypes of the ask2 gene were identified among the maize inbreds. Evolutionary relationship analysis performed among 22 orthologues grouped them into five major clusters. The number of exons varied from 7 to 17, with length varying from 12 to 495 bp among orthologues. ASK2 protein with 565 amino acids was predicted to be in homo-dimeric state with lysine and tartaric acid as binding ligands. Amino acid kinase and ACT domains were found to be conserved in maize and orthologues. The study depicted the presence of enough genetic diversity in ask2 gene in maize, and development of HMM can be accelerated through introgression of favourable allele of ask2 into the parental lines of elite hybrids using molecular breeding.

Dynamics of mRNA fate during light stress and recovery: from transcription to stability and translation

Transcript stability is an important determinant of its abundance and, consequently, translational output. Transcript destabilisation can be rapid and is well suited for modulating the cellular response. However, it is unclear the extent to which RNA stability is altered under changing environmental conditions in plants. We previously hypothesised that recovery-induced transcript destabilisation facilitated a phenomenon of rapid recovery gene downregulation (RRGD) in Arabidopsis thaliana (Arabidopsis) following light stress, based on mathematical calculations to account for ongoing transcription. Here, we test this hypothesis and investigate processes regulating transcript abundance and fate by quantifying changes in transcription, stability and translation before, during and after light stress. We adapt syringe infiltration to apply a transcriptional inhibitor to soil-grown plants in combination with stress treatments. Compared with measurements in juvenile plants and cell culture, we find reduced stability across a range of transcripts encoding proteins involved in RNA binding and processing. We also observe light-induced destabilisation of transcripts, followed by their stabilisation during recovery. We propose that this destabilisation facilitates RRGD, possibly in combination with transcriptional shut-off that was confirmed for HSP101, ROF1 and GOLS1. We also show that translation remains highly dynamic over the course of light stress and recovery, with a bias towards transcript-specific increases in ribosome association, independent of changes in total transcript abundance, after 30 min of light stress. Taken together, we provide evidence for the combinatorial regulation of transcription and stability that occurs to coordinate translation during light stress and recovery in Arabidopsis.

The nucleolar protein NOL12 is required for processing of large ribosomal subunit rRNA precursors in Arabidopsis

BACKGROUND

NOL12 5'-3' exoribonucleases, conserved among eukaryotes, play important roles in pre-rRNA processing, ribosome assembly and export. The most well-described yeast counterpart, Rrp17, is required for maturation of 5.8 and 25S rRNAs, whereas human hNOL12 is crucial for the separation of the large (LSU) and small (SSU) ribosome subunit rRNA precursors.

RESULTS

In this study we demonstrate that plant AtNOL12 is also involved in rRNA biogenesis, specifically in the processing of the LSU rRNA precursor, 27S pre-rRNA. Importantly, the absence of AtNOL12 alters the expression of many ribosomal protein and ribosome biogenesis genes. These changes could potentially exacerbate rRNA biogenesis defects, or, conversely, they might stem from the disturbed ribosome assembly caused by delayed pre-rRNA processing. Moreover, exposure of the nol12 mutant to stress factors, including heat and pathogen Pseudomonas syringae, enhances the observed molecular phenotypes, linking pre-rRNA processing to stress response pathways. The aberrant rRNA processing, dependent on AtNOL12, could impact ribosome function, as suggested by improved mutant resistance to ribosome-targeting antibiotics.

CONCLUSION

Despite extensive studies, the pre-rRNA processing pathway in plants remains insufficiently characterized. Our investigation reveals the involvement of AtNOL12 in the maturation of rRNA precursors, correlating this process to stress response in Arabidopsis. These findings contribute to a more comprehensive understanding of plant ribosome biogenesis.

Novel genomic regions on chromosome 5B controlling wheat powdery mildew seedling resistance under Egyptian conditions

Wheat powdery mildew (PM) causes significant yield losses worldwide. None of the Egyptian wheat cultivars was detected to be highly resistant to such a severe disease. Therefore, a diverse spring wheat panel was evaluated for PM seedling resistance using different Bgt conidiospores collected from Egyptian fields in two growing seasons. The evaluation was done in two separate experiments. Highly significant differences were found between the two experiments suggesting the presence of different isolates populations. Highly significant differences were found among the tested genotypes confirming the ability to improve PM resistance using the recent panel. Genome-wide association study (GWAS) was done for each experiment separately and a total of 71 significant markers located within 36 gene models were identified. The majority of these markers are located on chromosome 5B. Haplotype block analysis identified seven blocks containing the significant markers on chromosome 5B. Five gene models were identified on the short arm of the chromosome. Gene enrichment analysis identified five and seven pathways based on the biological process and molecular functions respectively for the detected gene models. All these pathways are associated with disease resistance in wheat. The genomic regions on 5B seem to be novel regions that are associated with PM resistance under Egyptian conditions. Selection of superior genotypes was done and Grecian genotypes seem to be a good source for improving PM resistance under Egyptian conditions.

A Protocol for the Automatic Construction of Highly Curated Genome-Scale Models of Human Metabolism

Genome-scale metabolic models (GEMs) have emerged as a tool to understand human metabolism from a holistic perspective with high relevance in the study of many diseases and in the metabolic engineering of human cell lines. GEM building relies on either automated processes that lack manual refinement and result in inaccurate models or manual curation, which is a time-consuming process that limits the continuous update of reliable GEMs. Here, we present a novel algorithm-aided protocol that overcomes these limitations and facilitates the continuous updating of highly curated GEMs. The algorithm enables the automatic curation and/or expansion of existing GEMs or generates a highly curated metabolic network based on current information retrieved from multiple databases in real time. This tool was applied to the latest reconstruction of human metabolism (Human1), generating a series of the human GEMs that improve and expand the reference model and generating the most extensive and comprehensive general reconstruction of human metabolism to date. The tool presented here goes beyond the current state of the art and paves the way for the automatic reconstruction of a highly curated, up-to-date GEM with high potential in computational biology as well as in multiple fields of biological science where metabolism is relevant.

Genome-wide identification and comparative analyses of key genes involved in C4 photosynthesis in five main gramineous crops

Compared to C₃ species, C₄ plants showed higher photosynthetic capacity as well as water and nitrogen use efficiency due to the presence of the C₄ photosynthetic pathway. Previous studies have shown that all genes required for the C₄ photosynthetic pathway exist in the genomes of C₃ species and are expressed. In this study, the genes encoding six key C₄ photosynthetic pathway enzymes (β-CA, PEPC, ME, MDH, RbcS, and PPDK) in the genomes of five important gramineous crops (C₄: maize, foxtail millet, and sorghum; C₃: rice and wheat) were systematically identified and compared. Based on sequence characteristics and evolutionary relationships, their C₄ functional gene copies were distinguished from non-photosynthetic functional gene copies. Furthermore, multiple sequence alignment revealed important sites affecting the activities of PEPC and RbcS between the C₃ and C₄ species. Comparisons of expression characteristics confirmed that the expression patterns of non-photosynthetic gene copies were relatively conserved among species, while C₄ gene copies in C₄ species acquired new tissue expression patterns during evolution. Additionally, multiple sequence features that may affect C₄ gene expression and subcellular localization were found in the coding and promoter regions. Our work emphasized the diversity of the evolution of different genes in the C₄ photosynthetic pathway and confirmed that the specific high expression in the leaf and appropriate intracellular distribution were the keys to the evolution of C₄ photosynthesis. The results of this study will help determine the evolutionary mechanism of the C₄ photosynthetic pathway in Gramineae and provide references for the transformation of C₄ photosynthetic pathways in wheat, rice, and other major C₃ cereal crops.

Na+/Cl- cotransporter 2 is not fish-specific and is widely found in amphibians, non-avian reptiles, and select mammals

Solute carrier 12 (Slc12) is a family of electroneutral cation-coupled chloride (Cl-) cotransporters. Na+/K+/2Cl- (Nkcc) and Na+/Cl- cotransporters (Ncc) belong to the Nkcc/Ncc subfamily. Human and mouse possess one gene for the Na+/Cl- cotransporter (ncc gene: slc12a3), whereas teleost fishes possess multiple ncc genes, slc12a3 (ncc1) and slc12a10 (ncc2), in addition to their species-specific paralogs. Amphibians and squamates have two ncc genes: slc12a3 (ncc1) and ncc3. However, the evolutionary relationship between slc12a10 and ncc3 remains unresolved, and the presence of slc12a10 (ncc2) in mammals has not been clarified. Synteny and phylogenetic analyses of vertebrate genome databases showed that ncc3 is the ortholog of slc12a10, and slc12a10 is present in most ray-finned fishes, coelacanths, amphibians, reptiles, and a few mammals (e.g., platypus and horse) but pseudogenized or deleted in birds, most mammals, and some ray-finned fishes (pufferfishes). This shows that slc12a10 is widely present among bony vertebrates and pseudogenized or deleted independently in multiple lineages. Notably, as compared with some fish that show varied slc12a10 tissue expression profile, spotted gar, African clawed frog, red-eared slider turtle, and horse express slc12a10 in the ovaries or premature gonads. In horse tissues, an unexpectedly large number of splicing variants for Slc12a10 have been cloned, many of which encode truncated forms of Slc12a10, suggesting that the functional constraints of horse slc12a10 are weakened, which may be in the process of becoming a pseudogene. Our results elaborate on the evolution of Nkcc/Ncc subfamily of Slc12 in vertebrates.NEW & NOTEWORTHY slc12a10 is not a fish-specific gene and is present in a few mammals (e.g., platypus and horse), non-avian reptiles, amphibians, but was pseudogenized or deleted in most mammals (e.g., human, mouse, cat, cow, and rhinoceros), birds, and some ray-finned fishes (pufferfishes).

A chromosome-level genome assembly for Erianthus fulvus provides insights into its biofuel potential and facilitates breeding for improvement of sugarcane

Erianthus produces substantial biomass, exhibits a good Brix value, and shows wide environmental adaptability, making it a potential biofuel plant. In contrast to closely related sorghum and sugarcane, Erianthus can grow in degraded soils, thus releasing pressure on agricultural lands used for biofuel production. However, the lack of genomic resources for Erianthus hinders its genetic improvement, thus limiting its potential for biofuel production. In the present study, we generated a chromosome-scale reference genome for Erianthus fulvus Nees. The genome size estimated by flow cytometry was 937 Mb, and the assembled genome size was 902 Mb, covering 96.26% of the estimated genome size. A total of 35 065 protein-coding genes were predicted, and 67.89% of the genome was found to be repetitive. A recent whole-genome duplication occurred approximately 74.10 million years ago in the E. fulvus genome. Phylogenetic analysis showed that E. fulvus is evolutionarily closer to S. spontaneum and diverged after S. bicolor. Three of the 10 chromosomes of E. fulvus formed through rearrangements of ancestral chromosomes. Phylogenetic reconstruction of the Saccharum complex revealed a polyphyletic origin of the complex and a sister relationship of E. fulvus with Saccharum sp., excluding S. arundinaceum. On the basis of the four amino acid residues that provide substrate specificity, the E. fulvus SWEET proteins were classified as mono- and disaccharide sugar transporters. Ortho-QTL genes identified for 10 biofuel-related traits may aid in the rapid screening of E. fulvus populations to enhance breeding programs for improved biofuel production. The results of this study provide valuable insights for breeding programs aimed at improving biofuel production in E. fulvus and enhancing sugarcane introgression programs.

A receptor-like kinase controls the amplitude of secondary cell wall synthesis in rice

Cell wall expansion is a key element in determining plant morphology and growth, and cell wall integrity changes are relayed to the cell to fine-tune growth responses. Here, we show that variations in the ectodomain of a cell wall-associated receptor-like kinase, WAK10, in temperate Oryza japonica accessions differentially amplify fluctuations in cell wall integrity to control rice stem height. Mutation in the WAK10 gene exhibited increased cell wall thickening in stem sclerenchyma and reduced cell expansion in the stem. Two WAK10 ectodomain variants bound pectic oligosaccharides with different affinities. The pectic oligosaccharide binding regulated WAK10 phosphorylation activity, the amplitude of secondary wall deposition, and ultimately, stem height. Rice population analyses revealed active enrichment of the short-stem WAK10 ectodomain alleles in japonica subspecies during domestication. Our study outlines not only a mechanism for how variations in ligand affinities of a receptor kinase control cell wall biosynthesis and plant growth, but it also provides breeding targets for new semi-dwarf rice cultivars.

Hepatitis C and HIV detection by blood RNA-sequencing in cohort of smokers

Detection of viruses by RNA and DNA sequencing has improved the understanding of the human virome. We sought to identify blood viral signatures through secondary use of RNA-sequencing (RNA-seq) data in a large study cohort. The ability to reveal undiagnosed infections with public health implications among study subjects with available sequencing data could enable epidemiologic surveys and may lead to diagnosis and therapeutic interventions, leveraging existing research data in a clinical context. We detected viral RNA in peripheral blood RNA-seq data from a COPD-enriched population of current and former smokers. Correlation between viral detection and both reported infections and relevant disease outcomes was evaluated. We identified Hepatitis C virus RNA in 228 subjects and HIV RNA in 30 subjects. Overall, we observed 31 viral species, including Epstein-Barr virus and Cytomegalovirus. We observed an enrichment of Hepatitis C and HIV infections among subjects reporting liver disease and HIV infections, respectively. Higher interferon expression scores were observed in the subjects with Hepatitis C and HIV infections. Through secondary use of RNA-seq from a cohort of current and former smokers, we detected peripheral blood viral signatures. We identified HIV and Hepatitis C virus (HCV), highlighting potential public health implications for the approach described this study. We observed correlations with reported infections, chronic infection outcomes and the host transcriptomic response, providing evidence to support the validity of the approach.

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

PEATmoss: A Gene Expression Atlas for Bryophytes

PEATmoss is an interactive gene expression atlas for bryophytes, which originally unified Physcomitrium patens RNA-seq and microarray expression data from multiple gene annotation versions. This atlas includes more than 100 experiments of P. patens, is expanding to host Anthoceros agrestis and Marchantia polymorpha, and aims to host data from more species in the future. PEATmoss has multiple visualization methods and tools for data downloading and is connected to the Physcomitrium patens Gene Model Lookup DB (PpGML DB), which links P. patens genes to annotations and resources from several databases and contains tools for gene version lookup and sequence and annotation extraction. Among the new features available in PEATmoss are dataset privacy control, multispecies menu, interactive color scale, co-expression network visualization, and replicate data downloading.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

An overview of online resources for intra-species detection of gene duplications

Gene duplication plays an important role in evolutionary mechanism, which can act as a new source of genetic material in genome evolution. However, detecting duplicate genes from genomic data can be challenging. Various bioinformatics resources have been developed to identify duplicate genes from single and/or multiple species. Here, we summarize the metrics used to measure sequence identity among gene duplicates within species, compare several computational approaches that have been used to predict gene duplicates, and review recent advancements of a Basic Local Alignment Search Tool (BLAST)-based web tool and database, allowing future researchers to easily identify intra-species gene duplications. This article is a quick reference guide for research tools used for detecting gene duplicates.

RSC and GRFs confer promoter directionality by restricting divergent noncoding transcription

The directionality of gene promoters-the ratio of protein-coding over divergent noncoding transcription-is highly variable. How promoter directionality is controlled remains poorly understood. Here, we show that the chromatin remodelling complex RSC and general regulatory factors (GRFs) dictate promoter directionality by attenuating divergent transcription relative to protein-coding transcription. At gene promoters that are highly directional, depletion of RSC leads to a relative increase in divergent noncoding transcription and thus to a decrease in promoter directionality. We find that RSC has a modest effect on nucleosome positioning upstream in promoters at the sites of divergent transcription. These promoters are also enriched for the binding of GRFs such as Reb1 and Abf1. Ectopic targeting of divergent transcription initiation sites with GRFs or the dCas9 DNA-binding protein suppresses divergent transcription. Our data suggest that RSC and GRFs play a pervasive role in limiting divergent transcription relative to coding direction transcription. We propose that any DNA-binding factor, when stably associated with cryptic transcription start sites, forms a barrier which represses divergent transcription, thereby promoting promoter directionality.

No evidence for widespread positive selection on double substitutions within codons in primates and yeasts

Nucleotide substitutions in protein-coding genes can be divided into synonymous (S) and non-synonymous (N) ones that alter amino acids (including nonsense mutations causing stop codons). The S substitutions are expected to have little effect on function. The N substitutions almost always are affected by strong purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases can modulate the deleterious effect of single N substitutions and, thus, could be subjected to the positive selection. This effect has been demonstrated for mutations in the serine codons, stop codons and double N substitutions in prokaryotes. In all abovementioned cases, a novel technique was applied that allows elucidating the effects of selection on double substitutions considering mutational biases. Here, we applied the same technique to study double N substitutions in eukaryotic lineages of primates and yeast. We identified markedly fewer cases of purifying selection relative to prokaryotes and no evidence of codon double substitutions under positive selection. This is consistent with previous studies of serine codons in primates and yeast. In general, the obtained results strongly suggest that there are major differences between studied pro- and eukaryotes; double substitutions in primates and yeasts largely reflect mutational biases and are not hallmarks of selection. This is especially important in the context of detection of positive selection in codons because it has been suggested that multiple mutations in codons cause false inferences of lineage-specific site positive selection. It is likely that this concern is applicable to previously studied prokaryotes but not to primates and yeasts where markedly fewer double substitutions are affected by positive selection.

Characterization of the CqCAMTA gene family reveals the role of CqCAMTA03 in drought tolerance

BACKGROUND

Calmodulin-binding transcription activators (CAMTAs) are relatively conserved calmodulin-binding transcription factors widely found in eukaryotes and play important roles in plant growth and stress response. CAMTA transcription factors have been identified in several plant species, but the family members and functions have not yet been identified and analyzed in quinoa.

RESULTS

In this study, we identified seven CAMTA genes across the whole quinoa genome and analyzed the expression patterns of CqCAMTAs in root and leaf tissues. Gene structure, protein domain, and phylogenetic analyses showed that the quinoa CAMTAs were structurally similar and clustered into the same three major groups as other plant CAMTAs. A large number of stress response-related cis-elements existed in the 2 kb promoter region upstream of the transcription start site of the CqCAMTA genes. qRT-PCR indicated that CqCAMTA genes were expressed differentially under PEG treatments in leaves, and responded to drought stress in leaves and roots. In particular, the CqCAMTA03 gene strongly responded to drought. The transient expression of CqCAMTA03-GFP fusion protein in the tobacco leaf showed that CqCAMTA03 was localized in the nucleus. In addition, transgenic Arabidopsis lines exhibited higher concentration levels of the antioxidant enzymes measured, including POD, SOD, and CAT, under drought conditions with very low levels of H₂O₂ and MDA. Moreover, relative water content and the degree of stomatal opening showed that the transgenic Arabidopsis lines were more tolerant of both stress factors as compared to their wild types.

CONCLUSION

In this study, the structures and functions of the CAMTA family in quinoa were systematically explored. Many CAMTAs may play vital roles in the regulation of organ development, growth, and responses to drought stress. The results of the present study serve as a basis for future functional studies on the quinoa CAMTA family.

Comparison of SWEET gene family between maize and foxtail millet through genomic, transcriptomic, and proteomic analyses

Foxtail millet [Setaria italica (L.) P. Beauv.] does not show high yield and biomass compared with maize (Zea mays L.) although it is a C₄ crop with the potential for high productivity. Because SWEET genes, which are important for sugar transport in plants, play critical roles in biomass production and seed filling in crops, genome-wide, transcriptomic, and proteomic comparison on SWEET gene family between these two species would provide some clues for unlocking this issue. In our study, 24 SWEET genes were identified in foxtail millet and maize. Sequence-based bioinformatics combined with gene expression analyses identified several candidate functional orthologs in these two species. A comparative analysis on expression characteristics of SWEET genes and proteins between maize and foxtail millet indicate that not only some critical major SWEET proteins show significant upregulation in maize compared with their orthologs in foxtail millet, but also there are more quantities of maize SWEET genes showing high expressions than that of foxtail millet genes, suggesting that compared with foxtail millet, maize possesses higher capacity of sugar transport, the crucial determinant for crop yield and biomass. These results provide a basis on revealing why foxtail millet exhibits low yield and biomass although it is a C₄ crop with the potential for high productivity.

24-nt phasiRNAs move from tapetal to meiotic cells in maize anthers

In maize, 24-nt phased, secondary small interfering RNAs (phasiRNAs) are abundant in meiotic stage anthers, but their distribution and functions are not precisely known. Using laser capture microdissection, we analyzed tapetal cells, meiocytes and other somatic cells at several stages of anther development to establish the timing of 24-PHAS precursor transcripts and the 24-nt phasiRNA products. By integrating RNA and small RNA profiling plus single-molecule and small RNA FISH (smFISH or sRNA-FISH) spatial detection, we demonstrate that the tapetum is the primary site of 24-PHAS precursor and Dcl5 transcripts and the resulting 24-nt phasiRNAs. Interestingly, 24-nt phasiRNAs accumulate in all cell types, with the highest levels in meiocytes, followed by tapetum. Our data support the conclusion that 24-nt phasiRNAs are mobile from tapetum to meiocytes and to other somatic cells. We discuss possible roles for 24-nt phasiRNAs in anther cell types.

Comprehensive analysis of the expression, prognostic significance, and function of FAM83 family members in breast cancer

BACKGROUND

The FAM83 family plays a key role in tumorigenesis and cancer progression. However, the role of the FAM83 family in the development of breast tumors is unclear to date. This report explores the expression, prognostic significance, and function of the FAM83 family members in breast cancer using public databases.

METHODS

UALCAN database was used to explore the expression of FAM83 family members in breast cancer. Furthermore, we validated the expression of FAM83 family members in twenty pairs of breast cancer and normal tissues by RT-PCR. Kaplan-Meier plotter database was used to explore the prognostic significance of FAM83 family members in breast cancer. GeneMANIA and DAVID databases were used for functional and pathway enrichment analysis of genes co-expressed with FAM83A, FAM83D, FAM83F, and FAM83G. MEXPRESS and UALCAN databases were used to analyze the level of DNA promoter methylation of FAM83A, FAM83D, FAM83F, and FAM83G in breast cancer. TIMER database was utilized to explore the relationships between immune cell infiltration and FAM83A, FAM83D, FAM83F, and FAM83G expression.

RESULTS

Among FAM83 family members, FAM83A, FAM83D, FAM83F, and FAM83G were higher expressed in breast cancer than in normal tissues. We also validated the significant high expression of FAM83A, FAM83D, FAM83F, and FAM83G mRNA in breast cancer than in normal samples. Their increased expression has an adverse prognostic effect on breast cancer patients. These genes co-expressed with FAM83A, FAM83D, FAM83F, and FAM83G might take part in cell proliferation, G2/M transition of the mitotic cell cycle, regulation of apoptosis process and other cancer-related biological processes. In addition, they were mainly enriched in the Hippo signaling pathway, Hedgehog signaling pathway, PI3K/AKT signaling pathway, and other cancer-related pathways. We also found that promoter DNA methylation might regulate the expression of FAM83A, FAM83D, FAM83F, and FAM83G mRNA in most CpG islands. At last, we found the expression of FAM83A, FAM83D, FAM83F, and FAM83G mRNA was significantly related to immune cell infiltration.

CONCLUSIONS

FAM83A, FAM83D, FAM83F, and FAM83G were highly expressed in breast cancer tissues and had an adverse effect on the survival outcomes of breast cancer patients. Also, they were involved in breast cancer-related signal pathways. Therefore, they might serve as potential therapeutic targets for breast cancer clinical treatment.

Allelic Variation in Zmfatb Gene Defines Variability for Fatty Acids Composition Among Diverse Maize Genotypes

Edible oil with lower saturated fatty acids is desired for perceived quality and health benefits to humans and livestock. fatb gene encoding acyl-ACP thioesterase is a key player in the conversion of palmitic acid to oleic acid, thereby modifying the ratio of saturated to unsaturated fatty acids in maize kernels. The present investigation characterised the full-length sequence of the Zmfatb gene (4.63 kb) in two mutants (Zmfatb) and eight wild-types (ZmfatB) inbreds to study allelic variation, gene-based diversity, phylogenetic-relationship, protein-modelling, and molecular-docking to identify novel candidates for modification of fatty acid profile. Sequence alignment revealed wide genomic variability for Zmfatb among the inbreds; identified five novel SNPs and two InDels that clearly differentiated the wild-type and mutant genotypes. Gene-based diversity using 11-InDel markers categorised 48-diverse maize-inbreds into two-clusters. The majority of mutant and wild-type inbreds were grouped in separate clusters and led to the generation of 41 haplotypes. Genetic relationship of maize fatb gene with orthologues among 40 accessions of 12 oilseed-crops using both nucleotide and protein sequence clustered maize, soybean, sunflower, opium-poppy, Citrulus lanata, quinoa, and prunus species into one cluster; and brassica, camelina, and arabidopsis into the different cluster. The clustering pattern revealed that the plant oil with higher unsaturated fatty acids, particularly oleic, linoleic, and linolenic acids grouped together in one cluster and higher proportions of other fractions like arachidic, eicosenoic, and erucic acids grouped in another cluster. Physico-chemical properties highlighted more similarity between maize and 29 orthologue proteins, but orthologues were found to have better thermostability. Homology models have been developed for maize mutant and wild-type inbreds using Umbellularia californica (PDB ID: 5x04) as a template. Predicted protein models possessed optimum confidence-score and RMSD values and validated stability via., Ramachandran plots. Molecular docking indicated most of the interactions of protein-ligand were having similar binding-affinity due to the broader specificity of fatty acyl-ACP thioesterases and the presence of conserved-domains across crops. This is the first report on the comprehensive molecular characterisation of the fatb gene in maize and various orthologues. The information generated here provided new insights into the genetic diversity of fatb gene which can be utilised for the enhanced nutritive value of oil in the breeding programme.

Flavonoids from Scutellaria barbata D. Don exert antitumor activity in colorectal cancer through inhibited autophagy and promoted apoptosis via ATF4/sestrin2 pathway

PURPOSE

Scutellaria barbata D. Don (SB), mainly containing flavonoids, has been frequently used for cancer treatment. However, little research has investigated the antitumor activity of flavonoids from SB (FSB). The current study aimed to assess the antitumor effect of TFSB and elucidate the probable underlying mechanism in vivo and in vitro.

STUDY DESIGN

FSB was prepared, and its chemical composition was characterized by HPLC-MS. Colorectal HCT116 cells were treated with various concentration of FSB. The viability, proliferation, apoptosis, migration, and autophagy of HCT116 cells were studied, as were further confirmed in tumor xenografts.

METHODS

Cell viability and proliferation were respectively examined by MTT and EdU staining. ROS was determined with DCFH-DA, and cell apoptosis was detected using flow cytometry. Transwell and wound-healing assays were performed to evaluate cell migration. Immunofluorescence was employed to evaluate sestrin2 and ATF4 level. The protein expressions of p-AMPK, p-ULK1, p-mTOR, 4E-BP1, LC3-I/II, cleaved-caspase-3, Bax, and bcl-2 were investigated by western blot. ATF4 was overexpressed in experiments to explore the role of ATF4/sestrin2 pathway in FSB-mediated efficacy.

RESULTS

FSB clearly reduced the cell viability, promoted ROS generation, and induced apoptosis in HCT116 cells by down-regulated Bcl-2, and increased cleaved-caspase-3 and Bax. Furthermore, FSB significantly inhibited migration of colorectal cells in a dose-dependent manner. Further mechanistic study indicated that FSB upregulated p-mTOR protein level, and reduced p-AMPK, p-ULK1, p-mTOR, p-4E-BP1 and LC3-I/II expression, which were major autophagy-related genes. In addition, FSB could cause downregulation of endogenous mTOR inhibitor sestrin2 and ATF4 expression. Transient overexpression of ATF4 resulted in mTOR and sestrin2 inhibition, and significantly compromised the effects of FSB on apoptosis and autophagy in HCT116 cells.

CONCLUSION

Our results reveal, for the first time, that FSB exerts antitumor activity through autophagy inhibition and apoptosis induction via ATF4/sestrin2 pathway in colorectal cancer cells. Scutellaria barbata D. Don may have great potential in the application for the prevention and treatment of human colorectal cancer.

Genomics and Informatics, Conjoined Tools Vital for Understanding and Protecting Plant Health

Genomics' impact on crop production continuously expands. The number of sequenced plant and microbial species and strains representing diverse populations of individual species rapidly increases thanks to the advent of next-generation sequencing technologies. Their genomic blueprints revealed candidate genes involved in various functions and processes crucial for crop health and helped in understanding how the sequenced organisms have evolved at the genome level. Functional genomics quickly translates these blueprints into a detailed mechanistic understanding of how such functions and processes work and are regulated; this understanding guides and empowers efforts to protect crops from diverse biotic and abiotic threats. Metagenome analyses help identify candidate microbes crucial for crop health and uncover how microbial communities associated with crop production respond to environmental conditions and cultural practices, presenting opportunities to enhance crop health by judiciously configuring microbial communities. Efficient conversion of disparate types of massive genomics data into actionable knowledge requires a robust informatics infrastructure supporting data preservation, analysis, and sharing. This review starts with an overview of how genomics came about and has quickly transformed life science. We illuminate how genomics and informatics can be applied to investigate various crop health-related problems using selected studies. We end the review by noting why community empowerment via crowdsourcing is crucial to harnessing genomics to protect global food and nutrition security without continuously expanding the environmental footprint of crop production.

A cascade of bHLH-regulated pathways programs maize anther development

The spatiotemporal development of somatic tissues of the anther lobe is necessary for successful fertile pollen production. This process is mediated by many transcription factors acting through complex, multi-layered networks. Here, our analysis of functional knockout mutants of interacting basic helix-loop-helix genes Ms23, Ms32, basic helix-loop-helix 122 (bHLH122), and bHLH51 in maize (Zea mays) established that male fertility requires all four genes, expressed sequentially in the tapetum (TP). Not only do they regulate each other, but also they encode proteins that form heterodimers that act collaboratively to guide many cellular processes at specific developmental stages. MS23 is confirmed to be the master factor, as the ms23 mutant showed the earliest developmental defect, cytologically visible in the TP, with the most drastic alterations in premeiotic gene expression observed in ms23 anthers. Notably, the male-sterile ms23, ms32, and bhlh122-1 mutants lack 24-nt phased secondary small interfering RNAs (phasiRNAs) and the precursor transcripts from the corresponding 24-PHAS loci, while the bhlh51-1 mutant has wild-type levels of both precursors and small RNA products. Multiple lines of evidence suggest that 24-nt phasiRNA biogenesis primarily occurs downstream of MS23 and MS32, both of which directly activate Dcl5 and are required for most 24-PHAS transcription, with bHLH122 playing a distinct role in 24-PHAS transcription.

Isolation and Molecular Characterisation of TtDro1A and TtDro1B Genes from Triticum turgidum Subspecies durum and turgidum, Study of Their Influences on Seedling Root Angles

Durum wheat (Triticum turgidum, 2n = 4x = AABB) includes several subspecies with differential characteristics in their root system architecture (RSA). Subspecies durum has longer and more vertical roots, while subspecies turgidum has smaller and shallower roots. The homeologous genes TtDro1A and TtDro1B of both subspecies have been identified and found to differ in their sizes, sequences and the proteins they encode. To determine whether there is a relationship between the level of expression of these two genes and the angle adopted by the roots of durum wheat seedlings, their expressions has been studied by RT-qPCR, both in the primary seminal root and in the other seminal roots. The results of the analyses showed that the TtDro1A gene is expressed 1.4 times more in the primary seminal root than in the other seminal roots. Furthermore, this gene is expressed 2.49 to 8.76 times more than TtDro1B depending on root type (primary or seminal) and subspecies. There are positive correlations between the expression ratio of both genes (TtDro1A/TtDro1B) and the mean of all root angles, the most vertical root angle and the most horizontal root angle of the seedlings. The higher the expression of TtDro1B gene, the lower the root growth angles.

Reduced reproductive success is associated with selective constraint on human genes

Genome-wide sequencing of human populations has revealed substantial variation among genes in the intensity of purifying selection acting on damaging genetic variants¹. Although genes under the strongest selective constraint are highly enriched for associations with Mendelian disorders, most of these genes are not associated with disease and therefore the nature of the selection acting on them is not known². Here we show that genetic variants that damage these genes are associated with markedly reduced reproductive success, primarily owing to increased childlessness, with a stronger effect in males than in females. We present evidence that increased childlessness is probably mediated by genetically associated cognitive and behavioural traits, which may mean that male carriers are less likely to find reproductive partners. This reduction in reproductive success may account for 20% of purifying selection against heterozygous variants that ablate protein-coding genes. Although this genetic association may only account for a very minor fraction of the overall likelihood of being childless (less than 1%), especially when compared to more influential sociodemographic factors, it may influence how genes evolve over time.

Identification and characterization of stress responsive homeodomain leucine zipper transcription factors in Medicago truncatula

BACKGROUND

Homeodomain leucine zipper (HD-ZIP) transcription factors play roles in regulating plant development and responses to abiotic stresses; however, how HD-ZIP genes in Medicago truncatula are involved in abiotic stress response remains elusive.

METHODS AND RESULTS

The HD-ZIP I genes in Medicago truncatula were identified and characterized, and their expression patterns in different tissues and under different abiotic stresses were analyzed. A total of 15 Medicago truncatula HD-ZIP I genes were identified and a phylogenetic analysis of HD-ZIP I proteins in Arabidopsis thaliana and Medicago truncatula was conducted. Fifteen HD-ZIP I genes showed diverse tissue preferences. Among them, expressions of MtHB22 and MtHB51 were specially detected in vegetative buds. In addition, they responded to various abiotic stresses, including salinity and osmotic stress and abscisic acid (ABA). For instance, MtHB7 and MtHB12 expression levels were found to be positively associated with salt, osmotic stress and ABA in both shoots and roots, while MtHB13 and MtHB23 were negatively associated with these stresses in Medicago truncatula.

CONCLUSION

The HD-ZIP I genes in Medicago truncatula are evolutionarily conserved, but also exhibit gene duplication and gene loss events. Differential expression analysis of Medicago truncatula HD-ZIP I genes indicated their crucial roles in abiotic stress responses. Our genome-wide analysis of the HD-ZIP I transcription factor family in Medicago truncatula provided a valuable reference for further research.

Analysis of mRNA-derived siRNAs in mutants of mRNA maturation and surveillance pathways in Arabidopsis thaliana

Defects in RNA maturation and RNA decay factors may generate substrates for the RNA interference machinery. This phenomenon was observed in plants where mutations in some RNA-related factors lead to the production of RNA-quality control small interfering RNAs and several mutants show enhanced silencing of reporter transgenes. To assess the potential of RNAi activation on endogenous transcripts, we sequenced small RNAs from a set of Arabidopsis thaliana mutants with defects in various RNA metabolism pathways. We observed a global production of siRNAs caused by inefficient pre-mRNA cleavage and polyadenylation leading to read-through transcription into downstream antisense genes. In addition, in the lsm1a lsm1b double mutant, we identified NIA1, SMXL5, and several miRNA-targeted mRNAs as producing siRNAs, a group of transcripts suggested being especially sensitive to deficiencies in RNA metabolism. However, in most cases, RNA metabolism perturbations do not lead to the widespread production of siRNA derived from mRNA molecules. This observation is contrary to multiple studies based on reporter transgenes and suggests that only a very high accumulation of defective mRNA species caused by specific mutations or substantial RNA processing defects trigger RNAi pathways.

Chronic calcium signaling in IgE+ B cells limits plasma cell differentiation and survival

In contrast to other antibody isotypes, B cells switched to IgE respond transiently and do not give rise to long-lived plasma cells (PCs) or memory B cells. To better understand IgE-BCR-mediated control of IgE responses, we developed whole-genome CRISPR screening that enabled comparison of IgE⁺ and IgG1⁺ B cell requirements for proliferation, survival, and differentiation into PCs. IgE⁺ PCs exhibited dependency on the PI3K-mTOR axis that increased protein amounts of the transcription factor IRF4. In contrast, loss of components of the calcium-calcineurin-NFAT pathway promoted IgE⁺ PC differentiation. Mice bearing a B cell-specific deletion of calcineurin B1 exhibited increased production of IgE⁺ PCs. Mechanistically, sustained elevation of intracellular calcium in IgE⁺ PCs downstream of the IgE-BCR promoted BCL2L11-dependent apoptosis. Thus, chronic calcium signaling downstream of the IgE-BCR controls the self-limiting character of IgE responses and may be relevant to the accumulation of IgE-producing cells in allergic disease.

Arabidopsi s Spliceosome Factor SmD3 Modulates Immunity to Pseudomonas syringae Infection

SmD3 is a core component of the small nuclear ribonucleoprotein (snRNP) that is essential for pre-mRNA splicing. The role of Arabidopsis SmD3 in plant immunity was assessed by testing sensitivity of smd3a and smd3b mutants to Pseudomonas syringae pv. tomato (Pst) DC3000 infection and its pathogenesis effectors flagellin (flg22), EF-Tu (elf18) and coronatine (COR). Both smd3 mutants exhibited enhanced susceptibility to Pst accompanied by marked changes in the expression of key pathogenesis markers. mRNA levels of major biotic stress response factors were also altered upon treatment with Pseudomonas effectors. Our genome-wide transcriptome analysis of the smd3b-1 mutant infected with Pst, verified by northern and RT-qPCR, showed that lack of SmD3-b protein deregulates defense against Pst infection at the transcriptional and posttranscriptional levels including defects in splicing and an altered pattern of alternative splicing. Importantly, we show that SmD3-b dysfunction impairs mainly stomatal immunity as a result of defects in stomatal development. We propose that it is the malfunction of the stomata that is the primary cause of an altered mutant response to the pathogen. Other changes in the smd3b-1 mutant involved enhanced elf18- and flg22-induced callose deposition, reduction of flg22-triggered production of early ROS and boost of secondary ROS caused by Pst infection. Together, our data indicate that SmD3 contributes to the plant immune response possibly via regulation of mRNA splicing of key pathogenesis factors.

Identification of WRKY transcription factors responding to abiotic stresses in Brassica napus L

A total of 278 BnWRKYs were identified and analyzed. Ectopic expression of BnWRKY149 and BnWRKY217 suggests that they function in the ABA signaling pathway. WRKY transcription factors play an important role in plant development, however, their function in Brassica napus L. abiotic stress response is still unclear. In this study, a total of 278 BnWRKY transcription factors were identified from the B. napus genome data, and they were subsequently distributed in three main groups. The protein motifs and classification of BnWRKY transcription factors were analyzed, and the locations of their corresponding encoding genes were mapped on the chromosomes of B. napus. Transcriptome analysis of rapeseed seedlings exposed to drought, salt, heat, cold and abscisic acid treatment revealed that 99 BnWRKYs responded to at least one of these stresses. The expression profiles of 12 BnWRKYs were examined with qPCR and the result coincided with RNA-seq analysis. Two genes of interest, BnWRKY149 and BnWRKY217 (homologs of AtWRKY40), were overexpressed in Arabidopsis, and the corresponding proteins were located to the nucleus. Transgene plants of BnWRKY149 and BnWRKY217 were less sensitive to ABA than Arabidopsis Col-0 plants, suggesting they might play important roles in the responses of rapeseed to abiotic stress.

Genomicus in 2022: comparative tools for thousands of genomes and reconstructed ancestors

Genomicus is a database and web-server dedicated to comparative genomics in eukaryotes. Its main functionality is to graphically represent the conservation of genomic blocks between multiple genomes, locally around a specific gene of interest or genome-wide through karyotype comparisons. Since 2010 and its first release, Genomicus has synchronized with 60 Ensembl releases and seen the addition of functions that have expanded the type of analyses that users can perform. Today, five public instances of Genomicus are supporting a total number of 1029 extant genomes and 621 ancestral reconstructions from all eukaryotes kingdoms available in Ensembl and Ensembl Genomes databases complemented with four additional instances specific to taxonomic groups of interest. New visualization and query tools are described in this manuscript. Genomicus is freely available at http://www.genomicus.bio.ens.psl.eu/genomicus.

Comparative analysis of RNA enrichment methods for preparation of Cryptococcus neoformans RNA sequencing libraries

RNA sequencing (RNA-Seq) experiments focused on gene expression involve removal of ribosomal RNA (rRNA) because it is the major RNA constituent of cells. This process, called RNA enrichment, is done primarily to reduce cost: without rRNA removal, deeper sequencing must be performed to compensate for the sequencing reads wasted on rRNA. The ideal RNA enrichment method removes all rRNA without affecting other RNA in the sample. We tested the performance of three RNA enrichment methods on RNA isolated from Cryptococcus neoformans, a fungal pathogen of humans. We find that the RNase H depletion method is more efficient in depleting rRNA and more specific in recapitulating non-rRNA levels present in unenriched controls than the commonly-used Poly(A) isolation method. The RNase H depletion method is also more effective than the Ribo-Zero depletion method as measured by rRNA depletion efficiency and recapitulation of protein-coding RNA levels present in unenriched controls, while the Ribo-Zero depletion method more closely recapitulates annotated non-coding RNA (ncRNA) levels. Finally, we leverage these data to accurately map the C. neoformans mitochondrial rRNA genes, and also demonstrate that RNA-Seq data generated with the RNase H and Ribo-Zero depletion methods can be used to explore novel C. neoformans long non-coding RNA genes.

Mitochondrial phosphate transporter and methyltransferase genes contribute to Fusarium head blight Type II disease resistance and grain development in wheat

Fusarium head blight (FHB) is an economically important disease of wheat that results in yield loss and grain contaminated with fungal mycotoxins that are harmful to human and animal health. Herein we characterised two wheat genes involved in the FHB response in wheat: a wheat mitochondrial phosphate transporter (TaMPT) and a methyltransferase (TaSAM). Wheat has three sub-genomes (A, B, and D) and gene expression studies demonstrated that TaMPT and TaSAM homoeologs were differentially expressed in response to FHB infection and the mycotoxigenic Fusarium virulence factor deoxynivalenol (DON) in FHB resistant wheat cv. CM82036 and susceptible cv. Remus. Virus-induced gene silencing (VIGS) of either TaMPT or TaSAM enhanced the susceptibility of cv. CM82036 to FHB disease, reducing disease spread (Type II disease resistance). VIGS of TaMPT and TaSAM significantly reduced grain number and grain weight. This indicates TaSAM and TaMPT genes also contribute to grain development in wheat and adds to the increasing body of evidence linking FHB resistance genes to grain development. Hence, Fusarium responsive genes TaSAM and TaMPT warrant further study to determine their potential to enhance both disease resistance and grain development in wheat.

LSTrAP-Kingdom: an automated pipeline to generate annotated gene expression atlases for kingdoms of life

MOTIVATION

There are now more than two million RNA sequencing experiments for plants, animals, bacteria and fungi publicly available, allowing us to study gene expression within and across species and kingdoms. However, the tools allowing the download, quality control and annotation of this data for more than one species at a time are currently missing.

RESULTS

To remedy this, we present the Large-Scale Transcriptomic Analysis Pipeline in Kingdom of Life (LSTrAP-Kingdom) pipeline, which we used to process 134 521 RNA-seq samples, achieving ∼12 000 processed samples per day. Our pipeline generated quality-controlled, annotated gene expression matrices that rival the manually curated gene expression data in identifying functionally related genes.

AVAILABILITY AND IMPLEMENTATION

LSTrAP-Kingdom is available from: https://github.com/wirriamm/plants-pipeline and is fully implemented in Python and Bash.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Nasal Immunization With Small Molecule Mast Cell Activators Enhance Immunity to Co-Administered Subunit Immunogens

Mast cell activators are a novel class of mucosal vaccine adjuvants. The polymeric compound, Compound 48/80 (C48/80), and cationic peptide, Mastoparan 7 (M7) are mast cell activators that provide adjuvant activity when administered by the nasal route. However, small molecule mast cell activators may be a more cost-efficient adjuvant alternative that is easily synthesized with high purity compared to M7 or C48/80. To identify novel mast cell activating compounds that could be evaluated for mucosal vaccine adjuvant activity, we employed high-throughput screening to assess over 55,000 small molecules for mast cell degranulation activity. Fifteen mast cell activating compounds were down-selected to five compounds based on in vitro immune activation activities including cytokine production and cellular cytotoxicity, synthesis feasibility, and selection for functional diversity. These small molecule mast cell activators were evaluated for in vivo adjuvant activity and induction of protective immunity against West Nile Virus infection in BALB/c mice when combined with West Nile Virus envelope domain III (EDIII) protein in a nasal vaccine. We found that three of the five mast cell activators, ST101036, ST048871, and R529877, evoked high levels of EDIII-specific antibody and conferred comparable levels of protection against WNV challenge. The level of protection provided by these small molecule mast cell activators was comparable to the protection evoked by M7 (67%) but markedly higher than the levels seen with mice immunized with EDIII alone (no adjuvant 33%). Thus, novel small molecule mast cell activators identified by high throughput screening are as efficacious as previously described mast cell activators when used as nasal vaccine adjuvants and represent next-generation mast cell activators for evaluation in mucosal vaccine studies.

Selection of transcripts related to low-temperature tolerance using RNA sequencing from F2 plants between japonica and indica rice (Oryza sativa L.) cultivars

At low temperatures (18°C), seedlings of an indica rice (Oryza sativa L.) cultivar Kasalath showed symptoms of chlorosis, although the leaves of a japonica cultivar Arroz da Terra remained green. In this study, transcripts related to the chlorophyll content of rice seedlings grown at 18°C were investigated using RNA-sequencing (RNA-Seq) data for F2 crosses between cultivars Arroz da Terra and Kasalath, as well as their parental cultivars. Differential expression analysis revealed that gene ontology terms related to 'photosynthesis' were significantly enriched in lowly expressed genes at 18°C than at 25°C in Kasalath. However, the gene ontology terms related to 'response to stress' were significantly enriched in highly expressed genes at 18°C than at 25°C in Kasalath. When the F2 plants were grown at 18°C, their chlorophyll contents varied. Transcripts with expression levels related to chlorophyll content were statistically selected using RNA-Seq data from 21 F2 plants. In regression models, frequently selected genes included four photosynthetic and two stress-responsive genes. The expression values of four photosynthetic and two stress-responsive genes in high-frequency selected genes were significantly correlated with chlorophyll content not only in plants analysed using RNA-Seq but also in 95 F2 plants.

lncExplore: a database of pan-cancer analysis and systematic functional annotation for lncRNAs from RNA-sequencing data

Over the past few years, with the rapid growth of deep-sequencing technology and the development of computational prediction algorithms, a large number of long non-coding RNAs (lncRNAs) have been identified in various types of human cancers. Therefore, it has become critical to determine how to properly annotate the potential function of lncRNAs from RNA-sequencing (RNA-seq) data and arrange the robust information and analysis into a useful system readily accessible by biological and clinical researchers. In order to produce a collective interpretation of lncRNA functions, it is necessary to integrate different types of data regarding the important functional diversity and regulatory role of these lncRNAs. In this study, we utilized transcriptomic sequencing data to systematically observe and identify lncRNAs and their potential functions from 5034 The Cancer Genome Atlas RNA-seq datasets covering 24 cancers. Then, we constructed the 'lncExplore' database that was developed to comprehensively integrate various types of genomic annotation data for collective interpretation. The distinctive features in our lncExplore database include (i) novel lncRNAs verified by both coding potential and translation efficiency score, (ii) pan-cancer analysis for studying the significantly aberrant expression across 24 human cancers, (iii) genomic annotation of lncRNAs, such as cis-regulatory information and gene ontology, (iv) observation of the regulatory roles as enhancer RNAs and competing endogenous RNAs and (v) the findings of the potential lncRNA biomarkers for the user-interested cancers by integrating clinical information and disease specificity score. The lncExplore database is to our knowledge the first public lncRNA annotation database providing cancer-specific lncRNA expression profiles for not only known but also novel lncRNAs, enhancer RNAs annotation and clinical analysis based on pan-cancer analysis. lncExplore provides a more complete pathway to highly efficient, novel and more comprehensive translation of laboratory discoveries into the clinical context and will assist in reinterpreting the biological regulatory function of lncRNAs in cancer research. Database URL http://lncexplore.bmi.nycu.edu.tw.

Signature of prognostic epithelial-mesenchymal transition related long noncoding RNAs (ERLs) in hepatocellular carcinoma

Reliable biomarkers are of great significance for the treatment and diagnosis of hepatocellular carcinoma (HCC). This study identified potential prognostic epithelial-mesenchymal transition related lncRNAs (ERLs) by the cancer genome atlas (TCGA) database and bioinformatics.The differential expression of long noncoding RNA (lncRNA) was obtained by analyzing the lncRNA data of 370 HCC samples in TCGA. Then, Pearson correlation analysis was carried out with EMT related genes (ERGs) from molecular signatures database. Combined with the univariate Cox expression analysis of the total survival rate of hepatocellular carcinoma (HCC) patients, the prognostic ERLs were obtained. Then use "step" function to select the optimal combination of constructing multivariate Cox expression model. The expression levels of ERLs in HCC samples were verified by real-time quantitative polymerase chain reaction.Finally, we identified 5 prognostic ERLs (AC023157.3, AC099850.3, AL031985.3, AL365203.2, CYTOR). The model showed that these prognostic markers were reliable independent predictors of risk factors (P value <.0001, hazard ratio [HR] = 2.400, 95% confidence interval [CI] = 1.667-3.454 for OS). In the time-dependent receiver operating characteristic analysis, this prognostic marker is a good predictor of HCC survival (area under the curve of 1 year, 2 years, 3 years, and 5 years are 0.754, 0.720, 0.704, and 0.662 respectively). We analyzed the correlation of clinical characteristics of these prognostic markers, and the results show that this prognostic marker is an independent factor that can predict the prognosis of HCC more accurately. In addition, by matching with the Molecular Signatures Database, we obtained 18 ERLs, and then constructed the HCC prognosis model and clinical feature correlation analysis using 5 prognostic ERLs. The results show that these prognostic markers have reliable independent predictive value. Bioinformatics analysis showed that these prognostic markers were involved in the regulation of EMT and related functions of tumor occurrence and migration.Five prognostic types of ERLs identified in this study can be used as potential biomarkers to predict the prognosis of HCC.

Localised expression of OsIAA29 suggests a key role for auxin in regulating development of the dorsal aleurone of early rice grains

Non-canonical AUX/IAA protein, OsIAA29, and ZmMPR-1 homologues, OsMRPLs, are part of an auxin-related signalling cascade operating in the dorsal aleurone during early rice grain development. Endosperm of rice and other cereals accumulates high concentrations of the predominant in planta auxin, indole-3-acetic acid (IAA) during early grain development. However, IAA signalling and function during endosperm development are poorly understood. Here, we report that OsYUC12 (an auxin biosynthesis gene) and OsIAA29 (encoding a non-canonical AUX/IAA) are both expressed exclusively in grains, reaching a maximum 5-6 days after pollination. OsYUC12 expression is localised in the aleurone, sub-aleurone and embryo, whereas OsIAA29 expression is restricted to a narrow strip in the dorsal aleurone, directly under the vascular bundle. Although rice has been reported to lack endosperm transfer cells (ETCs), this region of the aleurone is enriched with sugar transporters and is likely to play a key role in apoplastic nutrient transfer, analogous to ETCs in other cereals. OsIAA29 has orthologues only in grass species; expression of which is also specific to early grain development. OsYUC12 and OsIAA29 are temporally co-expressed with two genes (AL1 and OsPR602) previously linked to the development of dorsal aleurone or ETCs. Also up-regulated at the same time is a cluster of MYB-related genes (designated OsMRPLs) homologous to ZmMRP-1, which regulates maize ETC development. Wheat homologues of ZmMRP-1 are similarly expressed in ETCs. Although previous work has suggested that other cereals do not have orthologues of ZmMRP-1, our work suggests OsIAA29 and OsMRPLs and their homologues in other grasses are part of an auxin-regulated, conserved signalling network involved in the differentiation of cells with ETC-like function in developing cereal grains.

Reinvestigation of THOUSAND-GRAIN WEIGHT 6 grain weight genes in wheat and rice indicates a role in pollen development rather than regulation of auxin content in grains

Phylogenetic and expression analyses of grain weight genes TaTGW6 and OsTGW6 and investigation of substrate availability indicate TGW6 does not regulate auxin content of grains but may affect pollen development. The THOUSAND-GRAIN WEIGHT 6 genes (TaTGW6 and OsTGW6) are reported to result in larger grains of wheat and rice by reducing production of indole-3-acetic acid (IAA) in developing grains. However, a critical comparison of data on TaTGW6 and OsTGW6 with other reports on IAA synthesis in cereal grains requires that this hypothesis be reinvestigated. Here, we show that TaTGW6 and OsTGW6 are members of a large gene family that has undergone major, lineage-specific gene expansion. Wheat has nine genes, and rice three genes encoding proteins with more than 80% amino acid identity with TGW6, making it difficult to envisage how a single inactive allele could have a major effect on IAA levels in grains. In our study, we show that neither TaTGW6 nor OsTGW6 is expressed in developing grains. Instead, both genes and their close homologues are exclusively expressed in pre-emergent inflorescences; TaTGW6 is expressed particularly in microspores prior to mitosis. This evidence, combined with our observation that developing wheat grains have undetectable levels of ester IAA in comparison to free IAA and do not express an IAA-glucose synthase suggests that TaTGW6 and OsTGW6 do not regulate grain size via the hydrolysis of IAA-glucose. Instead, their similarity to rice strictosidine synthase-like (OsSTRL2) suggests they play a key role in pollen development.

Expression of key auxin biosynthesis genes correlates with auxin and starch content of developing wheat (Triticum aestivum) grains

The effect of auxin on wheat (Triticum aestivum L.) grain size is contentious. Additionally, the contributions to the IAA pool from de novo synthesis versus hydrolysis of IAA-glucose are unclear. Here, we describe the first comprehensive study of tryptophan aminotransferase and indole-3-pyruvate mono-oxygenase expression from 5 to 20 days after anthesis. A comparison of expression data with measurements of endogenous IAA via combined liquid chromatography-tandem mass spectrometry using heavy isotope labelled internal standards indicates that TaTAR2-B3, TaYUC9-A1, TaYUC9-B, TaYUC9-D1, TaYUC10-A and TaYUC10-D are primarily responsible for IAA production in developing grains. Furthermore, these genes are expressed specifically in developing grains, like those found in rice (Oryza sativa L.) and maize (Zea mays L.). Our results cast doubt on the proposed role of THOUSAND-GRAIN WEIGHT gene, TaTGW6, in promoting larger grain size via negative effects on grain IAA content. Work on this gene overlooked the contribution of IAA biosynthesis from tryptophan. Although IAA synthesis occurs primarily in the endosperm, we show the TaYUC9-1 group is also strongly expressed in the embryo. Within the endosperm, TaYUC9-1 expression is highest in aleurone and transfer cells, suggesting that IAA has a key role in differentiation of these tissues as has been proposed for other cereals.

Allelic variation in sugary1 gene affecting kernel sweetness among diverse-mutant and -wild-type maize inbreds

Sweet corn is popular worldwide as vegetable. Though large numbers of sugary1 (su1)-based sweet corn germplasm are available, allelic diversity in su1 gene encoding SU1 isoamylase among diverse maize inbreds has not been analyzed. Here, we characterized the su1 gene in maize and compared with allied species. The entire su1 gene (11,720 bp) was sequenced among six mutant (su1) and five wild (Su1) maize inbreds. Fifteen InDels of 2-45 bp were selected to develop markers for studying allelic diversity in su1 gene among 19 mutant- (su1) and 29 wild-type (Su1) inbreds. PIC ranged from 0.15 (SU-InDel7) to 0.37 (SU-InDel13). Major allele frequency varied from 0.52 to 0.90, while gene diversity ranged from 0.16 to 0.49. Phylogenetic tree categorized 48 maize inbreds in two clusters each for wild- type (Su1) and mutant (su1) types. 44 haplotypes of su1 were observed, with three haplotypes (Hap6, Hap22 and Hap29) sharing more than one genotype. Further, comparisons were made with 23 orthologues of su1 from 16 grasses and Arabidopsis. Maize possessed 15-19 exons in su1, while it was 11-24 exons among orthologues. Introns among the orthologues were longer (77-2206 bp) than maize (859-1718 bp). SU1 protein of maize and orthologues had conserved α-amylase and CBM_48 domains. The study also provided physicochemical properties and secondary structure of SU1 protein in maize and its orthologues. Phylogenetic analysis showed closer relationship of maize SU1 protein with P. hallii, S. bicolor and E. tef than Triticum sp. and Oryza sp. The study showed that presence of high allelic diversity in su1 gene which can be utilized in the sweet corn breeding program. This is the first report of comprehensive characterization of su1 gene and its allelic forms in diverse maize and related orthologues.

Differential responses of anthers of stress tolerant and sensitive wheat cultivars to high temperature stress

Transcriptomic analyses identified anther-expressed genes in wheat likely to contribute to heat tolerance and hence provide useful genetic markers. The genes included those involved in hormone biosynthesis, signal transduction, the heat shock response and anther development. Pollen development is particularly sensitive to high temperature heat stress. In wheat, heat-tolerant and heat-sensitive cultivars have been identified, although the underlying genetic causes for these differences are largely unknown. The effects of heat stress on the developing anthers of two heat-tolerant and two heat-sensitive wheat cultivars were examined in this study. Heat stress (35 °C) was found to disrupt pollen development in the two heat-sensitive wheat cultivars but had no visible effect on pollen or anther development in the two heat-tolerant cultivars. The sensitive anthers exhibited a range of developmental abnormalities including an increase in unfilled and clumped pollen grains, abnormal pollen walls and a decrease in pollen viability. This subsequently led to a greater reduction in grain yield in the sensitive cultivars following heat stress. Transcriptomic analyses of heat-stressed developing wheat anthers of the four cultivars identified a number of key genes which may contribute to heat stress tolerance during pollen development. Orthologs of some of these genes in Arabidopsis and rice are involved in regulation of the heat stress response and the synthesis of auxin, ethylene and gibberellin. These genes constitute candidate molecular markers for the breeding of heat-tolerant wheat lines.