Population genomic analyses from low-coverage RAD-Seq data: a case study on the non-model cucurbit bottle gourd

Genome-Wide Network Analysis of Above- and Below-Ground Co-growth in Populus euphratica

Tree growth is the consequence of developmental interactions between above- and below-ground compartments. However, a comprehensive view of the genetic architecture of growth as a cohesive whole is poorly understood. We propose a systems biology approach for mapping growth trajectories in genome-wide association studies viewing growth as a complex (phenotypic) system in which above- and below-ground components (or traits) interact with each other to mediate systems behavior. We further assume that trait-trait interactions are controlled by a genetic system composed of many different interactive genes and integrate the Lotka-Volterra predator-prey model to dissect phenotypic and genetic systems into pleiotropic and epistatic interaction components by which the detailed genetic mechanism of above- and below-ground co-growth can be charted. We apply the approach to analyze linkage mapping data of Populus euphratica, which is the only tree species that can grow in the desert, and characterize several loci that govern how above- and below-ground growth is cooperated or competed over development. We reconstruct multilayer and multiplex genetic interactome networks for the developmental trajectories of each trait and their developmental covariation. Many significant loci and epistatic effects detected can be annotated to candidate genes for growth and developmental processes. The results from our model may potentially be useful for marker-assisted selection and genetic editing in applied tree breeding programs. The model provides a general tool to characterize a complete picture of pleiotropic and epistatic genetic architecture in growth traits in forest trees and any other organisms.

A genome-wide association study of coat color in Chinese Rex rabbits

Coat color is an important phenotypic characteristic of the domestic rabbit (Oryctolagus cuniculus) and has specific economic importance in the Rex rabbit industry. Coat color varies considerably among different populations of rabbits, and several causal genes for this variation have been thoroughly studied. Nevertheless, the candidate genes affecting coat color variation in Chinese Rex rabbits remained to be investigated. In this study, we collected blood samples from 250 Chinese Rex rabbits with six different coat colors. We performed genome sequencing using a restriction site-associated DNA sequencing approach. A total of 91,546 single nucleotide polymorphisms (SNPs), evenly distributed among 21 autosomes, were identified. Genome-wide association studies (GWAS) were performed using a mixed linear model, in which the individual polygenic effect was fitted as a random effect. We detected a total of 24 significant SNPs that were located within a genomic region on chromosome 4 (OCU4). After re-fitting the most significant SNP (OCU4:13,434,448, p = 1.31e-12) as a covariate, another near-significant SNP (OCU4:11,344,946, p = 7.03e-07) was still present. Hence, we conclude that the 2.1-Mb genomic region located between these two significant SNPs is significantly associated with coat color in Chinese Rex rabbits. The well-studied coat-color-associated agouti signaling protein (ASIP) gene is located within this region. Furthermore, low genetic differentiation was also observed among the six coat color varieties. In conclusion, our results confirmed that ASIP is a putative causal gene affecting coat color variation in Chinese Rex rabbits.

Identification of Bottle Gourd (Lagenaria siceraria) OVATE Family Genes and Functional Characterization of LsOVATE1

The OVATE gene family is a class of conserved transcription factors that play significant roles in plant growth, development, and abiotic stress, and also affect fruit shape in vegetable crops. Bottle gourd (Lagenaria siceraria), commonly known as calabash or gourd, is an annual climber belonging to the Cucurbitaceae family. Studies on bottle gourd OVATE genes are limited. In this study, we performed genome-wide identification of the OVATE gene family in bottle gourd, and identified a total of 20 OVATE family genes. The identified genes were unevenly distributed across 11 bottle gourd chromosomes. We also analyzed the gene homology, amino acid sequence conservation, and three-dimensional protein structure (via prediction) of the 20 OVATE family genes. We used RNA-seq data to perform expression analysis, which found 20 OVATE family genes to be differentially expressed based on spatial and temporal characteristics, suggesting that they have varying functions in the growth and development of bottle gourd. In situ hybridization and subcellular localization analysis showed that the expression characteristics of the LsOVATE1 gene, located on chromosome 7 homologous to OVATE, is a candidate gene for affecting the fruit shape of bottle gourd. In addition, RT-qPCR data from bottle gourd roots, stems, leaves, and flowers showed different spatial expression of the LsOVATE1 gene. The ectopic expression of LsOVATE1 in tomato generated a phenotype with a distinct fruit shape and development. Transgenic-positive plants that overexpressed LsOVATE1 had cone-shaped fruit, calyx hypertrophy, petal degeneration, and petal retention after flowering. Our results indicate that LsOVATE1 could serve important roles in bottle gourd development and fruit shape determination, and provide a basis for future research into the function of LsOVATE1.

Assessment of genetic diversity of ninety-one bottle gourd [Lagenaria siceraria (Mol.) Standl.] genotypes from fourteen different agro-climatic zones of India using agro-morphological traits and SSR markers

BACKGROUND

The evaluation of the magnitude of genetic diversity present in the germplasm collection is prerequisite for bottle gourd improvement programme. The characterization of the diversity pattern of Indian bottle gourd accessions will facilitate the optimal use of genetic resources for breeding improved cultivars.

METHODS

In the present study, the magnitude of genetic diversity was evaluated in ninety-one genotypes of bottle gourd collected across fourteen different agro-climatic zones of India.

RESULTS

Significant variations were observed for all the studied ten quantitative and nine qualitative traits. The ninety-one genotypes were grouped into nine clusters based on cluster analysis of morphological characteristics. Eigen value from principal component analysis depicted first seven quantitative traits accounted for more than 97.5 cumulative percent of the total variations. The first two components accounted for 50 cumulative percent of the total variation, which signifies a high degree of correlation between the analyzed traits. Molecular diversity with the 40 SSR markers screened revealed 11 polymorphic markers in the genotypes studied. Population structure analysis divulged five populations, conforming to the Principal Coordinate Analysis. Molecular analysis revealed genetically diverse genotypes along with the morphologically divergent genotypes from the quantitative traits and highest inter-cluster distance would be the most appropriate parents for exploiting heterosis.

CONCLUSIONS

The results of this study will facilitate the optimal use of genetic resources for breeding improved cultivars of bottle gourd and the adoption of the identified superior genotypes directly by the breeders.

Ethylene plays a dual role in sex determination and fruit shape in cucurbits

Shapes of vegetables and fruits are the result of adaptive evolution and human selection. Modules controlling organ shape have been identified. However, little is known about signals coordinating organ development and shape. Here, we describe the characterization of a melon mutation rf1, leading to round fruit. Histological analysis of rf1 flower and fruits revealed fruit shape is determined at flower stage 8, after sex determination and before flower fertilization. Using positional cloning, we identified the causal gene as the monoecy sex determination gene CmACS7, and survey of melon germplasms showed strong association between fruit shape and sexual types. We show that CmACS7-mediated ethylene production in carpel primordia enhances cell expansion and represses cell division, leading to elongated fruit. Cell size is known to rise as a result of endoreduplication. At stage 8 and anthesis, we found no variation in ploidy levels between female and hermaphrodite flowers, ruling out endoreduplication as a factor in fruit shape determination. To pinpoint the gene networks controlling elongated versus round fruit phenotype, we analyzed the transcriptomes of laser capture microdissected carpels of wild-type and rf1 mutant. These high-resolution spatiotemporal gene expression dynamics revealed the implication of two regulatory modules. The first module implicates E2F-DP transcription factors, controlling cell elongation versus cell division. The second module implicates OVATE- and TRM5-related proteins, controlling cell division patterns. Our finding highlights the dual role of ethylene in the inhibition of the stamina development and the elongation of ovary and fruit in cucurbits.

