Human Chromosome Research Articles

Exploring 3D community inconsistency in human chromosome contact networks

Abstract Researchers have developed chromosome capture methods such as Hi-C to better understand DNA’s 3D folding in nuclei. The Hi-C method captures contact frequencies between DNA segment pairs across the genome. When analyzing Hi-C data sets, it is common to group these pairs using standard bioinformatics methods (e.g. PCA). Other approaches handle Hi-C data as weighted networks, where connected node pairs represent DNA segments in 3D proximity. In this representation, one can leverage community detection techniques developed in complex network theory to group nodes into mesoscale communities containing nodes with similar connection patterns. While there are several successful attempts to analyze Hi-C data in this way, it is common to report and study the most typical community structure. But in reality, there are often several valid candidates. Therefore, depending on algorithm design, different community detection methods focusing on slightly different connectivity features may have differing views on the ideal node groupings. In fact, even the same community detection method may yield different results if using a stochastic algorithm. This ambiguity is fundamental to community detection and shared by most complex networks whenever interactions span all scales in the network. This is known as community inconsistency. This paper explores this inconsistency of 3D communities in Hi-C data for all human chromosomes. We base our analysis on two inconsistency metrics, one local and one global, and quantify the network scales where the community separation is most variable. For example, we find that TADs are less reliable than A/B compartments and that nodes with highly variable node-community memberships are associated with open chromatin. Overall, our study provides a helpful framework for data-driven researchers and increases awareness of some inherent challenges when clustering Hi-C data into 3D communities.

Expert consensus on the standardization and clinical application of genotyping for human leukocyte antigens

The genes encoding human leukocyte antigens (HLA) are located in the region of 6p21.31 of the human chromosome, accounting for 1/3 000 of the human genome. The area is seen as a structurally complex, with the characteristics of multi-gene inheritance, high polymorphism, haplotype inheritance, linkage disequilibrium, racial diversity, and regional distribution. In recent years, with the continuous innovation of HLA genotyping technology, new issues and challenges have emerged in terms of the standardization of the technology platform, the implementation of the performance validation and verification plan for clinical testing projects, quality control of the entire testing process, analysis and interpretation of results, consultation and selection of donors, and clinical application. This consensus, in combination with the current status of HLA genotyping technology in clinical practice in China, systematically sorts out the key contents in the standardization and clinical application of the HLA genotyping technology platform. The first part proposes the basic requirements for the standardization of the HLA typing technology platform from the perspectives of personnel, instrumentation, reagents and consumables, experimental methods, environmental settings, and introduces the methods and performance indicators for confirming the performance of HLA genotyping with next-generation sequencing (NGS) technology. The second part elaborates on the key points of quality control during the entire process of HLA genotyping procedure, quality control of result reporting and releasing, and management of result interpretation and bioinformatics systems. The last part describes the clinical consultation and application promotion of HLA genotyping. After widely soliciting revisions and discussions from domestic experts in the field of HLA and transplantation, the consensus aims to promote the development and clinical application of the HLA gene in related medical field.

The Long Non-Coding RNA ANRIL in Cancers.

ANRIL (Antisense Noncoding RNA in the INK4 Locus), a long non-coding RNA encoded in the human chromosome 9p21 region, is a critical factor for regulating gene expression by interacting with multiple proteins and miRNAs. It has been found to play important roles in various cellular processes, including cell cycle control and proliferation. Dysregulation of ANRIL has been associated with several diseases like cancers and cardiovascular diseases, for instance. Understanding the oncogenic role of ANRIL and its potential as a diagnostic and prognostic biomarker in cancer is crucial. This review provides insights into the regulatory mechanisms and oncogenic significance of the 9p21 locus and ANRIL in cancer.

The emerging role of LncRNA AWPPH in multiple cancers: a review study.

Long non-coding RNAs (lncRNAs) are transcribed RNA molecules longer than 200 nucleotides in length that have no protein-coding potential. They are able to react with DNA, RNA, and protein. Hence they involve in regulating gene expression at the epigenetic, transcriptional, post-transcriptional, and translational levels. LncRNAs have been proven to play an important role in human malignancies and prognostic outcomes. In this review, we will comprehensively and functionally discuss the role of a novel identified lncRNA, namely lncRNA WAPPH located on human chromosome 2q13, in various cancers. Increasing research studies have shown that lncRNA AWPPH is deregulated in different malignancies, including breast cancer, gastric cancer, colorectal cancer, ovarian cancer, bladder cancer, leukemia, and others. LncRNA WAPPH serves as an oncogene in tumorigenesis and the development of cancer. Moreover, lncRNA AWPPH is involved in numerous biological processes of solid and blood cancers. Taken together, based on our scrutiny analysis, lncRNA AWPPH can be regarded as a putative biomarker for diagnosis or therapeutic target in human malignancies.

