Fine Mapping and Candidate Gene Analysis of Rice Grain Length QTL qGL9.1

Grain length (GL) is one of the crucial determinants of rice yield and quality. However, there is still a shortage of knowledge on the major genes controlling the inheritance of GL in japonica rice, which severely limits the improvement of japonica rice yields. Here, we systemically measured the GL of 667 F₂ and 1570 BC₃F₃ individuals derived from two cultivated rice cultivars, Pin20 and Songjing15, in order to identify the major genomic regions associated with GL. A novel major QTL, qGL9.1, was mapped on chromosome 9, which is associated with the GL, using whole-genome re-sequencing with bulked segregant analysis. Local QTL linkage analysis with F₂ and fine mapping with the recombinant plant revealed a 93-kb core region on qGL9.1 encoding 15 protein-coding genes. Only the expression level of LOC_Os09g26970 was significantly different between the two parents at different stages of grain development. Moreover, haplotype analysis revealed that the alleles of Pin20 contribute to the optimal GL (9.36 mm) and GL/W (3.31), suggesting that Pin20 is a cultivated species carrying the optimal GL variation of LOC_Os09g26970. Furthermore, a functional-type mutation (16398989-bp, G>A) located on an exon of LOC_Os09g26970 could be used as a molecular marker to distinguish between long and short grains. Our experiments identified LOC_Os09g26970 as a novel gene associated with GL in japonica rice. This result is expected to further the exploration of the genetic mechanism of rice GL and improve GL in rice japonica varieties by marker-assisted selection.

Genome-Wide Re-Sequencing Data Reveals the Population Structure and Selection Signatures of Tunchang Pigs in China

Tunchang pig is one population of Hainan pig in the Hainan Province of China, with the characteristics of delicious meat, strong adaptability, and high resistance to diseases. To explore the genetic diversity and population structure of Tunchang pigs and uncover their germplasm characteristics, 10 unrelated Tunchang pigs were re-sequenced using the Illumina NovaSeq 150 bp paired-end platform with an average depth of 10×. Sequencing data from 36 individuals of 7 other pig breeds (including 4 local Chinese pig breeds (5 Jinhua, 5 Meishan, 5 Rongchang, and 6 Wuzhishan), and 3 commonly used commercial pig breeds (5 Duorc, 5 Landrace, and 5 Large White)) were downloaded from the NCBI public database. After analysis of genetic diversity and population structure, it has been found that compared to commercial pigs, Tunchang pigs have higher genetic diversity and are genetically close to native Chinese breeds. Three methods, FST, θπ, and XP-EHH, were used to detect selection signals for three breeds of pigs: Tunchang, Duroc, and Landrace. A total of 2117 significantly selected regions and 201 candidate genes were screened. Gene enrichment analysis showed that candidate genes were mainly associated with good adaptability, disease resistance, and lipid metabolism traits. Finally, further screening was conducted to identify potential candidate genes related to phenotypic traits, including meat quality (SELENOV, CBR4, TNNT1, TNNT3, VPS13A, PLD3, SRFBP1, and SSPN), immune regulation (CD48, FBL, PTPRH, GNA14, LOX, SLAMF6, CALCOCO1, IRGC, and ZNF667), growth and development (SYT5, PRX, PPP1R12C, and SMG9), reproduction (LGALS13 and EPG5), vision (SLC9A8 and KCNV2), energy metabolism (ATP5G2), cell migration (EPS8L1), and olfaction (GRK3). In summary, our research results provide a genomic overview of the genetic variation, genetic diversity, and population structure of the Tunchang pig population, which will be valuable for breeding and conservation of Tunchang pigs in the future.

Identification of genetic loci conferring seed coat color based on a high-density map in soybean

Seed coat color is a typical evolutionary trait. Identification of the genetic loci that control seed coat color during the domestication of wild soybean could clarify the genetic variations between cultivated and wild soybean. We used 276 F₁₀ recombinant inbred lines (RILs) from the cross between a cultivated soybean (JY47) and a wild soybean (ZYD00321) as the materials to identify the quantitative trait loci (QTLs) for seed coat color. We constructed a high-density genetic map using re-sequencing technology. The average distance between adjacent markers was 0.31 cM on this map, comprising 9,083 bin markers. We identified two stable QTLs (qSC08 and qSC11) for seed coat color using this map, which, respectively, explained 21.933 and 26.934% of the phenotypic variation. Two candidate genes (CHS3C and CHS4A) in qSC08 were identified according to the parental re-sequencing data and gene function annotations. Five genes (LOC100786658, LOC100801691, LOC100806824, LOC100795475, and LOC100787559) were predicted in the novel QTL qSC11, which, according to gene function annotations, might control seed coat color. This result could facilitate the identification of beneficial genes from wild soybean and provide useful information to clarify the genetic variations for seed coat color in cultivated and wild soybean.

