TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions

RNA-Seq Data Analysis: A Practical Guide for Model and Non-Model Organisms

RNA sequencing (RNA-seq) has emerged as a powerful tool for assessing genome-wide gene expression, revolutionizing various fields of biology. However, analyzing large RNA-seq datasets can be challenging, especially for students or researchers lacking bioinformatics experience. To address these challenges, we present a comprehensive guide to provide step-by-step workflows for analyzing RNA-seq data, from raw reads to functional enrichment analysis, starting with considerations for experimental design. This is designed to aid students and researchers working with any organism, irrespective of whether an assembled genome is available. Within this guide, we employ various recognized bioinformatics tools to navigate the landscape of RNA-seq analysis and discuss the advantages and disadvantages of different tools for the same task. Our protocol focuses on clarity, reproducibility, and practicality to enable users to navigate the complexities of RNA-seq data analysis easily and gain valuable biological insights from the datasets. Additionally, all scripts and a sample dataset are available in a GitHub repository to facilitate the implementation of the analysis pipeline. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Analysis of data from a model plant with an available reference genome Basic Protocol 2: Gene ontology enrichment analysis Basic Protocol 3: De novo assembly of data from non-model plants.

Integrated multi-omics analysis reveals drought stress response mechanism in chickpea (Cicer arietinum L.)

Drought is one of the major constraints limiting chickpea productivity. To unravel complex mechanisms regulating drought response in chickpea, we generated transcriptomics, proteomics, and metabolomics datasets from root tissues of four contrasting drought-responsive chickpea genotypes: ICC 4958, JG 11, and JG 11+ (drought-tolerant), and ICC 1882 (drought-sensitive) under control and drought stress conditions. Integration of transcriptomics and proteomics data identified enriched hub proteins encoding isoflavone 4'-O-methyltransferase, UDP-d-glucose/UDP-d-galactose 4-epimerase, and delta-1-pyrroline-5-carboxylate synthetase. These proteins highlighted the involvement of pathways such as antibiotic biosynthesis, galactose metabolism, and isoflavonoid biosynthesis in activating drought stress response mechanisms. Subsequently, the integration of metabolomics data identified six metabolites (fructose, galactose, glucose, myoinositol, galactinol, and raffinose) that showed a significant correlation with galactose metabolism. Integration of root-omics data also revealed some key candidate genes underlying the drought-responsive "QTL-hotspot" region. These results provided key insights into complex molecular mechanisms underlying drought stress response in chickpea.

NQO1 regulates expression and alternative splicing of apoptotic genes associated with Alzheimer's disease in PC12 cells

PURPOSE

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by progressive memory loss and cognitive dysfunction. Quinone oxidoreductase 1 (NQO1) is an antioxidant enzyme that plays an important role in controlling cellular redox state, whose expression is altered in the brain tissues of AD patients. In addition to its traditional antioxidant effects, NQO1 also acts as a multifunctional RNA-binding protein involved in posttranscriptional regulation. Whether the RNA-binding activity of NQO1 influences AD pathology has not been investigated yet.

METHODS

The RNA-binding functions of NQO1 in rat pheochromocytoma (PC12) cells were investigated using siRNA knockdown followed by total RNA sequencing. Reverse transcription quantitative polymerase chain reaction was performed to explore the impact of NQO1 on the transcription and alternative splicing of apoptotic genes.

RESULTS

NQO1 knockdown led to a significant increase in cellular apoptosis. Genes involved in certain apoptosis pathways, such as positive regulation of apoptotic processes and mitogen-activated protein kinase signaling, were under global transcriptional and alternative splicing regulation. NQO1 regulated the transcription of apoptotic genes Cryab, Lgmn, Ngf, Apoe, Brd7, and Stat3, as well as the alternative splicing of apoptotic genes BIN1, Picalm, and Fyn.

CONCLUSION

Our findings suggest that NQO1 participates in the pathology of AD by regulating the expression and alternative splicing of the genes involved in apoptosis. These results extend our understanding of NQO1 in apoptotic pathways at the posttranscriptional level in AD.

Rfoot-seq: Transcriptomic RNase Footprinting for Mapping Stable RNA-Protein Complexes and Rapid Ribosome Profiling

