Significant efforts have been made to characterize viral diversity in bats from China. Many of these studies were prospective and focused mainly on Rhinolophus bats that could be related to zoonotic events. However, other species of bats that are part of ecosystems identified as virus diversity hotspots have not been studied in-depth. We analyzed the virome of a group of Myotis fimbriatus bats collected from the Yunnan Province during 2020. The virome of M. fimbriatus revealed the presence of families of pathogenic viruses such as Coronavirus, Astrovirus, Mastadenovirus, and Picornavirus, among others. The viral sequences identified in M. fimbriatus were characterized by significant divergence from other known viral sequences of bat origin. Complex phylogenetic landscapes implying a tendency of co-specificity and relationships with viruses from other mammals characterize these groups. The most prevalent and abundant virus in M. fimbriatus individuals was an alphacoronavirus. The genome of this virus shows evidence of recombination and is likely the product of ancestral host-switch. The close phylogenetic and ecological relationship of some species of the Myotis genus in China may have played an important role in the emergence of this alphacoronavirus.

Since the identification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the etiological agent of the current COVID-19 pandemic, a rapid and massive effort has been made to obtain the genomic sequences of this virus to monitor (in near real time) the phylodynamic and diversity of this new pathogen. However, less attention has been given to the assessment of intra-host diversity. RNA viruses such as SARS-CoV-2 inhabit the host as a population of variants called quasispecies. We studied the quasispecies diversity in four of the main SARS-CoV-2 genes (ORF1a, ORF1b, S and N genes), using a dataset consisting of 210 next-generation sequencing (NGS) samples collected between January and early April of 2020 in the State of Victoria, Australia. We found evidence of quasispecies diversity in 68% of the samples, 76% of which was nonsynonymous variants with a higher density in the spike (S) glycoprotein and ORF1a genes. About one-third of the nonsynonymous intra-host variants were shared among the samples, suggesting host-to-host transmission. Quasispecies diversity changed over time. Phylogenetic analysis showed that some of the intra-host single-nucleotide variants (iSNVs) were restricted to specific lineages, highlighting their potential importance in the epidemiology of this virus. A greater effort must be made to determine the magnitude of the genetic bottleneck during transmission and the epidemiological and/or evolutionary factors that may play a role in the changes in the diversity of quasispecies over time.

Coconut (Cocos nucifera) is the emblematic palm of tropical coastal areas all around the globe. It provides vital resources to millions of farmers. In an effort to better understand its evolutionary history and to develop genomic tools for its improvement, a sequence draft was recently released. Here, we present a dense linkage map (8402 SNPs) aiming to assemble the large genome of coconut (2.42 Gbp, 2n=32) into 16 pseudomolecules. As a result, 47% of the sequences (representing 77% of the genes) were assigned to 16 linkage groups and ordered. We observed segregation distortion in chromosome Cn15, which is a signature of strong selection among pollen grains, favouring the maternal allele. Comparing our results with the genome of the oil palm Elaeis guineensis allowed us to identify major events in the evolutionary history of palms. We find that coconut underwent a massive transposable element invasion in the last million years, which could be related to the fluctuations of sea level during the glaciations at Pleistocene that would have triggered a population bottleneck. Finally, to better understand the facultative halophyte trait of coconut, we conducted an RNA-seq experiment on leaves to identify key players of signaling pathways involved in salt stress response. Altogether, our findings represent a valuable resource for the coconut breeding community. Yaodong Yang et al. report a chromosome-scale genome assembly of the coconut and transcriptome analysis of coconut leaves. Their analysis reveals important steps along the evolutionary trajectory of palms and identifies genes that may play a key role in tolerance to salt stress.

