Systematic inference of copy-number genotypes from personal genome sequencing data reveals extensive olfactory receptor gene content diversity

Copy-number variations (CNVs) are widespread in the human genome, but comprehensive assignments of integer locus copy-numbers (i.e., copy-number genotypes) that, for example, enable discrimination of homozygous from heterozygous CNVs, have remained challenging. Here we present CopySeq, a novel computational approach with an underlying statistical framework that analyzes the depth-of-coverage of high-throughput DNA sequencing reads, and can incorporate paired-end and breakpoint junction analysis based CNV-analysis approaches, to infer locus copy-number genotypes. We benchmarked CopySeq by genotyping 500 chromosome 1 CNV regions in 150 personal genomes sequenced at low-coverage. The assessed copy-number genotypes were highly concordant with our performed qPCR experiments (Pearson correlation coefficient 0.94), and with the published results of two microarray platforms (95-99% concordance). We further demonstrated the utility of CopySeq for analyzing gene regions enriched for segmental duplications by comprehensively inferring copy-number genotypes in the CNV-enriched >800 olfactory receptor (OR) human gene and pseudogene loci. CopySeq revealed that OR loci display an extensive range of locus copy-numbers across individuals, with zero to two copies in some OR loci, and two to nine copies in others. Among genetic variants affecting OR loci we identified deleterious variants including CNVs and SNPs affecting ~15% and ~20% of the human OR gene repertoire, respectively, implying that genetic variants with a possible impact on smell perception are widespread. Finally, we found that for several OR loci the reference genome appears to represent a minor-frequency variant, implying a necessary revision of the OR repertoire for future functional studies. CopySeq can ascertain genomic structural variation in specific gene families as well as at a genome-wide scale, where it may enable the quantitative evaluation of CNVs in genome-wide association studies involving high-throughput sequencing.     

Functional constraint and small insertions and deletions in the ENCODE regions of the human genome

Background  

We describe the distribution of indels in the 44 Encyclopedia of DNA Elements (ENCODE) regions (about 1% of the human genome) and evaluate the potential contributions of small insertion and deletion polymorphisms (indels) to human genetic variation. We relate indels to known genomic annotation features and measures of evolutionary constraint.     

Genome sequencing reveals agronomically important loci in rice using MutMap

The majority of agronomic traits are controlled by multiple genes that cause minor phenotypic effects, making the identification of these genes difficult. Here we introduce MutMap, a method based on whole-genome resequencing of pooled DNA from a segregating population of plants that show a useful phenotype. In MutMap, a mutant is crossed directly to the original wild-type line and then selfed, allowing unequivocal segregation in second filial generation (F(2)) progeny of subtle phenotypic differences. This approach is particularly amenable to crop species because it minimizes the number of genetic crosses (n = 1 or 0) and mutant F(2) progeny that are required. We applied MutMap to seven mutants of a Japanese elite rice cultivar and identified the unique genomic positions most probable to harbor mutations causing pale green leaves and semidwarfism, an agronomically relevant trait. These results show that MutMap can accelerate the genetic improvement of rice and other crop plants.     

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.     

Comparative genome analysis reveals extensive conservation of genome organisation for Arabidopsis thaliana and Capsella rubella

Genome colinearity has been studied for two closely related diploid species of the Brassicaceae family, Arabidopsis thaliana and Capsella rubella. Markers mapping to chromosome 4 of A. thaliana were found on two linkage groups in Capsella and colinear segments spanning more than 10 cM were revealed. Detailed analysis of a 60 kbp region in A. thaliana and its counterpart in C. rubella showed virtually complete conservation of gene repertoire, order and orientation. The comparison of orthologous genes revealed very similar exon-intron structures and sequence identities of 90% or more were found for exon sequences. This extensive genome colinearity at the genetic and molecular level allows the efficient transfer of data from the well-studied A. thaliana genome to other species in the Brassicaceae family, substantially facilitating genome analysis studies for species of this family.     

Closing gaps in the human genome using sequencing by synthesis

