HairNet2: deep learning to quantify cotton leaf hairiness, a complex genetic and environmental trait.

BACKGROUND: Cotton accounts for 80% of the global natural fibre production. Its leaf hairiness affects insect resistance, fibre yield, and economic value. However, this phenotype is still qualitatively assessed by visually attributing a Genotype Hairiness Score (GHS) to a leaf/plant, or by using the HairNet deep-learning model which also outputs a GHS. Here, we introduce HairNet2, a quantitative deep-learning model which detects leaf hairs (trichomes) from images and outputs a segmentation mask and a Leaf Trichome Score (LTS). RESULTS: Trichomes of 1250 images were annotated (AnnCoT) and a combination of six Feature Extractor modules and five Segmentation modules were tested alongside a range of loss functions and data augmentation techniques. HairNet2 was further validated on the dataset used to build HairNet (CotLeaf-1), a similar dataset collected in two subsequent seasons (CotLeaf-2), and a dataset collected on two genetically diverse populations (CotLeaf-X). The main findings of this study are that (1) leaf number, environment and image position did not significantly affect results, (2) although GHS and LTS mostly correlated for individual GHS classes, results at the genotype level revealed a strong LTS heterogeneity within a given GHS class, (3) LTS correlated strongly with expert scoring of individual images. CONCLUSIONS: HairNet2 is the first quantitative and scalable deep-learning model able to measure leaf hairiness. Results obtained with HairNet2 concur with the qualitative values used by breeders at both extremes of the scale (GHS 1-2, and 5-5+), but interestingly suggest a reordering of genotypes with intermediate values (GHS 3-4+). Finely ranking mild phenotypes is a difficult task for humans. In addition to providing assistance with this task, HairNet2 opens the door to selecting plants with specific leaf hairiness characteristics which may be associated with other beneficial traits to deliver better varieties.

Genome-Wide Mining of CULLIN E3 Ubiquitin Ligase Genes from Uncaria rhynchophylla.

CULLIN (CUL) protein is a subtype of E3 ubiquitin ligase that is involved in a variety of biological processes and responses to stress in plants. In Uncaria rhynchophylla, the CUL gene family has not been identified and its role in plant development, stress response and secondary metabolite synthesis has not been studied. In this study, 12 UrCUL gene members all contained the typical N-terminal domain and C-terminal domain identified from the U. rhynchophylla genome and were classified into four subfamilies based on the phylogenetic relationship with CULs in Arabidopsis thaliana. They were unevenly distributed on eight chromosomes but had a similar structural composition in the same subfamily, indicating that they were relatively conserved and potentially had similar gene functions. An interspecific and intraspecific collinearity analysis showed that fragment duplication played an important role in the evolution of the CUL gene family. The analysis of the cis-acting elements suggests that the UrCULs may play an important role in various biological processes, including the abscisic acid (ABA) response. To investigate this hypothesis, we treated the roots of U. rhynchophylla tissue-cultured seedlings with ABA. The expression pattern analysis showed that all the UrCUL genes were widely expressed in roots with various expression patterns. The co-expression association analysis of the UrCULs and key enzyme genes in the terpenoid indole alkaloid (TIA) synthesis pathway revealed the complex expression patterns of 12 UrCUL genes and some key TIA enzyme genes, especially UrCUL1, UrCUL1-likeA, UrCUL2-likeA and UrCUL2-likeB, which might be involved in the biosynthesis of TIAs. The results showed that the UrCULs were involved in the response to ABA hormones, providing important information for elucidating the function of UrCULs in U. rhynchophylla. The mining of UrCULs in the whole genome of U. rhynchophylla provided new information for understanding the CUL gene and its function in plant secondary metabolites, growth and development.

Genome-wide identification of GATA transcription factor family and the effect of different light quality on the accumulation of terpenoid indole alkaloids in Uncaria rhynchophylla.