BACKGROUND:The ANRS12286/MOBIDIP trial showed that boosted protease inhibitor (bPI) plus lamivudine dual therapy was superior to bPI monotherapy as maintenance treatment in subjects with a history of M184V mutation. OBJECTIVES:We aimed to deep analyse the detection of M184V/I variants at time of switch and at the time of virological failure (VF). METHODS:Ultra-deep sequencing (UDS) was performed on proviral HIV-DNA at inclusion among 265 patients enrolled in the ANRS 12026/MOBIDIP trial, and on plasma from 31 patients experiencing VF. The proportion of M184V/I variants was described and the association between the M184V/I mutation at 1% of threshold and VF was explored with logistic regression models. RESULTS:M184V and I mutations were detected in HIV-DNA for 173/252 (69%) and 31/252 (12%) of participants, respectively. Longer duration of first-line treatment, higher plasma viral load at first-line treatment failure and higher baseline HIV-DNA load were associated with the archived M184V. M184I mutation was always associated with a STOP codon, suggesting defective virus. The 48 week estimated probability of remaining free from VF was comparable with or without the M184V/I mutation for dual therapy. At failure, M184V and major PI mutations were detected in 1/17 and 5/15 patients in the bPI arm and in 2/2 and 0/3 in the bPI+lamivudine arm, respectively. CONCLUSIONS:Using UDS evidenced that archiving of M184V in HIV-DNA is heterogeneous despite past historical M184V in 96% of cases. The antiviral efficacy of lamivudine-based dual therapy regimens is mainly due to the residual lamivudine activity.

Objectives: Human adenovirus (HAdV) infections are associated with a high morbidity and mortality in transplant patients requiring the use of antiviral treatments. Brincidofovir (BCV), a cytidine analog, inhibits HAdV replication through viral DNA elongation termination and likely through other mechanisms. To elucidate if BCV regulates cellular antiviral pathways, we analyzed its impact on HAdV-infected and non-HAdV-infected lung epithelial cells. Methods: We assessed the cellular and viral transcriptome of A549 cells infected and non-infected with HAdV C5 and treated or non-treated with BCV by RNAseq after 72 h. Results: BCV treatment of HAdV infected cells resulted in a profound decrease of viral transcription associated with a relative overexpression of the early genes E1A and E4 and of the late gene L1. BCV had also a profound impact on A549 cells' transcriptome. Ontologic analysis revealed an effect of BCV on several pathways known to interact with adenovirus replication as mTor signalling and Wnt pathways. A549 cells treated with BCV demonstrated a significant inhibition of the biological function of "viral replication" including 25 dysregulated genes involved in inflammation pathways. Conclusion: We demonstrated that BCV alters viral gene expression and promotes the expression of antiviral cellular pathways in A549 cells. These results provide new insights how to interfere with cellular pathways to control HAdV infections.

We used next-generation sequencing to evaluate the quantity and genetic diversity of the HIV envelope gene in various compartments in eight patients with acute infection. Plasma (PL) and seminal fluid (SF) were available for all patients, whole blood (WB) for seven, non-spermatozoid cells (NSC) for four, and saliva (SAL) for three. Median HIV-1 RNA was 6.2 log 10 copies/mL [IQR: 5.5–6.95] in PL, 4.9 log 10 copies/mL [IQR: 4.25–5.29] in SF, and 4.9 log 10 copies/mL [IQR: 4.46–5.09] in SAL. Median HIV-1 DNA was 4.1 log 10 copies/10 6 PBMCs [IQR: 3.15–4.15] in WB and 2.6 log 10 copies /10 6 Cells [IQR: 2.23–2.75] in NSC. The median overall diversity per patient varied from 0.0005 to 0.0232, suggesting very low diversity, confirmed by the clonal aspect of most of the phylogenetic trees. One single haplotype was present in all compartments for five patients in the earliest stage of infection. Evidence of higher diversity was established for two patients in PL and WB, suggesting compartmentalization. Our study shows low diversity of the env gene in the first stages of infection followed by the rapid establishment of cellular reservoirs of the virus. Such clonality could be exploited in the search for early patient-specific therapeutic solutions.

The myrtle rust disease, caused by the fungus Austropuccinia psidii, infects a wide range of host species within the Myrtaceae family worldwide. Since its first report in 2013 in New Caledonia, it was found on various types of native environments where Myrtaceae are the dominant or codominant species, as well as in several commercial nurseries. It is now considered as a significant threat to ecosystems biodiversity and Myrtaceae-related economy. The use of predictive molecular markers for resistance against myrtle rust is currently the most cost-effective and ecological approach to control the disease. Such an approach for neo Caledonian endemic Myrtaceae species was not possible because of the lack of genomic resources. The recent advancement in new generation sequencing technologies accompanied with relevant bioinformatics tools now provide new research opportunity for work in non-model organism at the transcriptomic level. The present study focuses on transcriptome analysis on three Myrtaceae species endemic to New Caledonia (Arillastrum gummiferum, Syzygium longifolium and Tristaniopsis glauca) that display contrasting responses to the pathogen (non-infected vs infected). Differential gene expression (DGE) and variant calling analysis were conducted on each species. We combined a dual approach by using 1) the annotated reference genome of a related Myrtaceae species (Eucalyptus grandis) and 2) a de novo transcriptomes of each species.

