Altitude adaptation in Tibetans caused by introgression of Denisovan-like DNA.

As modern humans migrated out of Africa, they encountered many new environmental conditions, including greater temperature extremes, different pathogens and higher altitudes. These diverse environments are likely to have acted as agents of natural selection and to have led to local adaptations. One of the most celebrated examples in humans is the adaptation of Tibetans to the hypoxic environment of the high-altitude Tibetan plateau. A hypoxia pathway gene, EPAS1, was previously identified as having the most extreme signature of positive selection in Tibetans, and was shown to be associated with differences in haemoglobin concentration at high altitude. Re-sequencing the region around EPAS1 in 40 Tibetan and 40 Han individuals, we find that this gene has a highly unusual haplotype structure that can only be convincingly explained by introgression of DNA from Denisovan or Denisovan-related individuals into humans. Scanning a larger set of worldwide populations, we find that the selected haplotype is only found in Denisovans and in Tibetans, and at very low frequency among Han Chinese. Furthermore, the length of the haplotype, and the fact that it is not found in any other populations, makes it unlikely that the haplotype sharing between Tibetans and Denisovans was caused by incomplete ancestral lineage sorting rather than introgression. Our findings illustrate that admixture with other hominin species has provided genetic variation that helped humans to adapt to new environments.

Introgression maintains the genetic integrity of the mating-type determining chromosome of the fungus Neurospora tetrasperma.

Genome evolution is driven by a complex interplay of factors, including selection, recombination, and introgression. The regions determining sexual identity are particularly dynamic parts of eukaryotic genomes that are prone to molecular degeneration associated with suppressed recombination. In the fungus Neurospora tetrasperma, it has been proposed that this molecular degeneration is counteracted by the introgression of nondegenerated DNA from closely related species. In this study, we used comparative and population genomic analyses of 92 genomes from eight phylogenetically and reproductively isolated lineages of N. tetrasperma, and its three closest relatives, to investigate the factors shaping the evolutionary history of the genomes.We found that suppressed recombination extends across at least 6 Mbp (∼ 63%) of the mating-type (mat) chromosome in N. tetrasperma and is associated with decreased genetic diversity, which is likely the result primarily of selection at linked sites. Furthermore, analyses of molecular evolution revealed an increased mutational load in this region, relative to recombining regions. However, comparative genomic and phylogenetic analyses indicate that the mat chromosomes are temporarily regenerated via introgression from sister species; six of eight lineages show introgression into one of their mat chromosomes, with multiple Neurospora species acting as donors. The introgressed tracts have been fixed within lineages, suggesting that they confer an adaptive advantage in natural populations, and our analyses support the presence of selective sweeps in at least one lineage. Thus, these data strongly support the previously hypothesized role of introgression as a mechanism for the maintenance of mating-type determining chromosomal regions.

The role of China in the global spread of the current cholera pandemic.

Epidemics and pandemics of cholera, a severe diarrheal disease, have occurred since the early 19th century and waves of epidemic disease continue today. Cholera epidemics are caused by individual, genetically monomorphic lineages of Vibrio cholerae: the ongoing seventh pandemic, which has spread globally since 1961, is associated with lineage L2 of biotype El Tor. Previous genomic studies of the epidemiology of the seventh pandemic identified three successive sub-lineages within L2, designated waves 1 to 3, which spread globally from the Bay of Bengal on multiple occasions. However, these studies did not include samples from China, which also experienced multiple epidemics of cholera in recent decades. We sequenced the genomes of 71 strains isolated in China between 1961 and 2010, as well as eight from other sources, and compared them with 181 published genomes. The results indicated that outbreaks in China between 1960 and 1990 were associated with wave 1 whereas later outbreaks were associated with wave 2. However, the previously defined waves overlapped temporally, and are an inadequate representation of the shape of the global genealogy. We therefore suggest replacing them by a series of tightly delineated clades. Between 1960 and 1990 multiple such clades were imported into China, underwent further microevolution there and then spread to other countries. China was thus both a sink and source during the pandemic spread of V. cholerae, and needs to be included in reconstructions of the global patterns of spread of cholera.

Genome sequence and analysis of the tuber crop potato.