Ribosome profiling isolates ribosome-protected fragments for sequencing and is a valuable method for studying different aspects of RNA translation. However, conventional protocols require millions of input cells and time-consuming steps to isolate translating ribosome complexes using ultracentrifugation or immunoprecipitation. These limitations have prevented their application to rare physiological samples. To address these technical barriers, we developed an RNase footprinting approach named Rfoot-seq to map stable transcriptomic RNA-protein complexes that allows rapid ribosome profiling using low-input samples (Li, Yang, Stroup, Wang, & Ji, 2022). In this assay, we treat a cell lysate with concentrated RNase without complex crosslinking and retained only RNA footprints associated with stable complexes for sequencing. The footprints in coding regions represent ribosome-protected fragments and can be used to study cytosolic and mitochondrial translation simultaneously. Rfoot-seq achieves comparable results to conventional ribosome profiling to quantify ribosome occupancy and works robustly for various cultured cells and primary tissue samples. Moreover, Rfoot-seq maps RNA fragments associated with stable non-ribosomal RNA-protein complexes in noncoding domains of small noncoding RNAs and some long noncoding RNAs. Taken together, Rfoot-seq opens an avenue to quantify transcriptomic translation and characterize functional noncoding RNA domains using low-input samples. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Harvesting and lysing adherent cells Alternate Protocol 1: Harvesting and lysing suspension cells Alternate Protocol 2: Harvesting and lysing primary tissue samples Basic Protocol 2: RNase treatment and footprint purification for low-input samples Alternate Protocol 3: RNase treatment and footprint purification for ultra-low-input samples Basic Protocol 3: Library preparation for high-throughput sequencing Support Protocol: Preparation of dsDNA markers for library size selection Basic Protocol 4: Data analysis and quality control after sequencing.

iGDP: An integrated genome decontamination pipeline for wild ciliated microeukaryotes

