Resource For Functional Genomics Research Articles

Comprehensive mapping and modelling of the rice regulome landscape unveils the regulatory architecture underlying complex traits

Unraveling the regulatory mechanisms that govern complex traits is pivotal for advancing crop improvement. Here we present a comprehensive regulome atlas for rice (Oryza sativa), charting the chromatin accessibility across 23 distinct tissues from three representative varieties. Our study uncovers 117,176 unique open chromatin regions (OCRs), accounting for ~15% of the rice genome, a notably higher proportion compared to previous reports in plants. Integrating RNA-seq data from matched tissues, we confidently predict 59,075 OCR-to-gene links, with enhancers constituting 69.54% of these associations, including many known enhancer-to-gene links. Leveraging this resource, we re-evaluate genome-wide association study results and discover a previously unknown function of OsbZIP06 in seed germination, which we subsequently confirm through experimental validation. We optimize deep learning models to decode regulatory grammar, achieving robust modeling of tissue-specific chromatin accessibility. This approach allows to predict cross-variety regulatory dynamics from genomic sequences, shedding light on the genetic underpinnings of cis-regulatory divergence and morphological disparities between varieties. Overall, our study establishes a foundational resource for rice functional genomics and precision molecular breeding, providing valuable insights into regulatory mechanisms governing complex traits.

De novo transcriptomes of cave and surface isopod crustaceans: insights from 11 species across three suborders

Isopods are a diverse group of crustaceans, that inhabit various environments, including terrestrial, freshwater, and marine, both on the surface and in the underground. The biological mechanisms underlying their wide range of adaptations to diverse ecological niches remain elusive. In order to unravel the molecular basis of their adaptability, we generated a comprehensive RNAseq dataset comprising 11 isopod species belonging to the three different suborders: freshwater Asellota, marine, brackish and freshwater Sphaeromatidea, and terrestrial Oniscidea, with representatives from families Asellidae, Sphaeromatidae, and Trichoniscidae, respectively. Representatives of each family were collected from both cave and surface environments, representing at least three independent cave colonization events. Three biological replicates were sequenced from each species to ensure data robustness. The 11 high-quality RNAseq datasets will serve as a valuable resource for understanding cave-specific adaptations, comparative and functional genomics, ecological annotation as well as aid in conservation efforts of these non-model organisms. Importantly, transcriptomes of eight featured species have been made publicly accessible for the first time.

Advancements of Breeding and Genomics in Wheat (Triticum aestivum L.): Enhancing Yield and Nutritional Value for Sustainable Agriculture

Advancements in wheat breeding and genomics presently explores the genomic interventions driving focusing on quantitative trait loci (QTL) mapping, marker-assisted selection (MAS) and genomic selection (GS). QTL mapping emerges as a pivotal method for pinpointing markers linked with desirable traits, facilitating MAS. Furthermore, genomic selection (GS) holds immense potential for crop improvement. It also delves into the current landscape of MAS and explores various prospects of GS for wheat biofortification. Looking ahead, accelerated mapping studies combined with MAS and GS schemes are poised to further enhance wheat breeding efficiency. Dense molecular maps and a large set of ESTs (Expressed Sequence Tags) have enabled genome-wide identification of gene-rich and gene-poor regions, as well as QTL, including eQTL (Expression quantitative trait loci). Additionally, markers associated with major economic traits have facilitated MAS programs in some countries and enabled map-based cloning of several genes/QTL. Resources for functional genomics, such as TILLING and RNA interference (RNAi), alongside emerging approaches like epigenetics and association mapping, are further enriching wheat genomics research. In this review, we initially present cutting-edge genome-editing technologies in crop plants, with a specific focus on wheat, addressing both functional genomics and genetic enhancement. We subsequently delineate the utilization of additional technologies, including GWAS, high-throughput genotyping and phenotyping, speed breeding, and synthetic biology, within the context of wheat breeding. We assert that integrating genome editing with other molecular breeding strategies will significantly expedite the genetic enhancement of wheat, thus contributing to sustainable global production.

Transcriptome-wide expression landscape and starch synthesis pathway co-expression network in sorghum.