Uncaria rhynchophylla is an evergreen vine plant, belonging to the Rubiaceae family, that is rich in terpenoid indole alkaloids (TIAs) that have therapeutic effects on hypertension and Alzheimer's disease. GATA transcription factors (TF) are a class of transcription regulators that participate in the light response regulation, chlorophyll synthesis, and metabolism, with the capability to bind to GATA cis-acting elements in the promoter region of target genes. Currently the charactertics of GATA TFs in U. rhynchophylla and how different light qualities affect the expression of GATA and key enzyme genes, thereby affecting the changes in U. rhynchophylla alkaloids have not been investigated. In this study, 25 UrGATA genes belonging to four subgroups were identified based on genome-wide analysis. Intraspecific collinearity analysis revealed that only segmental duplications were identified among the UrGATA gene family. Collinearity analysis of GATA genes between U. rhynchophylla and four representative plant species, Arabidopsis thaliana, Oryza sativa, Coffea Canephora, and Catharanthus roseus was also performed. U. rhynchophylla seedlings grown in either red lights or under reduced light intensity had altered TIAs content after 21 days. Gene expression analysis reveal a complex pattern of expression from the 25 UrGATA genes as well as a number of key TIA enzyme genes. UrGATA7 and UrGATA8 were found to have similar expression profiles to key enzyme TIA genes in response to altered light treatments, implying that they may be involved in the regulation TIA content. In this research, we comprehensively analyzed the UrGATA TFs, and offered insight into the involvement of UrGATA TFs from U. rhynchophylla in TIAs biosynthesis.

Genetic Mapping and Characterization of Verticillium Wilt Resistance in a Recombinant Inbred Population of Upland Cotton.

Verticillium wilt (VW) is an important and widespread disease of cotton and once established is long-lived and difficult to manage. In Australia, the non-defoliating pathotype of Verticillium dahliae is the most common, and extremely virulent. Breeding cotton varieties with increased VW resistance is the most economical and effective method of controlling this disease and is greatly aided by understanding the genetics of resistance. This study aimed to investigate VW resistance in 240 F7 recombinant inbred lines (RIL) derived from a cross between MCU-5, which has good resistance, and Siokra 1-4, which is susceptible. Using a controlled environment bioassay, we found that resistance based on plant survival or shoot biomass was complex but with major contributions from chromosomes D03 and D09, with genomic prediction analysis estimating a prediction accuracy of 0.73 based on survival scores compared to 0.36 for shoot biomass. Transcriptome analysis of MCU-5 and Siokra 1-4 roots uninfected or infected with V. dahliae revealed that the two cultivars displayed very different root transcriptomes and responded differently to V. dahliae infection. Ninety-nine differentially expressed genes were located in the two mapped resistance regions and so are potential candidates for further identifying the genes responsible for VW resistance.

Reference Genes Screening and Gene Expression Patterns Analysis Involved in Gelsenicine Biosynthesis under Different Hormone Treatments in Gelsemium elegans.

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) is an accurate method for quantifying gene expression levels. Choosing appropriate reference genes to normalize the data is essential for reducing errors. Gelsemium elegans is a highly poisonous but important medicinal plant used for analgesic and anti-swelling purposes. Gelsenicine is one of the vital active ingredients, and its biosynthesis pathway remains to be determined. In this study, G. elegans leaf tissue with and without the application of one of four hormones (SA, MeJA, ETH, and ABA) known to affect gelsenicine synthesis, was analyzed using ten candidate reference genes. The gene stability was evaluated using GeNorm, NormFinder, BestKeeper, ∆CT, and RefFinder. The results showed that the optimal stable reference genes varied among the different treatments and that at least two reference genes were required for accurate quantification. The expression patterns of 15 genes related to the gelsenicine upstream biosynthesis pathway was determined by RT-qPCR using the relevant reference genes identified. Three genes 8-HGO, LAMT, and STR, were found to have a strong correlation with the amount of gelsenicine measured in the different samples. This research is the first study to examine the reference genes of G. elegans under different hormone treatments and will be useful for future molecular analyses of this medically important plant species.

Evaluation of Reference Genes for Normalizing RT-qPCR and Analysis of the Expression Patterns of WRKY1 Transcription Factor and Rhynchophylline Biosynthesis-Related Genes in Uncaria rhynchophylla.

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

Optimization of Isolation and Transformation of Protoplasts from Uncaria rhynchophylla and Its Application to Transient Gene Expression Analysis.

