A revisited history of cacao domestication in pre-Columbian times revealed by archaeogenomic approaches.

Humans have a long history of transporting and trading plants, contributing to the evolution of domesticated plants. Theobroma cacao originated in the Neotropics from South America. However, little is known about its domestication and use in these regions. In this study, ceramic residues from a large sample of pre-Columbian cultures from South and Central America were analyzed using archaeogenomic and biochemical approaches. Here we show, for the first time, the widespread use of cacao in South America out of its native Amazonian area of origin, extending back 5000 years, likely supported by cultural interactions between the Amazon and the Pacific coast. We observed that strong genetic mixing between geographically distant cacao populations occurred as early as the middle Holocene, in South America, driven by humans, favoring the adaptation of T. cacao to new environments. This complex history of cacao domestication is the basis of today's cacao tree populations and its knowledge can help us better manage their genetic resources.

Target enrichment sequencing coupled with GWAS identifies MdPRX10 as a candidate gene in the control of budbreak in apple.

The timing of floral budbreak in apple has a significant effect on fruit production and quality. Budbreak occurs as a result of a complex molecular mechanism that relies on accurate integration of external environmental cues, principally temperature. In the pursuit of understanding this mechanism, especially with respect to aiding adaptation to climate change, a QTL at the top of linkage group (LG) 9 has been identified by many studies on budbreak, but the genes underlying it remain elusive. Here, together with a dessert apple core collection of 239 cultivars, we used a targeted capture sequencing approach to increase SNP resolution in apple orthologues of known or suspected A. thaliana flowering time-related genes, as well as approximately 200 genes within the LG9 QTL interval. This increased the 275 223 SNP Axiom® Apple 480 K array dataset by an additional 40 857 markers. Robust GWAS analyses identified MdPRX10, a peroxidase superfamily gene, as a strong candidate that demonstrated a dormancy-related expression pattern and down-regulation in response to chilling. In-silico analyses also predicted the residue change resulting from the SNP allele associated with late budbreak could alter protein conformation and likely function. Late budbreak cultivars homozygous for this SNP allele also showed significantly up-regulated expression of C-REPEAT BINDING FACTOR (CBF) genes, which are involved in cold tolerance and perception, compared to reference cultivars, such as Gala. Taken together, these results indicate a role for MdPRX10 in budbreak, potentially via redox-mediated signaling and CBF gene regulation. Moving forward, this provides a focus for developing our understanding of the effects of temperature on flowering time and how redox processes may influence integration of external cues in dormancy pathways.

Custom-made meniscus biofabrication.

Reconstructing the meniscus is not currently possible due to its intricate and heterogeneous structure. In this forum, we first discuss the shortcomings of current clinical strategies in meniscus repair. Then, we describe a new promising cell-based, ink-free 3D biofabrication technology to produce tailor-made large-scale functional menisci.

DNA metabarcoding suggests adaptive seasonal variation of individual trophic traits in a critically endangered fish.

Dietary studies are critical for understanding foraging strategies and have important applications in conservation and habitat management. We applied a robust metabarcoding protocol to characterize the diet of the critically endangered freshwater fish Zingel asper (the Rhone streber). We conducted modelling and simulation analyses to identify and characterize some of the drivers of individual trophic trait variation in this species. We found that population density and ontogeny had minor effects on the trophic niche of Z. asper. Instead, our results suggest that the majority of trophic niche variation was driven by seasonal variation in ecological opportunity. The total trophic niche width of Z. asper seasonally expanded to include a broader range of prey. Furthermore, null model simulations revealed that the increase of between-individual variation in autumn indicates that Z. asper become more opportunistic relative to summer and spring, rather than being associated with a seasonal specialization of individuals. Overall, our results suggest an adaptive variation of individual trophic traits in Z. asper: the species mainly consumes a few ephemeropteran taxa (Baetis fuscatus and Ecdyonurus) but seems to be capable of adapting its foraging strategy to maintain its body condition. This study illustrates how metabarcoding data obtained from faeces can be validated and combined with individual-based modelling and simulation approaches to explore inter- and intrapopulational individual trophic traits variation and to test hypotheses in the conventional analytic framework of trophic ecology.

