The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

The human OPA1delTTAG mutation induces adult onset and progressive auditory neuropathy in mice.

Dominant optic atrophy (DOA) is one of the most prevalent forms of hereditary optic neuropathies and is mainly caused by heterozygous variants in OPA1, encoding a mitochondrial dynamin-related large GTPase. The clinical spectrum of DOA has been extended to a wide variety of syndromic presentations, called DOAplus, including deafness as the main secondary symptom associated to vision impairment. To date, the pathophysiological mechanisms underlying the deafness in DOA remain unknown. To gain insights into the process leading to hearing impairment, we have analyzed the Opa1delTTAG mouse model that recapitulates the DOAplus syndrome through complementary approaches combining morpho-physiology, biochemistry, and cellular and molecular biology. We found that Opa1delTTAG mutation leads an adult-onset progressive auditory neuropathy in mice, as attested by the auditory brainstem response threshold shift over time. However, the mutant mice harbored larger otoacoustic emissions in comparison to wild-type littermates, whereas the endocochlear potential, which is a proxy for the functional state of the stria vascularis, was comparable between both genotypes. Ultrastructural examination of the mutant mice revealed a selective loss of sensory inner hair cells, together with a progressive degeneration of the axons and myelin sheaths of the afferent terminals of the spiral ganglion neurons, supporting an auditory neuropathy spectrum disorder (ANSD). Molecular assessment of cochlea demonstrated a reduction of Opa1 mRNA level by greater than 40%, supporting haploinsufficiency as the disease mechanism. In addition, we evidenced an early increase in Sirtuin 3 level and in Beclin1 activity, and subsequently an age-related mtDNA depletion, increased oxidative stress, mitophagy as well as an impaired autophagic flux. Together, these results support a novel role for OPA1 in the maintenance of inner hair cells and auditory neural structures, addressing new challenges for the exploration and treatment of OPA1-linked ANSD in patients.

The evolutionary history of three Baracoffea species from western Madagascar revealed by chloroplast and nuclear genomes.

The wild species of the Coffea genus present a very wide morphological, genetic, and biochemical diversity. Wild species are recognized more resistant to diseases, pests, and environmental variations than the two species currently cultivated worldwide: C. arabica (Arabica) and C. canephora (Robusta). Consequently, wild species are now considered as a crucial resource for adapting cultivated coffee trees to climate change. Within the Coffea genus, 79 wild species are native to the Indian Ocean islands of Comoros, Mayotte, Mauritius, Réunion and Madagascar, out of a total of 141 taxa worldwide. Among them, a group of 9 species called "Baracoffea" are particularly atypical in their morphology and adaptation to the sandy soils of the dry deciduous forests of western Madagascar. Here, we have attempted to shed light on the evolutionary history of three Baracoffea species: C. ambongensis, C. boinensis and C. bissetiae by analyzing their chloroplast and nuclear genomes. We assembled the complete chloroplast genomes de novo and extracted 28,800 SNP (Single Nucleotide Polymorphism) markers from the nuclear genomes. These data were used for phylogenetic analysis of Baracoffea with Coffea species from Madagascar and Africa. Our new data support the monophyletic origin of Baracoffea within the Coffea of Madagascar, but also reveal a divergence with a sister clade of four species: C. augagneurii, C. ratsimamangae, C. pervilleana and C. Mcphersonii (also called C. vohemarensis), belonging to the Subterminal botanical series and living in dry or humid forests of northern Madagascar. Based on a bioclimatic analysis, our work suggests that Baracoffea may have diverged from a group of Malagasy Coffea from northern Madagascar and adapted to the specific dry climate and low rainfall of western Madagascar. The genomic data generated in the course of this work will contribute to the understanding of the adaptation mechanisms of these particularly singular species.

Genome mining of metabolic gene clusters in the Rubiaceae family.

