A retrospective analysis of preemptive pharmacogenomic testing in 22,918 individuals from China.

BACKGROUND: Pharmacogenomics (PGx) examines the influence of genetic variation on drug responses. With more and more Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines published, PGx is gradually shifting from the reactive testing of single gene toward the preemptive testing of multiple genes. But the profile of PGx genes, especially for the intra-country diversity, is not well understood in China. METHODS: We retrospectively collected preemptive PGx testing data of 22,918 participants from 20 provinces of China, analyzed frequencies of alleles, genotypes and phenotypes of pharmacogenes, predicted drug responses for each participant, and performed comparisons between different provinces. RESULTS AND CONCLUSION: After analyzing 15 pharmacogenes from CPIC guidelines of 31 drugs, we found that 99.97% of individuals may have an atypical response to at least one drug; the participants carry actionable genotypes leading to atypical dosage recommendation for a median of eight drugs. Over 99% of the participants were recommended a decreased warfarin dose based on genetic factors. There were 20 drugs with high-risk ratios from 0.18% to 58.25%, in which clopidogrel showed the highest high-risk ratio. In addition, the high-risk ratio of rasburicase in GUANGDONG (risk ratio (RR) = 13.17, 95%CI:4.06-33.22, p < 0.001) and GUANGXI (RR = 23.44, 95%CI:8.83-52.85, p < 0.001) were significantly higher than that in all provinces. Furthermore, the diversity we observed among 20 provinces suggests that preemptive PGx testing in different geographical regions in China may need to pay more attention to specific genes. These results emphasize the importance of preemptive PGx testing and provide essential evidence for promoting clinical implementation in China.

Detecting mitochondrial mutations associated with aminoglycoside ototoxicity by noninvasive prenatal testing.

OBJECTIVES: Numerous diseases and disorders are associated with mitochondrial DNA (mtDNA) mutations, among which m.1555A > G and m.1494C > T mutations in the 12 S ribosomal RNA gene contribute to aminoglycoside-induced and nonsyndromic hearing loss worldwide. METHODS: A total of 76,842 qualified non-invasive prenatal (NIPT) samples were subjected to mtDNA mutation and haplogroup analysis. RESULTS: We detected 181 m.1555A > G and m.1494C > T mutations, 151 of which were subsequently sequenced for full-length mitochondrial genome verification. The positive predictive values for the m.1555A > G and m.1494C > T mutations were 90.78% and 90.00%, respectively, a performance comparable to that attained with newborn hearing screening. Furthermore, mitochondrial haplogroup analysis revealed that the 12 S rRNA 1555A > G mutation was enriched in sub-haplotype D5[p = 0, OR = 4.6706(2.81-7.78)]. CONCLUSIONS: Our findings indicate that the non-invasive prenatal testing of cell-free DNA obtained from maternal plasma can successfully detect m.1555A > G and m.1494C > T mutations.

Molecular diagnose of a large hearing loss population from China by targeted genome sequencing.

Hereditary hearing loss is genetically heterogeneous, with diverse clinical manifestations. Here we performed targeted genome sequencing of 227 hearing loss related genes in 1027 patients with bilateral hearing loss and 520 healthy volunteers with normal hearing to comprehensively identify the molecular etiology of hereditary hearing loss in a large cohort from China. We obtained a diagnostic rate of 57.25% (588/1027) for the patients, while 4.67% (48/1027) of the patients were identified with uncertain diagnoses. Of the implicated 35 hearing loss genes, three common genes, including SLC26A4(278/588), GJB2(207/588), MT-RNR1(19/588), accounted for 85.54% (503/588) of the diagnosed cases, while 32 uncommon hearing loss genes, including MYO15A, MITF, OTOF, POU3F4, PTPN11, etc. accounted for the remaining diagnostic rate of 14.46% (85/588). Apart from Pendred syndrome, other eight types of syndromic hearing loss were also identified. Of the 64 uncertain significant variants and 244 pathogenic/likely pathogenic variants identified in the patients, 129 novel variants were also detected. Thus, the molecular etiology presented with high heterogeneity with the leading causes to be SLC26A4 and GJB2 genes in the Chinese hearing loss population. It's urgent to develop a database of the ethnicity-matched healthy population as well as to perform functional studies for further classification of uncertain significant variants.