Protoplast-based engineering has become an important tool for basic plant molecular biology research and developing genome-edited crops. Uncaria rhynchophylla is a traditional Chinese medicinal plant with a variety of pharmaceutically important indole alkaloids. In this study, an optimized protocol for U. rhynchophylla protoplast isolation, purification, and transient gene expression was developed. The best protoplast separation protocol was found to be 0.8 M D-mannitol, 1.25% Cellulase R-10, and 0.6% Macerozyme R-10 enzymolysis for 5 h at 26 °C in the dark with constant oscillation at 40 rpm/min. The protoplast yield was as high as 1.5 × 107 protoplasts/g fresh weight, and the survival rate of protoplasts was greater than 90%. Furthermore, polyethylene glycol (PEG)-mediated transient transformation of U. rhynchophylla protoplasts was investigated by optimizing different crucial factors affecting transfection efficiency, including plasmid DNA amount, PEG concentration, and transfection duration. The U. rhynchophylla protoplast transfection rate was highest (71%) when protoplasts were transfected overnight at 24 °C with the 40 µg of plasmid DNA for 40 min in a solution containing 40% PEG. This highly efficient protoplast-based transient expression system was used for subcellular localization of transcription factor UrWRKY37. Finally, a dual-luciferase assay was used to detect a transcription factor promoter interaction by co-expressing UrWRKY37 with a UrTDC-promoter reporter plasmid. Taken together, our optimized protocols provide a foundation for future molecular studies of gene function and expression in U. rhynchophylla.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Siraitia siamensis.

WRKY transcription factors, as the largest gene family in higher plants, play an important role in various biological processes including growth and development, regulation of secondary metabolites, and stress response. In this study, we performed genome-wide identification and analysis of WRKY transcription factors in S. siamensis. A total of 59 SsWRKY genes were identified that were distributed on all 14 chromosomes, and these were classified into three major groups based on phylogenetic relationships. Each of these groups had similar conserved motifs and gene structures. We compared all the S. siamensis SsWRKY genes with WRKY genes identified from three diverse plant species, and the results implied that segmental duplication and tandem duplication play an important roles in the evolution processes of the WRKY gene family. Promoter region analysis revealed that SsWRKY genes included many cis-acting elements related to plant growth and development, phytohormone response, and both abiotic and biotic stress. Expression profiles originating from the transcriptome database showed expression patterns of these SsWRKY genes in four different tissues and revealed that most genes are expressed in plant roots. Fifteen SsWRKY genes with low-temperature response motifs were surveyed for their gene expression under cold stress, showing that most genes displayed continuous up-regulation during cold treatment. Our study provides a foundation for further study on the function and regulatory mechanism of the SsWRKY gene family.

Characterization and Genetic Mapping of Black Root Rot Resistance in Gossypium arboreum L.

Black root rot (BRR) is an economically important disease of cotton and other crops, especially in cooler regions with short growing seasons. Symptoms include black discoloration of the roots, reduced number of lateral roots and stunted or slow plant growth. The cultivated tetraploid Gossypium species are susceptible to BRR. Resistance to BRR was identified in G. arboreum accession BM13H and is associated with reduced and restricted hyphal growth and less sporulation. Transcriptome analysis indicates that BM13H responds to infection at early time points 2- and 3-days post-inoculation, but by day 5, few differentially expressed genes are observed between infected and uninfected roots. Inheritance of BM13H resistance to BRR was evaluated in an F6 recombinant inbred population and shows a single semi-dominant locus conferring resistance that was fine mapped to a region on chromosome 1, containing ten genes including five putative resistance-like genes.

The Research Progress of Taxol in Taxus.

BACKGROUND: Taxus is a valuable woody species with important medicinal value. The bark of Taxus can produce taxol, a natural antineoplastic drug that is widely used in the treatment of breast, ovarian and lung cancers. However, the low content of taxol in the bark of Taxus can not meet the growing clinical demands, so the current research aims at finding ways to increase taxol production. OBJECTIVE: In this review, the research progress of taxol including the factors affecting the taxol content, biosynthesis pathway of taxol, production of taxol in vitro and the application of multi-omics approaches in Taxus as well as future research prospects will be discussed. RESULTS: The taxol content is not only dependent on the species, age and tissues but is also affected by light, moisture levels, temperature, soil fertility and microbes. Most of the enzymes in the taxol biosynthesis pathway have been identified and characterized. Total chemical synthesis, semi-synthesis, plant cell culture and biosynthesis in endophytic fungi have been explored to product taxol. Multi-omics have been used to study Taxus and taxol. CONCLUSION: Further efforts in the identification of unknown enzymes in the taxol biosynthesis pathway, establishment of the genetic transformation system in Taxus and the regulatory mechanism of taxol biosynthesis and Taxus cell growth will play a significant role in improving the yield of taxol in Taxus cells and plants.