Identification and Validation of a Core Single-Nucleotide Polymorphism Marker Set for Genetic Diversity Assessment, Fingerprinting Identification, and Core Collection Development in Bottle Gourd

Germplasm collections are indispensable resources for the mining of important genes and variety improvement. To preserve and utilize germplasm collections in bottle gourd, we identified and validated a highly informative core single-nucleotide polymorphism (SNP) marker set from 1,100 SNPs. This marker set consisted of 22 uniformly distributed core SNPs with abundant polymorphisms, which were established to have strong representativeness and discriminatory power based on analyses of 206 bottle gourd germplasm collections and a multiparent advanced generation inter-cross (MAGIC) population. The core SNP markers were used to assess genetic diversity and population structure, and to fingerprint important accessions, which could provide an optimized procedure for seed authentication. Furthermore, using the core SNP marker set, we developed an accessible core population of 150 accessions that represents 100% of the genetic variation in bottle gourds. This core population will make an important contribution to the preservation and utilization of bottle gourd germplasm collections, cultivar identification, and marker-assisted breeding.

Long-read genome assembly and genetic architecture of fruit shape in the bottle gourd

The bottle gourd (Lagenaria siceraria, Cucurbitaceae) is an important horticultural crop exhibiting tremendous diversity in fruit shape. The genetic architecture of fruit shape variation in this species remains unknown. We assembled a long-read-based, high-quality reference genome (ZAAS_Lsic_2.0) with a contig N50 value over 390-fold greater than the existing reference genomes. We then focused on dissection of fruit shape using a one-step geometric morphometrics-based functional mapping approach. We identified 11 quantitative trait loci (QTLs) responsible for fruit shape (fsQTLs), reconstructed their visible effects and revealed syntenic relationships of bottle gourd fsQTLs with 12 fsQTLs previously reported in cucumber, melon or watermelon. Homologs of several well-known and newly identified fruit shape genes, including SUN, OFP, AP2 and auxin transporters, were comapped with bottle gourd QTLs.

Progress on Genetic Improvement and Analysis of Bottle Gourd [Lagenaria siceraria (Molina) Standl.] for Agronomic Traits, Nutrient Compositions, and Stress Tolerance: A Review

Bottle gourd [Lagenaria siceraria (Molina) Standl.] is an important multi-purpose cucurbit crop grown for its leaf, fruit, and seed. It is widely cultivated and used for human consumption in sub-Saharan Africa (SSA) providing vital human nutrition and serving as food security crop. There is wide genetic variation among bottle gourd genetic resources in Africa for diverse qualitative and quantitative attributes for effective variety design, product development, and marketing. However, the crop is under- researched and -utilized, and improved varieties are yet to be developed and commercialized in the region. Therefore, the objective of this review is to provide the progress on bottle gourd genetic improvement and genetic analysis targeting agronomic and horticultural attributes, nutritional composition, biotic, and abiotic stress tolerance to guide current and future cultivar development, germplasm access, and conservation in SSA. The first section of the paper presents progress on breeding of bottle gourd for horticultural traits, agronomic performance, nutritional and anti-nutritional composition, and biotic and abiotic stress tolerance. This is followed by important highlights on key genetic resources of cultivated and wild bottle gourd for demand driven breeding. Lastly, the review summaries advances in bottle gourd genomics, genetic engineering and genome editing. Information presented in this paper should aid bottle gourd breeders and agronomists to develop and deploy new generation and promising varieties with farmer- and market -preferred attributes.

Novel Genomic Regions of Fusarium Wilt Resistance in Bottle Gourd [Lagenaria siceraria (Mol.) Standl.] Discovered in Genome-Wide Association Study

Fusarium wilt (FW) is a typical soil-borne disease that seriously affects the yield and fruit quality of bottle gourd. Thus, to improve resistance to FW in bottle gourd, the genetic mechanism underlying FW resistance needs to be explored. In this study, we performed a genome-wide association study (GWAS) based on 5,330 single-nucleotide polymorphisms (SNPs) and 89 bottle gourd accessions. The GWAS results revealed a total of 10 SNPs (P ≤ 0.01, -log₁₀ P ≥ 2.0) significantly associated with FW resistance that were detected in at least two environments (2019DI, 2020DI, and the average across the 2 years); these SNPs were located on chromosomes 1, 2, 3, 4, 8, and 9. Linkage disequilibrium (LD) block structure analysis predicted three potential candidate genes for FW resistance. Genes HG_GLEAN_10001030 and HG_GLEAN_10001042 were within the range of the mean LD block of the marker BGReSe_14202; gene HG_GLEAN_10011803 was 280 kb upstream of the marker BGReSe_00818. Real-time quantitative PCR (qRT-PCR) analysis showed that HG_GLEAN_10011803 was significantly up-regulated in FW-infected plants of YD-4, Yin-10, and Hanbi; HG_GLEAN_10001030 and HG_GLEAN_10001042 were specifically up-regulated in FW-infected plants of YD-4. Therefore, gene HG_GLEAN_10011803 is likely the major effect candidate gene for resistance against FW in bottle gourd. This work provides scientific evidence for the exploration of candidate gene and development of functional markers in FW-resistant bottle gourd breeding programs.

QTL-Seq identifies quantitative trait loci of relative electrical conductivity associated with heat tolerance in bottle gourd (Lagenaria siceraria)

Heat is a major abiotic stress that seriously affects watermelon (Citrullus lanatus) production. However, its effects may be mitigated through grafting watermelon to heat tolerant bottle gourd (Lagenaria siceraria) rootstocks. Understanding the genetic basis of heat tolerance and development of reliable DNA markers to indirectly select for the trait are necessary in breeding for new varieties with heat tolerance. The objectives of this study were to investigate the inheritance of heat tolerance and identify molecular markers associated with heat tolerance in bottle gourd. A segregating F₂ population was developed from a cross between two heat tolerant and sensitive inbred lines. The population was phenotyped for relative electrical conductivity (REC) upon high temperature treatment which was used as an indicator for heat tolerance. QTL-seq was performed to identify regions associated with heat tolerance. We found that REC-based heat tolerance in this population exhibited recessive inheritance. Seven heat-tolerant quantitative trait loci (qHT1.1, qHT2.1, qHT2.2, qHT5.1, qHT6.1, qHT7.1, and qHT8.1) were identified with qHT2.1 being a promising major-effect QTL. In the qHT2.1 region, we identified three non-synonymous SNPs that were potentially associated with heat tolerance. These SNPs were located in the genes that may play roles in pollen sterility, intracellular transport, and signal recognition. Association of the three SNPs with heat tolerance was verified in segregating F₂ populations, which could be candidate markers for marker assisted selection for heat tolerance in bottle gourd. The qHT2.1 region is an important finding that may be used for fine mapping and discovery of novel genes associated with heat tolerance in bottle gourd.