The APOE-TOMM40 Humanized Mouse Model: Characterization of Age, Sex, and PolyT Variant Effects on Gene Expression.

The human chromosome 19q13.32 is a gene rich region and has been associated with multiple phenotypes, including late onset Alzheimer's disease (LOAD) and other age-related conditions. Here we developed the first humanized mouse model that contains the entire TOMM40 and APOE genes with all intronic and intergenic sequences including the upstream and downstream regions. Thus, the mouse model carries the human TOMM40 and APOE genes and their intact regulatory sequences. We generated the APOE-TOMM40 humanized mouse model in which the entire mouse region was replaced with the human (h)APOE-TOMM40 loci including their upstream and downstream flanking regulatory sequences using recombineering technologies. We then measured the expression of the human TOMM40 and APOE genes in the mice brain, liver, and spleen tissues using TaqMan based mRNA expression assays. We investigated the effects of the '523' polyT genotype (S/S or VL/VL), sex, and age on the human TOMM40- and APOE-mRNAs expression levels using our new humanized mouse model. The analysis revealed tissue specific and shared effects of the '523' polyT genotype, sex, and age on the regulation of the human TOMM40 and APOE genes. Noteworthy, the regulatory effect of the '523' polyT genotype was observed for all studied organs. The model offers new opportunities for basic science, translational, and preclinical drug discovery studies focused on the APOE genomic region in relation to LOAD and other conditions in adulthood.

Roles of increased NUCKS1 expression in endometriosis

BackgroundEndometriosis is still a difficult problem for women. The Nuclear Ubiquitous Casein and cyclin-dependent Kinase Substrate 1 (NUCKS1) gene is located on human chromosome 1q32.1. It encodes the NUCKS1 protein, a 27 kDa nuclear DNA binding protein that plays an important role in cell growth and proliferation. NUCKS1 plays an important role in the development of many diseases. However, its role in endometriosis is unclear.MethodsEctopic endometrial tissues and normal tissue specimens were collected, and the expression of NUCKS1, NF-κB and PI3K was detected by RT-qPCR and immunohistochemistry. Inhibition of NUCKS1 in hEM15A cells, study the changes in cell viability, apoptosis, migration and protein expression by CCK8 assay, flow cytometry, wound-healing assay, western blot and ELISA techniques. The comparison of differences between the two groups was implemented using unpaired sample t test or Mann-whitney U test. One-way analysis of variance or Kruskal-wallis test was used for comparisons among the three groups.Results(1) NUCKS1 is highly expressed in endometriosis tissues. (2) Inhibition of NUCKS1 decreases cell viability and capability of migration, and increases apoptosis in endometriosis cells. (3) Expressions of NF-κB and PI3K are increased in endometriosis tissues, and inhibition of NUCKS1 decreases the expression levels of PI3K and NF-κB in endometriosis cells. (4) Inhibition of NUCKS1 decreases the expression of VEGF.Conclusion(1) NUCKS1 is overexpressed in endometriosis, and inhibition of NUCKS1 inhibits cell viability and capability of migration, and increases apoptosis. (2) NUCKS1 promotes the progress of endometriosis through activating PI3K and NF-κB pathways, and VEFG is also involved in this process.

Understanding the genetic mechanisms and cognitive impairments in Down syndrome: towards a holistic approach.

The most common genetic cause of intellectual disability is Down syndrome (DS), trisomy 21. It commonly results from three copies of human chromosome 21 (HC21). There are no mutations or deletions involved in DS. Instead, the phenotype is caused by altered transcription of the genes on HC21. These transcriptional variations are responsible for a myriad of symptoms affecting every organ system. A very debilitating aspect of DS is intellectual disability (ID). Although tremendous advances have been made to try and understand the underlying mechanisms of ID, there is a lack of a unified, holistic view to defining the cause and managing the cognitive impairments. In this literature review, we discuss the mechanisms of neuronal over-inhibition, abnormal morphology, and other genetic factors in contributing to the development of ID in DS patients and to gain a holistic understanding of ID in DS patients. We also highlight potential therapeutic approaches to improve the quality of life of DS patients.