The Spatial Distribution and Genetic Diversity of the Soybean Cyst Nematode, Heterodera glycines, in China: It Is Time to Take Measures to Control Soybean Cyst Nematode

The continuous evolution and spread of virulent forms of the soybean cyst nematode (SCN) driven by the environment and anthropogenic intervention is a serious threat to the soybean production worldwide, including China. Especially in China, the implemented measures to control SCN are insufficient for sustainable agricultural development yet. We summarized our knowledge about the spread and spatial distribution of SCN in China and the virulence diversity in the main soybean growing areas. To reveal the genetic relatedness and diversity of SCN populations, we re-sequenced 53 SCN genomes from the Huang-Huai Valleys, one of the two main soybean growing areas in China. We identified spreading patterns linked to the local agroecosystems and topographies. Moreover, we disclosed the first evidence for the selection of complex virulence in the field even under low selection pressure in an example from North Shanxi. SCN is present in all soybean growing areas in China but SCN susceptible cultivars are still largely grown indicating that SCN-related damage and financial loss have not received the attention they deserve yet. To prevent increasing yield losses and to improve the acceptance of resistant cultivars by the growers, we emphasized that it is time to accelerate SCN resistance breeding, planting resistant cultivars to a larger extent, and to support farmers to implement a wider crop rotation for sustainable development of the soybean production in China.

Whole-genome re-sequencing association study on yearling wool traits in Chinese fine-wool sheep

To investigate single nucleotide polymorphism (SNP) loci associated with yearling wool traits of fine-wool sheep for optimizing marker-assisted selection and dissection of the genetic architecture of wool traits, we conducted a genome-wide association study (GWAS) based on the fixed and random model circulating probability unification (FarmCPU) for yearling staple length (YSL), yearling mean fiber diameter (YFD), yearling greasy fleece weight (YGFW), and yearling clean fleece rate (YCFR) by using the whole-genome re-sequenced data (totaling 577 sheep) from the following four fine-wool sheep breeds in China: Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Qinghai fine-wool sheep (QHS), and Aohan fine-wool sheep (AHS). A total of 16 SNPs were detected above the genome-wise significant threshold (P = 5.45E-09), and 79 SNPs were located above the suggestive significance threshold (P = 5.00E-07) from the GWAS results. For YFD and YGFW traits, 7 and 9 SNPs reached the genome-wise significance thresholds, whereas 10 and 12 SNPs reached the suggestive significance threshold, respectively. For YSL and YCFR traits, none of the SNPs reached the genome-wise significance thresholds, whereas 57 SNPs exceeded the suggestive significance threshold. We recorded 14 genes located at the region of ±50-kb near the genome-wise significant SNPs and 59 genes located at the region of ±50-kb near the suggestive significant SNPs. Meanwhile, we used the Average Information Restricted Maximum likelihood algorithm (AI-REML) in the “HIBLUP” package to estimate the heritability and variance components of the four desired yearling wool traits. The estimated heritability values (h²) of YSL, YFD, YGFW, and YCFR were 0.6208, 0.7460, 0.6758, and 0.5559, respectively. We noted that the genetic parameters in this study can be used for fine-wool sheep breeding. The newly detected significant SNPs and the newly identified candidate genes in this study would enhance our understanding of yearling wool formation, and significant SNPs can be applied to genome selection in fine-wool sheep breeding.

[Genetic diversity of protopine-6-hydroxylase in three medicinal Papaver plants]

Sanguinarine is the main active component of the Papaver plants, and protopine-6-hydroxylase(P6 H), involved in the sanguinarine biosynthetic pathway, can oxidize protopine to 6-hydroxyprotopine. The investigation on the diversity of P6 H genes in the medicinal Papaver plants contributes to the acquirement of P6 H with high activity to increase the biosynthesis of sanguinarine. Five P6 H genes in P. somniferum, P. orientale, and P. rhoeas were discovered based on the re-sequencing data of the Papaver species, followed by bioinformatics analysis. With the elongation factor 1α(EF-1α), which exhibits stable expression in the root and stem, as the internal reference gene, the transcription levels of P6H genes in roots and stems of the Papaver plants were detected by real-time fluorescent quantitative PCR. As indicated by the re-sequencing results, there were two genotypes of P6H in P. somniferum and P. orientale, respectively, and only one in P. rhoeas. The bioinformatics analysis showed that the P6 H proteins of the three Papaver plants contained the conserved domain cl12078, which is the characteristic of p450 supergene family, and transmembrane regions. The existence of signal peptide remained verification. Real-time fluorescent quantitative PCR results revealed that the transcription level of P6 H in roots of P. somniferum was about 1.44 times of that in stems(α=0.05). The present study confirmed genetic diversity of P6 H in the three medicinal Papaver plants, which lays a basis for the research on the biosynthesis pathway and mechanism of sanguinarine in Papaver species.

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

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

