Development of High-Density SNP Markers and Their Application in Evaluating Genetic Diversity and Population Structure in Elaeis guineensis

Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.

Genotyping-by-sequencing informs conservation of Andean palms sources of non-timber forest products

Conservation and sustainable management of lineages providing non-timber forest products are imperative under the current global biodiversity loss. Most non-timber forest species, however, lack genomic studies that characterize their intraspecific variation and evolutionary history, which inform species' conservation practices. Contrary to many lineages in the Andean biodiversity hotspot that exhibit high diversification, the genus Parajubaea (Arecaceae) has only three species despite the genus' origin 22 million years ago. Two of the three palm species, P. torallyi and P. sunkha, are non-timber forest species endemic to the Andes of Bolivia and are listed as IUCN endangered. The third species, P. cocoides, is a vulnerable species with unknown wild populations. We investigated the evolutionary relationships of Parajubaea species and the genetic diversity and structure of wild Bolivian populations. Sequencing of five low-copy nuclear genes (3753 bp) challenged the hypothesis that P. cocoides is a cultigen that originated from the wild Bolivian species. We further obtained up to 15,134 de novo single-nucleotide polymorphism markers by genotyping-by-sequencing of 194 wild Parajubaea individuals. Our total DNA sequencing effort rejected the taxonomic separation of the two Bolivian species. As expected for narrow endemic species, we observed low genetic diversity, but no inbreeding signal. We found three genetic clusters shaped by geographic distance, which we use to propose three management units. Different percentages of missing genotypic data did not impact the genetic structure of populations. We use the management units to recommend in situ conservation by creating new protected areas, and ex situ conservation through seed collection.

ZIP transporters-regulated Zn2+ homeostasis: A novel determinant of human diseases

As an essential trace element for organisms, zinc participates in various physiological processes, such as RNA transcription, DNA replication, cell proliferation, and cell differentiation. The destruction of zinc homeostasis is associated with various diseases. Zinc homeostasis is controlled by the cooperative action of zinc transporter proteins that are responsible for the influx and efflux of zinc. Zinc transporter proteins are mainly categorized into two families: Zrt/Irt-like protein (SLC39A/ZIP) family and zinc transporter (SLC30A/ZNT) family. ZIP transporters contain 14 members, namely ZIP1-14, which can be further divided into four subfamilies. Currently, ZIP transporters-regulated zinc homeostasis is one of the research hotspots. Cumulative evidence suggests that ZIP transporters-regulated zinc homeostasis may cause physiological dysfunction and contribute to the onset and progression of diverse diseases, such as cancers, neurological diseases, and cardiovascular diseases. In this review, we initially discuss the structure and distribution of ZIP transporters. Furthermore, we comprehensively review the latest research progress of ZIP transporters-regulated zinc homeostasis in diseases, providing a new perspective into new therapeutic targets for treating related diseases.

Assessing genetic diversity in critically endangered Chieniodendron hainanense populations within fragmented habitats in Hainan

Habitat fragmentation has led to a reduction in the geographic distribution of species, making small populations vulnerable to extinction due to environmental, demographic, and genetic factors. The wild plant Chieniodendron hainanense, a species with extremely small populations, is currently facing endangerment and thus requires urgent conservation efforts. Understanding its genetic diversity is essential for uncovering the underlying mechanisms of its vulnerability and for developing effective conservation strategies. In our study, we analyzed 35 specimens from six different populations of C. hainanense using genotyping-by-sequencing (GBS) and single nucleotide polymorphism (SNP) methodologies. Our findings indicate that C. hainanense has limited genetic diversity. The observed heterozygosity across the populations ranged from 10.79 to 14.55%, with an average of 13.15%. We categorized the six populations of C. hainanense into two distinct groups: (1) Diaoluoshan and Baishaling, and (2) Wuzhishan, Huishan, Bawangling, and Jianfengling. The genetic differentiation among these populations was found to be relatively weak. The observed loss of diversity is likely a result of the effects of natural selection.

Genome-Wide SNP and Indel Discovery in Abaca (Musa textilis Née) and among Other Musa spp. for Abaca Genetic Resources Management