Oncogenic Ras and ΔNp63α cooperate to recruit immunosuppressive polymorphonuclear myeloid-derived suppressor cells in a mouse model of squamous cancer pathogenesis.

Amplification of human chromosome 3q26-29, which encodes oncoprotein ΔNp63 among other isoforms of the p63 family, is a feature common to squamous cell carcinomas (SCCs) of multiple tissue origins. Along with overexpression of ΔNp63, activation of the protooncogene, RAS, whether by overexpression or oncogenic mutation, is frequently observed in many cancers. In this study, analysis of transcriptome data from The Cancer Genome Atlas (TCGA) demonstrated that expression of TP63 mRNA, particularly ΔNp63 isoforms, and HRAS are significantly elevated in advanced squamous cell carcinomas of the head and neck (HNSCCs), suggesting pathological significance. However, how co-overexpressed ΔNp63 and HRAS affect the immunosuppressive tumor microenvironment (TME) is incompletely understood. Here, we established and characterized an immune competent mouse model using primary keratinocytes with retroviral-mediated overexpression of ΔNp63α and constitutively activated HRAS (v-rasHa G12R) to evaluate the role of these oncogenes in the immune TME. In this model, orthotopic grafting of wildtype syngeneic keratinocytes expressing both v-rasHa and elevated levels of ΔNp63α consistently yield carcinomas in syngeneic hosts, while cells expressing v-rasHa alone yield predominantly papillomas. We found that polymorphonuclear (PMN) myeloid cells, experimentally validated to be immunosuppressive and thus representing myeloid-derived suppressor cells (PMN-MDSCs), were significantly recruited into the TME of carcinomas arising early following orthotopic grafting of ΔNp63α/v-rasHa-expressing keratinocytes. ΔNp63α/v-rasHa-driven carcinomas expressed higher levels of chemokines implicated in recruitment of MDSCs compared to v-rasHa-initiated tumors, providing a heretofore undescribed link between ΔNp63α/HRAS-driven carcinomas and the development of an immunosuppressive TME. These results support the utilization of a genetic carcinogenesis model harboring specific genomic drivers of malignancy to study mechanisms underlying the development of local immunosuppression.

A variant in the LDL receptor-related protein encoding gene LRP4 underlying polydactyly and phalangeal synostosis in a family of Pakistani origin.

A family of Pakistani origin, segregating polydactyly, and phalangeal synostosis in an autosomal dominant manner, has been investigated and presented in the present report. Whole-exome sequencing (WES), followed by segregation analysis using Sanger sequencing, revealed a heterozygous missense variant [c.G1696A, p.(Gly566Ser)] in the LRP4 gene located on human chromosome 11p11.2. Homology protein modeling revealed the mutant Ser566 generated new interactions with at least four other amino acids and disrupted protein folding and function. Our findings demonstrated the first direct evidence of involvement of LRP4 in causing polydactyly and phalangeal synostosis in the same family. This study highlighted the importance of inclusion of LRP4 gene in screening individuals presenting polydactyly in hands and feet, and phalangeal synostosis in the same family.

Race-specific association of an IRGM risk allele with cytokine expression in human subjects

Immunity-related GTPase family M (IRGM), located on human chromosome 5q33.1, encodes a protein that promotes autophagy and suppresses the innate immune response. The minor allele of rs13361189 (−4299T>C), a single nucleotide polymorphism in the IRGM promoter, has been associated with several diseases, including Crohn’s disease and tuberculosis. Although patterns of linkage disequilibrium and minor allele frequency for this polymorphism differ dramatically between subjects of European and African descent, studies of rs13361189 have predominantly been conducted in Europeans and the mechanism of association is poorly understood. We recruited a cohort of 68 individuals (30 White, 34 African American, 4 other race) with varying rs13361189 genotypes and assessed a panel of immune response measures including whole blood cytokine induction following ex vivo stimulation with Toll-like Receptor ligands. Minor allele carriers were found to have increased serum immunoglobulin M, C-reactive protein, and circulating CD8+ T cells. No differences in whole blood cytokines were observed between minor allele carriers and non-carriers in the overall study population; however, minor allele status was associated with increased induction of a subset of cytokines among African American subjects, and decreased induction among White subjects. These findings underline the importance of broad racial inclusion in genetic studies of immunity.

Emergence of the Synucleins.