The gene expression landscape across different tissues and developmental stages reflects their biological functions and evolutionary patterns. Integrative and comprehensive analyses of all transcriptomic data in an organism are instrumental to obtaining a comprehensive picture of gene expression landscape. Such studies are still very limited in sorghum, which limits the discovery of the genetic basis underlying complex agricultural traits in sorghum. We characterized the genome-wide expression landscape for sorghum using 873 RNA-sequencing (RNA-seq) datasets representing 19 tissues. Our integrative analysis of these RNA-seq data provides the most comprehensive transcriptomic atlas for sorghum, which will be valuable for the sorghum research community for functional characterizations of sorghum genes. Based on the transcriptome atlas, we identified 595 housekeeping genes (HKGs) and 2080 tissue-specific expression genes (TEGs) for the 19 tissues. We identified different gene features between HKGs and TEGs, and we found that HKGs have experienced stronger selective constraints than TEGs. Furthermore, we built a transcriptome-wide co-expression network (TW-CEN) comprising 35 modules with each module enriched in specific Gene Ontology terms. High-connectivity genes in TW-CEN tend to express at high levels while undergoing intensive selective pressure. We also built global and seed-preferential co-expression networks of starch synthesis pathways, which indicated that photosynthesis and microtubule-based movement play important roles in starch synthesis. The global transcriptome atlas of sorghum generated by this study provides an important functional genomics resource for trait discovery and insight into starch synthesis regulation in sorghum.

Expansion and improvement of ChinaMu by MuT-seq and chromosome-level assembly of the Mu-starter genome.

ChinaMu is the largest sequence-indexed Mutator (Mu) transposon insertional library in maize (Zea mays). In this study, we made significant improvements to the size and quality of the ChinaMu library. We developed a new Mu-tag isolation method Mu-Tn5-seq (MuT-seq). Compared to the previous method used by ChinaMu, MuT-seq recovered 1/3 more germinal insertions, while requiring only about 1/14 of the sequencing volume and 1/5 of the experimental time. Using MuT-seq, we identified 113,879 germinal insertions from 3,168 Mu-active F1 families. We also assembled a high-quality genome for the Mu-active line Mu-starter, which harbors the initial active MuDR element and was used as the pollen donor for the mutation population. Using the Mu-starter genome, we recovered 33,662 (15.6%) additional germinal insertions in 3,244 (7.4%) genes in the Mu-starter line. The Mu-starter genome also improved the assignment of 117,689 (54.5%) germinal insertions. The newly upgraded ChinaMu dataset currently contains 215,889 high-quality germinal insertions. These insertions cover 32,224 pan-genes in the Mu-starter and B73Ref5 genomes, including 23,006 (80.4%) core genes shared by the two genomes. As a test model, we investigated Mu insertions in the pentatricopeptide repeat (PPR) superfamily, discovering insertions for 92% (449/487) of PPR genes in ChinaMu, demonstrating the usefulness of ChinaMu as a functional genomics resource for maize.

Development of EST-SSR markers in Bergenia ciliata using de novo transcriptome sequencing.

Bergenia ciliata (Haw.) Sternb. is an important herb predominantly found in the Indian Himalayan Region. It is widely used in medicines, healthcare systems, cosmetics, fodder, and ornamental purposes. The Illumina sequencing and de novo transcriptome assembly were carried out in B. ciliata to develop and identify simple sequence repeat markers. A total of 18 226 simple sequence repeats (SSRs) were identified wherein di-nucleotides were found to be abundant (47.88%), followed by mono-nucleotide (35.03%) and tri-nucleotide (15.88%) repeats. A total of 11 839 EST-SSR primers were designed, of which 96 primer pairs were commercially synthesized. Finally, 17 primer pairs revealed clear, distinct polymorphic bands, and these primers were validated with 40 diverse B. ciliata accessions. The present study revealed moderate level of genetic diversity (Ho=0.389, He=0.542, and PIC=0.513). Furthermore, the transcriptome data and EST-SSR markers generated during the present investigation could be an important genetic resource for functional genomics, population studies, and conservation genetics of the genus Bergenia.

Rice Genotypes and DNA Markers: A Review

