Whole-genome sequencing and variant discovery in C. elegans

Massively parallel sequencing instruments enable rapid and inexpensive DNA sequence data production. Because these instruments are new, their data require characterization with respect to accuracy and utility. To address this, we sequenced a Caernohabditis elegans N2 Bristol strain isolate using the Solexa Sequence Analyzer, and compared the reads to the reference genome to characterize the data and to evaluate coverage and representation. Massively parallel sequencing facilitates strain-to-reference comparison for genome-wide sequence variant discovery. Owing to the short-read-length sequences produced, we developed a revised approach to determine the regions of the genome to which short reads could be uniquely mapped. We then aligned Solexa reads from C. elegans strain CB4858 to the reference, and screened for single-nucleotide polymorphisms (SNPs) and small indels. This study demonstrates the utility of massively parallel short read sequencing for whole genome resequencing and for accurate discovery of genome-wide polymorphisms.     

Comparisons with Caenorhabditis (approximately 100 Mb) and Drosophila (approximately 175 Mb) using flow cytometry show genome size in Arabidopsis to be approximately 157 Mb and thus approximately 25% larger than the Arabidopsis genome initiative estimate of approximately 125 Mb

Recent genome sequencing papers have given genome sizes of 180 Mb for Drosophila melanogaster Iso-1 and 125 Mb for Arabidopsis thaliana Columbia. The former agrees with early cytochemical estimates, but numerous cytometric estimates of around 170 Mb imply that a genome size of 125 Mb for arabidopsis is an underestimate. In this study, nuclei of species pairs were compared directly using flow cytometry. Co-run Columbia and Iso-1 female gave a 2C peak for arabidopsis only approx. 15 % below that for drosophila, and 16C endopolyploid Columbia nuclei had approx. 15 % more DNA than 2C chicken nuclei (with >2280 Mb). Caenorhabditis elegans Bristol N2 (genome size approx. 100 Mb) co-run with Columbia or Iso-1 gave a 2C peak for drosophila approx. 75 % above that for 2C C. elegans, and a 2C peak for arabidopsis approx. 57 % above that for C. elegans. This confirms that 1C in drosophila is approx. 175 Mb and, combined with other evidence, leads us to conclude that the genome size of arabidopsis is not approx. 125 Mb, but probably approx. 157 Mb. It is likely that the discrepancy represents extra repeated sequences in unsequenced gaps in heterochromatic regions. Complete sequencing of the arabidopsis genome until no gaps remain at telomeres, nucleolar organizing regions or centromeres is still needed to provide the first precise angiosperm C-value as a benchmark calibration standard for plant genomes, and to ensure that no genes have been missed in arabidopsis, especially in centromeric regions, which are clearly larger than once imagined.     

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest.Here,a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing,and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples.Overall,wild apples displayed higher levels of genetic diversity than cultivated apples.Linkage disequilibrium(LD) decays were observed quite rapidly in cultivated and wild apples,with an r~2-value below 0.2 at 440 and 280 bp,respectively.Moreover,bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples.Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality,thus suggesting that fruit quality had probably undergone selection during apple domestication.The genome of the earliest cultivated apple in China,Nai,was highly similar to that of Malus sieversii,and contained a small portion of genetic material from other wild apple species.This suggested that introgression could have been an important driving force during initial domestication of apple.These findings will facilitate future breeding and genetic dissection of complex traits in apple.     

Heading date 1 (Hd1), an ortholog of Arabidopsis CONSTANS, is a possible target of human selection during domestication to diversify flowering times of cultivated rice

During the domestication of rice (Oryza sativa L.), diversification of flowering time was important in expanding the areas of cultivation. Rice is a facultative short day (SD) plant and requires certain periods of dark to induce flowering. Heading date 1 (Hd1), a regulator of the florigen gene Hd3a, is one of the main factors used to generate diversity in flowering. Loss-of-function alleles of Hd1 are common in cultivated rice and cause the diversity of flowering time. However, it is unclear how these functional nucleotide polymorphisms of Hd1 accumulated in the course of evolution. Nucleotide polymorphisms within Hd1 and Hd3a were analyzed in 38 accessions of ancestral wild rice Oryza rufipogon and compared with those of cultivated rice. In contrast to cultivated rice, no nucleotide changes affecting Hd1 function were found in 38 accessions of wild rice ancestors. No functional changes were found in Hd3a in either cultivated or ancestral rice. A phylogenetic analysis indicated that evolution of the Hd1 alleles may have occurred independently in cultivars descended from various accessions of ancestral rice. The non-functional Hd1 alleles found in cultivated rice may be selected during domestication, because they were not found or very rare in wild ancestral rice. In contrast with Hd3a, which has been highly conserved, Hd1 may have undergone human selection to diversify the flowering times of rice during domestication or the early stage of the cultivation period.     

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.     

