Genome annotation: From human genetics to biodiversity genomics

Scaling metabolic model reconstruction up to the pan-genome level: A systematic review and prospective applications to photosynthetic organisms

Advances in genomics technologies have generated large data sets that provide tremendous insights into the genetic diversity of taxonomic groups. However, it remains challenging to pinpoint the effect of genetic diversity on different traits without performing resource-intensive phenotyping experiments. Pan-genome-scale metabolic models (panGEMs) extend traditional genome-scale metabolic models by considering the entire reaction repertoire that enables the prediction and comparison of metabolic capabilities within a taxonomic group. Here, we systematically review the state-of-the-art methodologies for constructing panGEMs, focusing on used tools, databases, experimental datasets, and orthology relationships. We highlight the unique advantages of panGEMs compared to single-species GEMs in predicting metabolic phenotypes and in guiding the experimental validation of genome annotations. In addition, we emphasize the disparity between the available (pan-)genomic data on photosynthetic organisms and their under-representation in current (pan)GEMs. Finally, we propose a perspective for tackling the reconstruction of panGEMs for photosynthetic eukaryotes that can help advance our understanding of the metabolic diversity in this taxonomic group.

De novo whole-genome assembly of the critically endangered southern muriqui (Brachyteles arachnoides)

The southern muriqui (Brachyteles arachnoides) is one of the 2 species of muriquis (genus Brachyteles), the largest body-sized nonhuman primate from the Neotropics. Deforestation and illegal hunting have led to a continuing decline in the muriqui population, leading to their current classification as critically endangered. The lack of a reference genome for the genus Brachyteles prevents scientists from taking full advantage of genomic tools to improve their conservation status. This study reports the first whole-genome assemblies of the genus Brachyteles, using DNA from 2 zoo-housed southern muriqui females. We performed sequencing with Oxford Nanopore Technologies' PromethION 2 Solo using a native DNA library preparation to preserve DNA modifications. We used Flye to assemble genomes for each individual. The best final assembly was 2.6 Gb, in 319 contigs, with an N50 of 58.8 Mb and an L50 of 17. BUSCO completeness score for this assembly was 99.5%. The assembly of the second individual had similar quality, with a length of 2.6 Gb, 759 contigs, an N50 of 47.9 Mb, an L50 of 18, and a BUSCO completeness score of 99.04%. Both assemblies had <1% duplicates, missing, or fragments. Gene model mapper detected 24,353 protein-coding genes, and repetitive elements accounted for 46% of the genome. We also reported the mitogenome, which had 16,562 bp over 37 genes, and global methylation of CpG sites, which revealed a mean of 80% methylation. Our study provides a high-quality reference genome assembly for the southern muriqui, expanding the tools that can be used to aid in their conservation efforts.

Long-read transcriptomics of a diverse human cohort reveals widespread ancestry bias in gene annotation

Accurate gene annotations are fundamental for interpreting genetic variation, cellular function, and disease mechanisms. However, current human gene annotations are largely derived from transcriptomic data of individuals with European ancestry, introducing potential biases that remain uncharacterized. Here, we generate over 800 million full-length reads with long-read RNA-seq in 43 lymphoblastoid cell line samples from eight genetically-diverse human populations and build a cross-ancestry gene annotation. We show that transcripts from non-European samples are underrepresented in reference gene annotations, leading to systematic biases in allele-specific transcript usage analyses. Furthermore, we show that personal genome assemblies enhance transcript discovery compared to the generic GRCh38 reference assembly, even though genomic regions unique to each individual are heavily depleted of genes. These findings underscore the urgent need for a more inclusive gene annotation framework that accurately represents global transcriptome diversity.

The BioGenome Portal: a web-based platform for biodiversity genomics data management

Biodiversity genomics projects are underway with the aim of sequencing the genomes of all eukaryotic species on Earth. Here we describe the BioGenome Portal, a web-based application to facilitate organization and access to the data produced by biodiversity genomics projects. The portal integrates user-generated data with data deposited in public repositories. The portal generates sequence status reports that can be eventually ingested by designated metadata tracking systems, facilitating the coordination task of these systems. The portal is open-source and fully customizable. It can be deployed at any site with minimum effort, contributing to the democratization of biodiversity genomics projects. We illustrate the features of the BioGenome Portal through a number of specific instances. One such instance is being used as the reference portal for the Catalan Initiative for the Earth Biogenome Project, a regional project aiming to sequencing the genomes of the species of the Catalan linguistic area.

Hookworm genes encoding intestinal excreted-secreted proteins are transcriptionally upregulated in response to the host's immune system

