Population genomic analyses of protected incense trees Aquilaria sinensis reveal the existence of genetically distinct subpopulations

The incense tree Aquilaria sinensis (Thymelaeaceae) can produce agarwood with commercial values and is now under threat from illegal exploitation in Hong Kong, impairing the local population and biodiversity. Together with other species of Aquilaria, it is listed in the CITES Appendix II, which strictly regulates its international trade. To understand the population structure of A. sinensis and to make relevant conservation measures, we have sequenced 346 individuals collected in Hong Kong and southern mainland China. Population genomic analyses including principal component analysis, neighbor-joining tree construction, ADMIXTURE, and hierarchical pairwise-FST analyses suggested that genetically distinct populations are contained in certain areas. Genomic scan analyses further detected single-nucleotide polymorphism (SNP) outliers related to plant defense, including the CYP71BE gene cluster. In addition to the population analyses, we have developed a modified hexadecyltrimethyl-ammonium bromide (CTAB) DNA extraction protocol for obtaining DNA from agarwood samples in this study, and resequencing of DNA extracted from two agarwood samples using this method allows us to successfully map to the sample corresponding localities in the phylogenetic tree. To sum up, this study suggested that there is a genetically distinct subpopulation of incense tree in Hong Kong that would require special conservation measures and established a foundation for future conservation measures.

[1]  S. Banu,et al.  Genome-wide investigation of Cytochrome P450 superfamily of Aquilaria agallocha: association with terpenoids and phenylpropanoids biosynthesis , 2022, bioRxiv.

[2]  The Darwin,et al.  Sequence locally, think globally: The Darwin Tree of Life Project , 2022, Proceedings of the National Academy of Sciences.

[3]  A. Tillmar,et al.  Getting the conclusive lead with investigative genetic genealogy - A successful case study of a 16 year old double murder in Sweden. , 2021, Forensic science international. Genetics.

[4]  Yong Kang.  Molecular identification of Aquilaria species with distribution records in China using DNA barcode technology , 2021, Mitochondrial DNA. Part B, Resources.

[5]  P. Haddrill.  Developments in forensic DNA analysis , 2021, Emerging topics in life sciences.

[6]  Shuhong Zhao,et al.  rMVP: A Memory-efficient, Visualization-enhanced, and Parallel-accelerated Tool for Genome-wide Association Study , 2020, bioRxiv.

[7]  Margaret H. Frank,et al.  TBtools - an integrative toolkit developed for interactive analyses of big biological data. , 2020, Molecular plant.

[8]  Zhangming Wang,et al.  Using genetic markers to identify the origin of illegally traded agarwood-producing Aquilaria sinensis trees , 2020, Global Ecology and Conservation.

[9]  D. T. Lau,et al.  Chromosomal‐level reference genome of the incense tree Aquilaria sinensis , 2020, bioRxiv.

[10]  P. Cui,et al.  Genome sequence of the agarwood tree Aquilaria sinensis (Lour.) Spreng: the first chromosome-level draft genome in the Thymelaeceae family , 2020, GigaScience.

[11]  Rozi Mohamed,et al.  Cross-amplification of microsatellite markers across agarwood-producing species of the Aquilarieae tribe (Thymelaeaceae) , 2020, 3 Biotech.

[12]  Jennifer Lu,et al.  Improved metagenomic analysis with Kraken 2 , 2019, Genome Biology.

[13]  Chen,et al.  Challenges in the Conservation of an Over-Harvested Plant Species with High Socioeconomic Values , 2019, Sustainability.

[14]  M. Lercher,et al.  Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees , 2019, Nucleic Acids Res..

[15]  J. Flowers,et al.  Origins and geographic diversification of African rice (Oryza glaberrima) , 2018, bioRxiv.

[16]  Rozi Mohamed,et al.  Genome Size, Molecular Phylogeny, and Evolutionary History of the Tribe Aquilarieae (Thymelaeaceae), the Natural Source of Agarwood , 2018, Front. Plant Sci..

[17]  Juan Carlos Castilla-Rubio,et al.  Earth BioGenome Project: Sequencing life for the future of life , 2018, Proceedings of the National Academy of Sciences.

[18]  Petri Auvinen,et al.  Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch , 2017, Nature Genetics.

[19]  Bernardo J. Clavijo,et al.  Genome sequence and genetic diversity of European ash trees , 2016, Nature.

[20]  David S. Wishart,et al.  Heatmapper: web-enabled heat mapping for all , 2016, Nucleic Acids Res..

[21]  Rozi Mohamed,et al.  DNA Barcoding of the Endangered Aquilaria (Thymelaeaceae) and Its Application in Species Authentication of Agarwood Products Traded in the Market , 2016, PloS one.

[22]  Zhian N. Kamvar,et al.  Novel R tools for analysis of genome-wide population genetic data with emphasis on clonality , 2015, Front. Genet..

[23]  Long-fang O. Chen,et al.  Identification of cucurbitacins and assembly of a draft genome for Aquilaria agallocha , 2014, BMC Genomics.

[24]  Yafang Yin,et al.  DNA barcoding for identification of the endangered species Aquilaria sinensis: comparison of data from heated or aged wood samples , 2014 .

[25]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[26]  M. Fumagalli,et al.  Assessing the Effect of Sequencing Depth and Sample Size in Population Genetics Inferences , 2013, PloS one.

[27]  Heng Li.  Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

[28]  Cheng Lu,et al.  Genetic Diversity and Differentiation of Aquilaria sinensis (Lour.) Gilg Revealed by ISSR and SRAP Markers , 2012 .

[29]  J. Bohlmann,et al.  Slow but not low: genomic comparisons reveal slower evolutionary rate and higher dN/dS in conifers compared to angiosperms , 2012, BMC Evolutionary Biology.

[30]  D. Higgins,et al.  Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega , 2011, Molecular systems biology.

[31]  Sean R. Eddy,et al.  Accelerated Profile HMM Searches , 2011, PLoS Comput. Biol..

[32]  R. Durbin,et al.  Inference of Human Population History From Whole Genome Sequence of A Single Individual , 2011, Nature.

[33]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[34]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[35]  Paramvir S. Dehal,et al.  FastTree 2 – Approximately Maximum-Likelihood Trees for Large Alignments , 2010, PloS one.

[36]  D. Nelson.  The Cytochrome P450 Homepage , 2009, Human Genomics.

[37]  David H. Alexander,et al.  Fast model-based estimation of ancestry in unrelated individuals. , 2009, Genome research.

[38]  O. Gaggiotti,et al.  A Genome-Scan Method to Identify Selected Loci Appropriate for Both Dominant and Codominant Markers: A Bayesian Perspective , 2008, Genetics.

[39]  Anne-Béatrice Dufour,et al.  The ade4 Package: Implementing the Duality Diagram for Ecologists , 2007 .

[40]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[41]  J. Goudet.  HIERFSTAT , a package for R to compute and test hierarchical F -statistics , 2005 .

[42]  张志权.  Aquilaria sinensis , 1982, CABI Compendium.

[43]  Yafang Yin,et al.  Wood Resources, Identification, and Utilization of Agarwood in China , 2016 .

[44]  C.Y. Jim.  Cross-border itinerant poaching of agarwood in Hong Kong's peri-urban forests , 2015 .

[45]  Z. Shu.  Extraction of Genomic DNA and Establishment of Siliver-stained AFLP System of Aquilaria sinensis(Lour.)Spreng , 2007 .