Identification and validation of QTLs for tuber quality related traits in greater yam (Dioscorea alata L.).

Two Dioscorea alata populations were generated by hand pollination between contrasted diploid genitors. Population A (74F × Kabusa) was composed of 121 progenies while population B (74F × 14M) involved 193 progenies. These two populations were assessed over two consecutive years regarding important tuber quality traits. Analysis of variance showed that the genotype had the greatest influence on the phenotypic scores. Also for some traits, effect of the year_replicate was strong. The heritabilities of most traits were high. Based on these data and a reference high-density genetic map of greater yam, a total of 34 quantitative trait loci (QTLs) were detected on 8 of the 20 yam chromosomes. They corresponded to five of each of the following traits: tuber size, shape regularity, tubercular roots, skin texture, tuber flesh oxidation, six for oxidation ratio and three for flesh colour. The fraction of total phenotypic variance attributable to a single QTL ranged from 11.1 to 43.5%. We detected significant correlations between traits and QTL colocalizations that were consistent with these correlations. A majority of QTLs (62%) were found on linkage group LG16, indicating that this chromosome could play a major role in genetic control of the investigated traits. In addition, an inversion involving this chromosome was detected in the Kabusa male. Nine QTLs were validated on a diversity panel, including three for tuber size, three for shape regularity, two for skin texture and one for tubercular roots. The approximate physical localization of validated QTLs allowed the identification of various candidates genes. The validated QTLs should be useful for breeding programs using marker-assisted selection to improve yam tuber quality.

Mechanical control of tissue shape: Cell-extrinsic and -intrinsic mechanisms join forces to regulate morphogenesis.

During vertebrate development, cells must proliferate, move, and differentiate to form complex shapes. Elucidating the mechanisms underlying the molecular and cellular processes involved in tissue morphogenesis is essential to understanding developmental programmes. Mechanical stimuli act as a major contributor of morphogenetic processes and impact on cell behaviours to regulate tissue shape and size. Specifically, cell extrinsic physical forces are translated into biochemical signals within cells, through the process of mechanotransduction, activating multiple mechanosensitive pathways and defining cell behaviours. Physical forces generated by tissue mechanics and the extracellular matrix are crucial to orchestrate tissue patterning and cell fate specification. At the cell scale, the actomyosin network generates the cellular tension behind the tissue mechanics involved in building tissue. Thus, understanding the role of physical forces during morphogenetic processes requires the consideration of the contribution of cell intrinsic and cell extrinsic influences. The recent development of multidisciplinary approaches, as well as major advances in genetics, microscopy, and force-probing tools, have been key to push this field forward. With this review, we aim to discuss recent work on how tissue shape can be controlled by mechanical forces by focusing specifically on vertebrate organogenesis. We consider the influences of mechanical forces by discussing the cell-intrinsic forces (such as cell tension and proliferation) and cell-extrinsic forces (such as substrate stiffness and flow forces). We review recently described processes supporting the role of intratissue force generation and propagation in the context of shape emergence. Lastly, we discuss the emerging role of tissue-scale changes in tissue material properties, extrinsic forces, and shear stress on shape establishment.

Extracellular mechanical forces drive endocardial cell volume decrease during zebrafish cardiac valve morphogenesis.

Organ morphogenesis involves dynamic changes of tissue properties while cells adapt to their mechanical environment through mechanosensitive pathways. How mechanical cues influence cell behaviors during morphogenesis remains unclear. Here, we studied the formation of the zebrafish atrioventricular canal (AVC) where cardiac valves develop. We show that the AVC forms within a zone of tissue convergence associated with the increased activation of the actomyosin meshwork and cell-orientation changes. We demonstrate that tissue convergence occurs with a reduction of cell volume triggered by mechanical forces and the mechanosensitive channel TRPP2/TRPV4. Finally, we show that the extracellular matrix component hyaluronic acid controls cell volume changes. Together, our data suggest that multiple force-sensitive signaling pathways converge to modulate cell volume. We conclude that cell volume reduction is a key cellular feature activated by mechanotransduction during cardiovascular morphogenesis. This work further identifies how mechanical forces and extracellular matrix influence tissue remodeling in developing organs.