This study establishes the origin and evolutionary history of the synuclein genes. A combination of phylogenetic analyses of the synucleins from twenty-two model species, characterization of local synteny similarities among humans, sharks and lampreys, and statistical comparisons among lamprey and human chromosomes, provides conclusive evidence for the current diversity of synuclein genes arising from the whole-genome duplications (WGDs) that occurred in vertebrates. An ancestral synuclein gene was duplicated in a first WGD, predating the diversification of all living vertebrates. The two resulting genes are still present in agnathan vertebrates. The second WGD, specific to the gnathostome lineage, led to the emergence of the three classical synuclein genes, SNCA, SNCB and SNCG, which are present in all jawed vertebrate lineages. Additional WGDs have added new genes in both agnathans and gnathostomes, while some gene losses have occurred in particular species. The emergence of synucleins through WGDs prevented these genes from experiencing dosage effects, thus avoiding the potential detrimental effects associated with individual duplications of genes that encode proteins prone to aggregation. Additional insights into the structural and functional features of synucleins are gained through the analysis of the highly divergent synuclein proteins present in chondrichthyans and agnathans.

Chaining for accurate alignment of erroneous long reads to acyclic variation graphs.

Aligning reads to a variation graph is a standard task in pangenomics, with downstream applications such as improving variant calling. While the vg toolkit [Garrison et al. (Variation graph toolkit improves read mapping by representing genetic variation in the reference. Nat Biotechnol 2018;36:875-9)] is a popular aligner of short reads, GraphAligner [Rautiainen and Marschall (GraphAligner: rapid and versatile sequence-to-graph alignment. Genome Biol 2020;21:253-28)] is the state-of-the-art aligner of erroneous long reads. GraphAligner works by finding candidate read occurrences based on individually extending the best seeds of the read in the variation graph. However, a more principled approach recognized in the community is to co-linearly chain multiple seeds. We present a new algorithm to co-linearly chain a set of seeds in a string labeled acyclic graph, together with the first efficient implementation of such a co-linear chaining algorithm into a new aligner of erroneous long reads to acyclic variation graphs, GraphChainer. We run experiments aligning real and simulated PacBio CLR reads with average error rates 15% and 5%. Compared to GraphAligner, GraphChainer aligns 12-17% more reads, and 21-28% more total read length, on real PacBio CLR reads from human chromosomes 1, 22, and the whole human pangenome. On both simulated and real data, GraphChainer aligns between 95% and 99% of all reads, and of total read length. We also show that minigraph [Li et al. (The design and construction of reference pangenome graphs with minigraph. Genome Biol 2020;21:265-19.)] and minichain [Chandra and Jain (Sequence to graph alignment using gap-sensitive co-linear chaining. In: Proceedings of the 27th Annual International Conference on Research in Computational Molecular Biology (RECOMB 2023). Springer, 2023, 58-73.)] obtain an accuracy of <60% on this setting. GraphChainer is freely available at https://github.com/algbio/GraphChainer. The datasets and evaluation pipeline can be reached from the previous address.

Concentration of inverted repeats along human DNA.

This work aims to describe the observed enrichment of inverted repeats in the human genome; and to identify and describe, with detailed length profiles, the regions with significant and relevant enriched occurrence of inverted repeats. The enrichment is assessed and tested with a recently proposed measure (z-scores based measure). We simulate a genome using an order 7 Markov model trained with the data from the real genome. The simulated genome is used to establish the critical values which are used as decision thresholds to identify the regions with significant enriched concentrations. Several human genome regions are highly enriched in the occurrence of inverted repeats. This is observed in all the human chromosomes. The distribution of inverted repeat lengths varies along the genome. The majority of the regions with severely exaggerated enrichment contain mainly short length inverted repeats. There are also regions with regular peaks along the inverted repeats lengths distribution (periodic regularities) and other regions with exaggerated enrichment for long lengths (less frequent). However, adjacent regions tend to have similar distributions.

Phenotypic features of dystrophin gene knockout pigs harboring a human artificial chromosome containing the entire dystrophin gene