Correlation of low fetal fraction of cell-free DNA at the early second-trimester and pregnancy complications related to placental dysfunction in twin pregnancy.

Introduction: This study aimed to determine the correlation between fetal fraction (FF) of cell-free DNA (cf-DNA) and pregnancy complications related to placental dysfunction in Twin Pregnancy. Methods: This retrospective cohort study analyzed twin pregnant women who underwent non-invasive prenatal testing (NIPT) at 12+0-26+6 weeks of gestation from April 2017 to April 2021. Low fetal fraction (LFF) was defined individually as less than the 25th, 10th, 5th, and 2.5th percentile among all fetal fractions in the cohort. Primary outcomes included gestational hypertension (GH), preeclampsia (PE), gestational diabetes mellitus (GDM), and small for gestational age (SGA). Logistic regression analysis was used to assess the relationship between LFF and pregnancy complications. Results: A total of 500 twin pregnancies (male-male twins, 245; female-female twins, 255) were included in this study. In LFF group (FF < 25th percentiles), maternal BMI was significantly higher than FF > 75th percentiles (23.6 kg/m2 vs. 21.3 kg/m2; P < 0.001). The risk of SGA increased gradually from FF < 25th percentiles [adjusted odds ratio (OR), 1.71; 95% confidence interval (CI), 1.07-2.99; P = 0.016] to FF < 2.5th percentiles (adjusted OR, 4.44; 95% CI,1.33-14.82; P < 0.015). In addition, the risks of SGA in both fetuses were higher than the risks of at least one fetus SGA in LFF group. LFF had no correlation with GH, PE, and GDM in twin pregnancy. Conclusion: LFF has a strong association with increased risk of SGA in twin pregnancy. Moreover, FF of cf-DNA may provide a new idea for the early screening of diseases related to placental dysfunction in twin pregnancy.

Impact of STAT1 polymorphisms on crizotinib-induced hepatotoxicity in ALK-positive non-small cell lung cancer patients.

PURPOSE: Crizotinib is the first-line small molecule tyrosine kinase inhibitor for ALK-positive non-small cell lung cancer. In this study, a retrospective pharmacogenomics investigation was conducted to explore the relationship between genes related to RTK downstream signaling pathways and crizotinib-induced hepatic toxicity in ALK-positive NSCLC patients. METHODS: The variable importance analysis of random forest algorithm was applied to identify the significant features which contribute to the crizotinib sensitivity in Cancer Cell Line Encyclopedia (CCLE) database. The KEGG and reactome pathway enrichment analysis were conducted with EnrichR. The differential expression genes were identified with R package DESeq2 in CCLE liver derived cell lines between crizotinib sensitive and resistant groups. From 2012 to 2015, 42 NSCLC patients were enrolled in this study. 90 polymorphisms were genotyped using the Sequenom Massarray system. Sequencing of HGFR (c-Met) genes was carried out on the Ion Torrent Proton. RESULTS: In total, 66.7% NSCLC patients suffered from crizotinib-induced liver toxicity and 11.9% progressed to severe hepatic toxicity. The features with the top importance from classification and regression random forest model were enriched in RTK downstream signaling pathways (JAK/STAT, RAS/RAF/MAPK, PI3K/AKT pathways) and immune system-related pathways. Collagen family genes (COL1A1, COL1A2, COL6A1, COL5A1) and other extracellular matrix protein (TNC, TAGLN, TENM2, EDIL3, VCAN, CNN1, SH3BP4, TAGLN), which were closely related to MAPK-ERK signaling pathways, were significantly enriched in crizotinib resistant cell lines. In multiple logistic regression, STAT1 rs10208033 (T > C) was significantly associated with crizotinib-induced liver toxicity (OR = 6.733, 95% CI 1.406-32.24, P = 0.017). Compared with non-CC, OR is 5.5 (95% CI 1.219-24.81, P = 0.027) for STAT1 rs10208033 CC genotype to develop crizotinib-induced liver toxicity. Further cell viability test in human fetal hepatocyte line, L-02, reveals that the STAT1 inhibitor might protect hepatocyte cells from the toxicity caused by crizotinib. CONCLUSION: Polymorphism of rs10208033 is a potential biomarker for predicting crizotinib-induced hepatotoxicity. These results suggest that STAT1 plays an important role in crizotinib-induced hepatotoxicity. Further studies are needed to confirm our finding and understand the underlying mechanisms.

