Ongoing transposition in cell culture reveals the phylogeny of diverse Drosophila S2 sublines

Cultured cells are widely used in molecular biology despite poor understanding of how cell line genomes change in vitro over time. Previous work has shown that Drosophila cultured cells have a higher transposable element content than whole flies, but whether this increase in transposable element content resulted from an initial burst of transposition during cell line establishment or ongoing transposition in cell culture remains unclear. Here, we sequenced the genomes of 25 sublines of Drosophila S2 cells and show that transposable element insertions provide abundant markers for the phylogenetic reconstruction of diverse sublines in a model animal cell culture system. DNA copy number evolution across S2 sublines revealed dramatically different patterns of genome organization that support the overall evolutionary history reconstructed using transposable element insertions. Analysis of transposable element insertion site occupancy and ancestral states support a model of ongoing transposition dominated by episodic activity of a small number of retrotransposon families. Our work demonstrates that substantial genome evolution occurs during long-term Drosophila cell culture, which may impact the reproducibility of experiments that do not control for subline identity.     

Screening of duplicated loci reveals hidden divergence patterns in a complex salmonid genome

A whole‐genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid‐origin lineages. However, little is known about general consequences of a WGD because gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid‐origin species. We demonstrate a new method that enables genome‐wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus‐specific copy number into account. We apply this method to RAD sequence data from different ecotypes of a polyploid‐origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences of WGDs and local segmental gene duplications.     

基因组拷贝数变异及其突变机理与人类疾病

拷贝数变异(Copy number variation,CNV)是由基因组发生重排而导致的,一般指长度为1 kb以上的基因组大片段的拷贝数增加或者减少,主要表现为亚显微水平的缺失和重复。CNV是基因组结构变异(Structural variation,SV)的重要组成部分。CNV位点的突变率远高于SNP(Single nucleotide polymorphism),是人类疾病的重要致病因素之一。目前,用来进行全基因组范围的CNV研究的方法有:基于芯片的比较基因组杂交技术(array-based comparative genomic hybridization,aCGH)、SNP分型芯片技术和新一代测序技术。CNV的形成机制有多种,并可分为DNA重组和DNA错误复制两大类。CNV可以导致呈孟德尔遗传的单基因病与罕见疾病,同时与复杂疾病也相关。其致病的可能机制有基因剂量效应、基因断裂、基因融合和位置效应等。对CNV的深入研究,可以使我们对人类基因组的构成、个体间的遗传差异、以及遗传致病因素有新的认识。     

Bayesian modeling of skewed X inactivation in genetically diverse mice identifies a novel Xce allele associated with copy number changes

Female mammals are functional mosaics of their parental X-linked gene expression due to X chromosome inactivation (XCI). This process inactivates one copy of the X chromosome in each cell during embryogenesis and that state is maintained clonally through mitosis. In mice, the choice of which parental X chromosome remains active is determined by the X chromosome controlling element (Xce), which has been mapped to a 176-kb candidate interval. A series of functional Xce alleles has been characterized or inferred for classical inbred strains based on biased, or skewed, inactivation of the parental X chromosomes in crosses between strains. To further explore the function structure basis and location of the Xce, we measured allele-specific expression of X-linked genes in a large population of F1 females generated from Collaborative Cross (CC) strains. Using published sequence data and applying a Bayesian “Pólya urn” model of XCI skew, we report two major findings. First, inter-individual variability in XCI suggests mouse epiblasts contain on average 20–30 cells contributing to brain. Second, CC founder strain NOD/ShiLtJ has a novel and unique functional allele, Xce^g, that is the weakest in the Xce allelic series. Despite phylogenetic analysis confirming that NOD/ShiLtJ carries a haplotype almost identical to the well-characterized C57BL/6J (Xce^b), we observed unexpected patterns of XCI skewing in females carrying the NOD/ShiLtJ haplotype within the Xce. Copy number variation is common at the Xce locus and we conclude that the observed allelic series is a product of independent and recurring duplications shared between weak Xce alleles.     

