Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequencing data annotation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Predicting promoters in phage genomes using PhagePromoter

Abstract Summary The growing interest in phages as antibacterial agents has led to an increase in the number of sequenced phage genomes, increasing the need for intuitive bioinformatics tools for performing genome annotation. The identification of phage promoters is indeed the most difficult step of this process. Due to the lack of online tools for phage promoter prediction, we developed PhagePromoter, a tool for locating promoters in phage genomes, using machine learning methods. This is the first online tool for predicting promoters that uses phage promoter data and the first to identify both host and phage promoters with different motifs. Availability and implementation This tool was integrated in the Galaxy framework and it is available online at: https://bit.ly/2Dfebfv. Supplementary information Supplementary data are available at Bioinformatics online.

Sampaio, Marta↗
DOE OSTI2022

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

Abstract Background Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity. Results COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer ( Ostrinia nubilalis ). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity. Conclusions The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Deeplasmid: deep learning accurately separates plasmids from bacterial chromosomes

Plasmids are mobile genetic elements that play a key role in microbial ecology and evolution by mediating horizontal transfer of important genes, such as antimicrobial resistance genes. Many microbial genomes have been sequenced by short read sequencers and have resulted in a mix of contigs that derive from plasmids or chromosomes. New tools that accurately identify plasmids are needed to elucidate new plasmid-borne genes of high biological importance. We have developed Deeplasmid, a deep learning tool for distinguishing plasmids from bacterial chromosomes based on the DNA sequence and its encoded biological data. It requires as input only assembled sequences generated by any sequencing platform and assembly algorithm and its runtime scales linearly with the number of assembled sequences. Deeplasmid achieves an AUC–ROC of over 89%, and it was more accurate than five other plasmid classification methods. Finally, as a proof of concept, we used Deeplasmid to predict new plasmids in the fish pathogen Yersinia ruckeri ATCC 29473 that has no annotated plasmids. Deeplasmid predicted with high reliability that a long assembled contig is part of a plasmid. Using long read sequencing we indeed validated the existence of a 102 kb long plasmid, demonstrating Deeplasmid's ability to detect novel plasmids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

FUSARIUM-ID v.3.0: An Updated, Downloadable Resource for Fusarium Species Identification

