Transposon activity, local duplications and propagation of structural variants across haplotypes drive the evolution of the Drosophila S2 cell line

BACKGROUND

Immortalized cell lines are widely used model systems whose genomes are often highly rearranged and polyploid. However, their genome structure is seldom deciphered and is thus not accounted for during analyses. We therefore used linked short- and long-read sequencing to perform haplotype-level reconstruction of the genome of a Drosophila melanogaster cell line (S2-DRSC) with a complex genome structure.

RESULTS

Using a custom implementation (that is designed to use ultra-long reads in complex genomes with nested rearrangements) to call structural variants (SVs), we found that the most common SV was repetitive sequence insertion or deletion (> 80% of SVs), with Gypsy retrotransposon insertions dominating. The second most common SV was local sequence duplication. SNPs and other SVs were rarer, but several large chromosomal translocations and mitochondrial genome insertions were observed. Haplotypes were highly similar at the nucleotide level but structurally very different. Insertion SVs existed at various haplotype frequencies and were unlinked on chromosomes, demonstrating that haplotypes have different structures and suggesting the existence of a mechanism that allows SVs to propagate across haplotypes. Finally, using public short-read data, we found that transposable element insertions and local duplications are common in other D. melanogaster cell lines.

CONCLUSIONS

The S2-DRSC cell line evolved through retrotransposon activity and vast local sequence duplications, that we hypothesize were the products of DNA re-replication events. Additionally, mutations can propagate across haplotypes (possibly explained by mitotic recombination), which enables fine-tuning of mutational impact and prevents accumulation of deleterious events, an inherent problem of clonal reproduction. We conclude that traditional linear homozygous genome representation conceals the complexity when dealing with rearranged and heterozygous clonal cells.

Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis

cis-Prenyltransferases (cis-PTs) constitute a large family of enzymes conserved during evolution and present in all domains of life. cis-PTs catalyze the cis-1,4-polymerization of isoprene units to generate isoprenoids with carbon skeletons varying from C10 (neryl pyrophosphate) to C > 10,000 (natural rubber). Though the previously reported CPTs in Hevea are designated based on sequence variations, their classification was done mostly by phylogenetic analysis using a mixture of partial as well as full length sequences often excluding the UTRs. In this context an attempt was made to reclassify the CPTs strictly based on their sequence similarity and distinguish the members putatively associated with rubber biosynthesis from the others. Extensive computational analysis was carried out on CPT sequences obtained from public resources and whole genome assemblies of Hevea. Based on the results from BLAST analysis, multiple sequence alignments of protein, nucleotide and untranslated regions, open reading frame analysis, gene prediction analysis and sequence length variations, we conclude that there exists mainly three CPTs namely RubCPT1, RubCPT2 and RubCPT3 putatively associated with rubber biosynthesis in Hevea brasiliensis. The rest were categorised as variants of dehydrodolichyl diphosphate synthase (DHDDS) involved in the synthesis of dolichols having short chain isoprenoids. Analysis of the sequence structure of the most highly expressed RubCPT1 in latex revealed the allele richness and diversity of this important variant prevailing in the popular rubber clones. Haplotypes consisting of SNPs with high degree of heterozygosity were also identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor towards the generation of allelic diversity rather than point mutations. Alternatively, gene expression analysis indicated the possibility of association between specific haplotypes and RubCPT1 expression in Hevea clones which may have downstream impact up to the level of rubber production. The conclusions from this study may pave way for the identification and better understanding of CPTs directly involved with natural rubber biosynthesis in Hevea and the SNP data generated may aid in the development of molecular markers putatively associated with yield in rubber.

MMAB, a novel candidate gene to be screened in the molecular diagnosis of Mevalonate Kinase Deficiency

Mevalonate kinase deficiency (MKD) is an autosomal recessive inflammatory disease. Mutations in MVK gene are associated with MKD with modest genotype-phenotype correlation. In spite of recent guidelines indicating specific MVK mutations for the more severe form or the milder one, little is known about MVK variability within and between populations. The aim of this work is to provide supplementary information about MVK variability useful in the molecular diagnosis of MKD, as well as to unravel the presence of novel genes potentially involved as involved in the clinical heterogeneity of MKD phenotype. We used a population-based approach, coupled with Combined Annotation-Dependent Depletion (CADD) score, to analyze the level of genetic variability for common and putatively deleterious MVK variants. We also performed Exome screening with the Illumina Human Exome Bead Chip on 21 MKD patients to double-check our in silico findings. Haplotype block detection in different populations revealed the existence of two blocks in MVK; interestingly, the first haploblock comprises the promoter region shared with MMAB gene. Analyses of MMAB and MVK genetic variants in 21 MKD patients strengthen our observations showing a novel scenario in which the same mutations commonly associated with MKD are found coupled with different combination of MMAB rs7134594 SNP was already described as associated with HDL cholesterol level and present in the haploblock promoter region. The rs7134594 SNP is reported as an eQTL for MVK and MMAB. Hypothesizing the presence of genetic variants modulating the complex phenotypic spectrum of MKD, we suggest that future directions in screening for MKD pathogenic variants should focus both MMAB and MVK genes.