Quantitative Trait Locus Analysis in Squash (Cucurbita moschata) Based on Simple Sequence Repeat Markers and Restriction Site-Associated DNA Sequencing Analysis

Squash (Cucurbita moschata) displays wide morphological and genetic variations; however, limited information is available regarding the genetic loci of squash that control its agronomic traits. To obtain basic genetic information for C. moschata, an F2 population was prepared derived from a cross between the Vietnamese cultivar ‘Bí Hồ Lô TN 6 (TN 6)’ and the Japanese cultivar ‘Shishigatani’, and flowering and fruit traits were examined. Overall, the traits showed a continuous distribution in the F2 population, suggesting that they were quantitative traits. A linkage map was constructed based on simple sequence repeat and restriction site-associated DNA (RAD) markers to detect quantitative trait loci (QTLs). Twelve QTLs for flowering and fruit traits, as well as one phenotypic trait locus, were successfully localized on the map. The present QTLs explained the phenotypic variations at a moderate to relatively high level (16.0%–47.3%). RAD markers linked to the QTLs were converted to codominant cleaved amplified polymorphic sequence (CAPS) and derived CAPS markers for the easy detection of alleles. The information reported here provides useful information for understanding the genetics of Cucurbita and other cucurbit species, and for the selection of individuals with ideal traits during the breeding of Cucurbita vegetables.

Molecular confirmation of varietal status in bottle gourd (Lagenaria siceraria) using genotyping-by-sequencing

Bottle gourd (Lagenaria siceraria) is an important food, medicinal, and utilitarian crop with a large pan-tropical distribution. Two morphologically different types in the siceraria subspecies are sufficiently different to be considered as varieties, but they are assigned into different taxonomic ranks. Genotyping-by-sequencing (GBS) of 95 different accessions from the Nangui Abrogoua University collection was used to confirm the varietal status in bottle gourd. This analysis produced 22 575 single-nucleotide polymorphisms (SNPs). Cluster analyses conducted with 2250 (9.96%) SNPs distinctly separated hard-shelled from soft-shelled types. Analysis of 23 SNPs located in 11 genes coding for traits that differentiate the two types of gourds revealed that genes in the soft-shelled types had about 21% fewer SNPs than genes within hard-shelled gourds, but the latter had more non-synonymous SNPs. Cluster analyses conducted with the 23 SNPs fitted well with the structure defined by the 2250 SNPs, suggesting the implication of these SNPs in the varietal differentiation of bottle gourd. These nucleotide changes along with the genetic relationships between the accessions provide molecular proof supporting the status of two varieties. To prevent the confusion inherent in the use of synonyms and homonyms in bottle gourd, we suggest the terms hard-shelled and soft-shelled to designate, respectively, the varieties used as utensils and those grown for their edible seeds.

Helping decision making for reliable and cost-effective 2b-RAD sequencing and genotyping analyses in non-model species

High-throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b-RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b-RAD protocols on non-model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b-RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade-off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective-base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b-RAD protocols on non-model organisms with different genome sizes, helping decision-making for a reliable and cost-effective genotyping.

A universal pipeline for mobile mRNA detection and insights into heterografting advantages under chilling stress

Heterografting has long been used to enhance the chilling tolerance of temperature-sensitive crops, including watermelon, whose mechanism is known to involve bidirectional long-distance mRNA movements. Despite several studies reporting on mobile mRNA (mb-mRNA) profiles in plants, accurate identification of mb-mRNAs is challenging owing to an array of technical problems. Here, we developed a bioinformatical pipeline that took most of the known technical concerns into consideration and is considered to be a universal tool for mb-mRNA detection in heterografts. By applying this pipeline to a commercial watermelon-bottle gourd heterografting system, we detected 130 and 1144 mb-mRNAs upwardly and 167 and 1051 mb-mRNAs downwardly transmitted under normal and chilling-stress conditions, respectively. Quantitative real-time PCR indicated a high accuracy rate (88.2%) of mb-mRNA prediction with our pipeline. We further revealed that the mobility of mRNAs was not associated with their abundance. Functional annotation and classification implied that scions may convey the stress signal to the rootstock, subsequently triggering energy metabolism reprogramming and abscisic acid-mediated stress responses by upward movement of effective mRNAs, ultimately leading to enhanced chilling tolerance. This study provides a universal tool for mb-mRNA detection in plant heterografting systems and novel insights into heterografting advantages under chilling stress.

Detection of significant SNP associated with production and oil quality traits in interspecific oil palm hybrids using RARSeq

A RARSeq based Association mapping study was performed in a population of 104 Elaeis oleifera x E. guineensis hybrids of five origins with the aim of finding functional markers associated to six productive and 19 oil quality traits. For this purpose mRNA of each genotype was isolated and double stranded cDNA was synthesized. Following digestion with two restriction enzymes and adapter ligation, a size selected pool of barcoded amplicons was produced and sequenced using Illumina MiSeq. The obtained sequences were processed with a "snakemake" pipeline, filtered and missing values were imputed. For all traits except two significant effects of the origin was observed. Genetic diversity analyses revealed high variability within origins and an excess of heterozygosity in the population. Two GLM models with Q matrix or PCA matrix as covariates and two MLM models incorporating in addition a Kinship matrix were tested for genotype-phenotype associations using GAPIT software. Using unadjusted p values (< 0.01) 78 potential associations were detected involving 25 SNP and 20 traits. When applying FDR multiple testing with p <  0.05, 25 significant associations remained involving eight SNP and six quality traits. Four SNP were located in genes with a potential relevant biological meaning.

Genet assignment and population structure analysis in a clonal forest-floor herb, Cardamine leucantha, using RAD-seq

To study the genetic structure of clonal plant populations, genotyping and genet detection using genetic markers are necessary to assign ramets to corresponding genets. Assignment is difficult as it involves setting a robust threshold of genetic distance for genet distinction as neighbouring genets in a plant population are often genetically related. Here, we used restriction site-associated DNA sequencing (RAD-seq) for a rhizomatous clonal herb, Cardamine leucantha [Brassicaceae] to accurately determine genet structure in a natural population. We determined a draft genome sequence of this species for the first time, which resulted in 66 617 scaffolds with N₅₀ = 6086 bp and an estimated genome size of approximately 253 Mbp. Using genetic distances based on the RAD-seq analysis, we successfully distinguished ramets that belonged to distinct genets even from a half-sib family. We applied these methods to 372 samples of C. leucantha collected at 1-m interval grids within a 20 × 20 m plot in a natural population in Hokkaido, Japan. From these samples, we identified 61 genets with high inequality in terms of genet size and patchy distribution. Spatial autocorrelation analyses indicated significant aggregation within 7 and 4 m at ramet and genet levels, respectively. An analysis of parallel DNA microsatellite loci (simple sequence repeats) suggested that RAD-seq can provide data that allows robust genet assignment. It remains unclear whether the large genets identified here became dominant stochastically or deterministically. Precise identification of genets will assist further study and characterization of dominant genets.