Abaca (Musa textilis Née) is an economically important fiber crop in the Philippines. Its economic potential, however, is hampered by biotic and abiotic stresses, which are exacerbated by insufficient genomic resources for varietal identification vital for crop improvement. To address these gaps, this study aimed to discover genome-wide polymorphisms among abaca cultivars and other Musa species and analyze their potential as genetic marker resources. This was achieved through whole-genome Illumina resequencing of abaca cultivars and variant calling using BCFtools, followed by genetic diversity and phylogenetic analyses. A total of 20,590,381 high-quality single-nucleotide polymorphisms (SNP) and DNA insertions/deletions (InDels) were mined across 16 abaca cultivars. Filtering based on linkage disequilibrium (LD) yielded 130,768 SNPs and 13,620 InDels, accounting for 0.396 ± 0.106 and 0.431 ± 0.111 of gene diversity across these cultivars. LD-pruned polymorphisms across abaca, M. troglodytarum, M. acuminata and M. balbisiana enabled genetic differentiation within abaca and across the four Musa spp. Phylogenetic analysis revealed the registered varieties Abuab and Inosa to accumulate a significant number of mutations, eliciting further studies linking mutations to their advantageous phenotypes. Overall, this study pioneered in producing marker resources in abaca based on genome-wide polymorphisms vital for varietal authentication and comparative genotyping with the more studied Musa spp.

Amaranth Genomic Resource Database: an integrated database resource of Amaranth genes and genomics

Amaranth (Amaranthus L.) is native to Mexico and North America, where it was cultivated thousands of years ago, but now amaranth is grown worldwide. Amaranth is one of the most promising food crops with high nutritional value and belongs to the family Amaranthaceae. The high-quality genome assembly of cultivated amaranth species (A. hypochondriacus, A. cruentus) and wild/weedy species (A. tuberculatus, A. hybridus, and A. palmeri) has already been reported; therefore, we developed an Amaranth Genomic Resource Database (AGRDB) to provide access to all the genomic information such as genes, SSRs, SNPs, TFs, miRNAs, and transporters in one place. The AGRDB database contains functionally annotated gene information with their sequence details, genic as well as genomic SSRs with their three sets of primers, transcription factors classified into different families with their sequence information and annotation details, putative miRNAs with their family, sequences, and targeted gene details, transporter genes with their superfamily, trans-membrane domain details, and details of genic as well as nongenic SNPs with 3' and 5' flanking sequence information of five amaranth species. A database search can be performed using the gene ID, sequence ID, sequence motif, motif repeat, family name, annotation keyword, scaffold or chromosome numbers, etc. This resource also includes some useful tools, including JBrowse for the visualization of genes, SSRs, SNPs, and TFs on the respective amaranth genomes and BLAST search to perform a BLAST search of the user's query sequence against the amaranth genome as well as protein sequences. The AGRDB database will serve as a potential platform for genetic improvement and characterization of this futuristic crop. The AGRDB database will be accessible via the link: http://www.nbpgr.ernet.in:8080/AmaranthGRD/.

Genetic diversity and structure of the 4th cycle breeding population of Chinese fir (Cunninghamia lanceolata (lamb.) hook)

Studying population genetic structure and diversity is crucial for the marker-assisted selection and breeding of coniferous tree species. In this study, using RAD-seq technology, we developed 343,644 high-quality single nucleotide polymorphism (SNP) markers to resolve the genetic diversity and population genetic structure of 233 Chinese fir selected individuals from the 4^th cycle breeding program, representing different breeding generations and provenances. The genetic diversity of the 4^th cycle breeding population was high with nucleotide diversity (P_i ) of 0.003, and H_o and H_e of 0.215 and 0.233, respectively, indicating that the breeding population has a broad genetic base. The genetic differentiation level between the different breeding generations and different provenances was low (F_st < 0.05), with population structure analysis results dividing the 233 individuals into four subgroups. Each subgroup has a mixed branch with interpenetration and weak population structure, which might be related to breeding rather than provenance, with aggregation from the same source only being in the local branches. Our results provide a reference for further research on the marker-assisted selective breeding of Chinese fir and other coniferous trees.