Molecular evolution and functional characterisation of haplotypes of an important rubber biosynthesis gene in Hevea brasiliensis

Hydroxy-methylglutaryl coenzyme-A synthase (HMGS) is a rate-limiting enzyme in the cytoplasmic isoprenoid biosynthesis pathway leading to natural rubber production in Hevea brasiliensis (rubber). Analysis of the structural variants of this gene is imperative to understand their functional significance in rubber biosynthesis so that they can be properly utilised for ongoing crop improvement programmes in Hevea. We report here allele richness and diversity of the HMGS gene in selected popular rubber clones. Haplotypes consisting of single nucleotide polymorphisms (SNPs) from the coding and non-coding regions with a high degree of heterozygosity were identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor to the generation of allelic diversity, rather than point mutations. The evolutionarily conserved nature of some SNPs was identified by comparative DNA sequence analysis of HMGS orthologues from diverse taxa, demonstrating the molecular evolution of rubber biosynthesis genes in general. In silico three-dimensional structural studies highlighting the structural positioning of non-synonymous SNPs from different HMGS haplotypes revealed that the ligand-binding site on the enzyme remains impervious to the reported sequence variations. In contrast, gene expression results indicated the possibility of association between specific haplotypes and HMGS expression in Hevea clones, which may have a downstream impact up to the level of rubber production. Moreover, haplotype diversity of the HMGS gene and its putative association with gene expression can be the basis for further genetic association studies in rubber. Furthermore, the data also show the role of SNPs in the evolution of candidate genes coding for functional traits in plants.

A hidden Markov model for investigating recent positive selection through haplotype structure

Recent positive selection can increase the frequency of an advantageous mutant rapidly enough that a relatively long ancestral haplotype will be remained intact around it. We present a hidden Markov model (HMM) to identify such haplotype structures. With HMM identified haplotype structures, a population genetic model for the extent of ancestral haplotypes is then adopted for parameter inference of the selection intensity and the allele age. Simulations show that this method can detect selection under a wide range of conditions and has higher power than the existing frequency spectrum-based method. In addition, it provides good estimate of the selection coefficients and allele ages for strong selection. The method analyzes large data sets in a reasonable amount of running time. This method is applied to HapMap III data for a genome scan, and identifies a list of candidate regions putatively under recent positive selection. It is also applied to several genes known to be under recent positive selection, including the LCT, KITLG and TYRP1 genes in Northern Europeans, and OCA2 in East Asians, to estimate their allele ages and selection coefficients.

Genetic diversity and haplotype structure of 24 Y-chromosomal STR in Chinese Hui ethnic group and its genetic relationships with other populations

In the present study, 24 Y-chromosomal short tandem repeat (Y-STR) loci were analyzed in 115 unrelated Hui male individuals from Haiyuan county or Tongxin county, Ningxia Hui Autonomous Region, China, to evaluate the forensic application of the 24 STR loci and to analyze interpopulation differentiations by making comparisons between the Hui group data and previously published data of other 13 populations. A total of 115 different haplotypes were observed on these 24 Y-STR loci. The gene diversities ranged from 0.4049 (DYS437) to 0.9729 (DYS385a, b). The overall haplotype diversity was 1 at AGCU 24 Y-STR loci level, while the values were reduced to 0.999237, 0.996949, and 0.996644 at the Y-filer 17 loci, 11 Y-STR loci of extended haplotype and 9 Y-STR loci of minimal haplotype levels, respectively; whereas, haplotype diversity for additional 7 loci (not included in Y-filer 17 loci) was 0.995271. The pairwise FST , multidimensional scaling plot and neighbor-joining tree indicated the Hui group had the closest genetic relationship with Sala in the paternal lineage in the present study. In summary, the results in our study indicated the 24 Y-STRs had a high level of polymorphism in Hui group and hence could be a powerful tool for forensic application and population genetic study.