How "simple" methodological decisions affect interpretation of population structure based on reduced representation library DNA sequencing: A case study using the lake whitefish

Reduced representation (RRL) sequencing approaches (e.g., RADSeq, genotyping by sequencing) require decisions about how much to invest in genome coverage and sequencing depth, as well as choices of values for adjustable bioinformatics parameters. To empirically explore the importance of these “simple” methodological decisions, we generated two independent sequencing libraries for the same 142 individual lake whitefish (Coregonus clupeaformis) using a nextRAD RRL approach: (1) a larger number of loci at low sequencing depth based on a 9mer (library A); and (2) fewer loci at higher sequencing depth based on a 10mer (library B). The fish were selected from populations with different levels of expected genetic subdivision. Each library was analyzed using the STACKS pipeline followed by three types of population structure assessment (F_ST, DAPC and ADMIXTURE) with iterative increases in the stringency of sequencing depth and missing data requirements, as well as more specific a priori population maps. Library B was always able to resolve strong population differentiation in all three types of assessment regardless of the selected parameters, largely due to retention of more loci in analyses. In contrast, library A produced more variable results; increasing the minimum sequencing depth threshold (-m) resulted in a reduced number of retained loci, and therefore lost resolution at high -m values for F_ST and ADMIXTURE, but not DAPC. When detecting fine population differentiation, the population map influenced the number of loci and missing data, which generated artefacts in all downstream analyses tested. Similarly, when examining fine scale population subdivision, library B was robust to changing parameters but library A lost resolution depending on the parameter set. We used library B to examine actual subdivision in our study populations. All three types of analysis found complete subdivision among populations in Lake Huron, ON and Dore Lake, SK, Canada using 10,640 SNP loci. Weak population subdivision was detected in Lake Huron with fish from sites in the north-west, Search Bay, North Point and Hammond Bay, showing slight differentiation. Overall, we show that apparently simple decisions about library construction and bioinformatics parameters can have important impacts on the interpretation of population subdivision. Although potentially more costly on a per-locus basis, early investment in striking a balance between the number of loci and sequencing effort is well worth the reduced genomic coverage for population genetics studies. More conservative stringency settings on STACKS parameters lead to a final dataset that was more consistent and robust when examining both weak and strong population differentiation. Overall, we recommend that researchers approach “simple” methodological decisions with caution, especially when working on non-model species for the first time.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

Molecular Genetic Mapping of Two Complementary Genes Underpinning Fruit Bitterness in the Bottle Gourd (Lagenaria siceraria [Mol.] Standl.)

Fruit bitterness is a serious problem threatening the bottle gourd (Lagenaria siceraria [Mol.] Standl.) industry worldwide. Previous genetic studies indicated that fruit bitterness in the bottle gourd was controlled by a pair of complementary genes. In this study, based on two non-bitter landraces "Hangzhou Gourd" and "Puxian Gourd," each of which carries a single bitterness gene, and their derived segregation populations, we mapped the complementary genes causing fruit bitterness. Quantitative trait locus (QTL) scanning based on an F₂ population detected two QTLs, which was QBt.1 locating in a 17.62-cM interval on linkage group (LG)2 corresponding to a 1.6-Mb region on chromosome 6, and QBt.2 mapped to a 8.44-cM interval on LG9 corresponding to a 1.9-Mb region on chromosome 7. An advanced bulked segregant analysis (A-BSA) well validated the QTL mapping results. Sequence-based comparative analysis showed no syntenic relationship between QBt.1/QBt.2 and the known bitterness genes in cucumber, melon, and watermelon, suggesting that causal genes underlying QBt.1 and QBt.2 were not direct orthologs of the reported cucurbit bitterness genes. Our results shed light on the molecular genetic mechanisms underlying fruit bitterness in the bottle gourd and is useful to guide breeders to properly select parental lines to avoid the occurrence of bitter fruits in breeding programs.

A First Insight into a Draft Genome of Silver Sillago (Sillago sihama) via Genome Survey Sequencing

Sillago sihama has high economic value and is one of the most attractive aquaculture species in China. Despite its economic importance, studies of its genome have barely been performed. In this study, we conducted a first genomic survey of S. sihama using next-generation sequencing (NGS). In total, 45.063 Gb of high-quality sequence data were obtained. For the 17-mer frequency distribution, the genome size was estimated to be 508.50 Mb. The sequence repeat ratio was calculated to be 21.25%, and the heterozygosity ratio was 0.92%. Reads were assembled into 1,009,363 contigs, with a N50 length of 1362 bp, and then into 814,219 scaffolds, with a N50 length of 2173 bp. The average Guanine and Cytosine (GC) content was 45.04%. Dinucleotide repeats (56.55%) were the dominant form of simple sequence repeats (SSR).

Population genomic analysis of mango (Mangifera indica) suggests a complex history of domestication

Humans have domesticated diverse species from across the plant kingdom, yet much of our foundational knowledge of domestication has come from studies investigating relatively few of the most important annual food crops. Here, we examine the impacts of domestication on genetic diversity in a tropical perennial fruit species, mango (Mangifera indica). We used restriction site associated DNA sequencing to generate genomic single nucleotide polymorphism (SNP) data from 106 mango cultivars from seven geographical regions along with 52 samples of closely related species and unidentified cultivars to identify centers of mango genetic diversity and examine how post-domestication dispersal shaped the geographical distribution of diversity. We identify two gene pools of cultivated mango, representing Indian and Southeast Asian germplasm. We found no significant genetic bottleneck associated with the introduction of mango into new regions of the world. By contrast, we show that mango populations in introduced regions have elevated levels of diversity. Our results suggest that mango has a more complex history of domestication than previously supposed, perhaps including multiple domestication events, hybridization and regional selection. Our work has direct implications for mango breeding and genebank management, and also builds on recent efforts to understand how woody perennial crops respond to domestication.

Genome survey, high-resolution genetic linkage map construction, growth-related quantitative trait locus (QTL) identification and gene location in Scylla paramamosain

