Search

Identifying alleles associated with adaptation to new environments will advance our understanding of evolution from the molecular level. Previous studies have found that the Populus davidiana southwest population in East Asia has differentiated from other populations in the range. We aimed to evaluate the contributions of the ancestral-state bases (ASBs) relative to derived bases (DBs) in the local adaptation of P. davidiana in the Yunnan–Guizhou Plateau from a quantitative perspective using whole-genome re-sequencing data from 90 P. davidiana samples from three regions across the species range. Our results showed that the uplift of the Qinghai–Tibet Plateau during the Neogene and associated climate fluctuations during the Middle Pleistocene were likely an important factor in the early divergence of P. davidiana. Highly differentiated genomic regions between populations were inferred to have undergone strong linked natural selection, and ASBs are the chief means by which populations of P. davidiana adapt to novel environmental conditions; however, when adapting to regions with high environmental differences relative to the ancestral range, the proportion of DBs was significantly higher than that of background regions, as ASBs are insufficient to cope with these environments. Finally, a number of genes were identified in the outlier region.

BackgroundWe have previously developed a rice-based oral vaccine against cholera diarrhea, MucoRice-CTB. UsingAgrobacterium-mediated co-transformation, we produced the selection marker–free MucoRice-CTB line 51A, which has three copies of the cholera toxin B subunit (CTB) gene and two copies of an RNAi cassette inserted into the rice genome. We determined the sequence and location of the transgenes on rice chromosomes 3 and 12. The expression of alpha-amylase/trypsin inhibitor, a major allergen protein in rice, is lower in this line than in wild-type rice. Line 51A was self-pollinated for five generations to fix the transgenes, and the seeds of the sixth generation produced by T5 plants were defined as the master seed bank (MSB). T6 plants were grown from part of the MSB seeds and were self-pollinated to produce T7 seeds (next seed bank; NSB). NSB was examined and its whole genome and proteome were compared with those of MSB.ResultsWe re-sequenced the transgenes of NSB and MSB and confirmed the positions of the three CTB genes inserted into chromosomes 3 and 12. The DNA sequences of the transgenes were identical between NSB and MSB. Using whole-genome sequencing, we compared the genome sequences of three NSB with three MSB samples, and evaluated the effects of SNPs and genomic structural variants by clustering. No functionally important mutations (SNPs, translocations, deletions, or inversions of genic regions on chromosomes) between NSB and MSB samples were detected. Analysis of salt-soluble proteins from NSB and MSB samples by shot-gun MS/MS detected no considerable differences in protein abundance. No difference in the expression pattern of storage proteins and CTB in mature seeds of NSB and MSB was detected by immuno-fluorescence microscopy.ConclusionsAll analyses revealed no considerable differences between NSB and MSB samples. Therefore, NSB can be used to replace MSB in the near future.

Despite the global biodiversity of terrestrial gastropods and their ecological and economic importance, the genomic basis of ecological adaptation and speciation in land snail taxa is still largely unknown. Here, we combined whole-genome re-sequencing with population genomics to evaluate the historical demography and the speciation process of two closely related species of land snails from western Europe, Candidula unifasciata and C. rugosiuscula. Historical demographic analysis indicated fluctuations in the size of ancestral populations, probably driven by Pleistocene climatic fluctuations. Although the current population distributions of both species do not overlap, our approximate Bayesian computation model selection approach on several speciation scenarios suggested that gene flow has occurred throughout the divergence process until recently. Positively selected genes diverging early in the process were associated with intragenomic and cyto-nuclear incompatibilities, respectively, potentially fostering reproductive isolation as well as ecological divergence. Our results suggested that the speciation between species entails complex processes involving both gene flow and ecological speciation, and that further research based on whole-genome data can provide valuable understanding on species divergence. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum This work was funded by LOEWE-Centre for Translational Biodiversity Genomics (LOEWE-TBG). L.J.C. was funded by a Post-doctoral Fellowship awarded by the Department of Education, Universities and Research of the Basque Government (Ref.: POS-2018-1-0012)

Drought tolerance is one of the most important targets for wheat breeding. Previous population genetics studies have uncovered 20 large-effect quantitative trait loci (QTLs) that contribute to stress-responsive agronomic traits. Here, we identified 19,035,814 single nucleotide polymorphisms and 719,049 insertion/deletion variations in the genomes of two popular winter wheat cultivars, Lu-Mai 14 and Han-Xuan 10, using a whole-genome re-sequencing assay. There were 4972 loss-of-function mutations carried by protein-coding genes, such as CCA1/LHY, AGO1, ABI3/VP1, EIN3, TPP, and ARFs. We carried out a time-course abscisic acid (ABA)-treatment experiment and profiled 61,251 expressed genes in the roots using a strand-specific RNA sequencing approach. A large number of genes showed time-point specific and/or cultivar-preferential responsive expression patterns. Gene ontology enrichment analysis revealed that ABA-responsive genes were associated with stress-related functions. Among the 20 QTLs, we uncovered 306 expressed genes with high- and/or moderate-effect variations and 472 differentially expressed genes. Detailed analysis and verification of the homozygous genomic variations in the candidate genes encoding sulfotransferase, proteinase, kinase, nitrate transporter, and transcription factors suggested previously unexpected pathways associated with abiotic stress responses in wheat.