Identification of the Genes Involved in Anthocyanin Biosynthesis and Accumulation in Taxus chinensis.

Taxus chinensis is a precious woody species with significant economic value. Anthocyanin as flavonoid derivatives plays a crucial role in plant biology and human health. However, the genes involved in anthocyanin biosynthesis have not been identified in T. chinensis. In this study, twenty-five genes involved in anthocyanin biosynthesis were identified, including chalcone synthase, chalcone isomerase, flavanone 3-hydroxylase, anthocyanidin synthase, flavonoid 3'-hydroxylase, flavonoid 3',5'-hydroxylase, dihydroflavonol 4-reductase, anthocyanidin reductase, and leucoanthocyanidin reductase. The conserved domains and phylogenetic relationships of these genes were characterized. The expression levels of these genes in different tissues and different ages of xylem were investigated. Additionally, the anthocyanin accumulation in xylem of different ages of T. chinensis was measured. The results showed the anthocyanin accumulation was correlated with the expression levels of dihydroflavonol 4-reductase, anthocyanidin synthase, flavonoid 3'-hydroxylase, and flavonoid 3',5'-hydroxylase. Our results provide a basis for studying the regulation of the biosynthetic pathway for anthocyanins and wood color formation in T. chinensis.

A comparative metabolomics analysis of the components of heartwood and sapwood in Taxus chinensis (Pilger) Rehd.

Taxus chinensis is a well-known gymnosperm with great ornamental and medicinal value. Its purple red brown heartwood (HW) has many attributes such as straight texture, high density, mechanical strength, rich elasticity and corrosion resistance that is highly prized commercially. T. chinensis sapwood (SW), in comparison, lacks these important traits. At present, little is known about the differences of metabolites between the SW and HW in T. chinensis. Widely targeted metabolic profiling was performed to analyze the metabolic profiles of HW and SW in T. chinensis using Liquid Chromatography-Electrospray Ionization-Mass Spectrometry (LC-EI-MS). A total of 607 metabolites were detected in HW and SW. Among them, 146 metabolites were significantly higher, and 167 metabolites significantly lower, in HW as compared to SW. These differential metabolites were mainly involved in metabolic pathways and biosynthesis of secondary metabolites, such as flavonoids, flavone and flavonol, phenylpropanoids and antibiotics. Moreover, 71 flavonoids and isoflavones were found to be significantly different between HW and SW. Our results show the difference of components between the HW and SW, which has potential significance to further elucidate the mechanism of HW color formation. The results will provide insight into the metabolites associated with wood color formation and useful information for understanding the metabolites associated with wood quality.

Transcriptomic Analysis of Betula halophila in Response to Salt Stress.

Soil salinization is a matter of concern worldwide. It can eventually lead to the desertification of land and severely damage local agricultural production and the ecological environment. Betula halophila is a tree with high salt tolerance, so it is of importance to understand and discover the salt responsive genes of B. halophila for breeding salinity resistant varieties of trees. However, there is no report on the transcriptome in response to salt stress in B. halophila. Using Illumina sequencing platform, approximately 460 M raw reads were generated and assembled into 117,091 unigenes. Among these unigenes, 64,551 unigenes (55.12%) were annotated with gene descriptions, while the other 44.88% were unknown. 168 up-regulated genes and 351 down-regulated genes were identified, respectively. These Differentially Expressed Genes (DEGs) involved in multiple pathways including the Salt Overly Sensitive (SOS) pathway, ion transport and uptake, antioxidant enzyme, ABA signal pathway and so on. The gene ontology (GO) enrichments suggested that the DEGs were mainly involved in a plant-type cell wall organization biological process, cell wall cellular component, and structural constituent of cell wall molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment showed that the top-four enriched pathways were 'Fatty acid elongation', 'Ribosome', 'Sphingolipid metabolism' and 'Flavonoid biosynthesis'. The expression patterns of sixteen DEGs were analyzed by qRT-PCR to verify the RNA-seq data. Among them, the transcription factor AT-Hook Motif Nuclear Localized gene and dehydrins might play an important role in response to salt stress in B. halophila. Our results provide an important gene resource to breed salt tolerant plants and useful information for further elucidation of the molecular mechanism of salt tolerance in B. halophila.