Potato (Solanum tuberosum L.) is the world's most important non-grain food crop and is central to global food security. It is clonally propagated, highly heterozygous, autotetraploid, and suffers acute inbreeding depression. Here we use a homozygous doubled-monoploid potato clone to sequence and assemble 86% of the 844-megabase genome. We predict 39,031 protein-coding genes and present evidence for at least two genome duplication events indicative of a palaeopolyploid origin. As the first genome sequence of an asterid, the potato genome reveals 2,642 genes specific to this large angiosperm clade. We also sequenced a heterozygous diploid clone and show that gene presence/absence variants and other potentially deleterious mutations occur frequently and are a likely cause of inbreeding depression. Gene family expansion, tissue-specific expression and recruitment of genes to new pathways contributed to the evolution of tuber development. The potato genome sequence provides a platform for genetic improvement of this vital crop.

Integrated genome sequence and linkage map of physic nut (Jatropha curcas L.), a biodiesel plant.

The family Euphorbiaceae includes some of the most efficient biomass accumulators. Whole genome sequencing and the development of genetic maps of these species are important components in molecular breeding and genetic improvement. Here we report the draft genome of physic nut (Jatropha curcas L.), a biodiesel plant. The assembled genome has a total length of 320.5 Mbp and contains 27,172 putative protein-coding genes. We established a linkage map containing 1208 markers and anchored the genome assembly (81.7%) to this map to produce 11 pseudochromosomes. After gene family clustering, 15,268 families were identified, of which 13,887 existed in the castor bean genome. Analysis of the genome highlighted specific expansion and contraction of a number of gene families during the evolution of this species, including the ribosome-inactivating proteins and oil biosynthesis pathway enzymes. The genomic sequence and linkage map provide a valuable resource not only for fundamental and applied research on physic nut but also for evolutionary and comparative genomics analysis, particularly in the Euphorbiaceae.

Detection of virus in CSF from the cases with meningoencephalitis by next-generation sequencing.

We screened for viral DNA in cerebrospinal fluid samples using next-generation sequencing (NGS) technology to diagnose CNS viral infections. We collected CSF samples from four cases with clinically suspected viral meningoencephalitis. DNA extracted from the samples was analyzed with NGS, and the results were further validated using PCR. Herpes simplex virus 1 (HSV-1) was detected in the CSF of two patients, HSV-2 and human herpes virus type 3 (HHV-3, VZV) in the CSF of two other patients separately. The number of unique reads of the identified viral genes ranged from 144 to 44205 (93.51 to 99.57%). The coverage of identified viral genes ranged from 12 to 98% with a depth value of 1.1 to 35, respectively. The results were further confirmed using PCR in three cases. The clinical presentation and outcomes of these four cases were consistent with the diagnostic results of NGS. NGS of CSF samples can be used as a diagnostic assay for CNS viral infection. Its further application for "pan-viral" or even "pan-microbial" screening of CSF might influence the diagnosis of CNS infectious diseases.

The sequence and de novo assembly of the giant panda genome.

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences.

The growing world population and shrinkage of arable land demand yield improvement of rice, one of the most important staple crops. To elucidate the genetic basis of yield and uncover its associated loci in rice, we resequenced the core recombinant inbred lines of Liang-You-Pei-Jiu, the widely cultivated super hybrid rice, and constructed a high-resolution linkage map. We detected 43 yield-associated quantitative trait loci, of which 20 are unique. Based on the high-density physical map, the genome sequences of paternal variety 93-11 and maternal cultivar PA64s of Liang-You-Pei-Jiu were significantly improved. The large recombinant inbred line population combined with plentiful high-quality single nucleotide polymorphisms and insertions/deletions between parental genomes allowed us to fine-map two quantitative trait loci, qSN8 and qSPB1, and to identify days to heading8 and lax panicle1 as candidate genes, respectively. The quantitative trait locus qSN8 was further confirmed to be days to heading8 by a complementation test. Our study provided an ideal platform for molecular breeding by targeting and dissecting yield-associated loci in rice.

The diploid genome sequence of an Asian individual.

Here we present the first diploid genome sequence of an Asian individual. The genome was sequenced to 36-fold average coverage using massively parallel sequencing technology. We aligned the short reads onto the NCBI human reference genome to 99.97% coverage, and guided by the reference genome, we used uniquely mapped reads to assemble a high-quality consensus sequence for 92% of the Asian individual's genome. We identified approximately 3 million single-nucleotide polymorphisms (SNPs) inside this region, of which 13.6% were not in the dbSNP database. Genotyping analysis showed that SNP identification had high accuracy and consistency, indicating the high sequence quality of this assembly. We also carried out heterozygote phasing and haplotype prediction against HapMap CHB and JPT haplotypes (Chinese and Japanese, respectively), sequence comparison with the two available individual genomes (J. D. Watson and J. C. Venter), and structural variation identification. These variations were considered for their potential biological impact. Our sequence data and analyses demonstrate the potential usefulness of next-generation sequencing technologies for personal genomics.