Comparison of Read Mapping and Variant Calling Tools for the Analysis of Plant NGS Data

High-throughput sequencing technologies have rapidly developed during the past years and have become an essential tool in plant sciences. However, the analysis of genomic data remains challenging and relies mostly on the performance of automatic pipelines. Frequently applied pipelines involve the alignment of sequence reads against a reference sequence and the identification of sequence variants. Since most benchmarking studies of bioinformatics tools for this purpose have been conducted on human datasets, there is a lack of benchmarking studies in plant sciences. In this study, we evaluated the performance of 50 different variant calling pipelines, including five read mappers and ten variant callers, on six real plant datasets of the model organism Arabidopsis thaliana. Sets of variants were evaluated based on various parameters including sensitivity and specificity. We found that all investigated tools are suitable for analysis of NGS data in plant research. When looking at different performance metrics, BWA-MEM and Novoalign were the best mappers and GATK returned the best results in the variant calling step.

QTL and Candidate Gene Identification for Silique Length Based on High-Dense Genetic Map in Brassica napus L

Silique length (SL) is an important yield trait and positively correlates with seeds per silique and seed weight. In the present study, two double haploid (DH) populations, established from crosses Zhongshuang11 × R11 (ZR) and R1 × R2 (RR), containing 280 and 95 DH lines, respectively, were used to map quantitative trait loci (QTL) for SL. A high-dense genetic map from ZR population was constructed comprising 14,658 bins on 19 linkage groups, with map length of 2,198.85 cM and an average marker distance of 0.15 cM. Genetic linkage map from RR population was constructed by using 2,046 mapped markers anchored to 19 chromosomes with 2,217-cM map length and an average marker distance of 1.08 cM. Major QTL qSL_ZR_A09 and qSL_RR_A09b on A09 were identified from ZR and RR populations, respectively. Both QTL could be stably detected in four environments. QTL qSL_RR_A09b and qSL_ZR_A09 were located on 68.5-70.8 cM and 91.33-91.94 cM interval with R² values of 14.99-39.07% and 15.00-20.36% in RR and ZR populations, respectively. Based on the physical positions of single nucleotide polymorphism (SNP) markers flanking qSL_ZR_A09 and gene annotation in Arabidopsis, 26 genes were identified with SNP/Indel variation between parents and two genes (BnaA09g41180D and BnaA09g41380D) were selected as the candidate genes. Expression analysis further revealed BnaA09g41180D, encoding homologs of Arabidopsis fasciclin-like arabinogalactan proteins (FLA3), as the most promising candidate gene for qSL_ZR_A09. The QTL identification and candidate gene analysis will provide new insight into the genomic regions controlling SL in Brassica napus as well as candidate genes underlying the QTL.

Sequence Analysis of APOA5 Among the Kuwaiti Population Identifies Association of rs2072560, rs2266788, and rs662799 With TG and VLDL Levels

Common variants of Apolipoprotein A5 (APOA5) have been associated with lipid levels yet very few studies have reported full sequence data from various ethnic groups. The purpose of this study was to analyse the full APOA5 gene sequence to identify variants in 100 healthy Kuwaitis of Arab ethnicities and assess their association with variation in lipid levels in a cohort of 733 samples. Sanger method was used in the direct sequencing of the full 3.7 Kb APOA5 and multiple sequence alignment was used to identify variants. The complete APOA5 sequence in Kuwaiti Arabs has been deposited in GenBank (KJ401315). A total of 20 reported single nucleotide polymorphisms (SNPs) were identified. Two novel SNPs were also identified: a synonymous 2197G>A polymorphism at genomic position 116661525 and a 3′ UTR 3222 C>T polymorphism at genomic position 116660500 based on human genome assembly GRCh37/hg:19. Five SNPs along with the two novel SNPs were selected for validation in the cohort. Association of those SNPs with lipid levels was tested and minor alleles of three SNPs (rs2072560, rs2266788, and rs662799) were found significantly associated with TG and VLDL levels. This is the first study to report the full APOA5 sequence and SNPs in an Arab ethnic group. Analysis of the variants identified and comparison to other populations suggests a distinctive genetic component in Arabs. The positive association observed for rs2072560 and rs2266788 with TG and VLDL levels confirms their role in lipid metabolism.

A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo

The melon (Cucumis melo) genome and genetic maps with hundreds to thousands of single nucleotide polymorphism markers were recently released. However, a high-resolution genetic map was lacking. Gummy stem blight (Gsb) is a destructive disease responsible for considerable economic losses during melon production. We herein describe the development of an ultra-dense genetic map consisting of 12,932 recombination bin markers covering 1,818 cM, with an average distance of 0.17 cM between adjacent tags. A comparison of the genetic maps for melon, watermelon, and cucumber revealed chromosome-level syntenic relationships and recombination events among the three Cucurbitaceae species. Our genetic map was useful for re-anchoring the genome scaffolds of melon. More than 92% assembly was anchored to 12 pseudo-chromosomes and 90% of them were oriented. Furthermore, 1,135 recombination hotspots revealed an unbalanced recombination rate across the melon genome. Genetic analyses of the Gsb-resistant and -susceptible lines indicated the resistance phenotype is mediated by a single dominant gene. We identified Gsb-resistance gene candidates in a 108-kb region on pseudo-chromosome 4. Our findings verify the utility of an ultra-dense genetic map for mapping a gene of interest, and for identifying new disease resistant genes.

Genome Analysis Identified Novel Candidate Genes for Ascochyta Blight Resistance in Chickpea Using Whole Genome Re-sequencing Data

Ascochyta blight (AB) is a fungal disease that can significantly reduce chickpea production in Australia and other regions of the world. In this study, 69 chickpea genotypes were sequenced using whole genome re-sequencing (WGRS) methods. They included 48 Australian varieties differing in their resistance ranking to AB, 16 advanced breeding lines from the Australian chickpea breeding program, four landraces, and one accession representing the wild chickpea species Cicer reticulatum. More than 800,000 single nucleotide polymorphisms (SNPs) were identified. Population structure analysis revealed relatively narrow genetic diversity amongst recently released Australian varieties and two groups of varieties separated by the level of AB resistance. Several regions of the chickpea genome were under positive selection based on Tajima's D test. Both Fst genome- scan and genome-wide association studies (GWAS) identified a 100 kb region (AB4.1) on chromosome 4 that was significantly associated with AB resistance. The AB4.1 region co-located to a large QTL interval of 7 Mb∼30 Mb identified previously in three different mapping populations which were genotyped at relatively low density with SSR or SNP markers. The AB4.1 region was validated by GWAS in an additional collection of 132 advanced breeding lines from the Australian chickpea breeding program, genotyped with approximately 144,000 SNPs. The reduced level of nucleotide diversity and long extent of linkage disequilibrium also suggested the AB4.1 region may have gone through selective sweeps probably caused by selection of the AB resistance trait in breeding. In total, 12 predicted genes were located in the AB4.1 QTL region, including those annotated as: NBS-LRR receptor-like kinase, wall-associated kinase, zinc finger protein, and serine/threonine protein kinases. One significant SNP located in the conserved catalytic domain of a NBS-LRR receptor-like kinase led to amino acid substitution. Transcriptional analysis using qPCR showed that some predicted genes were significantly induced in resistant lines after inoculation compared to non-inoculated plants. This study demonstrates the power of combining WGRS data with relatively simple traits to rapidly develop "functional makers" for marker-assisted selection and genomic selection.

Multi-lineages of Shiikuwasha (Citrus depressa Hayata) evaluated by using whole chloroplast genome sequences and its bio-diversity in Okinawa, Japan

Shiikuwasha (Citrus depressa Hayata) is distributed from the South-west of the Japanese archipelago to Taiwan. In this study, re-sequencing against the orange (C. sinensis (L.) Osbeck) chloroplast genome was applied to one superior landrace of Shiikuwasha cultivated in Oku ward, Okinawa, Japan. The chloroplast genome of the landrace was estimated to comprise 160,118 bp, including 48 indels and 71 nucleotide substitutions against the reference genome. The presumptive chloroplast indels were confirmed by subsequent experiments, and these identified multiple maternal lineages among other landraces. Some of the orange SSR markers were available for genotyping of other superior landraces and were able to distinguish among them. These molecular markers were then applied for evaluation of genetic diversity among wild and cultivated Shiikuwasha accessions. Except for Oku ward, the cultivated populations were found to have lost their genetic diversity in comparison with wild populations. Groves in Oku ward maintained, or showed even higher genetic diversity than wild accessions in the surrounding areas by the force of villagers.

Whole genome re-sequencing reveals genome-wide variations among parental lines of 16 mapping populations in chickpea (Cicer arietinum L.)

BACKGROUND

Chickpea (Cicer arietinum L.) is the second most important grain legume cultivated by resource poor farmers in South Asia and Sub-Saharan Africa. In order to harness the untapped genetic potential available for chickpea improvement, we re-sequenced 35 chickpea genotypes representing parental lines of 16 mapping populations segregating for abiotic (drought, heat, salinity), biotic stresses (Fusarium wilt, Ascochyta blight, Botrytis grey mould, Helicoverpa armigera) and nutritionally important (protein content) traits using whole genome re-sequencing approach.