Background. A better understanding of HIV transmission dynamics among populations at high risk is important for development of prevention strategies. We determined HIV transmission networks from infected individuals enrolled in the pre-exposure prophylaxis (PrEP) IPERGAY trial in combination with the ANRS PRIMO and Montreal PHI cohorts to identify and characterize active clusters of transmission in this high-risk population. Methods. Genotypic resistance tests were performed on plasma samples from 31 IPERGAY participants. Reverse transcriptase sequences were analyzed in combination with unique HIV pol sequences from 1351 individuals enrolled in the PRIMO ANRS cohort (1999-2014) and 511 individuals enrolled in the Montreal PHI cohort (1996-2016). Network analyses were performed to infer putative relationships between all participants. Results. Overall, 1893 participants were included. Transmission network analyses revealed that 14 individuals (45.2%) from the IPERGAY trial were involved in 13 clusters sampled over a median period (interquartile range) of 2 (0.3-7.8) years, including 7 dyads and 6 larger clusters ranging from 4 to 28 individuals. When comparing characteristics between clustering individuals enrolled in the PRIMO cohort (n = 377) and in IPERGAY (n = 14), we found that IPERGAY participants had a higher viral load (5.93 vs 5.20 log 10 copies/mL, P = .032) and reported a higher number of partners in the last 2 months (P < .01). Conclusions. These results demonstrate high rates of HIV transmission clustering among young high-risk MSM enrolled in the IPERGAY trial. In-depth sampling of high-risk populations may help to uncover unobserved transmission intermediaries and improve prevention efforts that could be targeted to the most active clusters.