An overview of remote monitoring methods in biodiversity conservation

Conservation of biodiversity is critical for the coexistence of humans and the sustenance of other living organisms within the ecosystem. Identification and prioritization of specific regions to be conserved are impossible without proper information about the sites. Advanced monitoring agencies like the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES) had accredited that the sum total of species that are now threatened with extinction is higher than ever before in the past and are progressing toward extinct at an alarming rate. Besides this, the conceptualized global responses to these crises are still inadequate and entail drastic changes. Therefore, more sophisticated monitoring and conservation techniques are required which can simultaneously cover a larger surface area within a stipulated time frame and gather a large pool of data. Hence, this study is an overview of remote monitoring methods in biodiversity conservation via a survey of evidence-based reviews and related studies, wherein the description of the application of some technology for biodiversity conservation and monitoring is highlighted. Finally, the paper also describes various transformative smart technologies like artificial intelligence (AI) and/or machine learning algorithms for enhanced working efficiency of currently available techniques that will aid remote monitoring methods in biodiversity conservation.

Analysis of the genetic diversity and population structure of Monochasma savatieri Franch. ex Maxim using novel EST-SSR markers

Background

Monochasma savatieri Franch. ex Maxim is a medicinally valuable herb. However, the collection and protection of the wild germplasm resources of M. savatieri are still insufficient, and their genetic diversity and population structure have been poorly studied.

Results

We collected and examined 46 M. savatieri individuals from Fujian, Hunan, Jiangxi, and Zhejiang provinces for genetic diversity and population structure, using 33 newly developed expressed sequence tag-simple sequence repeat (EST-SSR) markers. Applying these markers, we detected a total of 208 alleles, with an average of 6.303 alleles per locus. The polymorphic information content varied from 0.138 to 0.884 (average: 0.668), indicating a high level of polymorphism. At the population level, there was a low degree of genetic diversity among populations (I = 0.535, He = 0.342), with Zhejiang individuals showing the highest genetic diversity among the four populations (Fst = 0.497), which indicated little gene flow within the M. savatieri populations (Nm = 0.253). Mantel test analysis revealed a significant positive correlation between geographical and genetic distance among populations (R² = 0.3304, p < 0.05), and structure and principal coordinate analyses supported classification of populations into three clusters, which was consistent with the findings of cluster analysis.

Conclusions

As a rare medicinal plants, the protection of M. savatieri does not look optimistic, and accordingly, protective efforts should be beefed up on the natural wild populations. This study provided novel tools and insights for designing effective collection and conservation strategies for M. savatieri.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08832-x.

CerealsDB: A Whistle-Stop Tour of an Open Access SNP Resource

The CerealsDB website, created by members of the Functional Genomics Group at the University of Bristol, provides access to a database containing SNP and genotyping data for hexaploid wheat and, to a lesser extent, its progenitors and several of its relatives. The site is principally aimed at plant breeders and research scientists who wish to obtain information regarding SNP markers; for example, obtain primers used for their identification or the sequences upon which they are based. The database underpinning the website contains circa one million putative varietal SNPs of which several hundreds of thousands have been experimentally validated on a range of common genotyping platforms. For each SNP marker, the site also hosts the allelic scores for thousands of elite wheat varieties, landrace cultivars, and wheat relatives. Tools are available to help negotiate and visualize the datasets. The website has been designed to be simple and straightforward to use and is completely open access.

Genetic Diversity and Population Structure Analysis of the USDA Olive Germplasm Using Genotyping-By-Sequencing (GBS)