四倍体野生稻快速驯化: 启动人类新农业文明

通过人工选择优良遗传变异,将野生植物驯化为栽培作物,以满足人类对食物的需求,是人类发展历史中的重要事件,推动了人类文明的持续发展。随着世界人口持续增加,耕地面积不断减少,灾害性天气频发,全球粮食安全问题日趋严峻。基于作物驯化的分子机理及重要农艺性状的遗传基础,结合高通量基因组测序和高效基因组编辑技术,从头驯化野生植物,...     

Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly

Cucumber, Cucumis sativus L. is the only taxon with 2n = 2x = 14 chromosomes in the genus Cucumis. It consists of two cross‐compatible botanical varieties: the cultivated C. sativus var. sativus and the wild C. sativus var. hardwickii. There is no consensus on the evolutionary relationship between the two taxa. Whole‐genome sequencing of the cucumber genome provides a new opportunity to advance our understanding of chromosome evolution and the domestication history of cucumber. In this study, a high‐density genetic map for cultivated cucumber was developed that contained 735 marker loci in seven linkage groups spanning 707.8 cM. Integration of genetic and physical maps resulted in a chromosome‐level draft genome assembly comprising 193 Mbp, or 53% of the 367 Mbp cucumber genome. Strategically selected markers from the genetic map and draft genome assembly were employed to screen for fosmid clones for use as probes in comparative fluorescence in situ hybridization analysis of pachytene chromosomes to investigate genetic differentiation between wild and cultivated cucumbers. Significant differences in the amount and distribution of heterochromatins, as well as chromosomal rearrangements, were uncovered between the two taxa. In particular, six inversions, five paracentric and one pericentric, were revealed in chromosomes 4, 5 and 7. Comparison of the order of fosmid loci on chromosome 7 of cultivated and wild cucumbers, and the syntenic melon chromosome I suggested that the paracentric inversion in this chromosome occurred during domestication of cucumber. The results support the sub‐species status of these two cucumber taxa, and suggest that C. sativus var. hardwickii is the progenitor of cultivated cucumber.     

Species-specific accumulation of interspersed sequences in genus Saccharum

The genus Saccharum consists of two wild and four cultivated species. Novel interspersed sequences were isolated from cultivated sugar cane S. officinarum. These sequences were accumulated in all four cultivated species and their wild ancestral species S. robustum, but were not detected in the other wild species S. spontaneum and the relative Erianthus arundinaceus. The species-specific accumulation of interspersed sequences would correlate to the domestication of sugar canes.     

Caenorhabditis elegans mutant allele identification by whole-genome sequencing

Identification of the molecular lesion in Caenorhabditis elegans mutants isolated through forward genetic screens usually involves time-consuming genetic mapping. We used Illumina deep sequencing technology to sequence a complete, mutant C. elegans genome and thus pinpointed a single-nucleotide mutation in the genome that affects a neuronal cell fate decision. This constitutes a proof-of-principle for using whole-genome sequencing to analyze C. elegans mutants.     

Somatic excision of transposable element Tc1 from the Bristol genome of Caenorhabditis elegans.

We investigated the ability of the transposable element Tc1 to excise from the genome of the nematode Caenorhabditis elegans var. Bristol N2. Our results show that in the standard lab strain (Bristol), Tc1 excision occurred at a high frequency, comparable to that seen in the closely related Bergerac strain BO. We examined excision in the following way. We used a unique sequence flanking probe (pCeh29) to investigate the excision of Tc1s situated in the same location in both strains. Evidence of high-frequency excision from the genomes of both strains was observed. The Tc1s used in the first approach, although present in the same location in both genomes, were not known to be identical. Thus, a second approach was taken, which involved the genetic manipulation of a BO variant, Tc1(Hin). The ability of this BO Tc1(Hin) to excise was retained after its introduction into the N2 genome. Thus, we conclude that excision of Tc1 from the Bristol genome occurs at a high frequency and is comparable to that of Tc1 excision from the Bergerac genome. We showed that many Tc1 elements in N2 were apparently functionally intact and were capable of somatic excision. Even so, N2 Tc1s were prevented from exhibiting the high level of heritable transposition displayed by BO elements. We suggest that Bristol Tc1 elements have the ability to transpose but that transposition is heavily repressed in the gonadal tissue.     