Ciliates are a large group of ubiquitous and highly diverse single-celled eukaryotes that play an essential role in the functioning of microbial food webs. However, their genomic diversity is far from clear due to the need to develop cultivation methods for most species, so most research is based on wild organisms that almost invariably contain contaminants. Here we establish an integrated Genome Decontamination Pipeline (iGDP) that combines homology search, telomere reads-assisted and clustering approaches to filter contaminated ciliate genome assemblies from wild specimens. We benchmarked the performance of iGDP using genomic data from a contaminated ciliate culture and the results showed that iGDP could recall 91.9% of the target sequences with 96.9% precision. We also used a synthetic dataset to offer guidelines for the application of iGDP in the removal of various groups of contaminants. Compared with several popular metagenome binning tools, iGDP could show better performance. To further validate the effectiveness of iGDP on real-world data, we applied it to decontaminate genome assemblies of three wild ciliate specimens and obtained their genomes with high quality comparable to that of previously well-studied model ciliate genomes. It is anticipated that the newly generated genomes and the established iGDP method will be valuable community resources for detailed studies on ciliate biodiversity, phylogeny, ecology and evolution. The pipeline (https://github.com/GWang2022/iGDP) can be implemented automatically to reduce manual filtering and classification and may be further developed to apply to other microeukaryotes.

scRNA-seq data analysis method to improve analysis performance

With the development of single-cell RNA sequencing technology (scRNA-seq), we have the ability to study biological questions at the level of the individual cell transcriptome. Nowadays, many analysis tools, specifically suitable for single-cell RNA sequencing data, have been developed. In this review, the currently commonly used scRNA-seq protocols are discussed. The upstream processing flow pipeline of scRNA-seq data, including goals and popular tools for reads mapping and expression quantification, quality control, normalization, imputation, and batch effect removal is also introduced. Finally, methods to evaluate these tools in both cellular and genetic dimensions, clustering and differential expression analysis are presented.

Epigenetic and pharmacological control of pigmentation via Bromodomain Protein 9 (BRD9)

Lineage-specific differentiation programs are activated by epigenetic changes in chromatin structure. Melanin-producing melanocytes maintain a gene expression program ensuring appropriate enzymatic conversion of metabolites into the pigment, melanin, and transfer to surrounding cells. During neuroectodermal development, SMARCA4 (BRG1), the catalytic subunit of SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodeling complexes, is essential for lineage specification. SMARCA4 is also required for development of multipotent neural crest precursors into melanoblasts, which differentiate into pigment-producing melanocytes. In addition to the catalytic domain, SMARCA4 and several SWI/SNF subunits contain bromodomains which are amenable to pharmacological inhibition. We investigated the effects of pharmacological inhibitors of SWI/SNF bromodomains on melanocyte differentiation. Strikingly, treatment of murine melanoblasts and human neonatal epidermal melanocytes with selected bromodomain inhibitors abrogated melanin synthesis and visible pigmentation. Using functional genomics, iBRD9, a small molecule selective for the bromodomain of BRD9 was found to repress pigmentation-specific gene expression. Depletion of BRD9 confirmed a requirement for expression of pigmentation genes in the differentiation program from melanoblasts into pigmented melanocytes and in melanoma cells. Chromatin immunoprecipitation assays showed that iBRD9 disrupts the occupancy of BRD9 and the catalytic subunit SMARCA4 at melanocyte-specific loci. These data indicate that BRD9 promotes melanocyte pigmentation whereas pharmacological inhibition of BRD9 is repressive.

Performing in spite of starvation: How Saccharomyces cerevisiae maintains robust growth when facing famine zones in industrial bioreactors

In fed-batch operated industrial bioreactors, glucose-limited feeding is commonly applied for optimal control of cell growth and product formation. Still, microbial cells such as yeasts and bacteria are frequently exposed to glucose starvation conditions in poorly mixed zones or far away from the feedstock inlet point. Despite its commonness, studies mimicking related stimuli are still underrepresented in scale-up/scale-down considerations. This may surprise as the transition from glucose limitation to starvation has the potential to provoke regulatory responses with negative consequences for production performance. In order to shed more light, we performed gene-expression analysis of Saccharomyces cerevisiae grown in intermittently fed chemostat cultures to study the effect of limitation-starvation transitions. The resulting glucose concentration gradient was representative for the commercial scale and compelled cells to tolerate about 76 s with sub-optimal substrate supply. Special attention was paid to the adaptation status of the population by discriminating between first time and repeated entry into the starvation regime. Unprepared cells reacted with a transiently reduced growth rate governed by the general stress response. Yeasts adapted to the dynamic environment by increasing internal growth capacities at the cost of rising maintenance demands by 2.7%. Evidence was found that multiple protein kinase A (PKA) and Snf1-mediated regulatory circuits were initiated and ramped down still keeping the cells in an adapted trade-off between growth optimization and down-regulation of stress response. From this finding, primary engineering guidelines are deduced to optimize both the production host's genetic background and the design of scale-down experiments.

Evolution of mitochondrial RNA editing in extant gymnosperms

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

Chronological attenuation of NPRA/PKG/AMPK signaling promotes vascular aging and elevates blood pressure

Hypertension is common in elderly population. We designed to search comprehensively for genes that are chronologically shifted in their expressions and to define their contributions to vascular aging and hypertension. RNA sequencing was conducted to search for senescence-shifted transcripts in human umbilical vein endothelial cells (HUVECs). Small interfering RNA (siRNA), small-molecule drugs, CRISPR/Cas9 techniques, and imaging were used to determine genes' function and contributions to age-related phenotypes of the endothelial cell and blood vessel. Of 25 genes enriched in the term of "regulation of blood pressure," NPRA was changed most significantly. The decreased NPRA expression was replicated in aortas of aged mice. The knockdown of NPRA promoted HUVEC senescence and it decreased expressions of protein kinase cGMP-dependent 1 (PKG), sirtuin 1 (SIRT1), and endothelial nitric oxide synthase (eNOS). Suppression of NPRA also decreased the phosphorylation of AMP-activated protein kinase (AMPK) as well as the ratio of oxidized nicotinamide adenine dinucleotide (NAD⁺ )/reduced nicotinamide adenine dinucleotide (NADH) but increased the production of reactive oxygen species (ROS). 8-Br-cGMP (analog of cGMP), or AICAR (AMPK activator), counteracted the observed changes in HUVECs. The Npr1^+/- mice presented an elevated systolic blood pressure and their vessels became insensitive to endothelial-dependent vasodilators. Further, vessels from Npr1^+/- mice increased Cdkn1a but decreased eNos expressions. These phenotypes were rescued by intravenously administrated 8-Br-cGMP and viral overexpression of human PKG, respectively. In conclusion, we demonstrate NPRA/PKG/AMPK as a novel and critical signaling axis in the modulation of endothelial cell senescence, vascular aging, and hypertension.

Evolutionary fates of gene-body methylation and its divergent association with gene expression in pigeonpea

Pigeonpea (Cajanus cajan L. Huth) is an agronomically important legume cultivated worldwide. In this study, we extensively analyzed gene-body methylation (GbM) patterns in pigeonpea. We found a bimodal distribution of CG and CHG methylation patterns. GbM features- slow evolution rate and increased length remained conserved. Genes with moderate CG body methylation showed highest expression where as highly-methylated genes showed lowest expression. Transposable element (TE)-related genes were methylated in multiple contexts and hence classified as C-methylated genes. A low expression among C-methylated genes was associated with transposons insertion in gene-body and upstream regulatory regions. The CG methylation patterns were found to be conserved in orthologs compared with non-CG methylation. By comparing methylation patterns between differentially methylated regions (DMRs) of the three genotypes, we found that variably methylated marks are less likely to target evolutionary conserved sequences. Finally, our analysis showed enrichment of nitrogen-related genes in GbM orthologs of legumes, which could be promising candidates for generating epialleles for crop improvement.

Pathogen profiling of Australian rabbits by metatranscriptomic sequencing

Australia is known for its long history of using biocontrol agents, such as myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV), to manage wild European rabbit populations. Interestingly, while undertaking RHDV surveillance of rabbits that were found dead, we observed that approximately 40% of samples were negative for RHDV. To investigate whether other infectious agents are responsible for killing rabbits in Australia, we subjected a subset of these RHDV-negative liver samples to metatranscriptomic sequencing. In addition, we investigated whether the host transcriptome data could provide additional differentiation between likely infectious versus non-infectious causes of death. We identified transcripts from several Clostridia species, Pasteurella multocida, Pseudomonas spp., and Eimeria stiedae, in liver samples of several rabbits that had died suddenly, all of which are known to infect rabbits and are capable of causing disease and mortality. In addition, we identified Hepatitis E virus and Cyniclomyces yeast in some samples, both of which are not usually associated with severe disease. In one-third of the sequenced total liver RNAs, no infectious agent could be identified. While metatranscriptomic sequencing cannot provide definitive evidence of causation, additional host transcriptome analysis provided further insights to distinguish between pathogenic microbes and commensals or environmental contaminants. Interestingly, three samples where no pathogen could be identified showed evidence of up-regulated host immune responses, while immune response pathways were not up-regulated when E. stiedae, Pseudomonas, or yeast were detected. In summary, although no new putative rabbit pathogens were identified, this study provides a robust workflow for future investigations into rabbit mortality events.