Genetic control of flowering in greater yam (Dioscorea alata L.).

BACKGROUND: Greater yam (Dioscorea alata L.) is a major tropical and subtropical staple crop cultivated for its starchy tubers. Breeding of this dioecious species is hampered by its erratic flowering, yet little is currently known on the genetic determinism of its sexual reproduction. RESULT: Here we used a genome-wide association approach and identified a major genetic barrier to reproduction in yam on chromosome 1, as represented by two candidate genes. A deleterious effect on male fitness could be hypothesized considering the involvement of these two genes in male reproduction and the low frequency of this non-flowering dominant allele within the male genepool. We also extended the hypothesis of a XX/XY sex-determination system located on chromosome 6 in D. alata to encompass most of the species diversity. Moreover, a kompetitive allele-specific PCR (KASPar) marker was designed and validated that enables accurate cultivar sex estimation. The reconstruction of chromosome 6 associated with the detection of highly putative structural variations confirmed the possible involvement of a major part of the chromosome. CONCLUSION: The findings of this study, combined with proper estimation of accession ploidy levels to avoid endosperm incompatibility issues, could facilitate the design of future promising parental combinations in D. alata breeding programs. Moreover, the discovery of this genetic barrier to reproduction opens new avenues for gaining insight into yam reproductive biology and diversification.

Development of Microsatellite Markers Using Next-Generation Sequencing.

Among the molecular markers used for plant genetic studies, microsatellite markers are easy to implement and can provide suitable codominant markers for molecular taxonomy.Here we describe a method to obtain microsatellite primers from genomic DNA using a next-generation sequencer.

Isolation and characterization of microsatellite markers for the threatened African endemic tree species Pterocarpus erinaceus Poir.

To study the genetic diversity and structure of the forest species Pterocarpus erinaceus Poir., seventeen polymorphic nuclear microsatellite markers were isolated and characterized, using next-generation sequencing. Three hundred and sixty-five (365) individuals were analyzed within fifteen (15) West African populations. The number of alleles for these loci varied from 4 to 30, and the heterozygosity varied from 0.23 to 0.82. The seventeen (17) primers designed here will allow characterizing the genetic diversity of this threaten species on its natural stands and to better understand the population differentiation mechanisms shaping it.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis.

BACKGROUND: Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea. METHODS: Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20-22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators. RESULTS: Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

Mechanotransduction in cardiovascular morphogenesis and tissue engineering.

Cardiovascular morphogenesis involves cell behavior and cell identity changes that are activated by mechanical forces associated with heart function. Recently, advances in in vivo imaging, methods to alter blood flow, and computational modelling have greatly advanced our understanding of how forces produced by heart contraction and blood flow impact different morphogenetic processes. Meanwhile, traditional genetic approaches have helped to elucidate how endothelial cells respond to forces at the cellular and molecular level. Here we discuss the principles of endothelial mechanosensitity and their interplay with cellular processes during cardiovascular morphogenesis. We then discuss their implications in the field of cardiovascular tissue engineering.

Mechanically activated piezo channels modulate outflow tract valve development through the Yap1 and Klf2-Notch signaling axis.

Mechanical forces are well known for modulating heart valve developmental programs. Yet, it is still unclear how genetic programs and mechanosensation interact during heart valve development. Here, we assessed the mechanosensitive pathways involved during zebrafish outflow tract (OFT) valve development in vivo. Our results show that the hippo effector Yap1, Klf2, and the Notch signaling pathway are all essential for OFT valve morphogenesis in response to mechanical forces, albeit active in different cell layers. Furthermore, we show that Piezo and TRP mechanosensitive channels are important factors modulating these pathways. In addition, live reporters reveal that Piezo controls Klf2 and Notch activity in the endothelium and Yap1 localization in the smooth muscle progenitors to coordinate OFT valve morphogenesis. Together, this work identifies a unique morphogenetic program during OFT valve formation and places Piezo as a central modulator of the cell response to forces in this process.