Olives are one of the most important fruit and woody oil trees cultivated in many parts of the world. Olive oil is a critical component of the Mediterranean diet due to its importance in heart health. Olives are believed to have been brought to the United States from the Mediterranean countries in the 18th century. Despite the increase in demand and production areas, only a few selected olive varieties are grown in most traditional or new growing regions in the US. By understanding the genetic background, new sources of genetic diversity can be incorporated into the olive breeding programs to develop regionally adapted varieties for the US market. This study aimed to explore the genetic diversity and population structure of 90 olive accessions from the USDA repository along with six popular varieties using genotyping-by-sequencing (GBS)-generated SNP markers. After quality filtering, 54,075 SNP markers were retained for the genetic diversity analysis. The average gene diversity (GD) and polymorphic information content (PIC) values of the SNPs were 0.244 and 0.206, respectively, indicating a moderate genetic diversity for the US olive germplasm evaluated in this study. The structure analysis showed that the USDA collection was distributed across seven subpopulations; 63% of the accessions were grouped into an identifiable subpopulation. The phylogenetic and principal coordinate analysis (PCoA) showed that the subpopulations did not align with the geographical origins or climatic zones. An analysis of the molecular variance revealed that the major genetic variation sources were within populations. These findings provide critical information for future olive breeding programs to select genetically distant parents and facilitate future gene identification using genome-wide association studies (GWAS) or a marker-assisted selection (MAS) to develop varieties suited to production in the US.

Evidence of altitudinal gradient modifying genomic and chemical diversity in populations of Lychnophora pinaster Mart

Lychnophora pinaster Mart. (Asteraceae) is endemic to the Brazilian Cerrado. It is distributed along the altitudinal gradient of the mountainous ranges of the state of Minas Gerais. This study aimed to evaluate the influence of altitude on the genetic diversity of L. pinaster populations and the effects of altitude and climatic factors on essential oil chemical composition. Essential oils from L. pinaster populations from the north (North 01, North 02, and North 03, 700-859 m) and the Metropolitan region of Belo Horizonte (MhBH 01 and MrBH 02, 1366-1498 m) were analyzed. SNP markers from L. pinaster in these regions and Campos das Vertentes (CV 01, CV 02, and CV 03, 1055-1292 m) were also analyzed. The main compounds in essential oils were 14-hydroxy-α-humulene (North 01 and North 03), cedr-8(15)-en-9-α-ol (North 02), 14-acetoxy-α-humulene (MrBH 01), and 4-oxo-15-nor-eudesman-11-ene (MrBH 02). Hierarchical cluster and heatmap analyses showed that the North and MrBH populations included five different groups, indicating the chemical composition of essential oils is distinct in each population. Furthermore, principal component analysis showed that higher altitudes (1366 m and 1498 m) in the MrBH influence the chemical composition of essential oils, and climatic factors determine the chemical composition in North region. The genetic diversity showed that most alleles are in Hardy-Weinberg equilibrium and imply high genetic variation and genetic polymorphisms between populations. Furthermore, the results of Mantel tests (R = 0.3861517; p = 0.04709529; R = 0.9423121; p = 0.02739726) also showed that higher altitude (>1360 m) shapes the genetic diversity at the MrBH. The genetic structure showed that higher altitudes (>1360 m) contribute to the structure of the MrBH populations, but not to North and CV populations. Therefore, the altitudinal ranges of Minas Gerais mountainous ranges determine the higher genetic and chemical diversity of L. pinaster populations.

Insights on Genetic Diversity, Population Structure, and Linkage Disequilibrium in Globally Diverse Coconut Accessions Using Genotyping-by-Sequencing

Genotyping-by-sequencing (GBS) has emerged as a cost-effective approach for genome-wide discovery of single-nucleotide polymorphism (SNP) markers and high-throughput genotyping. In this study, 96 coconut palms, representing 16 accessions from globally diverse origins, were genotyped using the GBS strategy. A total of 10,835 high-quality SNPs, which were identified after stringent filtering, were utilized to assess genetic diversity, population structure, and linkage disequilibrium (LD) analyses. The polymorphism information content (PIC) values of SNPs ranged from 0.1 to 0.4, with a large proportion of SNPs (8633 nos.; 79.7%) having a higher PIC in the range of 0.3-0.4. The genetic diversity analysis revealed the existence of a high level of variation in coconut accessions, with an average expected heterozygosity (H_e) value of 0.43. Unweighted neighbor-joining phylogenetic tree and Bayesian-based model population structure grouped coconut genotypes into four main clusters. The accessions are generally clustered based on their height (tall or dwarf), with a few accession clusterings based on geographical origins. Investigation of LD pattern in coconut indicated a relatively rapid LD decay with a short range (9 kb). The results obtained in this study will contribute to enhancing the capacity of coconut researchers to utilize genetic diversity for further genetic improvement. In addition, it would open up possibilities for performing genomic studies such as genome-wide association studies and genomic selection to accelerate the efficiency and speed of coconut genetic improvement.