Scylla paramamosain is one of the most economically important crabs in China. In this study, the first genome survey sequencing of this crab was performed, and the results revealed that the estimated genome size was 1.21 Gb with high heterozygosity (1.3%). Then, RAD technology was used to construct a high-resolution linkage map for this species. A total of 24,444 single nucleotide polymorphism (SNP) makers were grouped into 47 linkage groups. The total length of the linkage groups was 3087.53 cM with a markers interval of 0.92 cM. With the aid of transcriptome and genome scaffold data, 4,271 markers were linked to genes, including several important growth-related genes such as transforming growth factor-beta regulator I, immune related-gene C-type lectin and ecdysone pathway gene broad-complex-like protein. Further, 442 markers, representing 279 QTLs, associated with 24 traits were identified, and of these markers, 78 were linked to genes. Some interesting genes, such as dedicator of cytokinesis protein 3, tenascin-X and DNA helicase MCM8, were believed to have important relationship with specific traits and merit further exploration. The results of this study will accelerate the genetic improvement and genome sequencing analysis of the mud crab.

RAD-Seq Data Point to a Distinct Split in Liriodendron (Magnoliaceae) and Obvious East–West Genetic Divergence in L. chinense

Liriodendron is a Tertiary period relic tree genus with a typical East Asian and North American disjunction distribution pattern. As an angiosperm base group of trees, Liriodendron provides a valuable resource for the study of evolution processes. Here, we reconstruct the phylogeny and population genetic structure of Liriodendron based on the restriction site-associated DNA sequencing (RAD-Seq) of a wide collection of individuals from 16 populations. Our results reveal a clear phylogenetic break between L. chinense and L. tulipifera and obvious genetic divergence between the eastern and western populations of L. chinense, which are consistent with the patterns of geographical distributions. The phylogeographic history and long-term geographical isolation of the genus may be responsible for this pattern. Furthermore, a closer relationship was found between L. tulipifera and the eastern populations of L. chinense, indicating the ancient phylogeny of L. chinense in this area. The results of this study will aid in the development of scientific strategies for the conservation and utilization of the Liriodendron germplasm.

Linkage Disequilibrium Estimation in Low Coverage High-Throughput Sequencing Data

High-throughput sequencing methods provide a cost-effective approach for genotyping and are commonly used in population genetics studies. A drawback of these methods, however, is that sequencing and genotyping errors can arise...

High-throughput sequencing methods that multiplex a large number of individuals have provided a cost-effective approach for discovering genome-wide genetic variation in large populations. These sequencing methods are increasingly being utilized in population genetic studies across a diverse range of species. Two side-effects of these methods, however, are (1) sequencing errors and (2) heterozygous genotypes called as homozygous due to only one allele at a particular locus being sequenced, which occurs when the sequencing depth is insufficient. Both of these errors have a profound effect on the estimation of linkage disequilibrium (LD) and, if not taken into account, lead to inaccurate estimates. We developed a new likelihood method, GUS-LD, to estimate pairwise linkage disequilibrium using low coverage sequencing data that accounts for undercalled heterozygous genotypes and sequencing errors. Our findings show that accurate estimates were obtained using GUS-LD, whereas underestimation of LD results if no adjustment is made for the errors.

GourdBase: a genome-centered multi-omics database for the bottle gourd (Lagenaria siceraria), an economically important cucurbit crop

GourdBase is an integrative data platform for the bottle gourd to examine its multifarious intuitive morphology and annotated genome. GourdBase consists of six main modules that store and interlink multi-omic data: the genome (with transcriptomic data integrated) module, the phenome module, the markers/QTLs module, the maps (genetic, physical and comparative) module, the cultivars module, and the publications module. These modules provide access to various type of data including the annotated reference genome sequence, gene models, transcriptomic data from various tissues, physical and comparative genome maps, molecular markers in different types, phenotypic data for featuring traits including fruit shape and umami taste, and quantitative trait loci (QTLs) that underlie these traits. GourdBase is intuitive, user-friendly and interlinked and is designed to allow researchers, breeders and trained farmers to browse, search and fetch information on interests and assist in genomics-driven studies and breeding. The knowledge base and web interface can be accessed at http://www.gourdbase.cn/.

QTL delineation for five fiber quality traits based on an intra-specific Gossypium hirsutum L. recombinant inbred line population

Gossypium hirsutum L. is the most important fiber crop worldwide and contributes to more than 95% of global cotton production. Marker-assisted selection (MAS) is an effective approach for improving fiber quality, and quantitative trait loci (QTL) mapping of fiber quality traits is important for cotton breeding. In this study, a permanent intra-specific recombinant inbred line (RIL) population containing 137 families was used for fiber quality testing. Based on a previously reported high-density genetic map with an average marker distance of 0.63 cM, 186 additive QTLs were obtained for five fiber quality traits over five consecutive years, including 39 for fiber length (FL), 36 for fiber strength (FS), 50 for fiber uniformity (FU), 33 for micronaire (MC) and 28 for fiber elongation (FE). Three stable QTLs, qMC-A4-1, qMC-D2-3 and qFS-D9-1, were detected in four datasets, and another eight stable QTLs, qMC-A4-2, qMC-D11-2, qFU-A9-1, qFU-A10-4, qFS-D11-1, qFL-D9-2, qFL-D11-1 and qFE-A3-2, were detected in three datasets. The annotated genes in these 11 stable QTLs were collected, and these genes included many transcription factors with functions during fiber development. 33 QTL coincidence regions were found, and these involved nearly half of the total QTLs. Four chromosome regions containing at least 6 QTLs were promising for fine mapping. In addition, 41 pairs of epistatic QTLs (e-QTLs) were screened, including 6 for FL, 30 for FS, 2 for FU and 3 for MC. The identification of stable QTLs adds valuable information for further QTL fine mapping and gene positional cloning for fiber quality genetic detection and provides useful markers for further molecular breeding in enhancing fiber quality.

Translational Research: Exploring and Creating Genetic Diversity

The crop selection process has created a genetic bottleneck ultimately restricting breeding output. Wild relatives of major crops as well as the so-called 'neglected plant' species represent a reservoir of genetic diversity that remains underutilized. These species could be used as a tool to discover new alleles of agronomic interest or could be the target of breeding programs. Targeted induced local lesions in the genome (TILLING) can be used to translate in neglected crops what has been discovered in major crops and reciprocally. However, random mutagenesis, used in TILLING approaches, provides only a limited density of mutational events at a defined target locus. Alternatively, clustered regularly interspaced short palindromic repeats (CRISPR) associated 9 (Cas9) fused to a cytidine deaminase could serve as a localized mutagenic agent to produce high-density mutant populations. Artificial evolution is at hand.

Adding loci improves phylogeographic resolution in red mangroves despite increased missing data: comparing microsatellites and RAD-Seq and investigating loci filtering

The widespread adoption of RAD-Seq data in phylogeography means genealogical relationships previously evaluated using relatively few genetic markers can now be addressed with thousands of loci. One challenge, however, is that RAD-Seq generates complete genotypes for only a small subset of loci or individuals. Simulations indicate that loci with missing data can produce biased estimates of key population genetic parameters, although the influence of such biases in empirical studies is not well understood. Here we compare microsatellite data (8 loci) and RAD-Seq data (six datasets ranging from 239 to 25,198 loci) from red mangroves (Rhizophora mangle) in Florida to evaluate how different levels of data filtering influence phylogeographic inferences. For all datasets, we calculated population genetic statistics and evaluated population structure, and for RAD-Seq datasets, we additionally examined population structure using coalescence. We found higher F_ST using microsatellites, but that RAD-Seq-based estimates approached those based on microsatellites as more loci with more missing data were included. Analyses of RAD-Seq datasets resolved the classic Gulf-Atlantic coastal phylogeographic break, which was not significant in the microsatellite analyses. Applying multiple levels of filtering to RAD-Seq datasets can provide a more complete picture of potential biases in the data and elucidate subtle phylogeographic patterns.

