Construction of an anchoring SSR marker genetic linkage map and detection of a sex-linked region in two dioecious populations of red bayberry

SSR Marker Acquisition and Application from Transcriptome of Captive Chinese Forest Musk Deer (Moschus berezovskii)

Forest musk deer (Moschus berezovskii) is one of the most endangered medicinally important wild animals in the world. Forest musk deer farming is the main way of production of musk. However, the single provenance and lack of genetic information lead to reduced genetic diversity of forest musk deer. Therefore, more SSR markers need to be developed to identify forest musk deer germplasm. In this study, bone marrow derived mesenchymal cells were used to construct cDNA library for transcriptome sequencing. The datasets were de novo assembled and annotated. 9 polymorphic simple sequence repeat (SSR) markers were finally identified and used to detect population genetic diversity. 6.07 Gb clean data were generated using Illumina sequencing technology, and de novo assembled into 138,591 transcripts and 81,553 unigenes. 5,777 simple sequence repeats (SSRs) were identified, in which there were 578 repeating motif types, with mono-nucleotide and tri-nucleotides comprising 55.88% and 25.60%, respectively. 100 primer pairs were designed to validate amplification and polymorphism using DNA from fecal samples. 9 polymorphic SSRs were developed and used to detect population genetic diversity of 122 forest musk deer in 2 farms. The average number of alleles per locus varied from 4 to 15 (average = 8.3). The observed heterozygosity (HO) per locus ranged from 0.102 to 0.941, while the expected heterozygosity (HE) per locus was from 0.111 to 0.651. All loci deviated significantly from the Hardy-Weinberg equilibrium (p < 0.001). The polymorphism information content (PIC) of these loci varied from 0.108 to 0.619. 9 polymorphic SSR markers were developed in this research. These sites can be used for breeding planning and conservation of germplasm resources.

T2T reference genome assembly and genome-wide association study reveal the genetic basis of Chinese bayberry fruit quality

Chinese bayberry (Myrica rubra or Morella rubra; 2n = 16) produces fruit with a distinctive flavor, high nutritional, and economic value. However, previous versions of the bayberry genome lack sequence continuity. Moreover, to date, no large-scale germplasm resource association analysis has examined the allelic and genetic variations determining fruit quality traits. Therefore, in this study, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar 'Zaojia' using PacBio HiFi long reads. The resulting 292.60 Mb T2T genome, revealed 8 centromeric regions, 15 telomeres, and 28 345 genes. This represents a substantial improvement in the genome continuity and integrity of Chinese bayberry. Subsequently, we re-sequenced 173 accessions, identifying 6 649 674 single nucleotide polymorphisms (SNPs). Further, the phenotypic analyses of 29 fruit quality-related traits enabled a genome-wide association study (GWAS), which identified 1937 SNPs and 1039 genes significantly associated with 28 traits. An SNP cluster pertinent to fruit color was identified on Chr6: 3407532 to 5 153 151 bp region, harboring two MYB genes (MrChr6G07650 and MrChr6G07660), exhibiting differential expression in extreme phenotype transcriptomes, linked to anthocyanin synthesis. An adjacent, closely linked gene, MrChr6G07670 (MLP-like protein), harbored an exonic missense variant and was shown to increase anthocyanin production in tobacco leaves tenfold. This SNP cluster, potentially a quantitative trait locus (QTL), collectively regulates bayberry fruit color. In conclusion, our study presented a complete reference genome, uncovered a suite of allelic variations related to fruit-quality traits, and identified functional genes that could be harnessed to enhance fruit quality and breeding efficiency of bayberries.

An overview of the nutritional value, health properties, and future challenges of Chinese bayberry

Chinese bayberry (CB) is among the most popular and valuable fruits in China owing to its attractive color and unique sweet/sour taste. Recent studies have highlighted the nutritional value and health-related benefits of CB. CB has special biological characteristics of evergreen, special aroma, dioecious, nodulation, nitrogen fixation. Moreover, the fruits, leaves, and bark of CB plants harbor a number of bioactive compounds including proanthocyanidins, flavonoids, vitamin C, phenolic acids, and anthocyanins that have been linked to the anti-cancer, anti-oxidant, anti-inflammatory, anti-obesity, anti-diabetic, and neuroprotective properties and to the treatment of cardiovascular and cerebrovascular diseases. The CB fruits have been used to produce a range of products: beverages, foods, and washing supplies. Future CB-related product development is thus expected to further leverage the health-promoting potential of this valuable ecological resource. The present review provides an overview of the botanical characteristics, processing, nutritional value, health-related properties, and applications of CB in order to provide a foundation for further research and development.