Integration analyses of structural variations and differential gene expression associated with beef fatty acid profile in Nellore cattle

This study aimed to integrate analyses of structural variations and differentially expressed genes (DEGs) associated with the beef fatty acid (FA) profile in Nellore cattle. Copy numbers variation (CNV) detection was performed using the penncnv algorithm and CNVRuler software in 3794 genotyped animals through the High-Density Bovine BeadChip. In order to perform the genomic wide association study (GWAS), a total of 963 genotyped animals were selected to obtain the intramuscular lipid concentration and quantify the beef FA profile. A total of 48 animals belonging to the same farm and management lot were extracted from the 963 genotyped and phenotyped animals to carry out the transcriptomic and differentially expressed gene analyses. The GWAS with extreme groups of FA profiles was performed using a logistic model. A total of 43, 42, 66 and 35 significant CNV regions (p < 0.05) for saturated, monounsaturated, polyunsaturated and omega 3 and 6 fatty acids were identified respectively. The paired-end sequencing of 48 samples was performed using the Illumina HiSeq2500 platform. Real-time quantitative PCR was used to validate the DEGs identified by RNA-seq analysis. The results showed several DEGs associated with the FA profile of Longissimus thoracis, such as BSCL2 and SAMD8. Enriched terms as the cellular response to corticosteroid (GO:0071384) and glucocorticoid stimulus (GO:0071385) could be highlighted. The identification of structural variations harboring candidate genes for beef FA must contribute to the elucidation of the genetic basis that determines the beef FA composition of intramuscular fat in Nellore cattle. Our results will contribute to the identification of potential biomarkers for complex phenotypes, such as the FA profile, to improve the reliability of the genomic predictions including pre-selected variants using differentiated weighting in the genomic models.

Parallel meta-transcriptome analysis reveals degradation of plant secondary metabolites by beetles and their gut symbionts