Comparative genomics of Leishmania isolates from Brazil confirms the presence of Leishmania major in the Americas

Leishmania (Leishmania) major is an important agent of cutaneous leishmaniasis, having as a vector sandflies belonging to the genus Phlebotomus. Although this species has been described as restricted to the Old World, parasites similar to L. major have been isolated from South American patients who have never travelled abroad. These parasites were named “L. major-like”, and several studies have been carried out to characterise them biochemically, molecularly, and biologically. However, the phylogenetic origin of these isolates is still unknown. In the present study we characterised three L. major-like isolates, named BH49, BH121 and BH129, using comparative genomics approaches. We evaluated the presence of gene and segmental duplications/deletions and the presence of aneuploidies that could explain the differences in infectivity observed in the BH49 and BH121 isolates. All isolates presented a pattern of mosaic aneuploidy and gene copy number variation, which are common in the genus Leishmania. Virulence factors such as phosphatases and peptidases were found to have increased gene copy numbers in the infective isolate, which could explain the difference in infectivity previously observed between BH121 and BH49. Phylogenetic analyses revealed that BH49, BH121 and BH129 L. major-like grouped with L. major isolates, and suggest they were imported from the Old World in at least two independent events. We suggest that new epidemiological inquiries should also evaluate L. major infections in South America, to assess the epidemiological importance of this species in the New World.     

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species,white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.     

人microRNA相关的遗传变异与肿瘤

微RNA(microRNA,miRNA)是新发现的一类进化上高度保守的重要的转录后调控因子,通过调节基因的表达而参与调控细胞凋亡、增殖及分化等生理过程,同时与肿瘤等疾病的发生发展密切相关。近年来研究发现,miRNA、miRNA生物合成通路基因及miRNA的靶基因结合位点的遗传变异(例如单核苷酸多态性和拷贝数变异等)可影响miRNA调控功能的发挥,并产生显著的遗传学效应。文章主要综述了miRNA相关的遗传变异与肿瘤易感性和临床转归等的研究进展。     

Analysis of genes associated with retrotransposons in the rice genome

Retrotransposons comprise a significant fraction of the rice genome. Despite their prevalence, the effects of retrotransposon insertions are not well understood, especially with regard to how they affect the expression of genes. In this study, we identified one-sixth of rice genes as being associated with retrotransposons, with insertions either in the gene itself or within its putative promoter region. Among genes with insertions in the promoter region, the likelihood of the gene being expressed was shown to be directly proportional to the distance of the retrotransposon from the translation start site. In addition, retrotransposon insertions in the transcribed region of the gene were found to be positively correlated with the presence of alternative splicing forms. Furthermore, preferential association of retrotransposon insertions with genes in several functional classes was identified. Some of the retrotransposons that are part of full-length cDNA (fl-cDNA) contribute splice sites and give rise to novel exons. Several interesting trends concerning the effects of retrotransposon insertions on gene expression were identified. Taken together, our data suggests that retrotransposon association with genes have a role in gene regulation. The data presented in this study provides a foundation for experimental studies to determine the role of retrotransposons in gene regulation.     

<Emphasis Type="Italic">Fastbreak</Emphasis>: a tool for analysis and visualization of structural variations in genomic data

Genomic studies are now being undertaken on thousands of samples requiring new computational tools that can rapidly analyze data to identify clinically important features. Inferring structural variations in cancer genomes from mate-paired reads is a combinatorially difficult problem. We introduce Fastbreak, a fast and scalable toolkit that enables the analysis and visualization of large amounts of data from projects such as The Cancer Genome Atlas.     

Comparative genomic analysis of Acinetobacter baumannii clinical isolates reveals extensive genomic variation and diverse antibiotic resistance determinants

Background

Acinetobacter baumannii is an important nosocomial pathogen that poses a serious health threat to immune-compromised patients. Due to its rapid ability to develop multidrug resistance (MDR), A. baumannii has increasingly become a focus of attention worldwide. To better understand the genetic variation and antibiotic resistance mechanisms of this bacterium at the genomic level, we reported high-quality draft genome sequences of 8 clinical isolates with various sequence types and drug susceptibility profiles.