RESULTS

A total of 192.19 Gb data, generated on 35 genotypes of chickpea, comprising 973.13 million reads, with an average sequencing depth of ~10 X for each line. On an average 92.18 % reads from each genotype were aligned to the chickpea reference genome with 82.17 % coverage. A total of 2,058,566 unique single nucleotide polymorphisms (SNPs) and 292,588 Indels were detected while comparing with the reference chickpea genome. Highest number of SNPs were identified on the Ca4 pseudomolecule. In addition, copy number variations (CNVs) such as gene deletions and duplications were identified across the chickpea parental genotypes, which were minimum in PI 489777 (1 gene deletion) and maximum in JG 74 (1,497). A total of 164,856 line specific variations (144,888 SNPs and 19,968 Indels) with the highest percentage were identified in coding regions in ICC 1496 (21 %) followed by ICCV 97105 (12 %). Of 539 miscellaneous variations, 339, 138 and 62 were inter-chromosomal variations (CTX), intra-chromosomal variations (ITX) and inversions (INV) respectively.

CONCLUSION

Genome-wide SNPs, Indels, CNVs, PAVs, and miscellaneous variations identified in different mapping populations are a valuable resource in genetic research and helpful in locating genes/genomic segments responsible for economically important traits. Further, the genome-wide variations identified in the present study can be used for developing high density SNP arrays for genetics and breeding applications.

2' and 'Qinguan' Varieties of Apple (Malus domestica Borkh.)

Apple (Malus domestica Borkh.) is a commercially important fruit worldwide. Detailed information on genomic DNA polymorphisms, which are important for understanding phenotypic traits, is lacking for the apple. We re-sequenced two elite apple varieties, 'Nagafu No. 2' and 'Qinguan,' which have different characteristics. We identified many genomic variations, including 2,771,129 single nucleotide polymorphisms (SNPs), 82,663 structural variations (SVs), and 1,572,803 insertion/deletions (INDELs) in 'Nagafu No. 2' and 2,262,888 SNPs, 63,764 SVs, and 1,294,060 INDELs in 'Qinguan.' The 'SNP,' 'INDEL,' and 'SV' distributions were non-random, with variation-rich or -poor regions throughout the genomes. In 'Nagafu No. 2' and 'Qinguan' there were 171,520 and 147,090 non-synonymous SNPs spanning 23,111 and 21,400 genes, respectively; 3,963 and 3,196 SVs in 3,431 and 2,815 genes, respectively; and 1,834 and 1,451 INDELs in 1,681 and 1,345 genes, respectively. Genetic linkage maps of 190 flowering genes associated with multiple flowering pathways in 'Nagafu No. 2,' 'Qinguan,' and 'Golden Delicious,' identified complex regulatory mechanisms involved in floral induction, flower bud formation, and flowering characteristics, which might reflect the genetic variation of the flowering genes. Expression profiling of key flowering genes in buds and leaves suggested that the photoperiod and autonomous flowering pathways are major contributors to the different floral-associated traits between 'Nagafu No. 2' and 'Qinguan.' The genome variation data provided a foundation for the further exploration of apple diversity and gene-phenotype relationships, and for future research on molecular breeding to improve apple and related species.

Prediction of Genes Related to Positive Selection Using Whole-Genome Resequencing in Three Commercial Pig Breeds

Selective sweep can cause genetic differentiation across populations, which allows for the identification of possible causative regions/genes underlying important traits. The pig has experienced a long history of allele frequency changes through artificial selection in the domestication process. We obtained an average of 329,482,871 sequence reads for 24 pigs from three pig breeds: Yorkshire (n = 5), Landrace (n = 13), and Duroc (n = 6). An average read depth of 11.7 was obtained using whole-genome resequencing on an Illumina HiSeq2000 platform. In this study, cross-population extended haplotype homozygosity and cross-population composite likelihood ratio tests were implemented to detect genes experiencing positive selection for the genome-wide resequencing data generated from three commercial pig breeds. In our results, 26, 7, and 14 genes from Yorkshire, Landrace, and Duroc, respectively were detected by two kinds of statistical tests. Significant evidence for positive selection was identified on genes ST6GALNAC2 and EPHX1 in Yorkshire, PARK2 in Landrace, and BMP6, SLA-DQA1, and PRKG1 in Duroc.These genes are reportedly relevant to lactation, reproduction, meat quality, and growth traits. To understand how these single nucleotide polymorphisms (SNPs) related positive selection affect protein function, we analyzed the effect of non-synonymous SNPs. Three SNPs (rs324509622, rs80931851, and rs80937718) in the SLA-DQA1 gene were significant in the enrichment tests, indicating strong evidence for positive selection in Duroc. Our analyses identified genes under positive selection for lactation, reproduction, and meat-quality and growth traits in Yorkshire, Landrace, and Duroc, respectively.