Comparative Transcriptome Analysis Reveals Sex-Biased Expression of Hormone-Related Genes at an Early Stage of Sex Differentiation in Red Bayberry (Morella rubra)

The molecular mechanism of sex development and differentiation in the economically important dioecious fruit tree, red bayberry (Morella rubra), was revealed using next-generation transcriptome sequencing (NGS), and comparative analyses were used to identify differentially expressed genes (DEGs) in female and male flower buds. A total of 7,029 of these DEGs were identified at two early development stages. KEGG pathway enrichment analysis revealed that plant hormone signal transduction was significantly overrepresented, and 91 genes related to hormones were identified. An analysis of 7,029 DEGs revealed 161 hormone-related genes, with the 42 related to auxin and 26 related to ethylene being the most highly represented. A total of 62 genes were significantly up-regulated in females and 29 were in males, with 18 of them specifically expressed in females and 10 in males. A total of 415 transcription factors were identified, with 129 genes up-regulated in females and 53 in males. Moreover, 38 had female-specific expression and 18 had male-specific expression. Using weighted gene co-expression network analysis (WGCNA), two modules were found to be associated with sexual type. In the module coded light-green, there were five genes related to hormones, one to flower development and ten transcription factors with four genes specifically expressed in the males and four in females. The hub gene in the light-green module is MR0TCONS_00017483.1 (ACO), which is involved in ethylene biosynthesis and had male-specific expression. Among the transcription factors, three of the four male-specific expressed genes involved in flavonoid biosynthesis, including the MYB gene MR1TCONS_00020658.1 and two BHLH genes, MR6G001563.1 and MR8G020751.1, played important roles in male floral differentiation. In the dark-cyan module, six hormone-related genes, five transcription factors and three flower development genes were identified with the hub gene MR1G019545.1 (ETR1), which participates in the ethylene signaling pathway, and MR4G023618.1, which encodes the C3H zinc finger transcription factor. These results indicate that ethylene is the key hormone that interacts with other hormones and transcription factors to regulate sex differentiation in the red bayberry, which also provides new insights into the mechanism of sex determination and differentiation in the red bayberry.

MrTPS3 and MrTPS20 Are Responsible for β-Caryophyllene and α-Pinene Production, Respectively, in Red Bayberry (Morella rubra)

Red bayberry is a sweet, tart fruit native to China and grown widely in the south. The key organic compounds forming the distinctive aroma in red bayberry, are terpenoids, mainly β-caryophyllene and α-pinene. However, the key genes responsible for different terpenoids are still unknown. Here, transcriptome analysis on samples from four cultivars, during fruit development, with different terpenoid production, provided candidate genes for volatile organic compound (VOC) production. Terpene synthases (TPS) are key enzymes regulating terpenoid biosynthesis, and 34 TPS family members were identified in the red bayberry genome. MrTPS3 in chromosome 2 and MrTPS20 in chromosome 7 were identified as key genes regulating β-caryophyllene and α-pinene synthesis, respectively, by qRT-PCR. Subcellular localization and enzyme activity assay showed that MrTPS3 was responsible for β-caryophyllene (sesquiterpenes) production and MrTPS20 for α-pinene (monoterpenes). Notably, one amino acid substitution between dark color cultivars and light color cultivars resulted in the loss of function of MrTPS3, causing the different β-caryophyllene production. Our results lay the foundation to study volatile organic compounds (VOCs) in red bayberry and provide potential genes for molecular breeding.

AGAMOUS Gene as a New Sex-Identification Marker in Fig (Ficus carica L.) Is More Efficient Than RAN1