Species within Fusarium are of global agricultural, medical, and food/feed safety concern and have been extensively characterized. However, accurate identification of species is challenging and usually requires DNA sequence data. FUSARIUM-ID ( http://isolate.fusariumdb.org/blast.php ) is a publicly available database designed to support the identification of Fusarium species using sequences of multiple phylogenetically informative loci, especially the highly informative ~680-bp 5' portion of the translation elongation factor 1-alpha (TEF1) gene that has been adopted as the primary barcoding locus in the genus. However, FUSARIUM-ID v.1.0 and 2.0 had several limitations, including inconsistent metadata annotation for the archived sequences and poor representation of some species complexes and marker loci. Here, we present FUSARIUM-ID v.3.0, which provides the following improvements: (i) additional and updated annotation of metadata for isolates associated with each sequence, (ii) expanded taxon representation in the TEF1 sequence database, (iii) availability of the sequence database as a downloadable file to enable local BLAST queries, and (iv) a tutorial file for users to perform local BLAST searches using either freely available software, such as SequenceServer, BLAST+ executable in the command line, and Galaxy, or the proprietary Geneious software. FUSARIUM-ID will be updated on a regular basis by archiving sequences of TEF1 and other loci from newly identified species and greater in-depth sampling of currently recognized species.

Plant Sciences↗
DOE OSTI2016

Backbone chemical shift assignments and secondary structure analysis of the U1 protein from the Bas-Congo virus

The Bas-Congo virus (BASV) is the first rhabdovirus associated with a human outbreak of acute hemorrhagic fever. The single-stranded, negative-sense RNA genome of BASV contains the five core genes present in all rhabdoviral genomes plus an additional three genes, annotated U1, U2, and U3, with weak (<21%) sequence similarity only to a handful of genes observed in a few other rhabdoviral genomes. The function of the rhabdoviral U proteins is unknown, but, they are hypothesized to play a role in viral infection or replication. To better understand this unique family of proteins, a construct containing residues 27–203 of the 216-residue U1 protein (BASV-U1*) was prepared. By collecting data in 0.5 M urea it was possible to eliminate transient association enough to enable the assignment of most of the observable 1 H N , 1 Hα, 15 N, 13 Cα, 13 Cβ, and 13 C´ chemical shifts for BASV-U1* that will provide a foundation to study its solution properties. Here, the analyses of these chemical shifts along with 15 N-edited NOESY data enabled the identification of the elements of secondary structure present in BASV-U1*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Genome Sequence and Analysis of the Flavinogenic Yeast Candida membranifaciens IST 626

The ascomycetous yeast Candida membranifaciens has been isolated from diverse habitats, including humans, insects, and environmental sources, exhibiting a remarkable ability to use different carbon sources that include pentoses, melibiose, and inulin. In this study, we isolated four C. membranifaciens strains from soil and investigated their potential to overproduce riboflavin. C. membranifaciens IST 626 was found to produce the highest concentrations of riboflavin. The volumetric production of this vitamin was higher when C. membranifaciens IST 626 cells were cultured in a commercial medium without iron and when xylose was the available carbon source compared to the same basal medium with glucose. Supplementation of the growth medium with 2 g/L glycine favored the metabolization of xylose, leading to biomass increase and consequent enhancement of riboflavin volumetric production that reached 120 mg/L after 216 h of cultivation. To gain new insights into the molecular basis of riboflavin production and carbon source utilization in this species, the first annotated genome sequence of C. membranifaciens is reported in this article, as well as the result of a comparative genomic analysis with other relevant yeast species. A total of 5619 genes were predicted to be present in C. membranifaciens IST 626 genome sequence (11.5 Mbp). Among them are genes involved in riboflavin biosynthesis, iron homeostasis, and sugar uptake and metabolism. This work put forward C. membranifaciens IST 626 as a riboflavin overproducer and provides valuable molecular data for future development of superior producing strains capable of using the wide range of carbon sources, which is a characteristic trait of the species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A large sequenced mutant library – valuable reverse genetic resource that covers 98% of sorghum genes

SUMMARY Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum , a drought and heat‐tolerant C4 species, requires a vast, large‐scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large‐scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)‐induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis ‐regulatory elements. A comparison of the distribution of mutations in the large‐scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS‐induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5′ UTR, gene body, and 3′‐UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the K a / K s value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB ( https://www.depts.ttu.edu/igcast/sorbmutdb.php ) and SorghumBase ( https://sorghumbase.org/ ). This current large‐scale sequence‐indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Beyond microbial abundance: metadata integration enhances disease prediction in human microbiome studies

Multiple studies have highlighted the interaction of the human microbiome with physiological systems such as the gut, immune, liver, and skin, via key axes. Advances in sequencing technologies and high-performance computing have enabled the analysis of large-scale metagenomic data, facilitating the use of machine learning to predict disease likelihood from microbiome profiles. However, challenges such as compositionality, high dimensionality, sparsity, and limited sample sizes have hindered the development of actionable models. One strategy to improve these models is by incorporating key metadata from both the human host and sample collection/processing protocols. This remains challenging due to sparsity and inconsistency in metadata annotation and availability. In this paper, we introduce a machine learning-based pipeline for predicting human disease states by integrating host and protocol metadata with microbiome abundance profiles from 68 different studies, processed through a consistent pipeline. Our findings indicate that metadata can enhance machine learning predictions, particularly at higher taxonomic ranks like Kingdom and Phylum, though this effect diminishes at lower ranks. Our study leverages a large collection of microbiome datasets comprising 11,208 samples, therefore enhancing the robustness and statistical confidence of our findings. This work is a critical step toward utilizing microbiome and metadata for predicting diseases such as gastrointestinal infections, diabetes, cancer, and neurological disorders.

Mathematics and Computing↗
DOE OSTI2018

Integrative analysis of large scale transcriptome data draws a comprehensive landscape of Phaeodactylum tricornutum genome and evolutionary origin of diatoms

Diatoms are one of the most successful and ecologically important groups of eukaryotic phytoplankton in the modern ocean. Deciphering their genomes is a key step towards better understanding of their biological innovations, evolutionary origins, and ecological underpinnings. Here, we have used 90 RNA-Seq datasets from different growth conditions combined with published expressed sequence tags and protein sequences from multiple taxa to explore the genome of the model diatom Phaeodactylum tricornutum, and introduce 1,489 novel genes. The new annotation additionally permitted the discovery of extensive alternative splicing in diatoms, including intron retention and exon skipping, which increase the diversity of transcripts generated in changing environments. In addition, we have used up-to-date reference sequence libraries to dissect the taxonomic origins of diatom genes. We show that the P. tricornutum genome is enriched in lineage-specific genes, with up to 47% of the gene models present only possessing orthologues in other stramenopile groups. Finally, we have performed a comprehensive de novo annotation of repetitive elements showing novel classes of transposable elements such as SINE, MITE and TRIM/LARD. This work provides a solid foundation for future studies of diatom gene function, evolution and ecology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis -regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center (https://danio-code.zfin.org) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Methane utilization in Methylomicrobium alcaliphilum 20Z R : a systems approach

Biological methane utilization, one of the main sinks of the greenhouse gas in nature, represents an attractive platform for production of fuels and value-added chemicals. Despite the progress made in our understanding of the individual parts of methane utilization, our knowledge of how the whole-cell metabolic network is organized and coordinated is limited. Attractive growth and methane-conversion rates, a complete and expert-annotated genome sequence, as well as large enzymatic, 13 C-labeling, and transcriptomic datasets make Methylomicrobium alcaliphilum 20Z R an exceptional model system for investigating methane utilization networks. Here we present a comprehensive metabolic framework of methane and methanol utilization in M. alcaliphilum 20Z R . A set of novel metabolic reactions governing carbon distribution across central pathways in methanotrophic bacteria was predicted by in-silico simulations and confirmed by global non-targeted metabolomics and enzymatic evidences. Our data highlight the importance of substitution of ATP-linked steps with PPi-dependent reactions and support the presence of a carbon shunt from acetyl-CoA to the pentose-phosphate pathway and highly branched TCA cycle. The diverged TCA reactions promote balance between anabolic reactions and redox demands. As a result, the computational framework of C 1 -metabolism in methanotrophic bacteria can represent an efficient tool for metabolic engineering or ecosystem modeling.

09 BIOMASS FUELS↗
DOE OSTI2016

Cleavage and polyadenylation specificity factor 30: An RNA-binding zinc-finger protein with an unexpected 2Fe–2S cluster

Cleavage and polyadenylation specificity factor 30 (CPSF30) is a key protein involved in pre-mRNA processing. CPSF30 contains five Cys 3 His domains (annotated as “zinc-finger” domains). Using inductively coupled plasma mass spectrometry, X-ray absorption spectroscopy, and UV-visible spectroscopy, we report that CPSF30 is isolated with iron, in addition to zinc. Iron is present in CPSF30 as a 2Fe–2S cluster and uses one of the Cys 3 His domains; 2Fe–2S clusters with a Cys 3 His ligand set are rare and notably have also been identified in MitoNEET, a protein that was also annotated as a zinc finger. These findings support a role for iron in some zinc-finger proteins. Using electrophoretic mobility shift assays and fluorescence anisotropy, we report that CPSF30 selectively recognizes the AU-rich hexamer (AAUAAA) sequence present in pre-mRNA, providing the first molecular-based evidence to our knowledge for CPSF30/RNA binding. Removal of zinc, or both zinc and iron, abrogates binding, whereas removal of just iron significantly lessens binding. As a result, from these data we propose a model for RNA recognition that involves a metal-dependent cooperative binding mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A chromosome-level genome assembly of the Chinese cork oak (Quercus variabilis)

Quercus variabilis (Fagaceae) is an ecologically and economically important deciduous broadleaved tree species native to and widespread in East Asia. It is a valuable woody species and an indicator of local forest health, and occupies a dominant position in forest ecosystems in East Asia. However, genomic resources from Q. variabilis are still lacking. Here, we present a high-quality Q. variabilis genome generated by PacBio HiFi and Hi-C sequencing. The assembled genome size is 787 Mb, with a contig N50 of 26.04 Mb and scaffold N50 of 64.86 Mb, comprising 12 pseudo-chromosomes. The repetitive sequences constitute 67.6% of the genome, of which the majority are long terminal repeats, accounting for 46.62% of the genome. We used ab initio , RNA sequence-based and homology-based predictions to identify protein-coding genes. A total of 32,466 protein-coding genes were identified, of which 95.11% could be functionally annotated. Evolutionary analysis showed that Q. variabilis was more closely related to Q. suber than to Q. lobata or Q. robur. We found no evidence for species-specific whole genome duplications in Quercus after the species had diverged. This study provides the first genome assembly and the first gene annotation data for Q. variabilis. These resources will inform the design of further breeding strategies, and will be valuable in the study of genome editing and comparative genomics in oak species.

Han, Biao↗
DOE OSTI2023

Genome features of a novel hydrocarbonoclastic Chryseobacterium oranimense strain and its comparison to bacterial oil-degraders and to other C. oranimense strains

For the first time, we report the whole genome sequence of a hydrocarbonoclastic Chryseobacterium oranimense strain isolated from Trinidad and Tobago (COTT) and its genes involved in the biotransformation of hydrocarbons and xenobiotics through functional annotation. The assembly consisted of 11 contigs with 2,794 predicted protein-coding genes which included a diverse group of gene families involved in aliphatic and polycyclic hydrocarbon degradation. Comparative genomic analyses with 18 crude-oil degrading bacteria in addition to two C. oranimense strains not associated with oil were carried out. The data revealed important differences in terms of annotated genes involved in the hydrocarbon degradation process that may explain the molecular mechanisms of hydrocarbon and xenobiotic biotransformation. Notably, many gene families were expanded to explain COTT’s competitive ability to manage habitat-specific stressors. Gene-based evidence of the metabolic potential of COTT supports the application of indigenous microbes for the remediation of polluted terrestrial environments and provides a genomic resource for improving our understanding of how to optimize these characteristics for more effective bioremediation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Metagenome-assembled genomes provide insight into the metabolic potential during early production of Hydraulic Fracturing Test Site 2 in the Delaware Basin

Demand for natural gas continues to climb in the United States, having reached a record monthly high of 104.9 billion cubic feet per day (Bcf/d) in November 2023. Hydraulic fracturing, a technique used to extract natural gas and oil from deep underground reservoirs, involves injecting large volumes of fluid, proppant, and chemical additives into shale units. This is followed by a “shut-in” period, during which the fracture fluid remains pressurized in the well for several weeks. The microbial processes that occur within the reservoir during this shut-in period are not well understood; yet, these reactions may significantly impact the structural integrity and overall recovery of oil and gas from the well. To shed light on this critical phase, we conducted an analysis of both pre-shut-in material alongside production fluid collected throughout the initial production phase at the Hydraulic Fracturing Test Site 2 (HFTS 2) located in the prolific Wolfcamp formation within the Permian Delaware Basin of west Texas, USA. Specifically, we aimed to assess the microbial ecology and functional potential of the microbial community during this crucial time frame. Prior analysis of 16S rRNA sequencing data through the first 35 days of production revealed a strong selection for a Clostridia species corresponding to a significant decrease in microbial diversity. Here, we performed a metagenomic analysis of produced water sampled on Day 33 of production. This analysis yielded three high-quality metagenome-assembled genomes (MAGs), one of which was a Clostridia draft genome closely related to the recently classified Petromonas tenebris. This draft genome likely represents the dominant Clostridia species observed in our 16S rRNA profile. Annotation of the MAGs revealed the presence of genes involved in critical metabolic processes, including thiosulfate reduction, mixed acid fermentation, and biofilm formation. These findings suggest that this microbial community has the potential to contribute to well souring, biocorrosion, and biofouling within the reservoir. Our research provides unique insights into the early stages of production in one of the most prolific unconventional plays in the United States, with important implications for well management and energy recovery.

natural gas↗
DOE OSTI2020

A clinically and genomically annotated nerve sheath tumor biospecimen repository

Nerve sheath tumors occur as a heterogeneous group of neoplasms in patients with neurofibromatosis type 1 (NF1). The malignant form represents the most common cause of death in people with NF1, and even when benign, these tumors can result in significant disfigurement, neurologic dysfunction, and a range of profound symptoms. Lack of human tissue across the peripheral nerve tumors common in NF1 has been a major limitation in the development of new therapies. To address this unmet need, we have created an annotated collection of patient tumor samples, patient-derived cell lines, and patient-derived xenografts, and carried out high-throughput genomic and transcriptomic characterization to serve as a resource for further biologic and preclinical therapeutic studies. In this work, we release genomic and transcriptomic datasets comprised of 55 tumor samples derived from 23 individuals, complete with clinical annotation. All data are publicly available through the NF Data Portal and at http://synapse.org/jhubiobank.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Sporophyte Stage Genes Exhibit Stronger Selection Than Gametophyte Stage Genes in Haplodiplontic Giant Kelp

Macrocystis pyrifera (giant kelp), a haplodiplontic brown macroalga that alternates between a macroscopic diploid (sporophyte) and a microscopic haploid (gametophyte) phase, provides an ideal system to investigate how ploidy background affects the evolutionary history of a gene. In M. pyrifera , the same genome is subjected to different selective pressures and environments as it alternates between haploid and diploid life stages. We assembled M. pyrifera gene models using available expression data and validated 8,292 genes models using the model alga Ectocarpus siliculosus . Differential expression analysis identified gene models expressed in either or both the haploid and diploid life stages while functional annotation identified processes enriched in each stage. Genes expressed preferentially or exclusively in the gametophyte stage were found to have higher nucleotide diversity (π = 2.3 × 10 –3 and 2.8 × 10 –3 , respectively) than those for sporophytes (π = 1.1 × 10 –3 and 1 × 10 –3 , respectively). While gametophyte-biased genes show faster sequence evolution, the sequence evolution exhibits less signatures of adaptations when compared to sporophyte-biased genes. Our findings contrast the standing masking hypothesis, which predicts higher standing genetic variation at the sporophyte stage, and support the strength of expression theory, which posits that genes expressed more strongly are expected to evolve slower. We argue that the sporophyte stage undergoes more stringent selection compared with the gametophyte stage, which carries a heavy genetic load associated with broadcast spawning. Furthermore, using whole-genome sequencing, we confirm the strong population structure in wild M. pyrifera populations previously established using microsatellite markers, and estimate population genetic parameters, such as pairwise genetic diversity and Tajima’s D , important for conservation and domestication of M. pyrifera .

Molano, Gary↗