Clinical performance of non-invasive prenatal served as a first-tier screening test for trisomy 21, 18, 13 and sex chromosome aneuploidy in a pilot city in China.

BACKGROUND: Cell-free fetal DNA (cffDNA) has opened up new approaches for non-invasive prenatal testing (NIPT), and it is often used as the second-tier test for high-risk pregnant women in detecting trisomy (T) 21, T18, and T13 after serum biochemistry screening. This study aims to discuss the clinical performance of NIPT as an alternative first-tier screening test for pregnant women in detecting T21, T18, T13, and sex chromosome aneuploidies (SCAs) in China. METHODS: A total of 42,924 samples were recruited. The cell-free plasma DNA was directly sequenced. Each of the chromosome aneuploidies of PPV was analyzed. A total of 22 placental samples were acquired, including 14 FP and 8 TP samples. The placental verification of FP NIPT results was performed. RESULTS: Among 42,924 samples, 281 (0.65%) positive cases, including 87 of T21, 31 of T18, 22 of T13, and 141 of SCAs were detected. For the detection of T21, the positive predictive value (PPV) was 78.46%, for trisomy 18, 62.96%, for trisomy 13, 10.00%, for SCAs, 47.22% in the total samples. For trisomy 21, the PPV was 86.67%, for trisomy 18, 80.00%, for trisomy 13, 20.00%, for SCAs, 56.52% in advanced maternal age (AMA) women. The PPV of T21 increased with age. For T18, the PPV showed an overall upward trend. For T13 and SCAs, PPV was raised first and then lowered. Placental verification of false positive (FP) NIPT results confirmed confined placental mosaicism(CPM) was the reason for false positives. CONCLUSIONS: This study represents the first time that NIPT has been used as a first-tier screening test for fetal aneuploidies in a pilot city with large clinical samples in China. We propose that NIPT could replace serum biochemistry screening as a first-tier test.

Targeted bile acids and gut microbiome profiles reveal the hepato-protective effect of WZ tablet (Schisandra sphenanthera extract) against LCA-induced cholestasis.

Cholestasis is caused by the obstacle of bile formation or secretion and can develop into severe liver diseases. We previously reported the ethanol extract of Schisandra sphenanthera (Wuzhi tablet, WZ) can significantly protect against lithocholic acid (LCA)-induced intrahepatic cholestasis in mice, partially due to the activation of PXR pathway and promotion of liver regeneration. However, the effect of WZ on the bile acids profile and gut microbiome in cholestastic mice remain unknown. In this study, the effect of WZ against LCA-induced liver injury was evaluated and its effect on the bile acids metabolome and gut microbiome profiles in cholestastic mice was further investigated. Targeted metabolomics analysis was performed to examine the change of bile acids in the serum, liver, intestine and feces. The change of intestinal flora were detected by the genomics method. Targeted metabolomics analysis revealed that WZ enhanced the excretion of bile acids from serum and liver to intestine and feces. Genomics analysis of gut microbiome showed that WZ can reverse LCA-induced gut microbiome disorder to the normal level. In conclusion, WZ protects against LCA-induced cholestastic liver injury by reversing abnormal bile acids profiles and alteration of gut microbiome.