Transcriptome data from three endemic Myrtaceae species from New Caledonia displaying contrasting responses to myrtle rust (Austropuccinia psidii).

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.

A reference high-density genetic map of greater yam (Dioscorea alata L.).

KEY MESSAGE: This study generated the first high-density genetic map for D. alata based on genotyping-by-sequencing and provides new insight on sex determination in yam. Greater yam (Dioscorea alata L.) is a major staple food in tropical and subtropical areas. This study aimed to produce the first reference genetic map of this dioecious species using genotyping-by-sequencing. In this high-density map combining information of two F1 outcrossed populations, 20 linkage groups were resolved as expected and 1579 polymorphic markers were ordered. The consensus map length was 2613.5 cM with an average SNP interval of 1.68 cM. An XX/XY sex determination system was identified on LG6 via the study of sex ratio, homology of parental linkage groups and the identification of a major QTL for sex determination. Homology with the sequenced D. rotundata is described, and the median physical distance between SNPs was estimated at 139.1 kb. The effects of segregation distortion and the presence of heteromorphic sex chromosomes are discussed. This D. alata linkage map associated with the available genomic resources will facilitate quantitative trait mapping, marker-assisted selection and evolutionary studies in the important yet scarcely studied yam species.

The use and domestication of Theobroma cacao during the mid-Holocene in the upper Amazon.

Cacao (Theobroma cacao L.) is an important economic crop, yet studies of its domestication history and early uses are limited. Traditionally, cacao is thought to have been first domesticated in Mesoamerica. However, genomic research shows that T. cacao's greatest diversity is in the upper Amazon region of northwest South America, pointing to this region as its centre of origin. Here, we report cacao use identified by three independent lines of archaeological evidence-cacao starch grains, absorbed theobromine residues and ancient DNA-dating from approximately 5,300 years ago recovered from the Santa Ana-La Florida (SALF) site in southeast Ecuador. To our knowledge, these findings constitute the earliest evidence of T. cacao use in the Americas and the first unequivocal archaeological example of its pre-Columbian use in South America. They also reveal the upper Amazon region as the oldest centre of cacao domestication yet identified.

Genotyping by sequencing reveals the interspecific C. maxima / C. reticulata admixture along the genomes of modern citrus varieties of mandarins, tangors, tangelos, orangelos and grapefruits.

The mandarin horticultural group is an important component of world citrus production for the fresh fruit market. This group formerly classified as C. reticulata is highly polymorphic and recent molecular studies have suggested that numerous cultivated mandarins were introgressed by C. maxima (the pummelos). C. maxima and C. reticulata are also the ancestors of sweet and sour oranges, grapefruit, and therefore of all the "small citrus" modern varieties (mandarins, tangors, tangelos) derived from sexual hybridization between these horticultural groups. Recently, NGS technologies have greatly modified how plant evolution and genomic structure are analyzed, moving from phylogenetics to phylogenomics. The objective of this work was to develop a workflow for phylogenomic inference from Genotyping By Sequencing (GBS) data and to analyze the interspecific admixture along the nine citrus chromosomes for horticultural groups and recent varieties resulting from the combination of the C. reticulata and C. maxima gene pools. A GBS library was established from 55 citrus varieties, using the ApekI restriction enzyme and selective PCR to improve the read depth. Diagnostic polymorphisms (DPs) of C. reticulata/C. maxima differentiation were identified and used to decipher the phylogenomic structure of the 55 varieties. The GBS approach was powerful and revealed 30,289 SNPs and 8,794 Indels with 12.6% of missing data. 11,133 DPs were selected covering the nine chromosomes with a higher density in genic regions. GBS combined with the detection of DPs was powerful for deciphering the "phylogenomic karyotypes" of cultivars derived from admixture of the two ancestral species after a limited number of interspecific recombinations. All the mandarins, mandarin hybrids, tangelos and tangors analyzed displayed introgression of C. maxima in different parts of the genome. C. reticulata/C. maxima admixture should be a major component of the high phenotypic variability of this germplasm opening up the way for association studies based on phylogenomics.