Genomewide analysis of the lateral organ boundaries domain gene family in Eucalyptus grandis reveals members that differentially impact secondary growth.

Lateral Organ Boundaries Domain (LBD) proteins are plant-specific transcription factors playing crucial roles in growth and development. However, the function of LBD proteins in Eucalyptus grandis remains largely unexplored. In this study, LBD genes in E. grandis were identified and characterized using bioinformatics approaches. Gene expression patterns in various tissues and the transcriptional responses of EgLBDs to exogenous hormones were determined by qRT-PCR. Functions of the selected EgLBDs were studied by ectopically overexpressing in a hybrid poplar (Populus alba × Populus glandulosa). Expression levels of genes in the transgenic plants were investigated by RNA-seq. Our results showed that there were forty-six EgLBD members in the E. grandis genome and three EgLBDs displayed xylem- (EgLBD29) or phloem-preferential expression (EgLBD22 and EgLBD37). Confocal microscopy indicated that EgLBD22, EgLBD29 and EgLBD37 were localized to the nucleus. Furthermore, we found that EgLBD22, EgLBD29 and EgLBD37 were responsive to the treatments of indol-3-acetic acid and gibberellic acid. More importantly, we demonstrated EgLBDs exerted different influences on secondary growth. Namely, 35S::EgLBD37 led to significantly increased secondary xylem, 35S::EgLBD29 led to greatly increased phloem fibre production, and 35S::EgLBD22 showed no obvious effects. We revealed that key genes related to gibberellin, ethylene and auxin signalling pathway as well as cell expansion were significantly up- or down-regulated in transgenic plants. Our new findings suggest that LBD genes in E. grandis play important roles in secondary growth. This provides new mechanisms to increase wood or fibre production.

Genome-wide identification and characterization of the SPL gene family in Ziziphus jujuba.

SQUAMOSA Promoter-Binding Protein-Likes (SPLs) are plant specific transcription factors playing important roles in plant growth and development. The SPL gene family has been studied in various plant species; however, there is no report about SPLs in Zizyphus jujuba. In this study, we identified 18 putative ZjSPL genes in Z. jujuba using a genome-wide analysis. Sequence features, gene structures, conserved domains and motifs were analyzed. The phylogenetic relationships of SPLs in Z. jujuba and A. thaliana were revealed. A total of 5 pairs of ZjSPLs were identified, suggesting the importance of gene duplication in SPL gene expansion in Z. jujuba. In addition, 11 of the 18 ZjSPLs, belonging to G1, G2 and G5 subgroups, were found to be targets of miR156, suggesting the conservation of miR156-mediated posttranscriptional regulation in plants. Expression analysis revealed that eight ZjSPL genes were responsive to the infection of witches'-broom phytoplasma. Our results provide a basis for the further elucidation of the biological function of ZjSPLs and their regulation in witches'-broom disease.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array.

BACKGROUND: Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits. RESULTS: The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content. CONCLUSIONS: Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Enhancing Integrated Pest Management in GM Cotton Systems Using Host Plant Resistance.