The Rubiaceae plant family, comprising 3 subfamilies and over 13,000 species, is known for producing significant bioactive compounds such as caffeine and monoterpene indole alkaloids. Despite an increase in available genomes from the Rubiaceae family over the past decade, a systematic analysis of the metabolic gene clusters (MGCs) encoded by these genomes has been lacking. In this study, we aim to identify and analyze metabolic gene clusters within complete Rubiaceae genomes through a comparative analysis of eight species. Applying two bioinformatics pipelines, we identified 2372 candidate MGCs, organized into 549 gene cluster families (GCFs). To enhance the reliability of these findings, we developed coexpression networks and conducted orthology analyses. Using genomic data from Solanum lycopersicum (Solanaceae) for comparative purposes, we provided a detailed view of predicted metabolic enzymes, pathways, and coexpression networks. We bring some examples of MGCs and GCFs involved in biological pathways of terpenes, saccharides and alkaloids. Such insights lay the groundwork for discovering new compounds and associated MGCs within the Rubiaceae family, with potential implications in developing more robust crop species and expanding the understanding of plant metabolism. This large-scale exploration also provides a new perspective on the evolution and structure-function relationship of these clusters, offering opportunities for the highly efficient utilization of these unique metabolites. The outcome of this study contributes to a broader comprehension of the biosynthetic pathways, elucidating multiple aspects of specialized metabolism and offering innovative avenues for biotechnological applications.

Genomic object detection: An improved approach for transposable elements detection and classification using convolutional neural networks.

Analysis of eukaryotic genomes requires the detection and classification of transposable elements (TEs), a crucial but complex and time-consuming task. To improve the performance of tools that accomplish these tasks, Machine Learning approaches (ML) that leverage computer resources, such as GPUs (Graphical Processing Unit) and multiple CPU (Central Processing Unit) cores, have been adopted. However, until now, the use of ML techniques has mostly been limited to classification of TEs. Herein, a detection-classification strategy (named YORO) based on convolutional neural networks is adapted from computer vision (YOLO) to genomics. This approach enables the detection of genomic objects through the prediction of the position, length, and classification in large DNA sequences such as fully sequenced genomes. As a proof of concept, the internal protein-coding domains of LTR-retrotransposons are used to train the proposed neural network. Precision, recall, accuracy, F1-score, execution times and time ratios, as well as several graphical representations were used as metrics to measure performance. These promising results open the door for a new generation of Deep Learning tools for genomics. YORO architecture is available at https://github.com/simonorozcoarias/YORO.

Efficient homology-based annotation of transposable elements using minimizers.

Premise: Transposable elements (TEs) make up more than half of the genomes of complex plant species and can modulate the expression of neighboring genes, producing significant variability of agronomically relevant traits. The availability of long-read sequencing technologies allows the building of genome assemblies for plant species with large and complex genomes. Unfortunately, TE annotation currently represents a bottleneck in the annotation of genome assemblies. Methods and Results: We present a new functionality of the Next-Generation Sequencing Experience Platform (NGSEP) to perform efficient homology-based TE annotation. Sequences in a reference library are treated as long reads and mapped to an input genome assembly. A hierarchical annotation is then assigned by homology using the annotation of the reference library. We tested the performance of our algorithm on genome assemblies of different plant species, including Arabidopsis thaliana, Oryza sativa, Coffea humblotiana, and Triticum aestivum (bread wheat). Our algorithm outperforms traditional homology-based annotation tools in speed by a factor of three to >20, reducing the annotation time of the T. aestivum genome from months to hours, and recovering up to 80% of TEs annotated with RepeatMasker with a precision of up to 0.95. Conclusions: NGSEP allows rapid analysis of TEs, especially in very large and TE-rich plant genomes.

Metagenomics of African Empogona and Tricalysia (Rubiaceae) reveals the presence of leaf endophytes.