This review provides an overview of the progress made in the field of rice genotypes and DNA markers. Rice is a crucial food source globally, and enhancing its nutritional value and resilience to stresses is of significant interest. The availability of high-quality rice genome sequences and functional genomic resources has accelerated genetic research, leading to the identification of genes that influence yield, grain quality, and stress tolerance. Rice genotypes are classified based on various factors, such as their response to salt stress, grain type, ionomic profile, arsenic accumulation, and appearance. DNA markers, including SSRs, SNPs, RAPD, RFLP, and AFLP, are used to study genetic variations and traits inheritance. DNA marker analysis has applications in studying genetic diversity, improving abiotic stress tolerance, developing salt-resistant germplasm, enhancing grain quality, developing resilient cultivars, and increasing crop yield. The review also includes case studies from Iraq, where DNA markers have been used to analyse genetic diversity and relationships among rice varieties. Overall, DNA markers play a crucial role in advancing rice research and breeding programs for improved productivity and food security.

When and how can we predict adaptive responses to climate change?

Predicting if, when, and how populations can adapt to climate change constitutes one of the greatest challenges in science today. Here, we build from contributions to the special issue on evolutionary adaptation to climate change, a survey of its authors, and recent literature to explore the limits and opportunities for predicting adaptive responses to climate change. We outline what might be predictable now, in the future, and perhaps never even with our best efforts. More accurate predictions are expected for traits characterized by a well-understood mapping between genotypes and phenotypes and traits experiencing strong, direct selection due to climate change. A meta-analysis revealed an overall moderate trait heritability and evolvability in studies performed under future climate conditions but indicated no significant change between current and future climate conditions, suggesting neither more nor less genetic variation for adapting to future climates. Predicting population persistence and evolutionary rescue remains uncertain, especially for the many species without sufficient ecological data. Still, when polled, authors contributing to this special issue were relatively optimistic about our ability to predict future evolutionary responses to climate change. Predictions will improve as we expand efforts to understand diverse organisms, their ecology, and their adaptive potential. Advancements in functional genomic resources, especially their extension to non-model species and the union of evolutionary experiments and "omics," should also enhance predictions. Although predicting evolutionary responses to climate change remains challenging, even small advances will reduce the substantial uncertainties surrounding future evolutionary responses to climate change.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution.

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20Mb long with a contig N50 length of 29.26Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

A large-scale whole-exome sequencing mutant resource for functional genomics in wheat.

Hexaploid wheat (Triticum aestivum), a major staple crop, has a remarkably large genome of ~14.4 Gb (containing 106 913 high-confidence [HC] and 159 840 low-confidence [LC] genes in the Chinese Spring v2.1 reference genome), which poses a major challenge for functional genomics studies. To overcome this hurdle, we performed whole-exome sequencing to generate a nearly saturated wheat mutant database containing 18 025 209 mutations induced by ethyl methanesulfonate (EMS), carbon (C)-ion beams, or γ-ray mutagenesis. This database contains an average of 47.1 mutations per kb in each gene-coding sequence: the potential functional mutations were predicted to cover 96.7% of HC genes and 70.5% of LC genes. Comparative analysis of mutations induced by EMS, γ-rays, or C-ion beam irradiation revealed that γ-ray and C-ion beam mutagenesis induced a more diverse array of variations than EMS, including large-fragment deletions, small insertions/deletions, and various non-synonymous single nucleotide polymorphisms. As a test case, we combined mutation analysis with phenotypic screening and rapidly mapped the candidate gene responsible for the phenotype of a yellow-green leaf mutant to a 2.8-Mb chromosomal region. Furthermore, a proof-of-concept reverse genetics study revealed that mutations in gibberellic acid biosynthesis and signalling genes could be associated with negative impacts on plant height. Finally, we built a publically available database of these mutations with the corresponding germplasm (seed stock) repository to facilitate advanced functional genomics studies in wheat for the broad plant research community.

The utility of transcriptomics in the conservation of sensitive and economically important species

The utility of transcriptomics in the conservation of sensitive and economically important species

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

Genome-scale targeted mutagenesis in Brassica napus using a pooled CRISPR library.