Noninvasive prenatal detection of hemoglobin Bart hydrops fetalis via maternal plasma dispensed with parental haplotyping using the semiconductor sequencing platform.

BACKGROUND: Thalassemia is one of the most common monogenetic diseases in the south of China and Southeast Asia. Hemoglobin Bart's hydrops fetalis syndrome was caused by a homozygous Southeast Asian deletion (-/-) in the HBA gene. Few studies have proved the potential of screen for Bart's hydrops fetalis using fetal cell-free DNA. However, the number of cases is still relatively small. Clinical trials of large samples would be needed. OBJECTIVE: In this study, we aimed to develop a noninvasive method of target-captured sequencing and genotyping by the Bayesian method using cell-free fetal DNA to identify the fetal genotype in pregnant women who are at risk of having hemoglobin Bart hydrops fetalis in a large-scale study. STUDY DESIGN: In total, 192,173 couples from 30 hospitals were enrolled in our study and 878 couples were recruited, among whom both the pregnant women and their husbands were detected to be carriers of Southeast Asian type (-/αα) of α-thalassemia. Prenatal diagnosis was performed by chorionic villus sampling, amniocentesis, or cordocentesis using gap-polymerase chain reaction considered as the golden standard. RESULTS: As a result, we found that the sensitivity and specificity of our noninvasive method were 98.81% and 94.72%, respectively, in the training set as well as 100% and 99.31%, respectively, in the testing set. Moreover, our method could identify all of 885 maternal samples with the Southeast Asian carrier and 36 trisomy samples with 100% of sensitivity in T13, T18, and T21 and 99.89% (1 of 917) and 99.88% (1 of 888) of specificity in T18 and T21, respectively. CONCLUSION: Our method opens the possibility of early screening for maternal genotyping of α-thalassemia, fetal aneuploidies in chromosomes 13/18/21, and hemoglobin Bart hydrops fetalis detection in 1 tube of maternal plasma.

The Gene Mutation Spectrum of Breast Cancer Analyzed by Semiconductor Sequencing Platform.

To use the semiconductor sequencing platform (SSP) to analyze the gene mutate spectrum of breast cancer patients. We recruited 46 breast cancer patients, and detected the ER/PR/HER2 expression level of the tumor tissue by immunohistochemistry. In addition, combined with SSP technology, we detected 207 hot mutation regions of 50 breast cancer related genes with multiple polymerase chain reaction (PCR) technology. There were 8 cases of grade I, 18 cases of grade II, 20 cases of grade III in 46 breast cancer patients according to histological grade and 12 cases of ER/PR + HER2 +, 18 cases of ER/PR + HER2-, 13 cases of ER/PR - HER2 +, 20 cases of ER/PR - HER2- according to ER/PR/HER2 status classification. Moreover, we found that there were 33 gene locus mutations of 8 genes including AKT1, APC, BRAF, CDKN2A, KRAS, PTEN, PIK3CA and TP53, but difference was not statistically significant (P > 0.05) when compared these gene mutations (except for PIK3CA) in each groups according to the histological classification of breast cancer and the ER/PR/HER2 classification. PIK3CA mutation rate of grade I was obviously higher than that of grade II ~ III histological grading in breast cancer patients (P < 0.05). Based on our results, we drew a conclusion that the occurrence and development of breast cancer was a process involved multiple genes. Here, we found that PIK3CA played a role in the development of the early stage of breast cancer, which could provide clinical basis for treatment of breast cancer. Moreover, SSP technology could be an effective and sensitive method for detection of gene mutation spectrum in breast cancer.

A foetus with 18p11.32-q21.2 duplication and Xp22.33-p11.1 deletion derived from a maternal reciprocal translocation t(X;18)(q13;q21.3).