Results

We sequenced 7 MDR and 1 drug-sensitive clinical A. baumannii isolates and performed comparative genomic analysis of these draft genomes with 16 A. baumannii complete genomes from GenBank. We found a high degree of variation in A. baumannii, including single nucleotide polymorphisms (SNPs) and large DNA fragment variations in the AbaR-like resistance island (RI) regions, the prophage and the type VI secretion system (T6SS). In addition, we found several new AbaR-like RI regions with highly variable structures in our MDR strains. Interestingly, we found a novel genomic island (designated as GI_BJ4) in the drug-sensitive strain BJ4 carrying metal resistance genes instead of antibiotic resistance genes inserted into the position where AbaR-like RIs commonly reside in other A. baumannii strains. Furthermore, we showed that diverse antibiotic resistance determinants are present outside the RIs in A. baumannii, including antibiotic resistance-gene bearing integrons, the bla_OXA-23-containing transposon Tn2009, and chromosomal intrinsic antibiotic resistance genes.

Conclusions

Our comparative genomic analysis revealed that extensive genomic variation exists in the A. baumannii genome. Transposons, genomic islands and point mutations are the main contributors to the plasticity of the A. baumannii genome and play critical roles in facilitating the development of antibiotic resistance in the clinical isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1163) contains supplementary material, which is available to authorized users.     

Copy number variation profile-based genomic typing of premenstrual dysphoric disorder in Chinese

Premenstrual dysphoric disorder (PMDD) affects nearly 5% of women of reproductive age. Symptomatic heterogeneity, together with largely unknown genetics, has greatly hindered its effective treatment. In the present study, analysis of genomic sequencing-based copy number variations (CNVs) called from 100 kb white blood cell DNA sequence windows by means of semisupervized clustering led to the segregation of patient genomes into the D and V groups, which correlated with the depression and invasion clinical types, respectively, with 89.0% consistency. Application of diagnostic CNV features selected using the correlation-based machine learning method enabled the classification of the CNVs obtained into the D group, V group, total patient group, and control group with an average accuracy of 83.0%. The power of the diagnostic CNV features was 0.98 on average, suggesting that these CNV features could be used for the molecular diagnosis of the major clinical types of PMDD. This demonstrated concordance between the CNV profiles and clinical types of PMDD supported the validity of symptom-based diagnosis of PMDD for differentiating between its two major clinical types, as well as the predominantly genetic nature of PMDD with a host of overlaps between multiple susceptibility genes/pathways and the diagnostic CNV features as indicators of involvement in PMDD etiology.     

Whole genome amplification strategy for forensic genetic analysis using single or few cell equivalents of genomic DNA

Evidentiary items sometimes contain an insufficient quantity of DNA for routine forensic genetic analysis. These so-called low copy number DNA samples (< 100 pg of genomic DNA) often fall below the sensitivity limitations of routine DNA analysis methods. Theoretically, one way of making such intractable samples amenable to analysis would be to increase the number of starting genomes available for subsequent STR (short tandem repeat) analysis by a whole genome amplification strategy (WGA). Although numerous studies employing WGA have focused primarily on clinical applications, few in-depth studies have been conducted to evaluate the potential usefulness of these methods in forensic casework. After an initial evaluation of existing methods, a modified WGA strategy was developed that appears to have utility for low copy number forensic casework specimens. The method employs a slight, but important, modification of the "improved primer extension preamplification PCR" method (I-PEP-PCR), which we term mIPEP (modified-I-PEP-PCR). Complete autosomal STR and Y-STR (Y chromosome short tandem repeat) profiles were routinely obtained with 5 pg of template DNA, which is equivalent to 1-2 diploid cells. Remarkably, partial Y- and autosomal STR profiles were obtained from mIPEP-treated DNA recovered from bloodstains exposed to the outside environment for 1 year whereas non-mIPEP-treated samples did not produce profiles. STR profiles were obtained from contact DNA from single dermal ridge fingerprints when the DNA was subjected to prior mIPEP amplification but not when the mIPEP step was omitted.     