Population Genetic Diversity and Structure of an Endangered Salicaceae Species in Northeast China: Chosenia arbutifolia (Pall.) A. Skv

Chosenia arbutifolia (Pall.) A. Skv. is a unique and endangered species belonging to the Salicaceae family. It has great potential for ornamental and industrial use. However, human interference has led to a decrease in and fragmentation of its natural populations in the past two decades. To effectively evaluate, utilize, and conserve available resources, the genetic diversity and population structure of C. arbutifolia were analyzed in this study. A total of 142 individuals from ten provenances were sampled and sequenced. Moderate diversity was detected among these, with a mean expected heterozygosity and Shannon’s Wiener index of 0.3505 and 0.5258, respectively. The inbreeding coefficient was negative, indicating a significant excess of heterozygotes. The fixation index varied from 0.0068 to 0.3063, showing a varied genetic differentiation between populations. Analysis of molecular variance demonstrated that differentiation accounted for 82.23% of the total variation among individuals, while the remaining 17.77% variation was between populations. Furthermore, the results of population structure analysis indicated that the 142 individuals originated from three primitive groups. To provide genetic information and help design conservation and management strategies, landscape genomics analysis was performed by investigating loci associated with environmental variables. Eighteen SNP markers were associated with altitude and annual average temperature, of which five were ascribed with specific functions. In conclusion, the current study furthers the understanding of C. arbutifolia genetic architecture and provides insights for germplasm protection.

Assessment of genetic diversity and population structure of oil palm (Elaeis guineensis Jacq.) field genebank: A step towards molecular-assisted germplasm conservation

Oil palm (Elaeis guineensis) germplasm is exclusively maintained as ex situ living collections in the field for genetic conservation and evaluation. However, this is not for long term and the maintenance of field genebanks is expensive and challenging. Large area of land is required and the germplasms are exposed to extreme weather conditions and casualty from pests and diseases. By using 107 SSR markers, this study aimed to examine the genetic diversity and relatedness of 186 palms from a Nigerian-based oil palm germplasm and to identify core collection for conservation. On average, 8.67 alleles per SSR locus were scored with average effective number of alleles per population ranging from 1.96 to 3.34 and private alleles were detected in all populations. Mean expected heterozygosity was 0.576 ranging from 0.437 to 0.661 and the Wright's fixation index calculated was -0.110. Overall moderate genetic differentiation among populations was detected (mean pairwise population FST = 0.120, gene flow Nm = 1.117 and Nei's genetic distance = 0.466) and this was further confirmed by AMOVA analysis. UPGMA dendogram and Bayesian structure analysis concomitantly clustered the 12 populations into eight genetic groups. The best core collection assembled by Core Hunter ver. 3.2.1 consisted of 58 palms accounting for 31.2% of the original population, which was a smaller core set than using PowerCore 1.0. This core set attained perfect allelic coverage with good representation, high genetic distance between entries, and maintained genetic diversity and structure of the germplasm. This study reported the first molecular characterization and validation of core collections for oil palm field genebank. The established core collection via molecular approach, which captures maximum genetic diversity with minimum redundancy, would allow effective use of genetic resources for introgression and for sustainable oil palm germplasm conservation. The way forward to efficiently conserve the field genebanks into next generation without losing their diversity was further discussed.

SNP detection and population structure evaluation of Salix gordejevii Y. L. Chang et Skv. in Hunshandake Sandland, Inner Mongolia, China