The Chongming white goat (CM) is an indigenous goat breed exhibits unique traits that are adapted to the local environment and artificial selection. By performing whole-genome re-sequencing, we generated 14–20× coverage sequences from 10 domestic goat breeds to explore the genomic characteristics and selection signatures of the CM breed. We identified a total of 23,508,551 single-nucleotide polymorphisms (SNPs) and 2,830,800 insertion–deletion mutations (indels) after read mapping and variant calling. We further specifically identified 1.2% SNPs (271,713) and 0.9% indels (24,843) unique to the CM breed in comparison with the other nine goat breeds. Missense (SIFT < 0.05), frameshift, splice-site, start-loss, stop-loss, and stop-gain variants were identified in 183 protein-coding genes of the CM breed. Of the 183, 36 genes, including AP4E1, FSHR, COL11A2, and DYSF, are involved in phenotype ontology terms related to the nervous system, short stature, and skeletal muscle morphology. Moreover, based on genome-wide FST and pooled heterozygosity (Hp) calculation, we further identified selection signature genes between the CM and the other nine goat breeds. These genes are significantly associated with the nervous system (C2CD3, DNAJB13, UCP2, ZMYND11, CEP126, SCAPER, and TSHR), growth (UCP2, UCP3, TSHR, FGFR1, ERLIN2, and ZNF703), and coat color (KITLG, ASIP, AHCY, RALY, and MC1R). Our results suggest that the CM breed may be differentiated from other goat breeds in terms of nervous system owing to natural or artificial selection. The whole-genome analysis provides an improved understanding of genetic diversity and trait exploration for this indigenous goat breed.

Whole-genome sequences of four EMS (ethyl methanesulfonate)-induced eggplant mutants were analyzed to identify genome-wide mutations. In total, 173.01 GB of paired-end reads were obtained for four EMS-induced mutants and (WT) wild type and 1,076,010 SNPs (single nucleotide polymorphisms) and 183,421 indels were identified. The most common mutation type was C/G to T/A transitions followed by A/T to G/C transitions. The mean densities were one SNP per 1.3 to 2.6 Mb. The effect of mutations on gene function was annotated and only 7.2% were determined to be deleterious. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis showed 10 and 11 genes, which were nonsynonymous mutation or frameshift deletion in 48-5 and L6-5 involved in the anthocyanin biosynthesis or flavone and flavonol biosynthesis. QRT-PCR results showed that only the Sme2.5_06210.1_g00004.1, which was annotated as UFGT (Flavonoid galactosidase transferase), expression significantly decreased in the L6-5 mutant compared with the WT. Also, the Sme2.5_06210.1_g00004.1 expression was lower in the colorless eggplant compared with colorful eggplant in the natural eggplant cultivar. These results suggest that Sme2.5_06210.1_g00004.1 may play a key role in eggplant anthocyanin synthesis.

There is increasing evidence that several genes are associated with an increased risk of type 2 diabetes (T2D); genome-wide association investigations and whole-genome re-sequencing investigations offer a useful approach for the identification of genes involved in common human diseases. To further investigate which polymorphisms confer susceptibility to T2D, the present study screened for high-contribution susceptibility gene variants Chinese patients with T2D using whole-genome re-sequencing with DNA pooling. In total, 100 Chinese individuals with T2D and 100 healthy Chinese individuals were analyzed using whole-genome re-sequencing using DNA pooling. To minimize the likelihood of systematic bias in sampling, paired-end libraries with an insert size of 500 bp were prepared for in T2D in all samples, which were then subjected to whole-genome sequencing. Each library contained four lanes. The average sequencing depth was 35.70. In the present study, 1.36 GB of clean sequence data were generated, and the resulting calculated T2D genome consensus sequence covered 99.88% of the hg19 sequence. A total of 3,974,307 single nucleotide polymorphisms were identified, of which 99.88% were in the dbSNP database. The present study also found 642,189 insertions and deletions, 5,590 structure variants (SVs), 4,713 copy number variants (CNVs) and 13,049 single nucleotide variants. A total of 1,884 somatic CNVs and 74 somatic SVs were significantly different between the cases and controls. Therefore, the present study provided validation of whole-genome re-sequencing using the DNA pooling approach. It also generated a whole-genome re-sequencing genotype database for future investigations of T2D.