Background: Leaf symbiosis is a phenomenon in which host plants of Rubiaceae interact with bacterial endophytes within their leaves. To date, it has been found in around 650 species belonging to eight genera in four tribes; however, the true extent in Rubiaceae remains unknown. Our aim is to investigate the possible occurrence of leaf endophytes in the African plant genera Empogona and Tricalysia and, if present, to establish their identity. Methods: Total DNA was extracted from the leaves of four species of the Coffeeae tribe (Empogona congesta, Tricalysia hensii, T. lasiodelphys, and T. semidecidua) and sequenced. Bacterial reads were filtered out and assembled. Phylogenetic analysis of the endophytes was used to reveal their identity and their relationship with known symbionts. Results: All four species have non-nodulated leaf endophytes, which are identified as Caballeronia. The endophytes are distinct from each other but related to other nodulated and non-nodulated endophytes. An apparent phylogenetic or geographic pattern appears to be absent in endophytes or host plants. Caballeronia endophytes are present in the leaves of Empogona and Tricalysia, two genera not previously implicated in leaf symbiosis. This interaction is likely to be more widespread, and future discoveries are inevitable.

Selection signatures and population dynamics of transposable elements in lima bean.

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

InpactorDB: A Plant LTR Retrotransposon Reference Library.

LTR retrotransposons (LTR-RT) are major components of plant genomes. These transposable elements participate in the structure and evolution of genes and genomes through their mobility and their copy number amplification. For example, they are commonly used as evolutionary markers in genetic, genomic, and cytogenetic approaches. However, the plant research community is faced with the near absence of free availability of full-length, curated, and lineage-level classified LTR retrotransposon reference sequences. In this chapter, we will introduce InpactorDB, an LTR retrotransposon sequence database of 181 plant species representing 98 plant families for a total of 67,241 non-redundant elements. We will introduce how to use newly sequenced genomes to identify and classify LTR-RTs in a similar way with a standardized procedure using the Inpactor tool. InpactorDB is freely available at https://inpactordb.github.io .

Is periventricular heterotopia a useful endpoint for developmental thyroid hormone system disruption in mouse toxicity studies?

In rats, hypothyroidism during fetal and neonatal development can disrupt neuronal migration and induce the formation of periventricular heterotopia in the brain. However, it remains uncertain if heterotopia also manifest in mice after developmental hypothyroidism and whether they could be used as a toxicological endpoint to detect TH-mediated effects caused by TH system disrupting chemicals. Here, we performed a mouse study where we induced severe hypothyroidism by exposing pregnant mice (n = 3) to a very high dose of propylthiouracil (PTU) (1500 ppm) in the diet. This, to obtain best chances of detecting heterotopia. We found what appears to be very small heterotopia in 4 out of the 8 PTU-exposed pups. Although the incidence rate could suggest some utility for this endpoint, the small size of the ectopic neuronal clusters at maximum hypothyroidism excludes the utility of heterotopia in mouse toxicity studies aimed to detect TH system disrupting chemicals. On the other hand, parvalbumin expression was manifestly lower in the cortex of hypothyroid mouse offspring demonstrating that offspring TH-deficiency caused an effect on the developing brain. Based on overall results, we conclude that heterotopia formation in mice is not a useful toxicological endpoint for examining TH-mediated developmental neurotoxicity.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes.

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

HERV-K (HML-2) insertion polymorphisms in the 8q24.13 region and their potential etiological associations with acute myeloid leukemia.

Human endogenous retroviruses (HERVs) are LTR retrotransposons that are present in the human genome. Among them, members of the HERV-K (HML-2) group are suspected to play a role in the development of different types of cancer, including lung, ovarian, and prostate cancer, as well as leukemia. Acute myeloid leukemia (AML) is an important disease that causes 1% of cancer deaths in the United States and has a survival rate of 28.7%. Here, we describe a method for assessing the statistical association between HERV-K (HML-2) transposable element insertion polymorphisms (or TIPs) and AML, using whole-genome sequencing and read mapping using TIP_finder software. Our results suggest that 101 polymorphisms involving HERV-K (HML-2) elements were correlated with AML, with a percentage between 44.4 to 56.6%, most of which (70) were located in the region from 8q24.13 to 8q24.21. Moreover, it was found that the TRIB1, LRATD2, POU5F1B, MYC, PCAT1, PVT1, and CCDC26 genes could be displaced or fragmented by TIPs. Furthermore, a general method was devised to facilitate analysis of the correlation between transposable element insertions and specific diseases. Finally, although the relationship between HERV-K (HML-2) TIPs and AML remains unclear, the data reported in this study indicate a statistical correlation, as supported by the χ2 test with p-values < 0.05.