Switching to a new host plant is a driving force for divergence and speciation in herbivorous insects. This process of incorporating a novel host plant into the diet may require a number of adaptations in the insect herbivores that allow them to consume host plant tissue that may contain toxic secondary chemicals. As a result, herbivorous insects are predicted to have evolved efficient ways to detoxify major plant defenses and increase fitness by either relying on their own genomes or by recruiting other organisms such as microbial gut symbionts. In the present study we used parallel meta-transcriptomic analyses of Altica flea beetles and their gut symbionts to explore the contributions of beetle detoxification mechanisms versus detoxification by their gut consortium. We compared the gut meta-transcriptomes of two sympatric Altica species that feed exclusively on different host plant species as well as their F₁ hybrids that were fed one of the two host plant species. These comparisons revealed that gene expression patterns of Altica are dependent on both beetle species identity and diet. The community structure of gut symbionts was also dependent on the identity of the beetle species, and the gene expression patterns of the gut symbionts were significantly correlated with beetle species and plant diet. Some of the enriched genes identified in the beetles and gut symbionts are involved in the degradation of secondary metabolites produced by plants, suggesting that Altica flea beetles may use their gut microbiota to help them feed on and adapt to their host plants.

Experimental sexual selection reveals rapid evolutionary divergence in sex-specific transcriptomes and their interactions following mating

Post copulatory interactions between the sexes in internally fertilizing species elicits both sexual conflict and sexual selection. Macroevolutionary and comparative studies have linked these processes to rapid transcriptomic evolution in sex‐specific tissues and substantial transcriptomic post mating responses in females, patterns of which are altered when mating between reproductively isolated species. Here, we tested multiple predictions arising from sexual selection and conflict theory about the evolution of sex‐specific and tissue‐specific gene expression and the post mating response at the microevolutionary level. Following over 150 generations of experimental evolution under either reduced (enforced monogamy) or elevated (polyandry) sexual selection in Drosophila pseudoobscura, we found a substantial effect of sexual selection treatment on transcriptomic divergence in virgin male and female reproductive tissues (testes, male accessory glands, the female reproductive tract and ovaries). Sexual selection treatment also had a dominant effect on the post mating response, particularly in the female reproductive tract – the main arena for sexual conflict – compared to ovaries. This effect was asymmetric with monandry females typically showing more post mating responses than polyandry females, with enriched gene functions varying across treatments. The evolutionary history of the male partner had a larger effect on the post mating response of monandry females, but females from both sexual selection treatments showed unique patterns of gene expression and gene function when mating with males from the alternate treatment. Our microevolutionary results mostly confirm comparative macroevolutionary predictions on the role of sexual selection on transcriptomic divergence and altered gene regulation arising from divergent coevolutionary trajectories between sexual selection treatments.

Chromosome-Scale Genome Assembly Provides Insights into Speciation of Allotetraploid and Massive Biomass Accumulation of Elephant Grass (Pennisetum purpureum Schum.)

Elephant grass (Pennisetum purpureum Schum) is an important forage, biofuels and industrial plant widely distributed in tropical and subtropical areas globally. It is characterized with robust growth and high biomass. We sequenced its allopolyploid genome and assembled 2.07 Gb into A' and B sub-genomes of 14 chromosomes with scaffold N50 of 8.47 Mb, yielding a total of 77,139 genes. The allotetraploid speciation occurred approximately 15 MYA after the divergence between Setaria italica and Pennisetum glaucum, according to a phylogenetic analysis of Pennisetum species. Double whole-genome duplication (WGD) and polyploidization events resulted in large scale gene expansion, especially in the key steps of growth and biomass accumulation. Integrated transcriptome profiling revealed the functional divergence between sub-genomes A' and B. A' sub-genome mainly contributed to plant growth, development and photosynthesis, whereas the B sub-genome was primarily responsible for effective transportation and resistance to stimulation. Some key gene families related to cellulose biosynthesis were expanded and highly expressed in stems, which could explain the high cellulose content in elephant grass. Our findings provide deep insights into genetic evolution of elephant grass and will aid future biological research and breeding, even for other grasses in the family Poaceae.

Transcriptome analysis of messenger RNA and long noncoding RNA related to different developmental stages of tail adipose tissues of sunite sheep