BACKGROUND: Non-invasive prenatal testing (NIPT) evaluates circulating cell-free DNA (cfDNA) and has been widely applied, with highly accurate results for detecting foetal trisomies 21, 18 and 13. Recently, increasing attention has been paid to the clinical application of the non-invasive detection of foetal sub-chromosomal duplications and deletions beyond common aneuploidies. CASE PRESENTATION: A 32-year-old healthy pregnant woman was referred to the Medical Genetic Centre of Ganzhou Maternal and Child Health Care Hospital. As routine practice, ultrasound examination at a gestational age of 16 weeks showed that the foetus is normal. To avoid invasive prenatal diagnosis procedures, an NIPT was offered to further screen for common foetal chromosomal abnormalities. The result showed that there was an approximately 50.94 Mb duplication in p11.32-q21.2 of chromosome 18 and an approximately 58.46 Mb deletion in p22.33-p11.1 of chromosome X. In addition, the chromosome karyotypes of the parents and foetus were also analysed. Chromosome karyotype analysis results showed that foetal karyotype was 46,X,der(18), the maternal karyotype was 46,XX,t(X;18)(q13;q21.3), and the paternal karyotype revealed no obvious abnormality. CONCLUSION: In this case, we successfully detected a healthy pregnant woman with balanced translocation X;18(q13;q21.3) and described the foetal karyotype as 46,X,der(18)t(X;18)(q11;q21.1)mat. Our report illustrated these cases which present complex X;autosome balance translocation and X;autosome unbalance translocation which may contribute to severe clinical phenotypes. In addition, our report also proved that the interruption of genes in the Xq critical region is not only reason of primary infertility. Finally, we prompted that NIPT might play a role in the first trimester screening of sub-chromosomal rearrangement.

Noninvasive detection of fetal subchromosomal abnormalities by semiconductor sequencing of maternal plasma DNA.

Noninvasive prenatal testing (NIPT) using sequencing of fetal cell-free DNA from maternal plasma has enabled accurate prenatal diagnosis of aneuploidy and become increasingly accepted in clinical practice. We investigated whether NIPT using semiconductor sequencing platform (SSP) could reliably detect subchromosomal deletions/duplications in women carrying high-risk fetuses. We first showed that increasing concentration of abnormal DNA and sequencing depth improved detection. Subsequently, we analyzed plasma from 1,456 pregnant women to develop a method for estimating fetal DNA concentration based on the size distribution of DNA fragments. Finally, we collected plasma from 1,476 pregnant women with fetal structural abnormalities detected on ultrasound who also underwent an invasive diagnostic procedure. We used SSP of maternal plasma DNA to detect subchromosomal abnormalities and validated our results with array comparative genomic hybridization (aCGH). With 3.5 million reads, SSP detected 56 of 78 (71.8%) subchromosomal abnormalities detected by aCGH. With increased sequencing depth up to 10 million reads and restriction of the size of abnormalities to more than 1 Mb, sensitivity improved to 69 of 73 (94.5%). Of 55 false-positive samples, 35 were caused by deletions/duplications present in maternal DNA, indicating the necessity of a validation test to exclude maternal karyotype abnormalities. This study shows that detection of fetal subchromosomal abnormalities is a viable extension of NIPT based on SSP. Although we focused on the application of cell-free DNA sequencing for NIPT, we believe that this method has broader applications for genetic diagnosis, such as analysis of circulating tumor DNA for detection of cancer.

Complex evolutionary trajectories of sex chromosomes across bird taxa.

Sex-specific chromosomes, like the W of most female birds and the Y of male mammals, usually have lost most genes owing to a lack of recombination. We analyze newly available genomes of 17 bird species representing the avian phylogenetic range, and find that more than half of them do not have as fully degenerated W chromosomes as that of chicken. We show that avian sex chromosomes harbor tremendous diversity among species in their composition of pseudoautosomal regions and degree of Z/W differentiation. Punctuated events of shared or lineage-specific recombination suppression have produced a gradient of "evolutionary strata" along the Z chromosome, which initiates from the putative avian sex-determining gene DMRT1 and ends at the pseudoautosomal region. W-linked genes are subject to ongoing functional decay after recombination was suppressed, and the tempo of degeneration slows down in older strata. Overall, we unveil a complex history of avian sex chromosome evolution.

Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle.