Mammalian artificial chromosomes have enabled the introduction of extremely large amounts of genetic information into animal cells in an autonomously replicating, nonintegrating format. However, the evaluation of human artificial chromosomes (HACs) as novel tools for curing intractable hereditary disorders has been hindered by the limited efficacy of the delivery system. We generated dystrophin gene knockout (DMD-KO) pigs harboring the HAC bearing the entire human DMD via asomatic cell cloning procedure (DYS-HAC-cloned pig). Restored human dystrophin expression was confirmed by immunofluorescence staining in the skeletal muscle of the DYS-HAC-cloned pigs. Viability at the first month postpartum of the DYS-HAC-cloned pigs, including motor function in the hind leg and serum creatinine kinase level, was improved significantly when compared with that in the original DMD-KO pigs. However, decrease in systemic retention of the DYS-HAC vector and limited production of the DMD protein might have caused severe respiratory impairment with general prostration by 3months postpartum. The results demonstrate that the use of transchromosomic cloned pigs permitted a straightforward estimation of the efficacy of the DYS-HAC carried in affected tissues/organs in a large-animal disease model, providing novel insights into the therapeutic application of exogenous mammalian artificial chromosomes.

Overexpression of the autism candidate gene Cyfip1 pathologically enhances olivo-cerebellar signaling in mice.

Cyfip1, the gene encoding cytoplasmic FMR1 interacting protein 1, has been of interest as an autism candidate gene for years. A potential role in autism spectrum disorder (ASD) is suggested by its location on human chromosome 15q11-13, an instable region that gives rise to a variety of copy number variations associated with syndromic autism. In addition, the CYFIP1 protein acts as a binding partner to Fragile X Messenger Ribonucleoprotein (FMRP) in the regulation of translation initiation. Mutation of FMR1, the gene encoding FMRP, causes Fragile X syndrome, another form of syndromic autism. Here, in mice overexpressing CYFIP1, we study response properties of cerebellar Purkinje cells to activity of the climbing fiber input that originates from the inferior olive and provides an instructive signal in sensorimotor input analysis and plasticity. We find that CYFIP1 overexpression results in enhanced localization of the synaptic organizer neurexin 1 (NRXN1) at climbing fiber synaptic input sites on Purkinje cell primary dendrites and concomitant enhanced climbing fiber synaptic transmission (CF-EPSCs) measured using whole-cell patch-clamp recordings from Purkinje cells in vitro. Moreover, using two-photon measurements of GCaMP6f-encoded climbing fiber signals in Purkinje cells of intact mice, we observe enhanced responses to air puff stimuli applied to the whisker field. These findings resemble our previous phenotypic observations in a mouse model for the human 15q11-13 duplication, which does not extend to the Cyfip1 locus. Thus, our study demonstrates that CYFIP1 overexpression shares a limited set of olivo-cerebellar phenotypes as those resulting from an increased number of copies of non-overlapping genes located on chromosome 15q11-13.

PCP4 Promotes Alzheimer's Disease Pathogenesis by Affecting Amyloid-β Protein Precursor Processing.

Down syndrome (DS) is caused by an extra copy of all or part of chromosome 21. The patients with DS develop typical Alzheimer's disease (AD) neuropathology, indicating the role of genes on human chromosome 21 (HSA21) in the pathogenesis of AD. Purkinje cell protein 4 (PCP4), also known as brain-specific protein 19, is a critical gene located on HSA21. However, the role of PCP4 in DS and AD pathogenesis is not clear. To explore the role of PCP4 in amyloid-β protein precursor (AβPP) processing in AD. In this study, we investigated the role of PCP4 in AD progression in vitro and in vivo. In vitro experiments, we overexpressed PCP4 in human Swedish mutant AβPP stable expression or neural cell lines. In vitro experiments, APP23/PS45 double transgenic mice were selected and treated with AAV-PCP4. Multiple topics were detected by western blot, RT-PCR, immunohistochemical and behavioral test. We found that PCP4 expression was altered in AD. PCP4 was overexpressed in APP23/PS45 transgenic mice and PCP4 affected the processing of AβPP. The production of amyloid-β protein (Aβ) was also promoted by PCP4. The upregulation of endogenous AβPP expression and the downregulation of ADAM10 were due to the transcriptional regulation of PCP4. In addition, PCP4 increased Aβ deposition and neural plaque formation in the brain, and exuberated learning and memory impairment in transgenic AD model mice. Our finding reveals that PCP4 contributes to the pathogenesis of AD by affecting AβPP processing and suggests PCP4 as a novel therapeutic target for AD by targeting Aβ pathology.

Treatment of CHO cells with Taxol and reversine improves micronucleation and microcell-mediated chromosome transfer efficiency