The genetic background of clinical mastitis in Holstein-Friesian cattle

Mastitis is an inflammatory disease of the mammary gland, which has a significant economic impact and is an animal welfare concern. This work examined the association between single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) with the incidence of clinical mastitis (CM). Using information from 16 half-sib pairs of Holstein-Friesian cows (32 animals in total) we searched for genomic regions that differed between a healthy (no incidence of CM) and a mastitis-prone (multiple incidences of CM) half-sib. Three cows with average sequence depth of coverage below 10 were excluded, which left 13 half-sib pairs available for comparisons. In total, 191 CNV regions were identified, which were deleted in a mastitis-prone cow, but present in its healthy half-sib and overlapped in at least nine half-sib pairs. These regions overlapped with exons of 46 genes, among which APP (BTA1), FOXL2 (BTA1), SSFA2 (BTA2), OTUD3 (BTA2), ADORA2A (BTA17), TXNRD2 (BTA17) and NDUFS6 (BTA20) have been reported to influence CM. Moreover, two duplicated CNV regions present in nine healthy individuals and absent in their mastitis-affected half-sibs overlapped with exons of a cholinergic receptor nicotinic α 10 subunit on BTA15 and a novel gene (ENSBTAG00000008519) on BTA27. One CNV region deleted in nine mastitis-affected sibs overlapped with two neighbouring long non-coding RNA sequences located on BTA12. Single nucleotide polymorphisms with differential genotypes between a healthy and a mastitis-affected sib included 17 polymorphisms with alternate alleles in eight affected and healthy half-sib families. Three of these SNPs were located introns of genes: MET (BTA04), RNF122 (BTA27) and WRN (BTA27). In summary, structural polymorphisms in form of CNVs, putatively play a role in susceptibility to CM. Specifically, sequence deletions have a greater effect on reducing resistance against mastitis, than sequence duplications have on increasing resistance against the disease.     

Copy number variations and polymorphisms in HSP90AB1 and risk of systemic lupus erythematosus and efficacy of glucocorticoids

The aim of our study was to assess the associations of HSP90AB1 copy number variations (CNVs) with systemic lupus erythematosus (SLE) risk and glucocorticoids (GCs) efficacy, as well as the relationship between HSP90AB1 single‐nucleotide polymorphisms (SNPs) and GCs efficacy. HSP90AB1 CNVs and SLE risk were analysed in 519 patients and 538 controls. Patients treated with GCs were followed up for 12 weeks and were divided into sensitive and insensitive groups to investigate the effects of CNVs (419 patients) and SNPs (457 patients) on the efficacy of GCs. Health‐related quality of life (HRQoL) was also measured by SF‐36 at baseline and week 12 to explore the relationship between CNVs/SNPs and HRQoL improvements in Chinese SLE patients. Our results indicated a statistically significant association between HSP90AB1 CNVs and SLE (P_BH = 0.039), and this association was more pronounced in the female subgroup (P_BH = 0.039). However, we did not detect association of HSP90AB1 CNVs/SNPs with efficacy of GCs. But we found a marginal association between SNP rs13296 and improvement in Role‐emotional, while this association was not strong enough to survive in the multiple testing corrections. Collectively, our findings suggest that the copy number of HSP90AB1 is associated with SLE susceptibility. But copy number and polymorphisms of HSP90AB1 may not be associated with efficacy of GCs.     

Systematic analysis of <Emphasis Type="Italic">NPK1</Emphasis>-like genes in rice reveals a stress-inducible gene cluster co-localized with a quantitative trait locus of drought resistance