Coconut is an integral part of the livelihood of millions of farmers in the coastal tropical areas and is facing an important mutation. While its place as a major cheap oil crop on the international market is slowly declining, promising markets are emerging, such as coconut water, virgin coconut oil, and sugar among others. Coconut is regaining its tradition role of multi-usage crop, leading to new opportunities. It is however facing serious threats such as widespread lethal diseases and insect pests. To meet these challenges, coconut genetic improvement should renew its practices and objectives. But it is a difficult task due to long generation duration, low planting densities and low prolificacy. Advances in coconut genomic studies will improve its efficiency in several ways: Neutral markers allow broadening the genetic base of selection; QTLs based on mapping populations or on whole genome association studies (GWAS) will reduce the time and the areas needed to establish a breeding program. Comparative genomics and transcriptomics provide an in-depth understanding of metabolic pathways involved in production, product quality and adaptation to biotic and abiotic stress. A coconut genome draft was published recently and will be converted into a reference sequence thanks to high-density linkage mapping. Other sequencing efforts have been undertaken, whose combined results will provide a preliminary basis for characterizing coconut genetic diversity at the gene level to be completed by more resequencing. A revival of coconut genetic improvement will depend, among other elements, on high quality phenotyping, in conjunction with large SNP sets, transcriptomics and comparative genomics. (Resume d'auteur)

Le cocotier (Cocos nucifera L.) est une espece cotiere d'importance economique pour les zones tropicales, pour laquelle l'amelioration genetique se heurte a la faiblesse relative des ressources moleculaires a disposition des selectionneurs. La connaissance de son genome de ses variations peut beneficier des avancees en termes de ressources genomiques propre au cocotier, mais aussi des avancees sur le sequencage des genomes d'autres arecacees, dont le palmier a huile (Elaeis guinensis / oleifera) et le palmier dattier (Phoenix dactylifera), en lien avec les autres monocotyledones (recherche translationnelle). En retour cette connaissance peut contribuer a l'etude de l'histoire genetique des arecacees, notamment pour ce qui concerne l'evolution des familles et reseaux de genes d'interet economique comme la tolerance a la salinite. L'analyse bioinformatique developpee au cours de ces travaux a permis d'obtenir un proteome de reference du cocotier par l'amelioration de donnees transcriptomiques publiques en utilisant le genome du palmier a huile comme reference. Ces donnees ont permis une premiere analyse d'evolution de familles de genes impliquees notamment dans les stress biotiques. Accompagnant le developpement des premieres donnees de sequencage du genome du cocotier, l'etude comparative du cocotier et d'autres palmacees a permis de completer les hypotheses evolutives ayant conduit au genome actuel du cocotier et identifier des retentions de genes cles pour expliquer la biologie et l'adaptation de cette espece. Enfin, la disponibilite, sur deux genotypes de cocotier, de donnees transcriptomiques retracant leur reponse au stress salin a permis d'emettre des hypotheses sur ces mecanismes et leur variation au sein de l'espece.

The palms are a family of tropical origin and one of the main constituents of the ecosystems of these regions around the world. The two main species of palm represent different challenges: coconut (Cocos nucifera L.) is a source of multiple goods and services in tropical communities, while oil palm (Elaeis guineensis Jacq) is the main protagonist of the oil market. In this study, we present a workflow that exploits the comparative genomics between a target species (coconut) and a reference species (oil palm) to improve the transcriptomic data, providing a proteome useful to answer functional or evolutionary questions. This workflow reduces redundancy and fragmentation, two inherent problems of transcriptomic data, while preserving the functional representation of the target species. Our approach was validated in Arabidopsis thaliana using Arabidopsis lyrata and Capsella rubella as references species. This analysis showed the high sensitivity and specificity of our strategy, relatively independent of the reference proteome. The workflow increased the length of proteins products in A. thaliana by 13%, allowing, often, to recover 100% of the protein sequence length. In addition redundancy was reduced by a factor greater than 3. In coconut, the approach generated 29,366 proteins, 1,246 of these proteins deriving from new contigs obtained with the BRANCH software. The coconut proteome presented a functional profile similar to that observed in rice and an important number of metabolic pathways related to secondary metabolism. The new sequences found with BRANCH software were enriched in functions related to biotic stress. Our strategy can be used as a complementary step to de novo transcriptome assembly to get a representative proteome of a target species. The results of the current analysis are available on the website PalmComparomics (http://palm-comparomics.southgreen.fr/).

We describe here the reconstruction of the genome of the most recent common ancestor (MRCA) of modern monocots and eudicots, accounting for 95% of extant angiosperms, with its potential repertoire of 22,899 ancestral genes conserved in present-day crops. The MRCA provides a starting point for deciphering the reticulated evolutionary plasticity between species (rapidly versus slowly evolving lineages), subgenomes (pre- versus post-duplication blocks), genomic compartments (stable versus labile loci), genes (ancestral versus species-specific genes) and functions (gained versus lost ontologies), the key mutational forces driving the success of polyploidy in crops. The estimation of the timing of angiosperm evolution, based on MRCA genes, suggested that this group emerged 214 million years ago during the late Triassic era, before the oldest recorded fossil. Finally, the MRCA constitutes a unique resource for scientists to dissect major agronomic traits in translational genomics studies extending from model species to crops.

Abstract Coconut palm (Cocos nucifera,2n = 32), a member of genus Cocos and family Arecaceae (Palmaceae), is an important tropical fruit and oil crop. Currently, coconut palm is cultivated in 93 countries, including Central and South America, East and West Africa, Southeast Asia and the Pacific Islands, with a total growth area of more than 12 million hectares [1]. Coconut palm is generally classified into 2 main categories: “Tall” (flowering 8–10 years after planting) and “Dwarf” (flowering 4–6 years after planting), based on morphological characteristics and breeding habits. This Palmae species has a long growth period before reproductive years, which hinders conventional breeding progress. In spite of initial successes, improvements made by conventional breeding have been very slow. In the present study, we obtained de novo sequences of the Cocos nucifera genome: a major genomic resource that could be used to facilitate molecular breeding in Cocos nucifera and accelerate the breeding process in this important crop. A total of 419.67 gigabases (Gb) of raw reads were generated by the Illumina HiSeq 2000 platform using a series of paired-end and mate-pair libraries, covering the predicted Cocos nucifera genome length (2.42 Gb, variety “Hainan Tall”) to an estimated ×173.32 read depth. A total scaffold length of 2.20 Gb was generated (N50 = 418 Kb), representing 90.91% of the genome. The coconut genome was predicted to harbor 28 039 protein-coding genes, which is less than in Phoenix dactylifera (PDK30: 28 889), Phoenix dactylifera (DPV01: 41 660), and Elaeis guineensis (EG5: 34 802). BUSCO evaluation demonstrated that the obtained scaffold sequences covered 90.8% of the coconut genome and that the genome annotation was 74.1% complete. Genome annotation results revealed that 72.75% of the coconut genome consisted of transposable elements, of which long-terminal repeat retrotransposons elements (LTRs) accounted for the largest proportion (92.23%). Comparative analysis of the antiporter gene family and ion channel gene families between C. nucifera and Arabidopsis thaliana indicated that significant gene expansion may have occurred in the coconut involving Na+/H+ antiporter, carnitine/acylcarnitine translocase, potassium-dependent sodium-calcium exchanger, and potassium channel genes. Despite its agronomic importance, C. nucifera is still under-studied. In this report, we present a draft genome of C. nucifera and provide genomic information that will facilitate future functional genomics and molecular-assisted breeding in this crop species.

In 2015, the South Green Bioinformatics Platform http://www.southgreen.fr/ is a network of 35 bioinformaticians from five biology research institutes working with two High - Performance Computing Data Centres to develop and use new tools for NGS/ Omic analytics of tropical and Mediterranean crops under projects studying relationsh ip between genetic diversity, agronomic performance and response to selection. South Green is affiliated to the South regional centre of the French Institute of Bioinformatics (the French node of the European research infrastructure, ELIXIR). This communit y and the HPC data centres are all located in Montpellier, which facilitates close collaboration and significant pooling to best meet the biologists' demands of our research units. Since 2004, we developed web - based applications with both generic and in - ho use components, for databases, analysis workflows and web interfaces, in order to: manage genetic and phenotypic information ( e.g. TropGeneDB), analyse molecular markers and genetic diversity ( e.g. SNiPlay), assemble transcriptomes ( e.g. ESTtik) map RNA - Se q ( e.g. ARCAD), annotate and compare genomes ( e.g. GNPAnnot), reconstruct evolutionary history of gene families by phylogenomics ( e.g. GreenPhyl). We also participate to the analysis of numerous crop species, that requires computing and storage facilities as well as interoperable information systems, such as rice ( e.g. OryGenesDB), wheat, sorghum, sugarcane, banana (Banana Genome Hub), palms, yam, coffee (CGH), rubber, cacao (CocoaGenDB), cotton, apple, grapevine, olive, eucalyptus, cassava. To face the dat a deluge, we must increase our analytics capabilities. We document our operation at both, administrator/ developer and user/ scientist level, to provide high quality services and reproducible research. We pool into working groups on key themes such as GBS, at both, developer (extreme pair programming) and user (interdisciplinary knowledge exchange) level. We provide training sessions each year. Finally, we implemented several instances of the Galaxy workflow manager and encapsulated our tools. These instanc es serve as a catalyst for massive NGS analyses but it remains to increase storage capacity and improve data management plans. (Resume d'auteur)

Modern plant genomes are diploidized paleopolyploids. We revisited grass genome paleohistory in response to the diploidization process through a detailed investigation of the evolutionary fate of duplicated blocks. Ancestrally duplicated genes can be conserved, deleted, and shuffled, defining dominant (bias toward duplicate retention) and sensitive (bias toward duplicate erosion) chromosomal fragments. We propose a new grass genome paleohistory deriving from an ancestral karyotype structured in seven protochromosomes containing 16,464 protogenes and following evolutionary rules where 1) ancestral shared polyploidizations shaped conserved dominant (D) and sensitive (S) subgenomes, 2) subgenome dominance is revealed by both gene deletion and shuffling from the S blocks, 3) duplicate deletion/movement may have been mediated by single-/double-stranded illegitimate recombination mechanisms, 4) modern genomes arose through centromeric fusion of protochromosomes, leading to functional monocentric neochromosomes, 5) the fusion of two dominant blocks leads to supradominant neochromosomes (D + D = D) with higher ancestral gene retention compared with D + S = D (i.e., fusion of blocks with opposite sensitivity) or even S + S = S (i.e., fusion of two sensitive ancestral blocks). A new user-friendly online tool named “PlantSyntenyViewer,” available at http://urgi.versailles.inra.fr/synteny-cereal, presents the refined comparative genomics data.