Genome sequences of three Helicobacter pylori strains isolated from atrophic gastritis and gastric ulcer patients in China.

Helicobacter pylori is a bacterial pathogen which can lead to several human gastric diseases. Here we describe the genome sequences of three strains isolated from atrophic gastritis and gastric ulcers patients in China. The data will permit genomic characterization of traits that may contribute to various gastric diseases.

Revised genome sequence of Burkholderia thailandensis MSMB43 with improved annotation.

There is growing interest in discovery of novel bioactive natural products from Burkholderia thailandensis. Here we report a significantly improved genome sequence and reannotation of Burkholderia thailandensis MSMB43, which will facilitate the discovery of new natural products through genome mining and studies of the metabolic versatility of this bacterium.

Draft genome sequence of Pseudomonas aeruginosa strain ATCC 27853.

Pseudomonas aeruginosa is a common bacterium that can cause disease. The versatility of Pseudomonas aeruginosa enables the organism to infect damaged tissues or those with reduced immunity which cause inflammation and sepsis. Here we report the genome sequence of the strain ATCC 27853.

Complete genome sequence of Pelagibacterium halotolerans B2(T).

Pelagibacterium halotolerans B2(T) is a marine halotolerant bacterium that was isolated from a seawater sample collected from the East China Sea. Here, we present the complete genome sequence of the type strain P. halotolerans B2(T), which consists of one chromosome (3,944,837 bp; 61.4% G+C content) and one plasmid (4,050 bp; 56.1% G+C content). This is the first complete genome of a member of the Pelagibacterium genus.

Genome sequences of three species in the family Planctomycetaceae.

Most of the species in the family Planctomycetaceae are of interest for their eukaryotic-like cell structures and characteristics of resistance to extreme environments. Here, we report draft genome sequences of three aquatic parasitic species of this family, Singulisphaera acidiphila (DSM 18658T), Schlesneria paludicola (DSM 18645T), and Zavarzinella formosa (DSM 19928T).

Genome sequence of Halobiforma lacisalsi AJ5, an extremely halophilic archaeon which harbors a bop gene.

The draft genome sequence (4,398,155 bp, with 65.35% G+C content) of Halobiforma lacisalsi AJ5, an extremely halophilic archaeon isolated from a salt lake, is reported here. This is the first genome report for a species of the Halobiforma genus.

Noninvasive prenatal testing of α-thalassemia and ß-thalassemia through population-based parental haplotyping.

BACKGROUND: Noninvasive prenatal testing (NIPT) of recessive monogenic diseases depends heavily on knowing the correct parental haplotypes. However, the currently used family-based haplotyping method requires pedigrees, and molecular haplotyping is highly challenging due to its high cost, long turnaround time, and complexity. Here, we proposed a new two-step approach, population-based haplotyping-NIPT (PBH-NIPT), using α-thalassemia and ß-thalassemia as prototypes. METHODS: First, we deduced parental haplotypes with Beagle 4.0 with training on a large retrospective carrier screening dataset (4356 thalassemia carrier screening-positive cases). Second, we inferred fetal haplotypes using a parental haplotype-assisted hidden Markov model (HMM) and the Viterbi algorithm. RESULTS: With this approach, we enrolled 59 couples at risk of having a fetus with thalassemia and successfully inferred 94.1% (111/118) of fetal alleles. We confirmed these alleles by invasive prenatal diagnosis, with 99.1% (110/111) accuracy (95% CI, 95.1-100%). CONCLUSIONS: These results demonstrate that PBH-NIPT is a sensitive, fast, and inexpensive strategy for NIPT of thalassemia.

NGS-based spinal muscular atrophy carrier screening of 10,585 diverse couples in China: a pan-ethnic study.