Phosphorylation by protein kinase is a ubiquitous key mechanism in translating external stimuli such as drought stress. NPK1 is a mitogen-activated protein kinase kinase kinase identified in Nicotiana tabacum and plays important roles in cytokinesis and auxin signaling transduction and responses to multiple stresses. Here we report the evolution, structure, and comprehensive expression profile of 21 NPK1-like genes in rice (Oryza sativa L.). Phylogenetic analysis of NPK1-like sequences in rice (OsNPKL), Arabidopsis, and other plants reveals that NPK1-like genes could be classified into three subgroups. Three OsNPKL gene clusters, located on chromosome 1 (OsNPKL1, 2, 3, and 4), 5 (OsNPKL14 and 15), and 10 (OsNPKL19 and 20), respectively, were identified in the rice genome. These clustered genes, which most likely evolved by tandem gene duplication, belong to the same phylogenetic subgroup, with similar genomic structures and conserved motifs in the kinase domain, which is unique to this subgroup. Expression analysis of OsNPKL genes under abiotic stresses suggests that the stress-responsive genes are mainly from the same subgroup. Especially interesting is that all the clustered genes are induced by drought, salt, or cold stress, and a few members are very strongly induced by drought. Some of the clustered genes are also induced by abscisic acid. The gene cluster on chromosome 1 is co-located with a quantitative trait locus (QTL) related to drought resistance. Although the drought-induced expression levels of the four genes in the cluster show no difference between the two parents used for QTL mapping, sequence variation in coding regions of the genes between the parents has provided some clues for further functional characterization of this gene cluster in abiotic stress tolerance in rice.     

Comparative analysis of the 100 kb region containing the Pi-k(h) locus between indica and japonica rice lines

We have recently cloned a pathogen inducible blast resistance gene Pi-kh from the indica rice line Tetep using a positional cloning approach. In this study, we carried out structural organization analysis of the Pi-kh locus in both indica and japonica rice lines. A 100 kb region containing 50 kb upstream and 50 kb down- stream sequences flanking to the Pi-kh locus was selected for the investigation. A total of 16 genes in indica and 15 genes in japonica were predicted and anno- tated in this region. The average GC content of indica and japonica genes in this region was 53.15% and 49.3%, respectively. Both indica and japonica sequences were polymorphic for simple sequence repeats having mono-, di-, tri-, tetra-, and pentanucleotides. Sequence analysis of the specific blast resistant Pi-kh allele of Tetep and the susceptible Pi-kh allele of the japonica rice line Nipponbare showed differences in the number and distribution of motifs involved in phosphorylation, resulting in the resistance phenotype in Tetep.     

Functional analysis of HapMap SNPs

Genome-wide association studies (GWAS) have successfully identified many genetic variants associated with complex diseases and traits. However, functional consequence of genetic variants studied in GWAS is not yet fully investigated, which would hinder the application of GWAS. We therefore performed a systematic functional analysis of HapMap SNPs, which have been most commonly used as the reference panel for GWAS. Our study highlights several characteristics of HapMap SNPs and identifies subsets of genetic variants with interesting functional implication. The results show that HapMap SNPs have good coverage within RefSeq genes, especially within known disease-related genes. On the other hand, only a small percentage of SNPs are non-synonymous SNPs while many SNPs are actually located at gene deserts. Moreover, many functionally important variants are not yet still interrogated. A redesigned SNP reference panel with additional functionally important variants would be useful to identify disease-causal variants in the future genome-wide studies.     

Plasmodium yoelii: combinatorial expression of variants of the 235 kDa rhoptry antigen during infection

The 235 kDa rhoptry protein Py235 of Plasmodium yoelii, has been implicated in erythrocyte invasion by the merozoite forms of the parasite. Py235 is encoded by a large, highly polymorphic gene family, members of which appear to be differentially transcribed. However, it is not clear how many variants are expressed at the protein level during an infection cycle and whether or not these variants are expressed selectively or combinatorially. Certain monoclonal antibodies to Py235 have been shown to attenuate parasite virulence upon passive transfer into mice, suggesting that this antigen or its derivatives may be useful vaccine candidates. To provide a basis for this, we sought to identify those variants that are recognised by the host immune system, and to establish the pattern of expression of the antigen in mice during infection. Using Py235 monoclonal antibodies as probes, we isolated distinct antigenic variants from an expression library, suggesting that the antigen repertoire is potentially large and that different Py235 variants may be produced during infection. The implications of these observations are discussed with respect to the ability of a cloned parasite line to express distinct antigenic variants in vivo.