Single nucleotide polymorphisms (SNPs) are the most abundant and richest form of genomic polymorphism and, hence, are highly favorable markers for genetic map construction and genome-wide association studies. Based on the DNA specific-locus amplified fragment sequencing (SLAF-seq) for large-scale SNP detection, the genetic diversity and population structure of Salix gordejevii Y. L. Chang et Skv., a valuable sand-fixing shrub, was assessed in 199 accessions from 20 populations in Hunshandake Sandland of northern China. A total of 623.15 M reads resulted in 30.49 × sequencing depth on average and a mean Q30 of 95.70%, and 2,287,715 SNPs in 178,509 polymorphic SLAF tags were obtained. By discarding minor allele frequency > 0.05 and integrity > 0.8, a total of 93,600 SNPs were retained for population genetic analyses, which revealed that 199 individuals could be divided into six groups based on cross-validation errors. However, this grouping pattern did not match the geographical distribution, indicating that there is no apparent geographic barrier in the blank areas where S. gordejevii was not distributed in Hunshandake Sandland. In addition, the physical distance of linkage disequilibrium decay in the analyzed S. gordejevii individuals was 18.5 kb when r ² = 0.1. The linkage disequilibrium decay distances for different chromosomes varied from 4.6 kb (chromosome 16) to 37.8 kb (chromosome 3). The obtained SNPs offer suitable marker resources for further genetic and genomic studies and will benefit S. gordejevii breeding programs.

Analysis of genetic diversity and discrimination of Oil Palm DxP populations based on the origins of pisifera elite parents

A total of 251 Dura cross Pisifera (DxP) hybrid palms from six populations descending from six parental African Pisifera origins and involving 12 progenies were analyzed with 19 selected Simple Sequence Repeats (SSR) markers. A total of 110 alleles were produced, ranging from three to eight per SSR, with a mean of 5.8 alleles per SSR locus. Of these, 68.5% were considered shared alleles by more than one population and the remaining 31.5% were population specific alleles. They generated between six and 21 haplotypes in all populations, and depending on the SSR marker, between one and 10 haplotypes within populations. Various parameters for analyzing genetic variability, differentiation and genetic structure were computed using GenAlEx, Structure and Darwin software. The obtained results confirmed microsatellites as a robust, feasible and trustful method for obtaining DNA fingerprints, tracing the source of oil palm samples. With respect to the authenticity of materials or for solving legitimacy issues, accession belonging to each population by SSR markers could be distinguished, but additional SSR should be screened for improving this process.

Genomic diversity and genome-wide association analysis related to yield and fatty acid composition of wild American oil palm

Existing Elaeis guineensis cultivars lack sufficient genetic diversity due to extensive breeding. Harnessing variation in wild crop relatives is necessary to expand the breadth of agronomically valuable traits. Using RAD sequencing, we examine the natural diversity of wild American oil palm populations (Elaeis oleifera), a sister species of the cultivated Elaeis guineensis oil palm. We genotyped 192 wild E. oleifera palms collected from seven Latin American countries along with four cultivated E. guineensis palms. Honduras, Costa Rica, Panama and Colombia palms are panmictic and genetically similar. Genomic patterns of diversity suggest that these populations likely originated from the Amazon Basin. Despite evidence of a genetic bottleneck and high inbreeding observed in these populations, there is considerable genetic and phenotypic variation for agronomically valuable traits. Genome-wide association revealed several candidate genes associated with fatty acid composition along with vegetative and yield-related traits. These observations provide valuable insight into the geographic distribution of diversity, phenotypic variation and its genetic architecture that will guide choices of wild genotypes for crop improvement.

Comparative analysis and implications of the chloroplast genomes of three thistles (Carduus L., Asteraceae)

BACKGROUND

Carduus, commonly known as plumeless thistles, is a genus in the Asteraceae family that exhibits both medicinal value and invasive tendencies. However, the genomic data of Carduus (i.e., complete chloroplast genomes) have not been sequenced.

METHODS

We sequenced and assembled the chloroplast genome (cpDNA) sequences of three Carduus species using the Illumina Miseq sequencing system and Geneious Prime. Phylogenetic relationships between Carduus and related taxa were reconstructed using Maximum Likelihood and Bayesian Inference analyses. In addition, we used a single nucleotide polymorphism (SNP) in the protein coding region of the matK gene to develop molecular markers to distinguish C. crispus from C. acanthoides and C. tenuiflorus.