Fig is an ancient gynodioecious fruit tree with females for commercial fruit production and hermaphrodites (males) sometimes used as pollen providers. An early sex-identification method would improve breeding efficiency. Three AGAMOUS (AG) genes were recruited from the Ficus carica genome using AG sequences from Ficus microcarpa and Ficus hispida. FcAG was 5230 bp in length, with 7 exons and 6 introns, and a 744-bp coding sequence. The gene was present in both female and male fig genomes, with a 15-bp deletion in the 7th exon. The other two AG genes (FcAG2-Gall_Stamen and FcAG3-Gall_Stamen) were male-specific, without the 15-bp deletion (759-bp coding sequence), and were only expressed in the gall and stamen of the male fig fruit. Using the deletion as the forward primer (AG-Marker), male plants were very efficiently identified by the presence of a 146-bp PCR product. The previously reported fig male and female polymorphism gene RESPONSIVE-TO-ANTAGONIST1 (RAN1) was also cloned and compared between male and female plants. Fifteen SNPs were found in the 3015-bp protein-coding sequence. Among them, 12 SNPs were identified as having sex-differentiating capacity by checking the sequences of 27 known male and 24 known female cultivars. A RAN1-Marker of 608 bp, including 6 SNPs, was designed, and a PCR and sequencing-based method was verified with 352 fig seedlings from two hybrid populations. Our results confirmed that the newly established AG-Marker is as accurate as the RAN1-Marker, and provide new clues to understanding Ficus sex determination.

Construction of a High-Density Genetic Map and Identification of Leaf Trait-Related QTLs in Chinese Bayberry (Myrica rubra)

Chinese bayberry (Myrica rubra) is an economically important fruit tree that is grown in southern China. Owing to its over 10-year seedling period, the crossbreeding of bayberry is challenging. The characteristics of plant leaves are among the primary factors that control plant architecture and potential yields, making the analysis of leaf trait-related genetic factors crucial to the hybrid breeding of any plant. In the present study, molecular markers associated with leaf traits were identified via a whole-genome re-sequencing approach, and a genetic map was thereby constructed. In total, this effort yielded 902.11 Gb of raw data that led to the identification of 2,242,353 single nucleotide polymorphisms (SNPs) in 140 F1 individuals and parents (Myrica rubra cv. Biqizhong × Myrica rubra cv. 2012LXRM). The final genetic map ultimately incorporated 31,431 SNPs in eight linkage groups, spanning 1,351.85 cM. This map was then used to assemble and update previous scaffold genomic data at the chromosomal level. The genome size of M. rubra was thereby established to be 275.37 Mb, with 94.98% of sequences being assembled into eight pseudo-chromosomes. Additionally, 18 quantitative trait loci (QTLs) associated with nine leaf and growth-related traits were identified. Two QTL clusters were detected (the LG3 and LG5 clusters). Functional annotations further suggested two chlorophyll content-related candidate genes being identified in the LG5 cluster. Overall, this is the first study on the QTL mapping and identification of loci responsible for the regulation of leaf traits in M. rubra, offering an invaluable scientific for future marker-assisted selection breeding and candidate gene analyses.

Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants

Simple sequence repeats (SSRs) are one of the most important genetic markers and widely exist in most species. Here, we identified 249,822 SSRs from 3,951,919 genes in 112 plants. Then, we conducted a comprehensive analysis of these SSRs and constructed a plant SSR database (PSSRD). Interestingly, more SSRs were found in lower plants than in higher plants, showing that lower plants needed to adapt to early extreme environments. Four specific enriched functional terms in the lower plant Chlamydomonas reinhardtii were detected when it was compared with seven other higher plants. In addition, Guanylate_cyc existed in more genes of lower plants than of higher plants. In our PSSRD, we constructed an interactive plotting function in the chart interface, and users can easily view the detailed information of SSRs. All SSR information, including sequences, primers, and annotations, can be downloaded from our database. Moreover, we developed Web SSR Finder and Batch SSR Finder tools, which can be easily used for identifying SSRs. Our database was developed using PHP, HTML, JavaScript, and MySQL, which are freely available at http://www.pssrd.info/ . We conducted an analysis of the Myb gene families and flowering genes as two applications of the PSSRD. Further analysis indicated that whole-genome duplication and whole-genome triplication played a major role in the expansion of the Myb gene families. These SSR markers in our database will greatly facilitate comparative genomics and functional genomics studies in the future.

Genome-Wide Analysis of MYB Gene Family in Chinese Bayberry (Morella rubra) and Identification of Members Regulating Flavonoid Biosynthesis