The bottle gourd genome provides insights into Cucurbitaceae evolution and facilitates mapping of a Papaya ring-spot virus resistance locus

Bottle gourd (Lagenaria siceraria) is an important vegetable crop as well as a rootstock for other cucurbit crops. In this study, we report a high-quality 313.4-Mb genome sequence of a bottle gourd inbred line, USVL1VR-Ls, with a scaffold N50 of 8.7 Mb and the longest of 19.0 Mb. About 98.3% of the assembled scaffolds are anchored to the 11 pseudomolecules. Our comparative genomic analysis identifies chromosome-level syntenic relationships between bottle gourd and other cucurbits, as well as lineage-specific gene family expansions in bottle gourd. We reconstructed the genome of the most recent common ancestor of Cucurbitaceae, which revealed that the ancestral Cucurbitaceae karyotypes consisted of 12 protochromosomes with 18 534 protogenes. The 12 protochromosomes are largely retained in the modern melon genome, while have undergone different degrees of shuffling events in other investigated cucurbit genomes. The 11 bottle gourd chromosomes derive from the ancestral Cucurbitaceae karyotypes followed by 19 chromosomal fissions and 20 fusions. The bottle gourd genome sequence has facilitated the mapping of a dominant monogenic locus, Prs, conferring Papaya ring-spot virus (PRSV) resistance in bottle gourd, to a 317.8-kb region on chromosome 1. We have developed a cleaved amplified polymorphic sequence (CAPS) marker tightly linked to the Prs locus and demonstrated its potential application in marker-assisted selection of PRSV resistance in bottle gourd. This study provides insights into the paleohistory of Cucurbitaceae genome evolution, and the high-quality genome sequence of bottle gourd provides a useful resource for plant comparative genomics studies and cucurbit improvement.

Genome-wide RAD sequencing data provide unprecedented resolution of the phylogeny of temperate bamboos (Poaceae: Bambusoideae)

The temperate bamboos (tribe Arundinarieae, Poaceae) are strongly supported as monophyly in recent molecular studies, but taxonomic delineation and phylogenetic relationships within the tribe lack resolution. Here, we sampled 39 species (36 temperate bamboos and 3 outgroups) for restriction-site associated DNA sequencing (RAD-seq) with an emphasis on Phyllostachys clade and related clades. Using the largest data matrix for the bamboos to date, we were able to infer phylogenetic relationships with unparalleled resolution. The Phyllostachys, Shibataea, and Arundinaria clades defined from plastid phylogeny, were not supported as monophyletic group. However, the RAD-seq phylogeny largely agreed with the morphology-based taxonomy, with two clades having leptomorph rhizomes strongly supported as monophyletic group. We also explored two approaches, BWA-GATK (a mapping system) and Stacks (a grouping system), for differences in SNP calling and phylogeny inference. For the same level of missing data, the BWA-GATK pipeline produced much more SNPs in comparison with Stacks. Phylogenetic analyses of the largest data matrices from both pipelines, using concatenation and coalescent methods provided similar tree topologies, despite the presence of missing data. Our study demonstrates the utility of RAD-seq data for elucidating phylogenetic relationships between genera and higher taxonomic levels in this important but phylogenetically challenging group.

The genetic architecture of shoot-root covariation during seedling emergence of a desert tree, Populus euphratica

The coordination of shoots and roots is critical for plants to adapt to changing environments by fine-tuning energy production in leaves and the availability of water and nutrients from roots. To understand the genetic architecture of how these two organs covary during developmental ontogeny, we conducted a mapping experiment using Euphrates poplar (Populus euphratica), a so-called hero tree able to grow in the desert. We geminated intraspecific F₁ seeds of Euphrates Poplar individually in a tube to obtain a total of 370 seedlings, whose shoot and taproot lengths were measured repeatedly during the early stage of growth. By fitting a growth equation, we estimated asymptotic growth, relative growth rate, the timing of inflection point and duration of linear growth for both shoot and taproot growth. Treating these heterochronic parameters as phenotypes, a univariate mapping model detected 19 heterochronic quantitative trait loci (hQTLs), of which 15 mediate the forms of shoot growth and four mediate taproot growth. A bivariate mapping model identified 11 pleiotropic hQTLs that determine the covariation of shoot and taproot growth. Most QTLs detected reside within the region of candidate genes with various functions, thus confirming their roles in the biochemical processes underlying plant growth.

Discovery of genome-wideSNPs by RAD-seqand the genetic diversity of captive hog deer (Axis porcinus)

The hog deer (Axis porcinus) is a small deer whose natural habitat is the wet or moist tall grasslands in South and Southeast Asia. Wild populations have dramatically decreased in recent decades. While wild hog deer were recently acknowledged to be extinct in China, a few captive populations have been maintained. In the present study, we successfully employed the restriction-site-associated DNA sequencing (RAD-seq) technique to generate a genome-wide profile of single-nucleotide polymorphisms (SNPs) in the captive population of hog deer from Chengdu Zoo, China (N = 11). Up to 4.7 million clean reads per sample were sequenced, and 11,155 SNPs and 8,247 haplotypes were simultaneously observed within more than 80% of sequenced individuals. In this population, the mean frequency of major alleles at each polymorphism site was 0.7903±0.0014, and the average nucleotide diversity (π) and inbreeding coefficient (FIS) were 0.3031±0.0015 and -0.0302±0.0062, respectively. Additionally, the Euclidean distance-based multidimensional scaling method revealed that the pairwise genetic relatedness was evenly distributed. However, the results of homologous searching by short reads did not provide any meaningful explanation of the phylogenetic relationship of hog deer, which should be further investigated. In conclusion, our results revealed current state of genetic diversity in this captive population of hog deer.Furthermore, these genome-wide SNPs would be useful for guiding the mating schedule to avoid sharp increase of inbreeding coefficient.

Application of RAD Sequencing for Evaluating the Genetic Diversity of Domesticated Panax notoginseng (Araliaceae)

Panax notoginseng, a traditional Chinese medicinal plant, has been cultivated and domesticated for approximately 400 years, mainly in Yunnan and Guangxi, two provinces in southwest China. This species was named according to cultivated rather than wild individuals, and no wild populations had been found until now. The genetic resources available on farms are important for both breeding practices and resource conservation. In the present study, the recently developed technology RADseq, which is based on next-generation sequencing, was used to analyze the genetic variation and differentiation of P. notoginseng. The nucleotide diversity and heterozygosity results indicated that P. notoginseng had low genetic diversity at both the species and population levels. Almost no genetic differentiation has been detected, and all populations were genetically similar due to strong gene flow and insufficient splitting time. Although the genetic diversity of P. notoginseng was low at both species and population levels, several traditional plantations had relatively high genetic diversity, as revealed by the He and π values and by the private allele numbers. These valuable genetic resources should be protected as soon as possible to facilitate future breeding projects. The possible geographical origin of Sanqi domestication was discussed based on the results of the genetic diversity analysis.