Neuronal Blockade of Thyroid Hormone Signaling Increases Sensitivity to Diet-Induced Obesity in Adult Male Mice.

Thyroid hormone increases energy expenditure. Its action is mediated by TR, nuclear receptors present in peripheral tissues and in the central nervous system, particularly in hypothalamic neurons. Here, we address the importance of thyroid hormone signaling in neurons, in general for the regulation of energy expenditure. We generated mice devoid of functional TR in neurons using the Cre/LoxP system. In hypothalamus, which is the center for metabolic regulation, mutations were present in 20% to 42% of the neurons. Phenotyping was performed under physiological conditions that trigger adaptive thermogenesis: cold and high-fat diet (HFD) feeding. Mutant mice displayed impaired thermogenic potential in brown and inguinal white adipose tissues and were more prone to diet-induced obesity. They showed a decreased energy expenditure on chow diet and gained more weight on HFD. This higher sensitivity to obesity disappeared at thermoneutrality. Concomitantly, the AMPK pathway was activated in the ventromedial hypothalamus of the mutants as compared with the controls. In agreement, sympathetic nervous system (SNS) output, visualized by tyrosine hydroxylase expression, was lower in the brown adipose tissue of the mutants. In contrast, absence of TR signaling in the mutants did not affect their ability to respond to cold exposure. This study provides the first genetic evidence that thyroid hormone signaling exerts a significant influence in neurons to stimulate energy expenditure in some physiological context of adaptive thermogenesis. TR function in neurons to limit weight gain in response to HFD and this effect is associated with a potentiation of SNS output.

High nucleotide similarity of three Copia lineage LTR retrotransposons among plant genomes.

Transposable elements (TEs) are mobile elements found in the majority of eukaryotic genomes. TEs deeply impact the structure and evolution of chromosomes and can induce mutations affecting coding genes. In plants, the major group of TEs is long terminal repeat retrotransposons (LTR-RTs). They are classified into superfamilies (Gypsy, Copia) and subclassified into lineages. Horizontal transfer (HT), defined as the nonsexual transmission of genetic material between species, is a process allowing LTR-RTs to invade a new genome. Although this phenomenon was considered rare, recent studies demonstrate numerous transfers of LTR-RTs. This study aims to determine which LTR-RT lineages are shared with high similarity among 69 plant genomes. We identified and classified 88 450 LTR-RTs and determined 143 cases of high similarities between pairs of genomes. Most of them involved three Copia lineages (Oryco/Ivana, Retrofit/Ale, and Tork/Tar/Ikeros). A detailed analysis of three cases of high similarities involving Tork/Tar/Ikeros group shows an uneven distribution in the phylogeny of the elements and incongruence with between phylogenetic trees topologies, indicating they could be originated from HTs. Overall, our results suggest that LTR-RT Copia lineages share outstanding similarity between distant species and may likely be involved in HT mechanisms more frequent than initially estimated.

G-SAIP: Graphical Sequence Alignment Through Parallel Programming in the Post-Genomic Era.

A common task in bioinformatics is to compare DNA sequences to identify similarities between organisms at the sequence level. An approach to such comparison is the dot-plots, a 2-dimensional graphical representation to analyze DNA or protein alignments. Dot-plots alignment software existed before the sequencing revolution, and now there is an ongoing limitation when dealing with large-size sequences, resulting in very long execution times. High-Performance Computing (HPC) techniques have been successfully used in many applications to reduce computing times, but so far, very few applications for graphical sequence alignment using HPC have been reported. Here, we present G-SAIP (Graphical Sequence Alignment in Parallel), a software capable of spawning multiple distributed processes on CPUs, over a supercomputing infrastructure to speed up the execution time for dot-plot generation up to 1.68× compared with other current fastest tools, improve the efficiency for comparative structural genomic analysis, phylogenetics because the benefits of pairwise alignments for comparison between genomes, repetitive structure identification, and assembly quality checking.

Inpactor2: a software based on deep learning to identify and classify LTR-retrotransposons in plant genomes.