In this study, we performed a spinal muscular atrophy carrier screening investigation with NGS-based method. First, the validation for NGS-based method was implemented in 2255 samples using real-time PCR. The concordance between the NGS-based method and real-time PCR for the detection of SMA carrier and patient were up to 100%. Then, we applied this NGS-based method in 10,585 self-reported normal couples (34 Chinese ethnic groups from 5 provinces in South China) for SMA carrier screening. The overall carrier frequency was 1 in 73.8 (1.4%). It varied substantially between ethnic groups, highest in Dai ethnicity (4.3%), and no significant difference was found between five provinces. One couple was detected as carriers with an elevated risk of having an SMA affected baby. The distribution of SMN1:SMN2 genotype was also revealed in this study. Among the individuals with normal phenotype, the exon 7 copy-number ratio of SMN1 to SMN2 proved the gene conversion between them. With NGS-based method, we investigated SMA carrier status in Chinese population for the first time, and our results demonstrated that it is a promising alternative for SMA carrier screening and could provide data support and reference for future clinical application.

A genetic variation map for chicken with 2.8 million single-nucleotide polymorphisms.

We describe a genetic variation map for the chicken genome containing 2.8 million single-nucleotide polymorphisms (SNPs). This map is based on a comparison of the sequences of three domestic chicken breeds (a broiler, a layer and a Chinese silkie) with that of their wild ancestor, red jungle fowl. Subsequent experiments indicate that at least 90% of the variant sites are true SNPs, and at least 70% are common SNPs that segregate in many domestic breeds. Mean nucleotide diversity is about five SNPs per kilobase for almost every possible comparison between red jungle fowl and domestic lines, between two different domestic lines, and within domestic lines--in contrast to the notion that domestic animals are highly inbred relative to their wild ancestors. In fact, most of the SNPs originated before domestication, and there is little evidence of selective sweeps for adaptive alleles on length scales greater than 100 kilobases.

Noninvasive prenatal diagnosis of hemophilia A by a haplotype-based approach using cell-free fetal DNA.

Aim: We aimed to demonstrate noninvasive prenatal diagnosis (NIPD) of hemophilia A (HA) using a haplotype-based approach. Methods: Two families at risk for HA were recruited for this study. First, maternal haplotypes associated with pathogenic variants were constructed using the genotypes of the mothers and probands. Then, fetal haplotypes were deduced using a maternal haplotype-assisted hidden Markov model. Finally, the NIPD results were further confirmed by invasive prenatal diagnosis. Results: Two fetal genotypes were successfully inferred, with one normal fetus and one carrier fetus. The NIPD results were confirmed by invasive prenatal diagnosis, with a 100% consistency rate. Conclusion: Our test has been shown to be accurate and reliable. With further validation in a large patient cohort, this haplotype-based approach could be feasible for the NIPD of HA and other X-linked single-gene disorders.

Circulating tumor DNA detection is correlated to histologic types in patients with early-stage non-small-cell lung cancer.

OBJECTIVES: Circulating tumor DNA (ctDNA) testing in plasma in patients with non-small-cell lung cancer (NSCLC) has the potential to be a supplemental or surrogate tool for tissue biopsy. Detection of genomic abnormalities in ctDNA and their association with clinical characteristics in early-stage NSCLC need to be clarified. MATERIALS AND METHODS: Here, we comprehensively analyzed gene variations of 48 tumor tissues and 48 matched preoperative (pre-op) plasma and 25 postoperative (post-op) plasma from early-stage NSCLC patients using a targeted 546 genes capture-based next generation sequencing (NGS) assay. RESULTS: In early-stage NSCLC, the average mutation allele frequency (MAF) in pre-op plasma ctDNA was lower than that in tissue DNA (tDNA). The concordant gene variations between pre-op ctDNA and tDNA were difficult to detect. However, we found the tissue- pre-op plasma concordant ctDNA mutation detection ratio in lung squamous cell carcinoma (LUSC) was much higher than that in lung adenocarcinoma (LUAD). We also established a LUSC-LUAD classification model by a least absolute shrinkage and selection operator (LASSO) based approach to help separate LUAD from LUSC based on ctDNA profiling. This model included 14 gene mutations and extracted an accuracy of 89.2% in the training set and 91.5% in the testing set. Correlation analysis showed tDNA-ctDNA concordant ratio was related to histologic subtype, gene mutations and tumor size in early-stage NSCLC. CONCLUSION: This study suggests histology subtype and gene mutations could affect ctDNA detection in early-stage NSCLC. NGS-based ctDNA profile has the potential utility in LUSC-LUAD classification.