Application of whole genome re-sequencing data in the development of diagnostic DNA markers tightly linked to a disease-resistance locus for marker-assisted selection in lupin (Lupinus angustifolius)

BACKGROUND

Molecular marker-assisted breeding provides an efficient tool to develop improved crop varieties. A major challenge for the broad application of markers in marker-assisted selection is that the marker phenotypes must match plant phenotypes in a wide range of breeding germplasm. In this study, we used the legume crop species Lupinus angustifolius (lupin) to demonstrate the utility of whole genome sequencing and re-sequencing on the development of diagnostic markers for molecular plant breeding.

RESULTS

Nine lupin cultivars released in Australia from 1973 to 2007 were subjected to whole genome re-sequencing. The re-sequencing data together with the reference genome sequence data were used in marker development, which revealed 180,596 to 795,735 SNP markers from pairwise comparisons among the cultivars. A total of 207,887 markers were anchored on the lupin genetic linkage map. Marker mining obtained an average of 387 SNP markers and 87 InDel markers for each of the 24 genome sequence assembly scaffolds bearing markers linked to 11 genes of agronomic interest. Using the R gene PhtjR conferring resistance to phomopsis stem blight disease as a test case, we discovered 17 candidate diagnostic markers by genotyping and selecting markers on a genetic linkage map. A further 243 candidate diagnostic markers were discovered by marker mining on a scaffold bearing non-diagnostic markers linked to the PhtjR gene. Nine out from the ten tested candidate diagnostic markers were confirmed as truly diagnostic on a broad range of commercial cultivars. Markers developed using these strategies meet the requirements for broad application in molecular plant breeding.

CONCLUSIONS

We demonstrated that low-cost genome sequencing and re-sequencing data were sufficient and very effective in the development of diagnostic markers for marker-assisted selection. The strategies used in this study may be applied to any trait or plant species. Whole genome sequencing and re-sequencing provides a powerful tool to overcome current limitations in molecular plant breeding, which will enable plant breeders to precisely pyramid favourable genes to develop super crop varieties to meet future food demands.

Nonsyndromic cleft lip with or without cleft palate: Increased burden of rare variants within Gremlin-1, a component of the bone morphogenetic protein 4 pathway

BACKGROUND

The genes Gremlin-1 (GREM1) and Noggin (NOG) are components of the bone morphogenetic protein 4 pathway, which has been implicated in craniofacial development. Both genes map to recently identified susceptibility loci (chromosomal region 15q13, 17q22) for nonsyndromic cleft lip with or without cleft palate (nsCL/P). The aim of the present study was to determine whether rare variants in either gene are implicated in nsCL/P etiology.

METHODS

The complete coding regions, untranslated regions, and splice sites of GREM1 and NOG were sequenced in 96 nsCL/P patients and 96 controls of Central European ethnicity. Three burden and four nonburden tests were performed. Statistically significant results were followed up in a second case-control sample (n = 96, respectively). For rare variants observed in cases, segregation analyses were performed.

RESULTS

In NOG, four rare sequence variants (minor allele frequency < 1%) were identified. Here, burden and nonburden analyses generated nonsignificant results. In GREM1, 33 variants were identified, 15 of which were rare. Of these, five were novel. Significant p-values were generated in three nonburden analyses. Segregation analyses revealed incomplete penetrance for all variants investigated.

CONCLUSION

Our study did not provide support for NOG being the causal gene at 17q22. However, the observation of a significant excess of rare variants in GREM1 supports the hypothesis that this is the causal gene at chr. 15q13. Because no single causal variant was identified, future sequencing analyses of GREM1 should involve larger samples and the investigation of regulatory elements.

Controlling for contamination in re-sequencing studies with a reproducible web-based phylogenetic approach

Polymorphism discovery is a routine application of next-generation sequencing technology where multiple samples are sent to a service provider for library preparation, subsequent sequencing, and bioinformatic analyses. The decreasing cost and advances in multiplexing approaches have made it possible to analyze hundreds of samples at a reasonable cost. However, because of the manual steps involved in the initial processing of samples and handling of sequencing equipment, cross-contamination remains a significant challenge. It is especially problematic in cases where polymorphism frequencies do not adhere to diploid expectation, for example, heterogeneous tumor samples, organellar genomes, as well as during bacterial and viral sequencing. In these instances, low levels of contamination may be readily mistaken for polymorphisms, leading to false results. Here we describe practical steps designed to reliably detect contamination and uncover its origin, and also provide new, Galaxy-based, readily accessible computational tools and workflows for quality control. All results described in this report can be reproduced interactively on the web as described at http://usegalaxy.org/contamination.