High-density linkage map construction and QTL analysis for earliness-related traits in Gossypium hirsutum L

Background

Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping.

Results

Using restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1.

Conclusions

We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3269-y) contains supplementary material, which is available to authorized users.

Development of a universal and simplified ddRAD library preparation approach for SNP discovery and genotyping in angiosperm plants

BACKGROUND

The double digest restriction-site associated DNA sequencing technology (ddRAD-seq) is a reduced representation sequencing technology by sampling genome-wide enzyme loci developed on the basis of next-generation sequencing. ddRAD-seq has been widely applied to SNP marker development and genotyping on animals, especially on marine animals as the original ddRAD protocol is mainly built and trained based on animal data. However, wide application of ddRAD-seq technology in plant species has not been achieved so far. Here, we aim to develop an optimized ddRAD library preparation protocol be accessible to most angiosperm plant species without much startup pre-experiment and costs.

RESULTS

We first tested several combinations of enzymes by in silico analysis of 23 plant species covering 17 families of angiosperm and 1 family of bryophyta and found AvaII + MspI enzyme pair produced consistently higher number of fragments in a broad range of plant species. Then we removed two purifying and one quantifying steps of the original protocol, replaced expensive consumables and apparatuses by conventional experimental apparatuses. Besides, we shortened P1 adapter from 37 to 25 bp and designed a new barcode-adapter system containing 20 pairs of barcodes of varying length. This is an optimized ddRAD strategy for angiosperm plants that is economical, time-saving and requires little technical expertise or investment in laboratory equipment. We refer to this simplified protocol as MiddRAD and we demonstrated the utility and flexibility of our approach by resolving phylogenetic relationships of two genera of woody bamboos (Dendrocalamus and Phyllostachys). Overall our results provide empirical evidence for using this method on different model and non-model plants to produce consistent data.

CONCLUSIONS

As MiddRAD adopts an enzyme pair that works for a broad range of angiosperm plants, simplifies library constructing procedure and requires less DNA input, it will greatly facilitate designing a ddRAD project. Our optimization of this method may make ddRAD be widely used in fields of plant population genetics, phylogenetics, phylogeography and molecular breeding.

Application of a dense genetic map for assessment of genomic responses to selection and inbreeding in Heliothis virescens

Adaptation of pest species to laboratory conditions and selection for resistance to toxins in the laboratory are expected to cause inbreeding and genetic bottlenecks that reduce genetic variation. Heliothis virescens, a major cotton pest, has been colonized in the laboratory many times, and a few laboratory colonies have been selected for Bacillus thuringiensis (Bt) resistance. We developed 350-bp double-digest restriction-site associated DNA-sequencing (ddRAD-seq) molecular markers to examine and compare changes in genetic variation associated with laboratory adaptation, artificial selection and inbreeding in this nonmodel insect species. We found that allelic and nucleotide diversity declined dramatically in laboratory-reared H. virescens as compared with field-collected populations. The declines were primarily a result of the loss of low frequency alleles present in field-collected H. virescens. A further, albeit modest decline in genetic diversity was observed in a Bt-selected population. The greatest decline was seen in H. virescens that were sib-mated for 10 generations, in which more than 80% of loci were fixed for a single allele. To determine which regions of the genome were resistant to fixation in our sib-mated line, we generated a dense intraspecific linkage map containing three PCR-based and 659 ddRAD-seq markers. Markers that retained polymorphism were observed in small clusters spread over multiple linkage groups, but this clustering was not statistically significant. Overall, we have confirmed and extended the general expectations for reduced genetic diversity in laboratory colonies, provided tools for further genomic analyses and produced highly homozygous genomic DNA for future whole genome sequencing of H. virescens.

The report of my death was an exaggeration: A review for researchers using microsatellites in the 21st century

Microsatellites, or simple sequence repeats (SSRs), have long played a major role in genetic studies due to their typically high polymorphism. They have diverse applications, including genome mapping, forensics, ascertaining parentage, population and conservation genetics, identification of the parentage of polyploids, and phylogeography. We compare SSRs and newer methods, such as genotyping by sequencing (GBS) and restriction site associated DNA sequencing (RAD-Seq), and offer recommendations for researchers considering which genetic markers to use. We also review the variety of techniques currently used for identifying microsatellite loci and developing primers, with a particular focus on those that make use of next-generation sequencing (NGS). Additionally, we review software for microsatellite development and report on an experiment to assess the utility of currently available software for SSR development. Finally, we discuss the future of microsatellites and make recommendations for researchers preparing to use microsatellites. We argue that microsatellites still have an important place in the genomic age as they remain effective and cost-efficient markers.

High-resolution genetic linkage mapping, high-temperature tolerance and growth-related quantitative trait locus (QTL) identification in Marsupenaeus japonicus

The Kuruma prawn, Marsupenaeus japonicus, is one of the most promising marine invertebrates in the industry in Asia, Europe and Australia. However, the increasing global temperatures result in considerable economic losses in M. japonicus farming. In the present study, to select genetically improved animals for the sustainable development of the Kuruma prawn industry, a high-resolution genetic linkage map and quantitative trait locus (QTL) identification were performed using the RAD technology. The maternal map contained 5849 SNP markers and spanned 3127.23 cM, with an average marker interval of 0.535 cM. Instead, the paternal map contained 3927 SNP markers and spanned 3326.19 cM, with an average marker interval of 0.847 cM. The consensus map contained 9289 SNP markers and spanned 3610.90 cM, with an average marker interval of 0.388 cM and coverage of 99.06 % of the genome. The markers were grouped into 41 linkage groups in the maps. Significantly, negative correlation was detected between high-temperature tolerance (UTT) and body weight (BW). The QTL mapping revealed 129 significant QTL loci for UTT and four significant QTL loci for BW at the genome-wide significance threshold. Among these QTLs, 129 overlapped with linked SNPs, and the remaining four were located in regions between contiguous SNPs. They explained the total phenotypic variance ranging from 8.9 to 12.4 %. Because of a significantly negative correlation between growth and high-temperature tolerance, we demonstrate that this high-resolution linkage map and QTLs would be useful for further marker-assisted selection in the genetic improvement of M. japonicus.

Genomic signature of successful colonization of Eurasia by the allopolyploid shepherd's purse (Capsella bursa-pastoris)