Hookworms are intestinal parasitic nematodes that chronically infect ~500 million people, with reinfection common even after clearance by drugs. How infecting hookworms successfully overcome host protective mechanisms is unclear, but it may involve hookworm proteins that digest host tissues, or counteract the host's immune system, or both. To find such proteins in the zoonotic hookworm Ancylostoma ceylanicum, we identified hookworm genes encoding excreted-secreted (ES) proteins, hookworm genes preferentially expressed in the hookworm intestine, and hookworm genes whose transcription is stimulated by the host immune system. We collected ES proteins from adult hookworms harvested from hamsters; mass spectrometry identified 565 A. ceylanicum genes encoding ES proteins. We also used RNA-seq to identify A. ceylanicum genes expressed both in young adults (12 days post-infection) and in intestinal and non-intestinal tissues dissected from mature adults (19 days post-infection), with hamster hosts that either had normal immune systems or were immunosuppressed by dexamethasone. In adult A. ceylanicum, we observed 1,670 and 1,196 genes with intestine- and non-intestine-biased expression, respectively. Comparing hookworm gene activity in normal versus immunosuppressed hosts, we observed almost no changes of gene activity in 12-day young adults or non-intestinal 19-day adult tissues. However, in intestinal 19-day adult tissues, we observed 1,951 positively immunoregulated genes (upregulated at least two-fold in normal hosts versus immunosuppressed hosts), and 137 genes that were negatively immunoregulated. Thus, immunoregulation was observed primarily in mature adult hookworm intestine directly exposed to host blood; it may include hookworm genes activated in response to the host immune system in order to neutralize the host immune system. We observed 153 ES genes showing positive immunoregulation in 19-day adult intestine; of these genes, 69 had ES gene homologs in the closely related hookworm Ancylostoma caninum, 24 in the human hookworm Necator americanus, and 24 in the more distantly related strongylid parasite Haemonchus contortus. Such a mixture of rapidly evolving and conserved genes could comprise virulence factors enabling infection, provide new targets for drugs or vaccines against hookworm, and aid in developing therapies for autoimmune diseases.

GENCODE: massively expanding the lncRNA catalog through capture long-read RNA sequencing

Accurate and complete gene annotations are indispensable for understanding how genome sequences encode biological functions. For twenty years, the GENCODE consortium has developed reference annotations for the human and mouse genomes, becoming a foundation for biomedical and genomics communities worldwide. Nevertheless, collections of important yet poorly-understood gene classes like long non-coding RNAs (lncRNAs) remain incomplete and scattered across multiple, uncoordinated catalogs, slowing down progress in the field. To address these issues, GENCODE has undertaken the most comprehensive lncRNAs annotation effort to date. This is founded on the manual annotation of full-length targeted long-read sequencing, on matched embryonic and adult tissues, of orthologous regions in human and mouse. Altogether 17,931 novel human genes (140,268 novel transcripts) and 22,784 novel mouse genes (136,169 novel transcripts) have been added to the GENCODE catalog representing a 2-fold and 6-fold increase in transcripts, respectively - the greatest increase since the sequencing of the human genome. Novel gene annotations display evolutionary constraints, have well-formed promoter regions, and link to phenotype-associated genetic variants. They greatly enhance the functional interpretability of the human genome, as they help explain millions of previously-mapped "orphan" omics measurements corresponding to transcription start sites, chromatin modifications and transcription factor binding sites. Crucially, our targeted design assigned human-mouse orthologs at a rate beyond previous studies, tripling the number of human disease-associated lncRNAs with mouse orthologs. The expanded and enhanced GENCODE lncRNA annotations mark a critical step towards deciphering the human and mouse genomes.

Quest for Orthologs in the Era of Biodiversity Genomics

The era of biodiversity genomics is characterized by large-scale genome sequencing efforts that aim to represent each living taxon with an assembled genome. Generating knowledge from this wealth of data has not kept up with this pace. We here discuss major challenges to integrating these novel genomes into a comprehensive functional and evolutionary network spanning the tree of life. In summary, the expanding datasets create a need for scalable gene annotation methods. To trace gene function across species, new methods must seek to increase the resolution of ortholog analyses, e.g. by extending analyses to the protein domain level and by accounting for alternative splicing. Additionally, the scope of orthology prediction should be pushed beyond well-investigated proteomes. This demands the development of specialized methods for the identification of orthologs to short proteins and noncoding RNAs and for the functional characterization of novel gene families. Furthermore, protein structures predicted by machine learning are now readily available, but this new information is yet to be integrated with orthology-based analyses. Finally, an increasing focus should be placed on making orthology assignments adhere to the findable, accessible, interoperable, and reusable (FAIR) principles. This fosters green bioinformatics by avoiding redundant computations and helps integrating diverse scientific communities sharing the need for comparative genetics and genomics information. It should also help with communicating orthology-related concepts in a format that is accessible to the public, to counteract existing misinformation about evolution.