The recently constructed mutant libraries of diploid crops by the CRISPR-Cas9 system have provided abundant resources for functional genomics and crop breeding. However, because of the genome complexity, it is a big challenge to accomplish large-scale targeted mutagenesis in polyploid plants. Here, we demonstrate the feasibility of using a pooled CRISPR library to achieve genome-scale targeted editing in an allotetraploid crop of Brassica napus A total of 18,414 sgRNAs were designed to target 10,480 genes of interest, and afterward, 1104 regenerated transgenic plants harboring 1088 sgRNAs were obtained. Editing interrogation results revealed that 93 of the 178 genes were identified as mutated, thus representing an editing efficiency of 52.2%. Furthermore, we have discovered that Cas9-mediated DNA cleavages tend to occur at all the target sites guided by the same individual sgRNA, a novel finding in polyploid plants. Finally, we show the strong capability of reverse genetic screening for various traits with the postgenotyped plants. Several genes, which might dominate the fatty acid profile and seed oil content and have yet to be reported, were unveiled from the forward genetic studies. Our research provides valuable resources for functional genomics, elite crop breeding, and a good reference for high-throughput targeted mutagenesis in other polyploid plants.

Predicting brain-regional gene regulatory networks from multi-omics for Alzheimer's disease phenotypes and Covid-19 severity.

Neuroinflammation and immune dysregulation play a key role in Alzheimer's disease (AD) and are also associated with severe Covid-19 and neurological symptoms. Also, genome-wide association studies found many risk single nucleotide polymorphisms (SNPs) for AD and Covid-19. However, our understanding of underlying gene regulatory mechanisms from risk SNPs to AD, Covid-19 and phenotypes is still limited. To this end, we performed an integrative multi-omics analysis to predict gene regulatory networks for major brain regions from population data in AD. Our networks linked transcription factors (TFs) to TF binding sites (TFBSs) on regulatory elements to target genes. Comparative network analyses revealed cross-region-conserved and region-specific regulatory networks, in which many immunological genes are present. Furthermore, we identified a list of AD-Covid genes using our networks involving known and Covid-19 genes. Our machine learning analysis prioritized 36AD-Covid candidate genes for predicting Covid severity. Our independent validation analyses found that these genes outperform known genes for classifying Covid-19 severity and AD. Finally, we mapped genome-wide association study SNPs of AD and severe Covid that interrupt TFBSs on our regulatory networks, revealing potential mechanistic insights of those disease risk variants. Our analyses and results are open-source available, providing an AD-Covid functional genomic resource at the brain region level.

Chemical-Genetic Approaches for Exploring Mode of Action of Antifungal Compounds in the Fungal Pathogen Candida albicans.

Candida albicans is a prevalent fungal pathogen of humans that can cause both superficial and life-threatening disease, primarily in immunocompromised populations. Currently, antifungal drug classes available to treat fungal infections remain limited and the emergence of drug-resistant strains threatens antifungal efficacy, necessitating the discovery and development of additional therapeutics. The construction of the C. albicans double-barcoded heterozygous deletion collection (DBC) enables the rapid and systematic assessment of haploinsufficiency phenotypes in a pooled format. Specifically, this functional genomics resource can be used to identify heterozygous deletion mutants that are hypersensitive to compounds in order to define putative cellular targets and/or other modifiers of compound activity. Here, we describe protocols to characterize the mode of action of small molecules using the C. albicans DBC, including how to prepare compound-treated cultures, isolate genomic DNA, amplify strain-specific barcodes, and prepare DNA libraries for high-throughput sequencing. This technique provides a powerful approach to elucidate the compound mechanism of action in order to bolster the antifungal pipeline.

Utilizing Large Functional and Population Genomics Resources for CRISPR/Cas Perturbation Experiment Design.

Genome sequencing technologies have rapidly evolved in the past decades, enabling us to interpret the human genome through multiple perspectives, ranging from cross-species comparisons, naturally occurring variation in health and disease state to regulatory mechanisms.Although such perspectives are all informative to narrow down the list of genes or variants for perturbation experiments based on specific biological aims, utilizing multiple sources of information is often challenging in practice.In this chapter, we provide an overview of major large-scale functional and population genomics resources, followed by a practical example of selecting target variants for genetic perturbation experiments involving genome engineering techniques such as CRISPR/Cas.

Annotating genetic variants to target genes using H-MAGMA.