Polyploidization is a dominant feature of flowering plant evolution. However, detailed genomic analyses of the interpopulation diversification of polyploids following genome duplication are still in their infancy, mainly because of methodological limits, both in terms of sequencing and computational analyses. The shepherd's purse (Capsella bursa-pastoris) is one of the most common weed species in the world. It is highly self-fertilizing, and recent genomic data indicate that it is an allopolyploid, resulting from hybridization between the ancestors of the diploid species Capsella grandiflora and Capsella orientalis. Here, we investigated the genomic diversity of C. bursa-pastoris, its population structure and demographic history, following allopolyploidization in Eurasia. To that end, we genotyped 261 C. bursa-pastoris accessions spread across Europe, the Middle East and Asia, using genotyping-by-sequencing, leading to a total of 4274 SNPs after quality control. Bayesian clustering analyses revealed three distinct genetic clusters in Eurasia: one cluster grouping samples from Western Europe and Southeastern Siberia, the second one centred on Eastern Asia and the third one in the Middle East. Approximate Bayesian computation (ABC) supported the hypothesis that C. bursa-pastoris underwent a typical colonization history involving low gene flow among colonizing populations, likely starting from the Middle East towards Europe and followed by successive human-mediated expansions into Eastern Asia. Altogether, these findings bring new insights into the recent multistage colonization history of the allotetraploid C. bursa-pastoris and highlight ABC and genotyping-by-sequencing data as promising but still challenging tools to infer demographic histories of selfing allopolyploids.

RAD-seq data point to a northern origin of the arctic-alpine genus Cassiope (Ericaceae)

Many arctic-alpine plants display a highly disjunct distribution between the Arctic/Boreal regions and the southern Asian mountains. Two main hypotheses have been proposed to explain the origin of this biogeographic pattern: (1) south-to-north migration in the late Pliocene/early Pleistocene, and (2) north-to-south migration during the Miocene. The genus Cassiope is disjunctly distributed between the Arctic/Boreal regions and the Himalayan-Hengduan Mountains (HHM) and was selected to test these hypotheses. We constructed a fossil-calibrated phylogeny of Ericaceae using two plastid regions to estimate the crown group age of Cassiope, and used sequence data from thousands of loci produced by restriction site associated DNA sequencing (RAD-seq) to reconstruct the phylogeny of Cassiope. We also performed Bayesian divergence time analysis and biogeographic analysis. The Cassiope crown group was estimated to have originated in the Miocene, which predates the onset of Northern hemisphere glaciation. All HHM species formed a clade together with one eastern Siberian species, and this clade was sister to all other Arctic/Boreal species. This topology implies a northern origin of Cassiope, which is confirmed by our biogeographic analysis. Our results thus suggest that the ancient north-to-south migration hypothesis is most consistent with the origin of Cassiope.

Phylogeny of Morella rubra and Its Relatives (Myricaceae) and Genetic Resources of Chinese Bayberry Using RAD Sequencing

Phylogenetic relationships among Chinese species of Morella (Myricaceae) are unresolved. Here, we use restriction site-associated DNA sequencing (RAD-seq) to identify candidate loci that will help in determining phylogenetic relationships among Morella rubra, M. adenophora, M. nana and M. esculenta. Three methods for inferring phylogeny, maximum parsimony (MP), maximum likelihood (ML) and Bayesian concordance, were applied to data sets including as many as 4253 RAD loci with 8360 parsimony informative variable sites. All three methods significantly favored the topology of (((M. rubra, M. adenophora), M. nana), M. esculenta). Two species from North America (M. cerifera and M. pensylvanica) were placed as sister to the four Chinese species. According to BEAST analysis, we deduced speciation of M. rubra to be at about the Miocene-Pliocene boundary (5.28 Ma). Intraspecific divergence in M. rubra occurred in the late Pliocene (3.39 Ma). From pooled data, we assembled 29378, 21902 and 23552 de novo contigs with an average length of 229, 234 and 234 bp for M. rubra, M. nana and M. esculenta respectively. The contigs were used to investigate functional classification of RAD tags in a BLASTX search. Additionally, we identified 3808 unlinked SNP sites across the four populations of M. rubra and discovered genes associated with fruit ripening and senescence, fruit quality and disease/defense metabolism based on KEGG database.

De Novo Assembly of Bitter Gourd Transcriptomes: Gene Expression and Sequence Variations in Gynoecious and Monoecious Lines

Bitter gourd (Momordica charantia L.) is a nutritious vegetable crop of Asian origin, used as a medicinal herb in Indian and Chinese traditional medicine. Molecular breeding in bitter gourd is in its infancy, due to limited molecular resources, particularly on functional markers for traits such as gynoecy. We performed de novo transcriptome sequencing of bitter gourd using Illumina next-generation sequencer, from root, flower buds, stem and leaf samples of gynoecious line (Gy323) and a monoecious line (DRAR1). A total of 65,540 transcripts for Gy323 and 61,490 for DRAR1 were obtained. Comparisons revealed SNP and SSR variations between these lines and, identification of gene classes. Based on available transcripts we identified 80 WRKY transcription factors, several reported in responses to biotic and abiotic stresses; 56 ARF genes which play a pivotal role in auxin-regulated gene expression and development. The data presented will be useful in both functions studies and breeding programs in bitter gourd.

High-density genetic map construction and QTLs analysis of grain yield-related traits in sesame (Sesamum indicum L.) based on RAD-Seq techonology

BACKGROUND

Sesame (Sesamum indicum L., 2n = 26) is an important oilseed crop with an estimated genome size of 369 Mb. The genetic basis, including the number and locations of quantitative trait loci (QTLs) of sesame grain yield and quality remain poorly understood, due in part to the lack of reliable markers and genetic maps. Here we report on the construction of a hitherto most high-density genetic map of sesame using the restriction-site associated DNA sequencing (RAD-seq) combined with 89 PCR markers, and the identification of grain yield-related QTLs using a recombinant inbred line (RIL) population.

RESULT

In total, 3,769 single-nucleotide polymorphism (SNP) markers were identified from RAD-seq, and 89 polymorphic PCR markers were identified including 44 expressed sequence tag-simple sequence repeats (EST-SSRs), 10 genomic-SSRs and 35 Insertion-Deletion markers (InDels). The final map included 1,230 markers distributed on 14 linkage groups (LGs) and was 844.46 cM in length with an average of 0.69 cM between adjacent markers. Using this map and RIL population, we detected 13 QTLs on 7 LGs and 17 QTLs on 10 LGs for seven grain yield-related traits by the multiple interval mapping (MIM) and the mixed linear composite interval mapping (MCIM), respectively. Three major QTLs had been identified using MIM with R2 > 10.0% or MCIM with ha 2 > 5.0%. Two co-localized QTL groups were identified that partially explained the correlations among five yield-related traits.

CONCLUSION

Three thousand eight hundred and four pairs of new DNA markers including SNPs and InDels were developed by RAD-seq, and a so far most high-density genetic map was constructed based on these markers in combination with SSR markers. Several grain yield-related QTLs had been identified using this population and genetic map. We report here the first QTL mapping of yield-related traits with a high-density genetic map using a RIL population in sesame. Results of this study solidified the basis for studying important agricultural traits and implementing marker-assisted selection (MAS) toward genetic improvement in sesame.