RESULTS

The cpDNA sequences of C. crispus, C. acanthoides, and C. tenuiflorus ranged from 152,342 bp to 152,617 bp in length. Comparative genomic analysis revealed high conservation in terms of gene content (including 80 protein-coding, 30 tRNA, and four rRNA genes) and gene order within the three focal species and members of subfamily Carduoideae. Despite their high similarity, the three species differed with respect to the number and content of repeats in the chloroplast genome. Additionally, eight hotspot regions, including psbI-trnS_GCU, trnE_UUC-rpoB, trnR_UCU-trnG_UCC, psbC-trnS_UGA, trnT_UGU-trnL_UAA, psbT-psbN, petD-rpoA, and rpl16-rps3, were identified in the study species. Phylogenetic analyses inferred from 78 protein-coding and non-coding regions indicated that Carduus is polyphyletic, suggesting the need for additional studies to reconstruct relationships between thistles and related taxa. Based on a SNP in matK, we successfully developed a molecular marker and protocol for distinguishing C. crispus from the other two focal species. Our study provides preliminary chloroplast genome data for further studies on plastid genome evolution, phylogeny, and development of species-level markers in Carduus.

Genome-wide discovery and characterization of long noncoding RNAs in African oil palm (Elaeis guineensis Jacq.)

Long noncoding RNAs (lncRNAs) are an important class of genes and play important roles in a range of biological processes. However, few reports have described the identification of lncRNAs in oil palm. In this study, we applied strand specific RNA-seq with rRNA removal to identify 1,363 lncRNAs from the equally mixed tissues of oil palm spear leaf and six different developmental stages of mesocarp (8–24 weeks). Based on strand specific RNA-seq data and 18 released oil palm transcriptomes, we systematically characterized the expression patterns of lncRNA loci and their target genes. A total of 875 uniq target genes for natural antisense lncRNAs (NAT-lncRNA, 712), long intergenic noncoding RNAs (lincRNAs, 92), intronic-lncRNAs (33), and sense-lncRNAs (52) were predicted. A majority of lncRNA loci (77.8%–89.6%) had low expression in 18 transcriptomes, while only 89 lncRNA loci had medium to high expression in at least one transcriptome. Coexpression analysis between lncRNAs and their target genes indicated that 6% of lncRNAs had expression patterns positively correlated with those of target genes. Based on single nucleotide polymorphism (SNP) markers derived from our previous research, 6,882 SNPs were detected for lncRNAs and 28 SNPs belonging to 21 lncRNAs were associated with the variation of fatty acid contents. Moreover, seven lncRNAs showed expression patterns positively correlated expression pattern with those of genes in de novo fatty acid synthesis pathways. Our study identified a collection of lncRNAs for oil palm and provided clues for further research into lncRNAs that may regulate mesocarp development and lipid metabolism.

Reconstruction of the Largest Pedigree Network for Pear Cultivars and Evaluation of the Genetic Diversity of the USDA-ARS National Pyrus Collection

The USDA-ARS National Clonal Germplasm Repository (NCGR) in Corvallis, Oregon, maintains one of the world's largest and most diverse living Pyrus collection. A thorough genetic characterization of this germplasm will provide relevant information to optimize the conservation strategy of pear biodiversity, support the use of this germplasm in breeding, and increase our knowledge of Pyrus taxonomy, evolution, and domestication. In the last two decades simple sequence repeat (SSR) markers have been used at the NCGR for cultivar identification and small population structure analysis. However, the recent development of the Applied Biosystems Axiom Pear 70K Genotyping Array has allowed high-density single nucleotide polymorphism (SNP)-based genotyping of almost the entire collection. In this study, we have analyzed this rich dataset to discover new synonyms and mutants, identify putative labeling errors in the collection, reconstruct the largest pear cultivar pedigree and further elucidate the genetic diversity of Pyrus.

Genetic diversity and population structure of Quercus fabri Hance in China revealed by genotyping-by-sequencing