16S rRNA Amplicon Sequencing for Epidemiological Surveys of Bacteria in Wildlife.

The human impact on natural habitats is increasing the complexity of human-wildlife interactions and leading to the emergence of infectious diseases worldwide. Highly successful synanthropic wildlife species, such as rodents, will undoubtedly play an increasingly important role in transmitting zoonotic diseases. We investigated the potential for recent developments in 16S rRNA amplicon sequencing to facilitate the multiplexing of the large numbers of samples needed to improve our understanding of the risk of zoonotic disease transmission posed by urban rodents in West Africa. In addition to listing pathogenic bacteria in wild populations, as in other high-throughput sequencing (HTS) studies, our approach can estimate essential parameters for studies of zoonotic risk, such as prevalence and patterns of coinfection within individual hosts. However, the estimation of these parameters requires cleaning of the raw data to mitigate the biases generated by HTS methods. We present here an extensive review of these biases and of their consequences, and we propose a comprehensive trimming strategy for managing these biases. We demonstrated the application of this strategy using 711 commensal rodents, including 208 Mus musculus domesticus, 189 Rattus rattus, 93 Mastomys natalensis, and 221 Mastomys erythroleucus, collected from 24 villages in Senegal. Seven major genera of pathogenic bacteria were detected in their spleens: Borrelia, Bartonella, Mycoplasma, Ehrlichia, Rickettsia, Streptobacillus, and Orientia. Mycoplasma, Ehrlichia, Rickettsia, Streptobacillus, and Orientia have never before been detected in West African rodents. Bacterial prevalence ranged from 0% to 90% of individuals per site, depending on the bacterial taxon, rodent species, and site considered, and 26% of rodents displayed coinfection. The 16S rRNA amplicon sequencing strategy presented here has the advantage over other molecular surveillance tools of dealing with a large spectrum of bacterial pathogens without requiring assumptions about their presence in the samples. This approach is therefore particularly suitable to continuous pathogen surveillance in the context of disease-monitoring programs. IMPORTANCE Several recent public health crises have shown that the surveillance of zoonotic agents in wildlife is important to prevent pandemic risks. High-throughput sequencing (HTS) technologies are potentially useful for this surveillance, but rigorous experimental processes are required for the use of these effective tools in such epidemiological contexts. In particular, HTS introduces biases into the raw data set that might lead to incorrect interpretations. We describe here a procedure for cleaning data before estimating reliable biological parameters, such as positivity, prevalence, and coinfection, using 16S rRNA amplicon sequencing on an Illumina MiSeq platform. This procedure, applied to 711 rodents collected in West Africa, detected several zoonotic bacterial species, including some at high prevalence, despite their never before having been reported for West Africa. In the future, this approach could be adapted for the monitoring of other microbes such as protists, fungi, and even viruses.

Characterization of microsatellites identified by next-generation sequencing in the Neotropical tree Handroanthus billbergii (Bignoniaceae).

PREMISE OF THE STUDY: We developed microsatellite (simple sequence repeat [SSR]) markers in the Neotropical tree Handroanthus billbergii (Bignoniaceae), to be applied in assessment of genetic diversity in this species as a reference for inferring the impact of dry forest fragmentation in Ecuador. METHODS AND RESULTS: Using next-generation sequencing, we detected a total of 26,893 putative SSRs reported here. Using an ABI 3500xl sequencer, we identified and characterized a set of polymorphic markers in 23 individuals belonging to three populations of H. billbergii. CONCLUSIONS: We report a set of 30 useful SSR markers for H. billbergii and a large list of potential microsatellites for developing new markers for this or related species.

Development of microsatellite-enriched libraries.

Among the molecular markers used for plant genetic studies, microsatellite markers are easy to implement and can provide suitable codominant markers for molecular taxonomy. Here we describe a method to obtain enriched libraries in microsatellite loci from genomic DNA, using capture method with synthetic oligonucleotide probes.