Chinese bayberry (Morella rubra), the most economically important fruit tree in the Myricaceae family, is a rich source of natural flavonoids. Recently the Chinese bayberry genome has been sequenced, and this provides an opportunity to investigate the organization and evolutionary characteristics of MrMYB genes from a whole genome view. In the present study, we performed the genome-wide analysis of MYB genes in Chinese bayberry and identified 174 MrMYB transcription factors (TFs), including 122 R2R3-MYBs, 43 1R-MYBs, two 3R-MYBs, one 4R-MYB, and six atypical MYBs. Collinearity analysis indicated that both syntenic and tandem duplications contributed to expansion of the MrMYB gene family. Analysis of transcript levels revealed the distinct expression patterns of different MrMYB genes, and those which may play important roles in leaf and flower development. Through phylogenetic analysis and correlation analyses, nine MrMYB TFs were selected as candidates regulating flavonoid biosynthesis. By using dual-luciferase assays, MrMYB12 was shown to trans-activate the MrFLS1 promoter, and MrMYB39 and MrMYB58a trans-activated the MrLAR1 promoter. In addition, overexpression of 35S:MrMYB12 caused a significant increase in flavonol contents and induced the expression of NtCHS, NtF3H, and NtFLS in transgenic tobacco leaves and flowers and significantly reduced anthocyanin accumulation, resulting in pale-pink or pure white flowers. This indicates that MrMYB12 redirected the flux away from anthocyanin biosynthesis resulting in higher flavonol content. The present study provides valuable information for understanding the classification, gene and motif structure, evolution and predicted functions of the MrMYB gene family and identifies MYBs regulating different aspects of flavonoid biosynthesis in Chinese bayberry.

The de novo genome assembly of Tapiscia sinensis and the transcriptomic and developmental bases of androdioecy

Tapiscia sinensis (Tapisciaceae) possesses an unusual androdioecious breeding system that has attracted considerable interest from evolutionary biologists. Key aspects of T. sinensis biology, including its biogeography, genomics, and sex-linked genes, are unknown. Here, we report the first de novo assembly of the genome of T. sinensis. The genome size was 410 Mb, with 22,251 predicted genes. Based on whole-genome resequencing of 55 trees from 10 locations, an analysis of population genetic structure indicated that T. sinensis has fragmented into five lineages, with low intrapopulation genetic diversity and little gene flow among populations. By comparing whole-genome scans of male versus hermaphroditic pools, we identified 303 candidate sex-linked genes, 79 of which (25.9%) were located on scaffold 25. A 24-kb region was absent in hermaphroditic individuals, and five genes in that region, TsF-box4, TsF-box10, TsF-box13, TsSUT1, and TsSUT4, showed expression differences between mature male and hermaphroditic flowers. The results of this study shed light on the breeding system evolution and conservation genetics of the Tapisciaceae.

Genetic Diversity, Pedigree Relationships, and A Haplotype-Based DNA Fingerprinting System of Red Bayberry Cultivars

High throughput sequencing was used to reveal the distribution of whole-genome variations in cultivated Morella rubra (Sieb. et Zucc.). A total of 3,151,123 SNPs, 371,757 small indels, and 15,904 SVs were detected in 52 accessions. Verification by Sanger sequencing demonstrated that the positive rate of the SNPs was approximately 97.3%. Search for more genetic variations was expanded to 141 red bayberry accessions, most of which were cultivars, by sequencing 19 selected genomic segments (SEG1-19). The results showed that each segment harbored, on average, 7.8 alleles (haplotypes), a haplotype diversity of 0.42, and a polymorphic information content (PIC) of 0.40. Seventy-two different genotypes were identified from the 141 accessions, and statistical analysis showed that the accessions with duplicated genotypes were either somatic mutants or simply synonyms. Core set selection results showed that a minimum of 34 genotypes could already have covered all the alleles on the segments. A DNA fingerprinting system was developed for red bayberry, which used the diversity information of only 8 DNA segments yet still achieved a very high efficiency without losing robustness. No large clade was robustly supported by hierarchical clustering, and well-supported small clusters mainly included close relatives. These results should lead to an improved understanding of the genetic diversity of red bayberry and be valuable for future molecular breeding and variety protection.