Analysis of genetic diversity and population structure among Quercus fabri populations is essential for the conservation and utilization of Q. fabri resources. Here, the genetic diversity and structure of 158 individuals from 13 natural populations of Quercus fabri in China were analyzed using genotyping-by-sequencing (GBS). A total of 459,564 high-quality single nucleotide polymorphisms (SNPs) were obtained after filtration for subsequent analysis. Genetic structure analysis revealed that these individuals can be clustered into two groups and the structure can be explained mainly by the geographic barrier, showed gene introgression from coastal to inland areas and high mountains could significantly hinder the mutual introgression of genes. Genetic diversity analysis indicated that the individual differences within groups are greater than the differences between the two groups. These results will help us better understand the genetic backgrounds of Q. fabri.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

Construction of a high density linkage map in Oil Palm using SPET markers

A high-density genetic linkage map from a controlled cross of two oil palm (Elaeis guineensis) genotypes was constructed based on Single Primer Enrichment Technology (SPET) markers. A 5K panel of hybridization probes were used for this purpose which was derived from previously developed SNP primers in oil palm. Initially, 13,384 SNPs were detected which were reduced to 13,073 SNPs after filtering for only bi-allelic SNP. Around 75% of the markers were found to be monomorphic in the progeny, reducing the markers left for linkage mapping to 3,501. Using Lep-MAP3 software, a linkage map was constructed which contained initially 2,388 markers and had a total length of 1,370 cM. In many cases several adjacent SNP were located on the same locus, due to missing recombination events between them, leading to a total of 1,054 loci on the 16 LG. Nevertheless, the marker density of 1.74 markers per cM (0.57 cM/marker) should allow the detection of QTLs in the future. This study shows that cost efficient SPET markers are suitable for linkage map construction in oil palm and probably, also in other species.

Identifying Vitamin E Biosynthesis Genes in Elaeis guineensis by Genome-Wide Association Study

Elaeis guineensis is a tropical oil crop and has the highest oil yield per unit area. Palm oil has high palmitic acid content and is also rich in vitamins, including vitamin E. We conducted genome-wide association studies in a diversity panel of 161 E. guineensis accessions to identify single-nucleotide polymorphisms (SNPs) linked with vitamin E and validated candidate genes in these marker-associated intervals. Based on the SNPs reported in our previous research, 47 SNP markers were detected to be significantly associated with the variation of tocopherol and tocotrienol content at a cutoff P value of 6.3 × 10^-7. A total of 656 candidate genes in the flanking regions of the 47 SNPs were identified, followed by pathway enrichment analysis. Of these candidate genes, EgHGGT (homogentisate geranylgeranyl transferase) involved in the biosynthesis of tocotrienols had a higher expression level in the mesocarp compared to other tissues. Expression of the EgHGGT gene was positively correlated with the variation in α-tocotrienol content. Induced overexpression of the gene in Arabidopsis caused a significant increase in vitamin E content and production of α-tocotrienols compared to wild Arabidopsis.

Identification and Validation of Candidate Genes Involved in Fatty Acid Content in Oil Palm by Genome-Wide Association Analysis

Oil palm (Elaeis guineensis) is the highest yielding oil crop per unit area worldwide, but its oil is considered unhealthy for human consumption due to its high palmitic acid content (C16:0). In order to facilitate breeding for fatty acid content in oil palm, genome-wide association analysis (GWAS) was used to identify and validate single-nucleotide polymorphism (SNP) markers and underlying candidate genes associated with fatty acid content in a diversity panel of 200 oil palm individuals. A total of 1,261,501 SNP markers previously developed using SLAF-seq (specific locus amplified fragment sequencing) were used for GWAS. Based on this analysis, 62 SNP markers were significantly associated with fatty acid composition, and 223 candidate genes were identified in the flanking regions of these SNPs. We found one gene (acyl-ACP thioesterase B genes) that was involved in fatty acid biosynthesis and that was associated with high palmitic acid content in the mesocarp. Over-expression of this gene caused a significant increase in palmitic acid content. Our study provides key loci that can be used for breeding oil palm cultivars with low palmitic acid content.