Genetic sex determination by W and Z chromosomes has developed independently in different groups of organisms. To better understand the evolution of sex chromosomes and the plasticity of sex-determination mechanisms, we sequenced the whole genomes of a male (ZZ) and a female (ZW) half-smooth tongue sole (Cynoglossus semilaevis). In addition to insights into adaptation to a benthic lifestyle, we find that the sex chromosomes of these fish are derived from the same ancestral vertebrate protochromosome as the avian W and Z chromosomes. Notably, the same gene on the Z chromosome, dmrt1, which is the male-determining gene in birds, showed convergent evolution of features that are compatible with a similar function in tongue sole. Comparison of the relatively young tongue sole sex chromosomes with those of mammals and birds identified events that occurred during the early phase of sex-chromosome evolution. Pertinent to the current debate about heterogametic sex-chromosome decay, we find that massive gene loss occurred in the wake of sex-chromosome 'birth'.

Ultra-deep sequencing enables high-fidelity recovery of biodiversity for bulk arthropod samples without PCR amplification.

BACKGROUND: Next-generation-sequencing (NGS) technologies combined with a classic DNA barcoding approach have enabled fast and credible measurement for biodiversity of mixed environmental samples. However, the PCR amplification involved in nearly all existing NGS protocols inevitably introduces taxonomic biases. In the present study, we developed new Illumina pipelines without PCR amplifications to analyze terrestrial arthropod communities. RESULTS: Mitochondrial enrichment directly followed by Illumina shotgun sequencing, at an ultra-high sequence volume, enabled the recovery of Cytochrome c Oxidase subunit 1 (COI) barcode sequences, which allowed for the estimation of species composition at high fidelity for a terrestrial insect community. With 15.5 Gbp Illumina data, approximately 97% and 92% were detected out of the 37 input Operational Taxonomic Units (OTUs), whether the reference barcode library was used or not, respectively, while only 1 novel OTU was found for the latter. Additionally, relatively strong correlation between the sequencing volume and the total biomass was observed for species from the bulk sample, suggesting a potential solution to reveal relative abundance. CONCLUSIONS: The ability of the new Illumina PCR-free pipeline for DNA metabarcoding to detect small arthropod specimens and its tendency to avoid most, if not all, false positives suggests its great potential in biodiversity-related surveillance, such as in biomonitoring programs. However, further improvement for mitochondrial enrichment is likely needed for the application of the new pipeline in analyzing arthropod communities at higher diversity.

Stump-like mathematical model and computer simulation on dynamic separation of capillary zone electrophoresis with different sample injections.

In this paper, a novel mathematical model of stump-like peak was constructed for alternative computer simulation of capillary zone electrophoresis (CZE) by using Haarhoff-Van der Linde (HVL) function. Unlike a classical model of Gaussian peak, the developed model contains both Gaussian and plateau concentration distributions of analytes. Based on the model, the relevant computer software was developed and implemented in Borland Delphi 7 environment. The relevant results revealed that (i) the software could freely simulate the plateau and Gaussian distribution peaks; (ii) the simulator could simulate a dynamic process of CZE properly; and (iii) the program could display the final electropherogram of numerous analytes in CZE. We further conducted the relevant experiments and compared them with the simulations. The comparisons demonstrated the high agreement between the simulation results and the experiments as well as those cited from references. In addition, the software could calculate effective mobility, diffusion coefficient and concentration of analytes based on the physico-chemical parameters input by users. The developed model and software have evident significance for understanding of electrophoretic separation, conditional optimization and basic computation on physico-chemical parameters of analytes in CZE.

Whole-genome sequencing of Oryza brachyantha reveals mechanisms underlying Oryza genome evolution.