Cotton has lost many ancestral defensive traits against key invertebrate pests. This is suggested by the levels of resistance to some pests found in wild cotton genotypes as well as in cultivated landraces and is a result of domestication and a long history of targeted breeding for yield and fiber quality, along with the capacity to control pests with pesticides. Genetic modification (GM) allowed integration of toxins from a bacteria into cotton to control key Lepidopteran pests. Since the mid-1990s, use of GM cotton cultivars has greatly reduced the amount of pesticides used in many cotton systems. However, pests not controlled by the GM traits have usually emerged as problems, especially the sucking bug complex. Control of this complex with pesticides often causes a reduction in beneficial invertebrate populations, allowing other secondary pests to increase rapidly and require control. Control of both sucking bug complex and secondary pests is problematic due to the cost of pesticides and/or high risk of selecting for pesticide resistance. Deployment of host plant resistance (HPR) provides an opportunity to manage these issues in GM cotton systems. Cotton cultivars resistant to the sucking bug complex and/or secondary pests would require fewer pesticide applications, reducing costs and risks to beneficial invertebrate populations and pesticide resistance. Incorporation of HPR traits into elite cotton cultivars with high yield and fiber quality offers the potential to further reduce pesticide use and increase the durability of pest management in GM cotton systems. We review the challenges that the identification and use of HPR against invertebrate pests brings to cotton breeding. We explore sources of resistance to the sucking bug complex and secondary pests, the mechanisms that control them and the approaches to incorporate these defense traits to commercial cultivars.

Identification of novel and conserved microRNAs in Panax notoginseng roots by high-throughput sequencing.

BACKGROUND: MicroRNAs (miRNAs) are small, non-coding RNAs that are important regulators of gene expression, and play major roles in plant development and their response to the environment. Root extracts from Panax notoginseng contain triterpene saponins as their principal bioactive constituent, and demonstrate medicinal properties. To investigate the novel and conserved miRNAs in P. notoginseng, three small RNA libraries constructed from 1-, 2-, and 3-year-old roots in which root saponin levels vary underwent high-throughput sequencing. METHODS: P. notoginseng roots, purified from 1-, 2-, and 3-year-old roots, were extracted for RNA, respectively. Three small libraries were constructed and subjected to next generation sequencing. RESULTS: Sequencing of the three libraries generated 67,217,124 clean reads from P. notoginseng roots. A total of 316 conserved miRNAs (belonging to 67 miRNA families and one unclassified family) and 52 novel miRNAs were identified. MIR156 and MIR166 were the largest miRNA families, while miR156i and miR156g showed the highest abundance of miRNA species. Potential miRNA target genes were predicted and annotated using Cluster of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes. Comparing these miRNAs between root samples revealed 33 that were differentially expressed between 2- and 1-year-old roots (8 increased, 25 decreased), 27 differentially expressed between 3- and 1-year-old roots (7 increased, 20 decreased), and 29 differentially expressed between 3- and 2-year-old roots (8 increased, 21 decreased). Two significantly differentially expressed miRNAs and four miRNAs predicted to target genes involved in the terpenoid backbone biosynthesis pathway were selected and validated by quantitative reverse transcription PCR. Furthermore, the expression patterns of these six miRNAs were analyzed in P. notoginseng roots, stems, and leaves at different developmental stages. CONCLUSIONS: This study identified a large number of P. notoginseng miRNAs and their target genes, functional annotations, and gene expression patterns. It provides the first known miRNA profiles of the P. notoginseng root development cycle.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

Baseline survey of root-associated microbes of Taxus chinensis (Pilger) Rehd.

Taxol (paclitaxel) a diterpenoid is one of the most effective anticancer drugs identified. Biosynthesis of taxol was considered restricted to the Taxus genera until Stierle et al. discovered that an endophytic fungus isolated from Taxus brevifolia could independently synthesize taxol. Little is known about the mechanism of taxol biosynthesis in microbes, but it has been speculated that its biosynthesis may differ from plants. The microbiome from the roots of Taxus chinensis have been extensively investigated with culture-dependent methods to identify taxol synthesizing microbes, but not using culture independent methods.,Using bar-coded high-throughput sequencing in combination with a metagenomics approach, we surveyed the microbial diversity and gene composition of the root-associated microbiomefrom Taxus chinensis (Pilger) Rehd. High-throughput amplicon sequencing revealed 187 fungal OTUs which is higher than any previously reported fungal number identified with the culture-dependent method, suggesting that T. chinensis roots harbor novel and diverse fungi. Some operational taxonomic units (OTU) identified were identical to reported microbe strains possessing the ability to synthesis taxol and several genes previously associated with taxol biosynthesis were identified through metagenomics analysis.