Mobilization of quiet, endogenous Tc3 transposons of Caenorhabditis elegans by forced expression of Tc3 transposase.

The commonly studied Caenorhabditis elegans strain Bristol N2 contains approximately 15 copies per genome of the transposon Tc3. However, Tc3 is not active in Bristol N2. Tc3 contains one major open reading frame (Tc3A). We have fused this open reading frame to an inducible promoter and expressed it in a transgenic Bristol N2 line. Tc3A expression resulted in frequent excision and transposition of endogenous Tc3 elements. This shows that the Bristol N2 genome contains Tc3 transposons that are cis proficient for transposition, but are immobile because Tc3A is absent. We demonstrate that recombinant Tc3A binds specifically to the terminal nucleotides of the Tc3 inverted repeat, indicating that Tc3A is the Tc3 transposase. Activation of Tc3 transposition in vivo was accompanied by the appearance of extrachromosomal, linear copies of Tc3. These may be intermediates in Tc3 transposition.     

栽培茶树的驯化起源与传播

茶作为世界上最重要的饮品之一, 其栽培类型的驯化起源一直是人们关注的热点。本文总结了近年相关研究的进展, 讨论了存在的问题, 并对未来的研究方向提出建议。长江流域及以南地区分布有众多栽培茶树的野生近缘种, 特别集中于云南、贵州、广西等地; 一方面南方各族语言中“茶”发音的相似, 暗示了茶知识起源的单一性, 最可能起源于古代的巴蜀或云南, 另一方面遗传分析揭示栽培的茶存在多个起源中心, 即使Camellia sinensis (L.) O. Kuntze的几个栽培变种也可能起源于不同的地区; 文献记载, 茶的栽培中心曾经从西向东再向南迁移, 遗传多样性的变化也揭示了这一可能性, 但考古发现却提示最早的栽培茶可能出现在长江流域的最东部。我们推测在茶知识及栽培品种的传播过程中, 各地野生近缘植物的基因渗入栽培类型中, 或各地居民直接用当地野生茶培育出新的栽培茶类型, 从而导致遗传上的复杂性和语言上的一致性并存。茶树的祖先类型、起源地点、起源时间以及栽培品种的演变历程都还需要更为明确的证据, 未来应该以整个茶组植物为对象, 将茶文化、群体遗传学、谱系地理、人类学、气候变化、考古等多学科研究进行整合分析。     

The relative roles of three DNA repair pathways in preventing Caenorhabditis elegans mutation accumulation

Mutation is a central biological process whose rates and spectra are influenced by a variety of complex and interacting forces. Although DNA repair pathways are generally known to play key roles in maintaining genetic stability, much remains to be understood about the relative roles of different pathways in preventing the accumulation of mutations and the extent of heterogeneity in pathway-specific repair efficiencies across different genomic regions. In this study we examine mutation processes in base excision repair-deficient (nth-1) and nucleotide excision repair-deficient (xpa-1) Caenorhabditis elegans mutation-accumulation (MA) lines across 24 regions of the genome and compare our observations to previous data from mismatch repair-deficient (msh-2 and msh-6) and wild-type (N2) MA lines. Drastic variation in both average and locus-specific mutation rates, ranging two orders of magnitude for the latter, was detected among the four sets of repair-deficient MA lines. Our work provides critical insights into the relative roles of three DNA repair pathways in preventing C. elegans mutation accumulation and provides evidence for the presence of pathway-specific DNA repair territories in the C. elegans genome.     

Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice

Varying degrees of reduction of genetic diversity in crops relative to their wild progenitors occurred during the process of domestication. Such information, however, has not been available for the Asian cultivated rice (Oryza sativa) despite its importance as a staple food and a model organism. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationship and demographic history of O. sativa and its close relatives (Oryza rufipogon and Oryza nivara), we investigated nucleotide diversity data from 10 unlinked nuclear loci in species-wide samples of these species. The results indicated that O. rufipogon and O. nivara possessed comparable levels of nucleotide variation ((sil) = 0.0077 approximately 0.0095) compared with the relatives of other crops. In contrast, nucleotide diversity of O. sativa was as low as (sil) = 0.0024 and even lower ((sil) = 0.0021 for indica and 0.0011 for japonica), if we consider the 2 subspecies separately. Overall, only 20-10% of the diversity in the wild species was retained in 2 subspecies of the cultivated rice (indica and japonica), respectively. Because statistic tests did not reject the assumption of neutrality for all 10 loci, we further used coalescent to simulate bottlenecks under various lengths and population sizes to better understand the domestication process. Consistent with the dramatic reduction in nucleotide diversity, we detected a severe domestication bottleneck and demonstrated that the sequence diversity currently found in the rice genome could be explained by a founding population of 1,500 individuals if the initial domestication event occurred over a 3,000-year period. Phylogenetic analyses revealed close genetic relationships and ambiguous species boundary of O. rufipogon and O. nivara, providing additional evidence to treat them as 2 ecotypes of a single species. Lowest linkage disequilibrium (LD) was found in the perennial O. rufipogon where the r(2) value dropped to a negligible level within 400 bp, and the highest in the japonica rice where LD extended to the entirely sequenced region ( approximately 900 bp), implying that LD mapping by genome scans may not be feasible in wild rice due to the high density of markers needed.     

The molecular bases of cereal domestication and the history of rice

In this review, we discuss the development of molecular genetics and genomics that has allowed one to identify and characterize some of the key genes involved in cereal domestication. The list is far from being complete, but the first conclusion that can be drawn from the published works is that only a few loci have been the target of human selection in the first stages of the domestication process at the late neolithic. Mutations at these few loci have led to dramatic changes in plant morphology and phenology, transforming a wild into a cultivated plant. We also show that in the case of rice, for which the complete genome sequence is available, the development of new molecular markers based on retrotransposon insertion polymorphisms helped to resolve some of the questions regarding the origin of the domestication of the crop in Asia. To cite this article: O. Panaud, C. R. Biologies 332 (2009).     

CemaT1 is an active transposon within the Caenorhabditis elegans genome

The maT clade of transposons is a group of transposable elements intermediate in sequence and predicted protein structure to mariner and Tc transposons, with a distribution thus far limited to a few invertebrate species. In the nematode Caenorhabditis elegans, there are eight copies of CemaT1 that are predicted to encode a functional transposase, with five copies being >99% identical. We present evidence, based on searches of publicly available databases and on PCR-based mobility assays, that the CemaT1 transposase is expressed in C. elegans and that the CemaT transposons are capable of excising in both somatic and germline tissues. We also show that the frequency of CemaT1 excisions within the genome of the N2 strain of C. elegans is comparable to that of the Tc1 transposon. However, unlike Tc transposons in mutator strains of C. elegans, maT transposons do not exhibit increased frequencies of mobility, suggesting that maT is not regulated by the same factors that control Tc activity in these strains. Finally, we show that CemaT1 transposons are capable of precise transpositions as well as orientation inversions at some loci, and thereby become members of an increasing number of identified active transposons within the C. elegans genome.     

The mutational dynamics of the SARS-CoV-2 virus in serial passages in vitro

Since its outbreak in 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) keeps surprising the medical community by evolving diverse immune escape mutations in a rapid and effective manner. To gain deeper insight into mutation frequency and dynamics, we isolated ten ancestral strains of SARS-CoV-2 and performed consecutive serial incubation in ten replications in a suitable and common cell line and subsequently analysed them using RT-qPCR and whole genome sequencing. Along those lines we hoped to gain fundamental insights into the evolutionary capacity of SARS-CoV-2 in vitro. Our results identified a series of adaptive genetic changes, ranging from unique convergent substitutional mutations and hitherto undescribed insertions. The region coding for spike proved to be a mutational hotspot, evolving a number of mutational changes including the already known substitutions at positions S:484 and S:501. We discussed the evolution of all specific adaptations as well as possible reasons for the seemingly inhomogeneous potential of SARS-CoV-2 in the adaptation to cell culture. The combination of serial passage in vitro with whole genome sequencing uncovers the immense mutational potential of some SARS-CoV-2 strains. The observed genetic changes of SARS-CoV-2 in vitro could not be explained solely by selectively neutral mutations but possibly resulted from the action of directional selection accumulating favourable genetic changes in the evolving variants, along the path of increasing potency of the strain. Competition among a high number of quasi-species in the SARS-CoV-2 in vitro population gene pool may reinforce directional selection and boost the speed of evolutionary change.