Molecular responses of seaweeds to biotic interactions: A systematic review

Seaweed farming is the single largest aquaculture commodity with >30 million tonnes produced each year. Furthermore, the restoration of lost seaweed forests is gaining significant momentum, particularly for kelps in warming temperate areas. Whether in aquaculture settings, following restoration practices, or in the wild, all seaweeds undergo biotic interactions with a diverse range of co-occurring or cocultured organisms. To date, most research assessing such biotic interactions has focused on the response of the organism interacting with seaweeds, rather than on the seaweeds themselves. However, understanding how seaweeds respond to other organisms, particularly on a molecular scale, is crucial for optimizing outcomes of seaweed farming or restoration efforts and, potentially, also for the conservation of natural populations. In this systematic review, we assessed the molecular processes that seaweeds undergo during biotic interactions and propose priority areas for future research. Despite some insights into the response of seaweeds to biotic interactions, this review specifically highlights a lack of characterization of biomolecules involved in the response to chemical cues derived from interacting organisms (four studies in the last 20 years) and a predominant use of laboratory-based experiments conducted under controlled conditions. Additionally, this review reveals that studies targeting metabolites (70%) are more common than those examining the role of genes (22%) and proteins (8%). To effectively inform seaweed aquaculture efforts, it will be crucial to conduct larger scale experiments simulating natural environments. Also, employing a holistic approach targeting genes and proteins would be beneficial to complement the relatively well-established role of metabolites.

The Catalan initiative for the Earth BioGenome Project: contributing local data to global biodiversity genomics

The Catalan Initiative for the Earth BioGenome Project (CBP) is an EBP-affiliated project network aimed at sequencing the genome of the >40 000 eukaryotic species estimated to live in the Catalan-speaking territories (Catalan Linguistic Area, CLA). These territories represent a biodiversity hotspot. While covering less than 1% of Europe, they are home to about one fourth of all known European eukaryotic species. These include a high proportion of endemisms, many of which are threatened. This trend is likely to get worse as the effects of global change are expected to be particularly severe across the Mediterranean Basin, particularly in freshwater ecosystems and mountain areas. Following the EBP model, the CBP is a networked organization that has been able to engage many scientific and non-scientific partners. In the pilot phase, the genomes of 52 species are being sequenced. As a case study in biodiversity conservation, we highlight the genome of the Balearic shearwater Puffinus mauretanicus, sequenced under the CBP umbrella.

The promise and challenges of genomics-informed periodontal disease diagnoses

Recent advances in human genomics and the advent of molecular medicine have catapulted our ability to characterize human and health and disease. Scientists and healthcare practitioners can now leverage information on genetic variation and gene expression at the tissue or even individual cell level, and an enormous potential exists to refine diagnostic categories, assess risk in unaffected individuals, and optimize disease management among those affected. This review investigates the progress made in the domains of molecular medicine and genomics as they relate to periodontology. The review summarizes the current evidence of association between genomics and periodontal diseases, including the current state of knowledge that approximately a third of the population variance of periodontitis may be attributable to genetic variation and the management of several monogenic forms of the disease can be augmented by knowledge of the underlying genetic cause. Finally, the paper discusses the potential utility of polygenic risk scores and genetic testing for periodontitis diagnosis now and in the future, in light of applications that currently exist in other areas of medicine and healthcare.

BRD2 and BRD3 genes independently evolved RNA structures to control unproductive splicing

The mammalian BRD2 and BRD3 genes encode structurally related proteins from the bromodomain and extraterminal domain protein family. The expression of BRD2 is regulated by unproductive splicing upon inclusion of exon 3b, which is located in the region encoding a bromodomain. Bioinformatic analysis indicated that BRD2 exon 3b inclusion is controlled by a pair of conserved complementary regions (PCCR) located in the flanking introns. Furthermore, we identified a highly conserved element encoding a cryptic poison exon 5b and a previously unknown PCCR in the intron between exons 5 and 6 of BRD3, however, outside of the homologous bromodomain. Minigene mutagenesis and blockage of RNA structure by antisense oligonucleotides demonstrated that RNA structure controls the rate of inclusion of poison exons. The patterns of BRD2 and BRD3 expression and splicing show downregulation upon inclusion of poison exons, which become skipped in response to transcription elongation slowdown, further confirming a role of PCCRs in unproductive splicing regulation. We conclude that BRD2 and BRD3 independently acquired poison exons and RNA structures to dynamically control unproductive splicing. This study describes a convergent evolution of regulatory unproductive splicing mechanisms in these genes, providing implications for selective modulation of their expression in therapeutic applications.