The tail fat of sheep is the most typical deposited fat, and it can be widely used in human daily life, such as diet, cosmetics, and industrial raw materials. To understand the potential regulatory mechanism of different growth stages of tail fat in Sunite sheep, we performed high-throughput RNA sequencing to characterize the long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles of the sheep tail fat at the age of 6, 18, and 30 months. A total of 223 differentially expressed genes (DEGs) and 148 differentially expressed lncRNAs were found in the tail fat of 6-, 18-, and 30-month-old sheep. Based on functional analysis, we found that fat-related DEGs were mainly expressed at 6 months of age and gradually decreased at 18 and 30 months of age. The target gene prediction analysis shows that most of the lncRNAs target more than 20 mRNAs as their transregulators. Further, we obtained several fat-related differentially expressed target genes; these target genes interact with different differentially expressed lncRNAs at various ages and play an important role in the development of tail fat. Based on the DEGs and differentially expressed lncRNAs, we established three co-expression networks for each comparison group. Finally, we concluded that the development of the sheep tail fat is more active during the early stage of growth and gradually decreases with the increase in age. The mutual regulation of lncRNAs and mRNAs may play a key role in this complex biological process.

Can functional genomic diversity provide novel insights into mechanisms of community assembly? A pilot study from an invaded alpine streambed

An important focus of community ecology, including invasion biology, is to investigate functional trait diversity patterns to disentangle the effects of environmental and biotic interactions. However, a notable limitation is that studies usually rely on a small and easy-to-measure set of functional traits, which might not immediately reflect ongoing ecological responses to changing abiotic or biotic conditions, including those that occur at a molecular or physiological level. We explored the potential of using the diversity of expressed genes-functional genomic diversity (FGD)-to understand ecological dynamics of a recent and ongoing alpine invasion. We quantified FGD based on transcriptomic data measured for 26 plant species occurring along adjacent invaded and pristine streambeds. We used an RNA-seq approach to summarize the overall number of expressed transcripts and their annotations to functional categories, and contrasted this with functional trait diversity (FTD) measured from a suite of characters that have been traditionally considered in plant ecology. We found greater FGD and FTD in the invaded community, independent of differences in species richness. However, the magnitude of functional dispersion was greater from the perspective of FGD than from FTD. Comparing FGD between congeneric alien-native species pairs, we did not find many significant differences in the proportion of genes whose annotations matched functional categories. Still, native species with a greater relative abundance in the invaded community compared with the pristine tended to express a greater fraction of genes at significant levels in the invaded community, suggesting that changes in FGD may relate to shifts in community composition. Comparisons of diversity patterns from the community to the species level offer complementary insights into processes and mechanisms driving invasion dynamics. FGD has the potential to illuminate cryptic changes in ecological diversity, and we foresee promising avenues for future extensions across taxonomic levels and macro-ecosystems.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

Multiple Maize Reference Genomes Impact the Identification of Variants by Genome-Wide Association Study in a Diverse Inbred Panel

Use of a single reference genome for genome-wide association studies (GWAS) limits the gene space represented to that of a single accession. This limitation can complicate identification and characterization of genes located within presence-absence variations (PAVs). In this study, we present the draft de novo genome assembly of 'PHJ89', an 'Oh43'-type inbred line of maize ( L.). From three separate reference genome assemblies ('B73', 'PH207', and PHJ89) that represent the predominant germplasm groups of maize, we generated three separate whole-seedling gene expression profiles and single nucleotide polymorphism (SNP) matrices from a panel of 942 diverse inbred lines. We identified 34,447 (B73), 39,672 (PH207), and 37,436 (PHJ89) transcripts that are not present in the respective reference genome assemblies. Genome-wide association studies were conducted in the 942 inbred panel with both the SNP and expression data values to map (SCMV) resistance. Highlighting the impact of alternative reference genomes in gene discovery, the GWAS results for SCMV resistance with expression values as a surrogate measure of PAV resulted in robust detection of the physical location of a known resistance gene when the B73 reference that contains the gene was used, but not the PH207 reference. This study provides the valuable resource of the Oh43-type PHJ89 genome assembly as well as SNP and expression data for 942 individuals generated from three different reference genomes.

Hybrid Genome Assembly of a Major Quantitative Disease Resistance Locus in Soybean Toward Fusarium graminearum