Summary To map resistance genes for Fusarium wilt (FW) and sterility mosaic disease (SMD) in pigeonpea, sequencing‐based bulked segregant analysis (Seq‐BSA) was used. Resistant (R) and susceptible (S) bulks from the extreme recombinant inbred lines of ICPL 20096 × ICPL 332 were sequenced. Subsequently, SNP index was calculated between R‐ and S‐bulks with the help of draft genome sequence and reference‐guided assembly of ICPL 20096 (resistant parent). Seq‐BSA has provided seven candidate SNPs for FW and SMD resistance in pigeonpea. In parallel, four additional genotypes were re‐sequenced and their combined analysis with R‐ and S‐bulks has provided a total of 8362 nonsynonymous (ns) SNPs. Of 8362 nsSNPs, 60 were found within the 2‐Mb flanking regions of seven candidate SNPs identified through Seq‐BSA. Haplotype analysis narrowed down to eight nsSNPs in seven genes. These eight nsSNPs were further validated by re‐sequencing 11 genotypes that are resistant and susceptible to FW and SMD. This analysis revealed association of four candidate nsSNPs in four genes with FW resistance and four candidate nsSNPs in three genes with SMD resistance. Further, In silico protein analysis and expression profiling identified two most promising candidate genes namely C.cajan_01839 for SMD resistance and C.cajan_03203 for FW resistance. Identified candidate genomic regions/SNPs will be useful for genomics‐assisted breeding in pigeonpea.

With the aim of understanding relationship between genetic and phenotypic variations in cultivated tomato, single nucleotide polymorphism (SNP) markers covering the whole genome of cultivated tomato were developed and genome-wide association studies (GWAS) were performed. The whole genomes of six tomato lines were sequenced with the ABI-5500xl SOLiD sequencer. Sequence reads covering ~13.7× of the genome for each line were obtained, and mapped onto tomato reference genomes (SL2.40) to detect ~1.5 million SNP candidates. Of the identified SNPs, 1.5% were considered to confer gene functions. In the subsequent Illumina GoldenGate assay for 1536 SNPs, 1293 SNPs were successfully genotyped, and 1248 showed polymorphisms among 663 tomato accessions. The whole-genome linkage disequilibrium (LD) analysis detected highly biased LD decays between euchromatic (58 kb) and heterochromatic regions (13.8 Mb). Subsequent GWAS identified SNPs that were significantly associated with agronomical traits, with SNP loci located near genes that were previously reported as candidates for these traits. This study demonstrates that attractive loci can be identified by performing GWAS with a large number of SNPs obtained from re-sequencing analysis.

Heading date is one of the most important traits in rice breeding. It is governed by multiple genes, including known quantitative trait loci (QTLs). In general, almost all japonica cultivars, including Nipponbare, head early under short-day (SD) conditions, but some indica cultivars, including Kasalath, head late. To explain this difference, we identified QTLs controlling heading date under SD conditions. We used NILs, CSSLs, and BILs from a cross between Nipponbare and Kasalath, and evaluated days to heading (DTH) under SD conditions. No NILs or CSSLs showed late heading, but two BILs (BIL-55 and BIL-78) had almost the same DTH as Kasalath. We developed an F2 population from a cross between BIL-55 and Nipponbare and performed QTL analysis using SSR markers. The late-heading phenotype was controlled by two known genes and at least two novel QTLs on chromosomes 4 and 6, named qDTH4.5 and qDTH6.3. These QTLs were confirmed by QTL-seq. The QTLs and polymorphisms detected here will provide useful information for further genetic studies and breeding under SD conditions at lower latitudes.

Despite the global economic and ecological importance of forest trees, the genomic basis of differential adaptation and speciation in tree species is still poorly understood.Populus tremulaandP. tremuloidesare two of the most widespread tree species in the Northern Hemisphere. Using whole-genome re-sequencing data of 24P. tremulaand 22P. tremuloidesindividuals, we find that the two species diverged ∼2.2-3.1 million years ago, coinciding with the severing of the Bering land bridge and the onset of dramatic climatic oscillations during the Pleistocene. Both species have experienced substantial population expansions following long-term declines after species divergence. We detect widespread and heterogeneous genomic differentiation between species, and in accordance with the expectation of allopatric speciation, coalescent simulations suggest that neutral evolutionary processes can account for most of the observed patterns of genomic differentiation. However, there is an excess of regions exhibiting extreme differentiation relative to those expected under demographic simulations, which is indicative of the action of natural selection. Overall genetic differentiation is negatively associated with recombination rate in both species, providing strong support for a role of linked selection in generating the heterogeneous genomic landscape of differentiation between species. Finally, we identify a number of candidate regions and genes that may have been subject to positive and/or balancing selection during the speciation process.