The wild species of the genus Oryza contain a largely untapped reservoir of agronomically important genes for rice improvement. Here we report the 261-Mb de novo assembled genome sequence of Oryza brachyantha. Low activity of long-terminal repeat retrotransposons and massive internal deletions of ancient long-terminal repeat elements lead to the compact genome of Oryza brachyantha. We model 32,038 protein-coding genes in the Oryza brachyantha genome, of which only 70% are located in collinear positions in comparison with the rice genome. Analysing breakpoints of non-collinear genes suggests that double-strand break repair through non-homologous end joining has an important role in gene movement and erosion of collinearity in the Oryza genomes. Transition of euchromatin to heterochromatin in the rice genome is accompanied by segmental and tandem duplications, further expanded by transposable element insertions. The high-quality reference genome sequence of Oryza brachyantha provides an important resource for functional and evolutionary studies in the genus Oryza.

Insights into salt tolerance from the genome of Thellungiella salsuginea.

Thellungiella salsuginea, a close relative of Arabidopsis, represents an extremophile model for abiotic stress tolerance studies. We present the draft sequence of the T. salsuginea genome, assembled based on ~134-fold coverage to seven chromosomes with a coding capacity of at least 28,457 genes. This genome provides resources and evidence about the nature of defense mechanisms constituting the genetic basis underlying plant abiotic stress tolerance. Comparative genomics and experimental analyses identified genes related to cation transport, abscisic acid signaling, and wax production prominent in T. salsuginea as possible contributors to its success in stressful environments.

Deciphering neo-sex and B chromosome evolution by the draft genome of Drosophila albomicans.

BACKGROUND: Drosophila albomicans is a unique model organism for studying both sex chromosome and B chromosome evolution. A pair of its autosomes comprising roughly 40% of the whole genome has fused to the ancient X and Y chromosomes only about 0.12 million years ago, thereby creating the youngest and most gene-rich neo-sex system reported to date. This species also possesses recently derived B chromosomes that show non-Mendelian inheritance and significantly influence fertility. METHODS: We sequenced male flies with B chromosomes at 124.5-fold genome coverage using next-generation sequencing. To characterize neo-Y specific changes and B chromosome sequences, we also sequenced inbred female flies derived from the same strain but without B's at 28.5-fold. RESULTS: We assembled a female genome and placed 53% of the sequence and 85% of the annotated proteins into specific chromosomes, by comparison with the 12 Drosophila genomes. Despite its very recent origin, the non-recombining neo-Y chromosome shows various signs of degeneration, including a significant enrichment of non-functional genes compared to the neo-X, and an excess of tandem duplications relative to other chromosomes. We also characterized a B-chromosome linked scaffold that contains an actively transcribed unit and shows sequence similarity to the subcentromeric regions of both the ancient X and the neo-X chromosome. CONCLUSIONS: Our results provide novel insights into the very early stages of sex chromosome evolution and B chromosome origination, and suggest an unprecedented connection between the births of these two systems in D. albomicans.

Genome sequencing and comparison of two nonhuman primate animal models, the cynomolgus and Chinese rhesus macaques.

The nonhuman primates most commonly used in medical research are from the genus Macaca. To better understand the genetic differences between these animal models, we present high-quality draft genome sequences from two macaque species, the cynomolgus/crab-eating macaque and the Chinese rhesus macaque. Comparison with the previously sequenced Indian rhesus macaque reveals that all three macaques maintain abundant genetic heterogeneity, including millions of single-nucleotide substitutions and many insertions, deletions and gross chromosomal rearrangements. By assessing genetic regions with reduced variability, we identify genes in each macaque species that may have experienced positive selection. Genetic divergence patterns suggest that the cynomolgus macaque genome has been shaped by introgression after hybridization with the Chinese rhesus macaque. Macaque genes display a high degree of sequence similarity with human disease gene orthologs and drug targets. However, we identify several putatively dysfunctional genetic differences between the three macaque species, which may explain functional differences between them previously observed in clinical studies.

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.