Reduced representation approaches to interrogate genome diversity in large repetitive plant genomes

Technology and software improvements in the last decade now provide methodologies to access the genome sequence of not only a single accession, but also multiple accessions of plant species. This provides a means to interrogate species diversity at the genome level. Ample diversity among accessions in a collection of species can be found, including single-nucleotide polymorphisms, insertions and deletions, copy number variation and presence/absence variation. For species with small, non-repetitive rich genomes, re-sequencing of query accessions is robust, highly informative, and economically feasible. However, for species with moderate to large sized repetitive-rich genomes, technical and economic barriers prevent en masse genome re-sequencing of accessions. Multiple approaches to access a focused subset of loci in species with larger genomes have been developed, including reduced representation sequencing, exome capture and transcriptome sequencing. Collectively, these approaches have enabled interrogation of diversity on a genome scale for large plant genomes, including crop species important to worldwide food security.

Improvement of marker-based predictability of Apparent Amylose Content in japonica rice through GBSSI allele mining

BACKGROUND

Apparent Amylose Content (AAC), regulated by the Waxy gene, represents the key determinant of rice cooking properties. In occidental countries high AAC rice represents the most requested market class but the availability of molecular markers allowing specific selection of high AAC varieties is limited.

RESULTS

In this study, the effectiveness of available molecular markers in predicting AAC was evaluated in a collection of 127 rice accessions (125 japonica ssp. and 2 indica ssp.) characterized by AAC values from glutinous to 26%. The analyses highlighted the presence of several different allelic patterns identifiable by a few molecular markers, and two of them, i.e., the SNPs at intron1 and exon 6, were able to explain a maximum of 79.5% of AAC variation. However, the available molecular markers haplotypes did not provide tools for predicting accessions with AAC higher than 24.5%. To identify additional polymorphisms, the re-sequencing of the Waxy gene and 1kbp of the putative upstream regulatory region was performed in 21 genotypes representing all the AAC classes identified. Several previously un-characterized SNPs were identified and four of them were used to develop dCAPS markers.

CONCLUSIONS

The addition of the SNPs newly identified slightly increased the AAC explained variation and allowed the identification of a haplotype almost unequivocally associated to AAC higher than 24.5%. Haplotypes at the waxy locus were also associated to grain length and length/width (L/W) ratio. In particular, the SNP at the first intron, which identifies the Wxa and Wxb alleles, was associated with differences in the width of the grain, the L/W ratio and the length of the kernel, most likely as a result of human selection.

Small MAF genes variants and chronic myeloid leukemia

Chronic myeloid leukemia (CML) is one of the most frequent hematological neoplasia worldwide. The abnormal accumulation of reactive oxygen species may be an important factor in CML development. The transcription factor NRF2 can regulate the transcription of a battery of antioxidant and detoxificant genes after heterodimerizing with small-Maf proteins. Although the participation of NRF2 in the development of chronic degenerative diseases has been thoroughly studied, the role of small-Maf genes has not been documented. We have identified polymorphisms in the three MAF genes (F, G and K) and assessed their association with CML. Over 266 subjects with CML and 399 unrelated healthy donors have been studied. After sequencing each MAF gene by Sanger technology, we found 17 variants in MAFF gene, eight in MAFG and seven in MAFK. In the case-control study, the homozygote genotype CC for the rs9610915 SNP of MAFF was significantly associated with CML. The frequency of the ACC haplotype from MAFK was significantly lower than controls. After stratification by gender, the ACC and GTG haplotypes were associated only with males with CML. These novel data suggest an association between MAFF and MAFG and the development of CML.

Peeling back the evolutionary layers of molecular mechanisms responsive to exercise-stress in the skeletal muscle of the racing horse

The modern horse (Equus caballus) is the product of over 50 million yrs of evolution. The athletic abilities of the horse have been enhanced during the past 6000 yrs under domestication. Therefore, the horse serves as a valuable model to understand the physiology and molecular mechanisms of adaptive responses to exercise. The structure and function of skeletal muscle show remarkable plasticity to the physical and metabolic challenges following exercise. Here, we reveal an evolutionary layer of responsiveness to exercise-stress in the skeletal muscle of the racing horse. We analysed differentially expressed genes and their co-expression networks in a large-scale RNA-sequence dataset comparing expression before and after exercise. By estimating genome-wide dN/dS ratios using six mammalian genomes, and FST and iHS using re-sequencing data derived from 20 horses, we were able to peel back the evolutionary layers of adaptations to exercise-stress in the horse. We found that the oldest and thickest layer (dN/dS) consists of system-wide tissue and organ adaptations. We further find that, during the period of horse domestication, the older layer (FST) is mainly responsible for adaptations to inflammation and energy metabolism, and the most recent layer (iHS) for neurological system process, cell adhesion, and proteolysis.