The most recent release of the finished human genome contains 260 euchromatic gaps (excluding chromosome Y). Recent work has helped explain a large number of these unresolved regions as 'structural' in nature. Another class of gaps is likely to be refractory to clone-based approaches, and cannot be approached in ways previously described. We present an approach for closing these gaps using 454 sequencing. As a proof of principle, we closed all three remaining non-structural gaps in chromosome 15.     

Automated DNA sequencing and analysis of the human genome

In the past few years, striking advances have been made in automating DNA sequence analysis. Currently, efforts are underway to automate and improve DNA purification, mapping, and data processing procedures. The predictable advances in these technologies should soon place us in a position to sequence the entire human genome. The information derived from this project will have profound implications for basic biology and clinical medicine alike.     

人类基因组中巢式基因对的系统分析

基因组上一个基因的所有或大部分外显子位于另一基因内含子和UTR中被称为巢式基因(Nested gene)。巢式基因对(Nested gene pair)是由主基因(Host gene)和巢式基因组成的基因对。文章针对人类基因组上的顺式巢式基因对(简称为巢式基因对)进行全基因组鉴定,并通过比较巢式基因对在人和小鼠基因组上的保守性分布,研究巢式基因对的进化方式。保守性分析表明基因转座、序列原位突变和主基因转录起始/终止位点突变是巢式基因对进化的主要原因,其中主基因转录起始/终止位点突变是巢式基因对独特的一种进化方式。Gene On-tology分析揭示大部分巢式基因与主基因的产物在功能上没有相关性。     

Automated DNA sequencing and the analysis of the human genome

The Human Genome Initiative is a complex, multifaceted, international effort to establish a massive data base of map and sequence information for humans and other organisms. The success of this initiative is dependent upon the development of new technologies for the analysis of genomes. In this paper, an overview of the Human Genome Initiative is presented, and the current status of efforts to automate large-scale DNA sequence analysis is reviewed.     

Detecting structural variations in the human genome using next generation sequencing

Structural variations are widespread in the human genome and can serve as genetic markers in clinical and evolutionary studies. With the advances in the next-generation sequencing technology, recent methods allow for identification of structural variations with unprecedented resolution and accuracy. They also provide opportunities to discover variants that could not be detected on conventional microarray-based platforms, such as dosage-invariant chromosomal translocations and inversions. In this review, we will describe some of the sequencing-based algorithms for detection of structural variations and discuss the key issues in future development.     

Comparative hybridization reveals extensive genome variation in the AIDS-associated pathogen Cryptococcus neoformans

Background

Genome variability can have a profound influence on the virulence of pathogenic microbes. The availability of genome sequences for two strains of the AIDS-associated fungal pathogen Cryptococcus neoformans presented an opportunity to use comparative genome hybridization (CGH) to examine genome variability between strains of different mating type, molecular subtype, and ploidy.

Results

Initially, CGH was used to compare the approximately 100 kilobase MAT a and MATα mating-type regions in serotype A and D strains to establish the relationship between the Log2 ratios of hybridization signals and sequence identity. Subsequently, we compared the genomes of the environmental isolate NIH433 (MAT a) and the clinical isolate NIH12 (MATα) with a tiling array of the genome of the laboratory strain JEC21 derived from these strains. In this case, CGH identified putative recombination sites and the origins of specific segments of the JEC21 genome. Similarly, CGH analysis revealed marked variability in the genomes of strains representing the VNI, VNII, and VNB molecular subtypes of the A serotype, including disomy for chromosome 13 in two strains. Additionally, CGH identified differences in chromosome content between three strains with the hybrid AD serotype and revealed that chromosome 1 from the serotype A genome is preferentially retained in all three strains.

Conclusion

The genomes of serotypes A, D, and AD strains exhibit extensive variation that spans the range from small differences (such as regions of divergence, deletion, or amplification) to the unexpected disomy for chromosome 13 in haploid strains and preferential retention of specific chromosomes in naturally occurring diploids.