Microcell-mediated chromosome transfer is an attractive technique for transferring chromosomes from donor cells to recipient cells and has enabled the generation of cell lines and humanized animal models that contain megabase-sized gene(s). However, improvements in chromosomal transfer efficiency are still needed to accelerate the production of these cells and animals. The chromosomal transfer protocol consists of micronucleation, microcell formation, and fusion of donor cells with recipient cells. We found that the combination of Taxol (paclitaxel) and reversine rather than the conventional reagent colcemid resulted in highly efficient micronucleation and substantially improved chromosomal transfer efficiency from Chinese hamster ovary donor cells to HT1080 and NIH3T3 recipient cells by up to 18.3- and 4.9-fold, respectively. Furthermore, chromosome transfer efficiency to human induced pluripotent stem cells, which rarely occurred with colcemid, was also clearly improved after Taxol and reversine treatment. These results might be related to Taxol increasing the number of spindle poles, leading to multinucleation and delaying mitosis, and reversine inducing mitotic slippage and decreasing the duration of mitosis. Here, we demonstrated that an alternative optimized protocol improved chromosome transfer efficiency into various cell lines. These data advance chromosomal engineering technology and the use of human artificial chromosomes in genetic and regenerative medical research.

Surfactant Metabolism Dysfunction Type 3 (SMDP3)

Surfactant (surface-active-agent) is a compound of phospholipids and proteins which are synthesized and secreted into the alveolus by type II epithelial cells, where it functions to decrease surface tension, maintaining alveolar expansion, to facilitate pulmonary compliance. Surfactant proteins (SP)A, B, C, D represent around 8% of total components, but has vital role in optimizing rapid adsorption and spreading of phospholipids. ATP-binding cassette sub-family A member 3, protein that encoded by ABCA3 gene, which located in human chromosome 16p13.3, is synthesized in endoplasmic reticulum and migrated to lysosomal-derived organelles of alveolar type II cells, formally known as lamellar bodies. Once accumulate into the membrane, ABCA3 can directed surfactant phospholipid into the lumen of lamellar bodies and create tight packed of surfactant lipids and proteins. Mutation of the ATP-binding cassette transporter gene ABCA3 cause failure in lamellar body synthesis and result in decreased production of surfactant, along with respiratory distress syndrome, and fatal respiratory failure.

Deep-Learning-Based Human Chromosome Classification: Data Augmentation and Ensemble

Object classification is a crucial task in deep learning, which involves the identification and categorization of objects in images or videos. Although humans can easily recognize common objects, such as cars, animals, or plants, performing this task on a large scale can be time-consuming and error-prone. Therefore, automating this process using neural networks can save time and effort while achieving higher accuracy. Our study focuses on the classification step of human chromosome karyotyping, an important medical procedure that helps diagnose genetic disorders. Traditionally, this task is performed manually by expert cytologists, which is a time-consuming process that requires specialized medical skills. Therefore, automating it through deep learning can be immensely useful. To accomplish this, we implemented and adapted existing preprocessing and data augmentation techniques to prepare the chromosome images for classification. We used ResNet-50 convolutional neural network and Swin Transformer, coupled with an ensemble approach to classify the chromosomes, obtaining state-of-the-art performance in the tested dataset.

Random Integration Analysis of Recombinant Adeno-Associated Virus 6 Packaged in Sf9 Insect Cells

Recently, there have been growing concerns over the integration of recombinant adeno-associated virus (rAAV) used in gene therapy. Wild-type adeno-associated virus (AAV) site specifically integrates into AAVS1 site of human genome, while rAAV randomly integrates into host chromosomes at low frequencies. This research aims to study the random integration events of rAAV6-EGFP packaged in Sf9 insect cells. Baculo-Sf9 manufacturing platform has the advantages of high-density suspension culture of Sf9 insect cells and large-scale production of rAAV vectors. In this study, we used different doses of Baculo-Sf9 produced rAAV6-EGFP to transduce HEK293T cells and A549-implanted tumors in vitro and in vivo. Using flow cytometry and fluorescence microscopy, we studied their EGFP gene expression efficiencies and EGFP fluorescence intensities. Using inverse nested PCR and DNA sequencing, random integration sites of rAAV6-EGFP genome into human chromosomes were identified. In vitro results showed that gene expression efficiencies became stable after 20 days and random integration frequencies were 0.2‒4.2%. Both in vitro and in vivo results indicated that random integration of Baculo-Sf9 rAAV6 was dose-dependent. Sequencing results showed two random integration sites, which were on human chromosomes 8 and 12. The findings suggest that we should use as low dose of rAAV vector as possible for safe gene therapy.