SNP-based analysis of genetic diversity in anther-derived rice by whole genome sequencing

BACKGROUND

Anther culture has advantage to obtain a homozygous progeny by induced doubling of haploid chromosomes and to improve selection efficiency for invaluable agronomical traits. Therefore, anther culturing is widely utilized to breed new varieties and to induce genetic variations in several crops including rice. Genome sequencing technologies allow the detection of a massive number of DNA polymorphism such as SNPs and Indels between closely related cultivars. These DNA polymorphisms permit the rapid identification of genetic diversity among cultivars and genomic locations of heritable traits. To estimate sequence diversity derived from anther culturing, we performed whole-genome resequencing of five Korean rice accessions, including three anther culture lines (BLB, HY-04 and HY-08), their progenitor cultivar (Hwayeong), and an additional japonica cultivar (Dongjin).

RESULTS

A total of 1,165 × 106 raw reads were generated with over 58× coverage that detected 1,154,063 DNA polymorphisms between the Korean rice accessions and Nipponbare. We observed that in Hwayeong and its progenies, 0.64 SNP was found per one kb of Nipponbare genome, while Dongjin, bred by a conventional breeding method, had a lower number of SNPs (0.45 SNP/kb). Among 1,154,063 DNA polymorphisms, 29,269 non-synonymous SNPs located on 30,013 genes and these genes were functionally classified based on gene ontology (GO). We also analyzed line-specific SNPs which were estimated 1 ~ 3% of the total SNPs. The frequency of non-synonymous SNPs in each accession ranged from 26 SNPs in Hwayeong to 214 SNPs in HY-04.

CONCLUSIONS

The genetic difference we detected between the progenies derived from anther culture and their mother cultivar is due to somaclonal variation during tissue culture process, such as karyotype change, chromosome rearrangement, gene amplification and deletion, transposable element, and DNA methylation. Detection of genome-wide DNA polymorphisms by high-throughput sequencer enabled to identify sequence diversity derived from anther culturing and genomic locations of heritable traits. Furthermore, it will provide an invaluable resource to identify molecular markers and genes associated with diverse traits of agronomical importance.

Re-sequencing of mitochondrial genes in a standard rice cultivar Nipponbare

BACKGROUND

Genomic sequence of a rice cultivar Nipponbare has been often used as a reference sequence since a whole-genomic sequence was first determined in 2005 by the International Rice Genome Sequencing Project. As for mitochondrial genomic sequence of Nipponbare, two groups have deposited their sequences into DDBJ/EMBL/GenBank under the accession numbers BA000029 and DQ167400. However, there are 19 discrepancies in the nucleotide sequences of 7 genes between BA000029 and DQ167400.

FINDINGS

We performed PCR to amplify these 7 genes and to perform direct sequencing. Nucleotides of the discrepant sites were all identical to those in DQ167400.1. The sequence in BA000029.3 is thought to contain sequencing errors.

CONCLUSION

Nucleotide sequences of the mitochondrial genes in BA000029.3 need to be updated using the data in this study when used as a reference genome.

Sequencing technologies and genome sequencing

The high-throughput - next generation sequencing (HT-NGS) technologies are currently the hottest topic in the field of human and animals genomics researches, which can produce over 100 times more data compared to the most sophisticated capillary sequencers based on the Sanger method. With the ongoing developments of high throughput sequencing machines and advancement of modern bioinformatics tools at unprecedented pace, the target goal of sequencing individual genomes of living organism at a cost of $1,000 each is seemed to be realistically feasible in the near future. In the relatively short time frame since 2005, the HT-NGS technologies are revolutionizing the human and animal genome researches by analysis of chromatin immunoprecipitation coupled to DNA microarray (ChIP-chip) or sequencing (ChIP-seq), RNA sequencing (RNA-seq), whole genome genotyping, genome wide structural variation, de novo assembling and re-assembling of genome, mutation detection and carrier screening, detection of inherited disorders and complex human diseases, DNA library preparation, paired ends and genomic captures, sequencing of mitochondrial genome and personal genomics. In this review, we addressed the important features of HT-NGS like, first generation DNA sequencers, birth of HT-NGS, second generation HT-NGS platforms, third generation HT-NGS platforms: including single molecule Heliscope™, SMRT™ and RNAP sequencers, Nanopore, Archon Genomics X PRIZE foundation, comparison of second and third HT-NGS platforms, applications, advances and future perspectives of sequencing technologies on human and animal genome research.