LTR-retrotransposons are the most abundant repeat sequences in plant genomes and play an important role in evolution and biodiversity. Their characterization is of great importance to understand their dynamics. However, the identification and classification of these elements remains a challenge today. Moreover, current software can be relatively slow (from hours to days), sometimes involve a lot of manual work and do not reach satisfactory levels in terms of precision and sensitivity. Here we present Inpactor2, an accurate and fast application that creates LTR-retrotransposon reference libraries in a very short time. Inpactor2 takes an assembled genome as input and follows a hybrid approach (deep learning and structure-based) to detect elements, filter partial sequences and finally classify intact sequences into superfamilies and, as very few tools do, into lineages. This tool takes advantage of multi-core and GPU architectures to decrease execution times. Using the rice genome, Inpactor2 showed a run time of 5 minutes (faster than other tools) and has the best accuracy and F1-Score of the tools tested here, also having the second best accuracy and specificity only surpassed by EDTA, but achieving 28% higher sensitivity. For large genomes, Inpactor2 is up to seven times faster than other available bioinformatics tools.

Brown adipocytes local response to thyroid hormone is required for adaptive thermogenesis in adult male mice.

Thyroid hormone (T3) and its nuclear receptors (TR) are important regulators of energy expenditure and adaptive thermogenesis, notably through their action in the brown adipose tissue (BAT). However, T3 acts in many other peripheral and central tissues which are also involved in energy expenditure. The general picture of how T3 regulates BAT thermogenesis is currently not fully established, notably due to the absence of extensive omics analyses and the lack of specific mice model. Here, we first used transcriptome and cistrome analyses to establish the list of T3/TR direct target genes in brown adipocytes. We then developed a novel model of transgenic mice, in which T3 signaling is specifically suppressed in brown adipocytes at adult stage. We addressed the capacity of these mice to mount a thermogenic response when challenged by either a cold exposure or a high-fat diet, and analyzed the associated changes in BAT transcriptome. We conclude that T3 plays a crucial role in the thermogenic response of the BAT, controlling the expression of genes involved in lipid and glucose metabolism and regulating BAT proliferation. The resulting picture provides an unprecedented view on the pathways by which T3 activates energy expenditure through an efficient adaptive thermogenesis in the BAT.

An Atlas of Thyroid Hormone Receptors' Target Genes in Mouse Tissues.

We gathered available RNA-seq and ChIP-seq data in a single database to better characterize the target genes of thyroid hormone receptors in several cell types. This database can serve as a resource to analyze the mode of action of thyroid hormone (T3). Additionally, it is an easy-to-use and convenient tool to obtain information on specific genes regarding T3 regulation or to extract large gene lists of interest according to the users' criteria. Overall, this atlas is a unique compilation of recent sequencing data focusing on T3, its receptors, modes of action, targets and roles, which may benefit researchers within the field. A preliminary analysis indicates extensive variations in the repertoire of target genes where transcription is upregulated by chromatin-bound nuclear receptors. Although it has a major influence, chromatin accessibility is not the only parameter that determines the cellular selectivity of the hormonal response.

Bactris gasipaes Kunth var. gasipaes complete plastome and phylogenetic analysis.

Bactris gasipaes var. gasipaes (Arecaceae, Palmae) is an economically and socially important plant species for populations across tropical South and Central America. It has been domesticated from its wild variety, B. gasipaes var. chichagui, since pre-Columbian times. In this study, we sequenced the plastome of the cultivated variety, B. gasipaes Kunth var. gasipaes and compared it with the published plastome of the wild variety. The chloroplast sequence obtained was 156,580 bp. The cultivated chloroplast sequence was conserved compared to the wild type sequence with 99.8% of nucleotide identity. We did, however, identify multiple Single Nucleotide Variants (SNVs), insertions, microsatellites and a resolved region of missing nucleotides. A SNV in one of the core barcode markers (matK) was detected between the wild and cultivated accessions. Phylogenetic analysis was carried out across the Arecaceae family and compared to previous reports, resulting in an identical topology. This study is a step forward in understanding the genome evolution of this species.