Schwabe [teleomorph: Gibberella zeae (Schweintiz) Petch] has been identified as a pathogen of soybean [ (L.) Merr.] causing seed, seedling damping-off and root rot in North America. A major quantitative disease resistance locus (QDRL) that contributed 38.5% of the phenotypic variance toward in soybean was previously identified through mapping of a recombinant inbred line (RIL) population derived from a cross between 'Wyandot' and PI 567301B. This major QDRL mapped to chromosome 8 to a predicted 305 kb region harboring 36 genes. This locus maps near the locus for soybean cyst nematode (SCN) and the locus contributing to seed coat color. Long-read sequencing of the region was completed and variations in gene sequence and gene order compared with the 'Williams 82' reference were identified. Molecular markers were developed for genes within this region and mapped in the original population, slightly narrowing the region of interest. Analyses of the hybrid genome reassembly using three previously published bacterial artificial chromosome (BAC) sequences (BAC56G2, BAC104J7, and BAC77G7-a) combined with RNA-sequencing narrowed the region making candidate gene identification possible. The markers within this region may be used for marker-assisted selection (MAS). There were 10 differentially expressed genes between resistant and susceptible lines, with four of these candidates also located within the genomic interval defined by the flanking markers. These genes included an actin-related protein 2/3 complex subunit, an unknown protein, a hypothetical protein, and a chalcone synthase 3.

Natural variation for gene expression responses to abiotic stress in maize

Plants respond to abiotic stress through a variety of physiological, biochemical, and transcriptional mechanisms. Many genes exhibit altered levels of expression in response to abiotic stress, which requires concerted action of both cis- and trans-regulatory features. In order to study the variability in transcriptome response to abiotic stress, RNA sequencing was performed using 14-day-old maize seedlings of inbreds B73, Mo17, Oh43, PH207 and B37 under control, cold and heat conditions. Large numbers of genes that responded differentially to stress between parental inbred lines were identified. RNA sequencing was also performed on similar tissues of the F₁ hybrids produced by crossing B73 and each of the three other inbred lines. By evaluating allele-specific transcript abundance in the F₁ hybrids, we were able to measure the abundance of cis- and trans-regulatory variation between genotypes for both steady-state and stress-responsive expression differences. Although examples of trans-regulatory variation were observed, cis-regulatory variation was more common for both steady-state and stress-responsive expression differences. The genes with cis-allelic variation for response to cold or heat stress provided an opportunity to study the basis for regulatory diversity.

Comprehensive definition of genome features in Spirodela polyrhiza by high-depth physical mapping and short-read DNA sequencing strategies

Spirodela polyrhiza is a fast-growing aquatic monocot with highly reduced morphology, genome size and number of protein-coding genes. Considering these biological features of Spirodela and its basal position in the monocot lineage, understanding its genome architecture could shed light on plant adaptation and genome evolution. Like many draft genomes, however, the 158-Mb Spirodela genome sequence has not been resolved to chromosomes, and important genome characteristics have not been defined. Here we deployed rapid genome-wide physical maps combined with high-coverage short-read sequencing to resolve the 20 chromosomes of Spirodela and to empirically delineate its genome features. Our data revealed a dramatic reduction in the number of the rDNA repeat units in Spirodela to fewer than 100, which is even fewer than that reported for yeast. Consistent with its unique phylogenetic position, small RNA sequencing revealed 29 Spirodela-specific microRNA, with only two being shared with Elaeis guineensis (oil palm) and Musa balbisiana (banana). Combining DNA methylation data and small RNA sequencing enabled the accurate prediction of 20.5% long terminal repeats (LTRs) that doubled the previous estimate, and revealed a high Solo:Intact LTR ratio of 8.2. Interestingly, we found that Spirodela has the lowest global DNA methylation levels (9%) of any plant species tested. Taken together our results reveal a genome that has undergone reduction, likely through eliminating non-essential protein coding genes, rDNA and LTRs. In addition to delineating the genome features of this unique plant, the methodologies described and large-scale genome resources from this work will enable future evolutionary and functional studies of this basal monocot family.

Rascaf: Improving Genome Assembly with RNA Sequencing Data

Abundant but short second-generation sequencing reads make assembly difficult, leading to fragmented genomes and gene annotations. Gene structure information from RNA sequences can be used to improve the completeness and contiguity of an assembly, but bioinformatics methods have been lacking. Rascaf is a highly efficient tool leveraging long-range continuity information from intron spanning RNA sequencing (RNA-seq) read pairs to detect new contig connections. It determines a heaviest path in an exon block graph that simultaneously represents a gene and the underlying contig relationships. Rascaf is more accurate than its competitors, highly precise, and finds thousands of new verifiable connections in several draft Rosaceae genomes. Lightweight and practical, it can be readily incorporated into sequencing pipelines to improve an assembly and its gene annotations.

Candidate Gene Sequence Analyses toward Identifying -Type Resistance to