An outstanding goal in modern genomics is to systematically predict the functional outcome of noncoding variation associated with complex traits. To address this, we developed Hi-C-coupled multi-marker analysis of genomic annotation (H-MAGMA), which builds on traditional MAGMA-a gene-based analysis tool that assigns variants to their target genes-by incorporating 3D chromatin configuration in assigning variants to their putative target genes. Applying this approach, we identified key biological pathways implicated in a wide range of brain disorders and showed its utility in complementing other functional genomic resources such as expression quantitative trait loci-based variant annotation. Here, we provide a detailed protocol for generating the H-MAGMA variant-gene annotation file by using chromatin interaction data from the adult human brain. In addition, we provide an example of how H-MAGMA is run by using genome-wide association study summary statistics of Parkinson's disease. Lastly, we generated variant-gene annotation files for 28 tissues and cell types, with the hope of contributing a resource for researchers studying a broad set of complex genetic disorders. H-MAGMA can be performed in <2 h for any cell type in which Hi-C data are available.

Dynamics of alkannin/shikonin biosynthesis in response to jasmonate and salicylic acid in Lithospermum officinale

Alkannin/shikonin and their derivatives are specialised metabolites of high pharmaceutical and ecological importance exclusively produced in the periderm of members of the plant family Boraginaceae. Previous studies have shown that their biosynthesis is induced in response to methyl jasmonate but not salicylic acid, two phytohormones that play important roles in plant defence. However, mechanistic understanding of induction and non-induction remains largely unknown. In the present study, we generated the first comprehensive transcriptomic dataset and metabolite profiles of Lithospermum officinale plants treated with methyl jasmonate and salicylic acid to shed light on the underlying mechanisms. Our results highlight the diverse biological processes activated by both phytohormones and reveal the important regulatory role of the mevalonate pathway in alkannin/shikonin biosynthesis in L. officinale. Furthermore, by modelling a coexpression network, we uncovered structural and novel regulatory candidate genes connected to alkannin/shikonin biosynthesis. Besides providing new mechanistic insights into alkannin/shikonin biosynthesis, the generated methyl jasmonate and salicylic acid elicited expression profiles together with the coexpression networks serve as important functional genomic resources for the scientific community aiming at deepening the understanding of alkannin/shikonin biosynthesis.

Improved pea reference genome and pan-genome highlight genomic features and evolutionary characteristics

Complete and accurate reference genomes and annotations provide fundamental resources for functional genomics and crop breeding. Here we report a de novo assembly and annotation of a pea cultivar ZW6 with contig N50 of 8.98 Mb, which features a 243-fold increase in contig length and evident improvements in the continuity and quality of sequence in complex repeat regions compared with the existing one. Genome diversity of 118 cultivated and wild pea demonstrated that Pisum abyssinicum is a separate species different from P. fulvum and P. sativum within Pisum. Quantitative trait locus analyses uncovered two known Mendel’s genes related to stem length (Le/le) and seed shape (R/r) as well as some candidate genes for pod form studied by Mendel. A pan-genome of 116 pea accessions was constructed, and pan-genes preferred in P. abyssinicum and P. fulvum showed distinct functional enrichment, indicating the potential value of them as pea breeding resources in the future.

Genome sequencing reveals the evolution and pathogenic mechanisms of the wheat sharp eyespot pathogen Rhizoctonia cerealis

The necrotrophic fungus Rhizoctonia cerealis is the causal agent of devastating diseases of cereal crops including wheat (Triticum aestivum). We present a high-quality genome assembly of R. cerealis Rc207, a virulent strain causing wheat sharp eyespot. The assembly (56.36 Mb) is composed of 17.87% repeat sequences and 14,433 predicted protein-encoding genes. The Rc207 genome encodes a large and diverse set of genes involved in pathogenicity, especially rich in those encoding secreted proteins, carbohydrate-active enzymes (CAZymes), peptidases, nucleases, cytochrome P450, and secondary metabolism-associated enzymes. Most secretory protein-encoding genes, including CAZymes, peroxygenases, dehydrogenases, and cytochrome P450, were up-regulated during fungal infection of wheat. We identified 831 candidate secretory effectors and validated the functions of 10 up-regulated candidate effector proteins. Of them, nine were confirmed as necrotrophic pathogen’s effectors promoting fungal infection. Abundant potential mobile or plastic genomic regions rich in repeat sequences suggest their roles in fungal adaption and virulence-associated genomic evolution. This study provides valuable resources for further comparative and functional genomics on important fungal pathogens, and provides essential tools for development of effective disease control strategies.