is one of three genetic loci conferring strain-specific resistance to (SMV). The locus has been mapped to a 154-kb region on chromosome 14, containing a cluster of five nucleotide-binding leucine-rich repeat (NB-LRR) resistance genes. High sequence similarity between the candidate genes challenges fine mapping of the locus. Among the five, Glyma14g38533 showed the highest transcript abundance in 1 to 3 h of SMV-G7 inoculation. Comparative sequence analyses were conducted with the five candidate NB-LRR genes from susceptible (-type) soybean [ (L.) Merr.] cultivar Williams 82, resistant (-type) cultivar Hwangkeum, and resistant lines L29 and RRR. Sequence comparisons revealed that Glyma14g38533 had far more polymorphisms than the other candidate genes. Interestingly, Glyma14g38533 gene from -type lines exhibited 150 single-nucleotide polymorphism (SNP and six insertion-deletion (InDel) markers relative to -type line, Furthermore, the polymorphisms identified in three -type lines were highly conserved. Several polymorphisms were validated in 18 -type resistant and six -type susceptible lines and were found associated with their disease response. The majority of the polymorphisms were located in LRR domain encoding region, which is involved in pathogen recognition via protein-protein interactions. These findings associating Glyma14g38533 with -type resistance to SMV suggest it is the most likely candidate gene for .

Global Transcriptome Analysis Reveals Differences in Gene Expression Patterns Between Nonhyperhydric and Hyperhydric Peach Leaves

Hyperhydricity is a morphophysiological disorder of plants in tissue culture characterized morphologically by the presence of translucent, thick, curled, and fragile leaves as a result of excessive water intake. Since clonal propagation is a major in vitro technique for multiplying plants vegetatively, the emergence of hyperhydricity-related symptoms causes significant economic losses to agriculture and horticulture. Although numerous efforts have been hitherto devoted to the morphological and anatomical responses of plants to hyperhydricity, the underlying molecular mechanism remains largely unknown. Here, a genome-wide transcriptome analysis was performed to identify differentially expressed genes in hyperhydric and nonhyperhydric leaves of peach [ (L.) Batsch]. The RNA sequencing (RNA-Seq) analysis showed that the expression of >300 transcripts was altered between control and hyperhydric leaf cells. The top 30 differentially expressed transcripts (DETs) were related to the posttranscriptional regulators of organelle gene expression and photosynthesis, cellular elimination, plant cuticle development, and abiotic stress response processes. The expression of 10 DETs was also conformed by quantitative real-time polymerase chain reaction (RT-qPCR) in hyperhydric and nonhyperhydric leaves. As a complex biological process, hyperhydricity alters the expression of various transcripts including transcription factor (), RNA binding protein (pentatricopeptide, ), transporter protein (), and . Thus, this genome-wide transcriptome profiling study may help elucidate the molecular mechanism of hyperhydricity.

Transcriptome Profiling of Rust Resistance in Switchgrass Using RNA-Seq Analysis

Switchgrass rust caused by Puccinia emaculata is a major limiting factor for switchgrass (Panicum virgatum L.) production, especially in monoculture. Natural populations of switchgrass displayed diverse reactions to P. emaculata when evaluated in an Ardmore, OK, field. To identify the differentially expressed genes during the rust infection process and the mechanisms of switchgrass rust resistance, transcriptome analysis using RNA-Seq was conducted in two pseudo-F₁ parents ('PV281' and 'NFGA472'), and three moderately resistant and three susceptible progenies selected from a three-generation, four-founder switchgrass population (K5 × A4) × (AP13 × VS16). On average, 23.5 million reads per sample (leaf tissue was collected at 0, 24, and 60 h post-inoculation (hpi)) were obtained from paired-end (2 × 100 bp) sequencing on the Illumina HiSeq2000 platform. Mapping of the RNA-Seq reads to the switchgrass reference genome (AP13 ver. 1.1 assembly) constructed a total of 84,209 transcripts from 98,007 gene loci among all of the samples. Further analysis revealed that host defense-related genes, including the nucleotide binding site-leucine-rich repeat domain containing disease resistance gene analogs, play an important role in resistance to rust infection. Rust-induced gene (RIG) transcripts inherited across generations were identified. The rust-resistant gene transcripts can be a valuable resource for developing molecular markers for rust resistance. Furthermore, the rust-resistant genotypes and gene transcripts identified in this study can expedite rust-resistant cultivar development in switchgrass.