Duplicated Region Research Articles

An efficient CRISPR-Cas9 enrichment sequencing strategy for characterizing complex and highly duplicated genomic regions. A case study in the Prunus salicina LG3-MYB10 genes cluster

BackgroundGenome complexity is largely linked to diversification and crop innovation. Examples of regions with duplicated genes with relevant roles in agricultural traits are found in many crops. In both duplicated and non-duplicated genes, much of the variability in agronomic traits is caused by large as well as small and middle scale structural variants (SVs), which highlights the relevance of the identification and characterization of complex variability between genomes for plant breeding.ResultsHere we improve and demonstrate the use of CRISPR-Cas9 enrichment combined with long-read sequencing technology to resolve the MYB10 region in the linkage group 3 (LG3) of Japanese plum (Prunus salicina). This region, which has a length from 90 to 271 kb according to the P. salicina genomes available, is associated with fruit color variability in Prunus species. We demonstrate the high complexity of this region, with homology levels between Japanese plum varieties comparable to those between Prunus species. We cleaved MYB10 genes in five plum varieties using the Cas9 enzyme guided by a pool of crRNAs. The barcoded fragments were then pooled and sequenced in a single MinION Oxford Nanopore Technologies (ONT) run, yielding 194 Mb of sequence. The enrichment was confirmed by aligning the long reads to the plum reference genomes, with a mean read on-target value of 4.5% and a depth per sample of 11.9x. From the alignment, 3261 SNPs and 287 SVs were called and phased. A de novo assembly was constructed for each variety, which also allowed detection, at the haplotype level, of the variability in this region.ConclusionsCRISPR-Cas9 enrichment is a versatile and powerful tool for long-read targeted sequencing even on highly duplicated and/or polymorphic genomic regions, being especially useful when a reference genome is not available. Potential uses of this methodology as well as its limitations are further discussed.

SHFLD3 phenotypes caused by 17p13.3 triplication/ duplication encompassing Fingerin (BHLHA9) invariably

BackgroundSplit-hand/ foot malformation with long bone deficiency 3 (SHFLD3) is an extremely rare condition associated with duplications located on 17p13.3, which invariably encompasses the BHLHA9 gene. The disease inherits with variable expressivity and significant incomplete penetrance as high as 50%.ResultsWe have detected 17p13.3 locus one-allele triplication in a male proband from family 1 (F1.1), and duplication in a male proband from family 2 (F2.1) applying array comparative genomic hybridization (array CGH). The rearrangements mapped to the following chromosomal regions–arr[GRCh38] 17p13.3(960254–1291856)×4 in F1.1 and arr[GRCh38] 17p13.3(1227482–1302716)×3 in F2.1. The targeted quantitative PCR revealed that the 17p13.3 locus was also duplicated in the second affected member from family 2 (F2.2; brother of F2.1). In the next step, we performed segregation studies using quantitative PCR and revealed that F1.1 inherited the triplication from his healthy father—F1.2, whereas the locus was unremarkable in the mother of F2.1 & F2.2 and the healthy son of F2.1. However, the duplication was present in a healthy daughter of F2.2, an asymptomatic carrier. The breakpoint analysis allowed to define the exact size and span of the duplicated region in Family 2, i.e., 78,948 bp chr17:1225063–1304010 (HG38). Interestingly, all symptomatic carriers from both families presented with variable SHFLD3 phenotype. The involvement of secondary modifying locus could not be excluded, however, the Sanger sequencing screening of BHLHA9 entire coding sequence was unremarkable for both families.ConclusionsWe have shed light on the one-allele CNV triplication occurrence that should be considered when a higher probe (over duplication range) signal is noted. Second, all SHFLD3 patients were accurately described regarding infrequent limb phenotypes, which were highly variable even when familial. Of note, all symptomatic individuals were males. SHFLD3 still remains a mysterious ultra-rare disease and our findings do not answer crucial questions regarding the disease low penetrance, variable expression and heterogeneity. However, we have presented some clinical and molecular aspects that may be helpful in daily diagnostic routine, both dysmorphological and molecular assessment, of patients affected with SHFLD3.

A deletion/duplication in the Ligon lintless-2 locus induces siRNAs that inhibit cotton fiber cell elongation.

Most cultivated cotton (Gossypium hirsutum L.) varieties have two types of seed fibers: short fuzz fiber strongly adhered to the seed coat, and long lint fiber used in the textile industry. The Ligon lintless-2 (Li2) cotton mutant has a normal vegetative phenotype but produces very short lint fiber on the seeds. The Li2 mutation is controlled by a single dominant gene. We discovered a large structural rearrangement at the end of chromosome D13 in the Li2 mutant based on whole-genome sequencing and genetic mapping of segregating populations. The rearrangement contains a 177-kb deletion and a 221-kb duplication positioned as a tandem inverted repeat. The gene Gh_D13G2437 is located at the junction of the inverted repeat in the duplicated region. During transcription such structure spontaneously forms self-complementary hairpin RNA of Gh_D13G2437 followed by production of small interfering RNA (siRNA). Gh_D13G2437 encodes a Ran-Binding Protein 1 (RanBP1) that preferentially expresses during cotton fiber elongation. The abundance of siRNA produced from Gh_D13G2437 reciprocally corresponds with the abundance of highly homologous (68%-98% amino acid sequence identity) RanBP1 family transcripts during fiber elongation, resulting in a shorter fiber phenotype in the Li2. Overexpression of Gh_D13G2437 in the Li2 mutant recovered the long lint fiber phenotype. Taken together, our findings revealed that siRNA-induced silencing of a family of RanBP1s inhibit elongation of cotton fiber cells in the Li2 mutant.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae.

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication.

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

CSWS—Suppress confusion subregion in indoor image fingerprint localization

AbstractImage data can provide rich content information, which has attracted a lot of attention in the field of indoor fingerprint positioning. However, indoor image information from different locations is characterized by high content repetition, and these repetitive regions can cause the problem of insufficient differentiation of adjacent image fingerprints or even fingerprint misclassification. To solve this problem, this article proposes a confusion subregion weighted suppression strategy in the image fingerprint database. First, duplicate regions (confusion subregions) are extracted from the fingerprint database using an appropriate salient region detection method. Then the similarity of these duplicate regions in Euclidean space is defined, and the degree of influence of these regions on the distinguishability of the fingerprint database is measured. Finally, the suppression of these duplicate regions is achieved in the original image fingerprint database by introducing weighted suppression coefficients. Experiments show that the algorithm proposed in this article can achieve significant results in the fingerprint localization task of real indoor scenes and effectively improve the localization accuracy.

Evaluation of the Effect of Gene Duplication by Genome Editing on Drug Resistance in Plasmodium falciparum.

The emergence and spread of drug-resistant Plasmodium falciparum have compromised antimalarial efficacy and threatened the global malaria elimination campaign using artemisinin combination therapies. The impacts of amino acid substitutions in antimalarial drug resistance-associated genes on drug susceptibility have been investigated; however, the effects of amplification of those genes remain unexplored due to the lack of robust genetic approaches. Here, we generated transgenic P. falciparum parasites with an additional copy of a drug resistance-associated gene using the highly efficient CRISPR/Cas9 system and investigated their drug response. Insertion of a drug resistance-associated gene expression cassette in the genome resulted in a roughly twofold increase of mRNA levels of the target gene mdr1, which encodes multidrug resistance protein 1. The gene duplication event contributed to resistance to mefloquine, lumefantrine, and dihydroartemisinin, while the duplication of a genomic region encoding plasmepsin 2 and plasmepsin 3 did not affect resistance to antimalarial drugs, including piperaquine. We further demonstrated that mdr1 mRNA expression levels are strongly associated with mefloquine resistance in several field-derived P. falciparum lines with various genetic backgrounds. This study provides compelling evidence that mdr1 could serve as a molecular marker for the surveillance of mefloquine-resistant parasites. Long DNA integration into parasite genomes using the CRISPR/Cas9 system provides a useful tool for the evaluation of the effect of copy number variation on drug response.

Multiple Recurrent Copy Number Variations (CNVs) in Chromosome 22 Including 22q11.2 Associated with Autism Spectrum Disorder.

IntroductionAutism spectrum disorder (ASD) is a developmental disorder that can cause substantial social, communication, and behavioral challenges. Genetic factors play a significant role in ASD, where the risk of ASD has been increased for unclear reasons. Twin studies have shown important evidence of both genetic and environmental contributions in ASD, where the level of contribution of these factors has not been proven yet. It has been suggested that copy number variation (CNV) duplication and the deletion of many genes in chromosome 22 (Ch22) may have a strong association with ASD. This study screened the CNVs in Ch22 in autistic Saudi children and assessed the candidate gene in the CNVs region of Ch22 that is most associated with ASD.MethodsThis study included 15 autistic Saudi children as well as 4 healthy children as controls; DNA was extracted from samples and analyzed using array comparative genomic hybridization (aCGH) and DNA sequencing.ResultsThe aCGH detected (in only 6 autistic samples) deletion and duplication in many regions of Ch22, including some critical genes. Moreover, DNA sequencing determined a genetic mutation in the TBX1 gene sequence in autistic samples. This study, carried out using aCGH, found that six autistic patients had CNVs in Ch22, and DNA sequencing revealed mutations in the TBX1 gene in autistic samples but none in the control.ConclusionCNV deletion and the duplication of the TBX1 gene could be related to ASD; therefore, this gene needs more analysis in terms of expression levels.

Passive Framework of Sparse Region Duplication Detection from Digital Images

Currently, digital images are widely communicated by media using social media applications. The general public captures the digital images for preserving the family and personal memories and to share with their friends and family. Digital images have been used extensively in forensic science to present the digital images as proof in the court and law enforcement agencies, which present a loophole for the culprits to forge the digital image and change the proofs and evidence. Copy-move forgery (CMF) is among the most widely employed image manipulation methods. In this method, the area of the image is duplicated to some other part to modify its content by applying different postprocessing operations on images like blurring, color reduction, and scaling which is a challenging research problem in copy-move forgery detection (CMFD). In this paper, an efficient and effective CMFD method is presented to identify the single and multiple altered areas in an image in the presence of postprocessing operations. The proposed CMFD method divides the image into circular blocks. It computes a rotation-invariant feature vector from each circular block of the image by applying local intensity order pattern (LIOP) features. The computed feature vectors are then compared using Euclidean distance to locate the suspected image’s forged areas. The experimental results of the proposed CMFD method are reported on three standard datasets of the CMF, namely, CoMoFoD, KLTCI, and MICC-F220. The experimental analysis of the proposed CMFD method on these datasets indicates that it produces robust performance (detection accuracies of 97.29% on the CoMoFoD dataset, 98.53% on the KLTCI dataset, and 97.57% on the MICC-F220 dataset) as compared with state-of-the-art CMFD methods in terms of the standard performance evaluation parameters of the CMF.

Exploring the characteristics and most bothersome symptoms in MECP2 duplication syndrome to pave the path toward developing parent-oriented outcome measures.

BackgroundMECP2 Duplication Syndrome (MDS), resulting from the duplication of Xq28 region, including MECP2, is a rare disorder with a nascent understanding in clinical features and severity. Studies using antisense oligonucleotides revealed a broad phenotypic rescue in transgenic mice. With human clinical trials on the horizon, there is a need to develop clinical outcome measures for MDS.MethodsWe surveyed caregivers of MDS individuals to explore the frequency and severity of MDS clinical features, and identify the most meaningful symptoms/domains that need to be included in the outcome measure scales.ResultsA total of 101 responses were eligible for the survey. The top six most meaningful symptoms to caregivers in descending order included epilepsy, gross motor, fine motor, communication, infection, and constipation problems. Epilepsy was present in 58.4% of the subjects and 75% were drug‐resistant, Furthermore, ~12% required intensive care unit (ICU) admission. Infections were present in 55% of the subjects, and one‐fourth of them required ICU admission. Constipation was present in ~85% of the subjects and one‐third required enemas/suppositories.ConclusionOur study is one of the largest cohorts conducted on MDS individuals characterizing the frequency and severity of MDS symptoms. Additionally, these study results will contribute to establishing a foundation to develop parent‐reported outcomes in MDS.

Integrative Analysis of the GRAS Genes From Chinese White Pear (Pyrus bretschneideri): A Critical Role in Leaf Regeneration.

GRAS is a transcription regulator factor, which plays an important role in plant growth and development. Previous analyses found that several GRAS functions have been identified, such as axillary bud meristem formation, radial root elongation, gibberellin signaling, light signaling, and abiotic stress. The GRAS family has been comprehensively evaluated in several species. However, little finding is on the GRAS transcription factors (TFs) in Chinese white pear. In this study, 99 PbGRAS were systemically characterized and renamed PbGRAS1 to PbGRAS99 according to their chromosomal localizations. Phylogenetic analysis and structural features revealed that could be classified into eight subfamilies (LISCL, Ls, SHR, HAM, SCL, PAT, SCR, and DELLA). Further analysis of introns/exons and conserved motifs revealed that they are diverse and functionally differentiated in number and structure. Synteny analysis among Pyrus bretschenedri, Prunus mume, Prunus avium, Fragaria vesca, and Prunus persica showed that GRAS duplicated regions were more conserved. Dispersed duplication events are the most common mechanism and may play a crucial role in the expansion of the GRAS gene family. In addition, cis-acting elements of the PbGRAS gene were found in promoter regions associated with hormone and environmental stress responses. Notably, the expression pattern detected by qRT-PCR indicated that PbGRAS genes were differentially expressed under gibberellin (GA), abscisic acid (ABA), and auxin (IAA) conditions, which are responsive to abiotic stress. PbGRAS89 and PbGRAS99 were highly expressed at different stages of hormone treatment and may play important role in leaf development. Therefore, we selected PbGRAS89 and PbGRAS99 to clone and construct pCAMBIA1301-PbGRAS89, 99 and transferred them into Arabidopsis thaliana. Finally, we observed and compared the changes of overexpressed plants and wild-type plants during regeneration. This method was used to analyze their roles in leaf regeneration of Chinese white pear. In addition, we also constructed pCAMBIA1305-PbGRAS89, 99, and transferred them into onion cells to determine the subcellular localization. Subcellular localization experiments showed that PbGRAS89 and PbGRAS99 were localized in the nucleus. In summary, the results of this study indicate that PbGRAS89 and PbGRAS99 are mainly responsible for leaf regeneration of Chinese white pear, which plays a positive role in callus formation and provides rich resources for studying GRAS gene functions.

Characterization of the mitochondrial GlutamyltRNAGln amidotransferase (GatCAB) as a new model for mitochondrial translation disorders.

Mitochondrial OXPHOS diseases are provoked by dysfunction of the OXPHOS system, showing heterogeneity from a genetic, biochemical and clinical perspective. This complexity is partly due to the involvement of proteins encoded in two genomes located in different cellular compartments, nuclear and mitochondrial DNA (mtDNA). mtDNA encodes part of the mtRNAs translation machinery, including 22 tRNAs. All these are charged by their corresponding aminoacyl-tRNA synthetases (ARS2), except mt-Gln-tRNAGln whose synthesis is carried out through an indirect pathway. The mt-tRNAGln is charged with glutamic acid (Glu) by a non-discriminating mitochondrial glutamyl-tRNA synthetase (EARS2), then the Glu is converted to Gln by a Glutamyl-tRNAGln amidotransferase, using free glutamine as an amide donor, yielding Gln- tRNAGln. This mitochondrial amidotransferase activity lies in the GatCAB complex, which is formed by three subunits: GatA (QSRL1), GatB (GATB), and GatC (GATC). Patients with mutations in these genes presented with severe cardiomyopathy and lactic acidosis, which underscores the importance of the GatCAB complex as an essential component in the translational machinery of mitochondrial protein synthesis. We have generated Knockout (KO) lines for all subunits using the genomic editing system CRISPR/Cas9 in HEK293T cells. We measured the levels of the GatCAB subunits and mitochondrial proteins by western blot. We studied the oxygen consumption levels using a Clark electrode and detected the GatCAB complex by Blue Native-western blotting. In all cases, KO cells show reduced levels of mitochondrial proteins, decreased oxygen consumption and an absence of the GatCAB complex compared to wildtype cells. In the QRSL1 and GATC KO cells the levels of all subunits diminish, but in GATB KO cells the levels of subunits A and C remain stable. Possibly, QRSL1 and GATC could form a stable but nonfunctional dimer, to which GATB would join. GATB has two potential ATG codons at the start of its sequence separated by 12 nucleotides. From over 90 clones analyzed, we only obtained one KO, by generating an upstream ORF, and two cell lines with very low expression of GATB. The latter present a duplicated region at the beginning of GATB which changes the reading frame from the first ATG; suggesting the downstream ATG as an alternative translation start codon for those ribosome small subunits that do not assemble on the first one, yielding a GATB protein missing 4 amino acids. Characterization of the molecular mechanisms that lead to the synthesis of mt-Gln-tRNAGln and its role in the physiology of the cell could allow us to comprehend the pathological manifestation of the defects in these genes and propose possible therapeutic avenues.

Repeated translocation of a supergene underlying rapid sex chromosome turnover in Takifugu pufferfish

Recent studies have revealed a surprising diversity of sex chromosomes in vertebrates. However, the detailed mechanism of their turnover is still elusive. To understand this process, it is necessary to compare closely related species in terms of sex-determining genes and the chromosomes harboring them. Here, we explored the genus Takifugu, in which one strong candidate sex-determining gene, Amhr2, has been identified. To trace the processes involved in transitions in the sex-determination system in this genus, we studied 12 species and found that while the Amhr2 locus likely determines sex in the majority of Takifugu species, three species have acquired sex-determining loci at different chromosomal locations. Nevertheless, the generation of genome assemblies for the three species revealed that they share a portion of the male-specific supergene that contains a candidate sex-determining gene, GsdfY, along with genes that potentially play a role in male fitness. The shared supergene spans ∼100 kb and is flanked by two duplicated regions characterized by CACTA transposable elements. These results suggest that the shared supergene has taken over the role of sex-determining locus from Amhr2 in lineages leading to the three species, and repeated translocations of the supergene underlie the turnover of sex chromosomes in these lineages. These findings highlight the underestimated role of a mobile supergene in the turnover of sex chromosomes in vertebrates.

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale.

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNAIle, ND1, tRNALeu2, 16S, tRNAVal, 12S, tRNAPhe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

Identification and expression analysis of capsaicin biosynthesis pathway genes at genome level in Capsicum chinense

Capsicum chinense is one of the best ingredients for making chili sauces and extracting capsaicin. In this study, we identified 59 genes associated with shikimic acid, phenylpropane, unsaturated fatty acids, phenylalanine, valine, leucine, isoleucine and other metabolic pathways based on pepper genome. In the present study, 58 of the 59 genes involved in capsaicin synthesis were found to be unevenly distributed on 11 chromosomes, and several genes with the same number of exons and high homology were found to be clustered in a duplication region of the chromosome. The promoter regions of these genes were predicted and analyzed, and multiple gene promoters were found to contain MYB binding sites, which provided a reference for the subsequent exploration of MYB target genes. At the same time, the expression levels of genes related to capsaicin synthesis were analyzed based on different transcriptomes, and it was found that these genes showed significant differences in different plant development stages and different tissues. The results of this study will be helpful to further study the synthesis mechanism of capsaicin in Capsicum chinense and the functional relationship between genes.

Array Comparative Genomic Hybridisation and Droplet Digital PCR Uncover Recurrent Copy Number Variation of the TTN Segmental Duplication Region.

Intragenic segmental duplication regions are potential hotspots for recurrent copy number variation and possible pathogenic aberrations. Two large sarcomeric genes, nebulin and titin, both contain such segmental duplication regions. Using our custom Comparative Genomic Hybridisation array, we have previously shown that a gain or loss of more than one copy of the repeated block of the nebulin triplicate region constitutes a recessive pathogenic mutation. Using targeted array-CGH, similar copy number variants can be detected in the segmental duplication region of titin. Due to the limitations of the array-CGH methodology and the repetitiveness of the region, the exact copy numbers of the blocks could not be determined. Therefore, we developed complementary custom Droplet Digital PCR assays for the titin segmental duplication region to confirm true variation. Our combined methods show that the titin segmental duplication region is subject to recurrent copy number variation. Gains and losses were detected in samples from healthy individuals as well as in samples from patients with different muscle disorders. The copy number variation observed in our cohort is likely benign, but pathogenic copy number variants in the segmental duplication region of titin cannot be excluded. Further investigations are needed, however, this region should no longer be neglected in genetic analyses.

Polymorphism of Sex Determination Amongst Wild Populations Suggests its Rapid Turnover Within the Nile Tilapia Species.

Sex-determining regions have been identified in the Nile tilapia on linkage groups (LG) 1, 20 and 23, depending on the domesticated strains used. Sex determining studies on wild populations of this species are scarce. Previous work on two wild populations, from Lake Volta (Ghana) and from Lake Koka (Ethiopia), found the sex-determining region on LG23. These populations have a Y-specific tandem duplication containing two copies of the Anti-Müllerian Hormone amh gene (named amhY and amhΔY). Here, we performed a whole-genome short-reads analysis using male and female pools on a third wild population from Lake Hora (Ethiopia). We found no association of sex with LG23, and no duplication of the amh gene. Furthermore, we found no evidence of sex linkage on LG1 or on any other LGs. Long read whole genome sequencing of a male from each population confirmed the absence of a duplicated region on LG23 in the Lake Hora male. In contrast, long reads established the structure of the Y haplotype in Koka and Kpandu males and the order of the genes in the duplicated region. Phylogenies constructed on the nuclear and mitochondrial genomes, showed a closer relationship between the two Ethiopian populations compared to the Ghanaian population, implying an absence of the LG23Y sex-determination region in Lake Hora males. Our study supports the hypothesis that the amh region is not the sex-determining region in Hora males. The absence of the Y amh duplication in the Lake Hora population reflects a rapid change in sex determination within Nile tilapia populations. The genetic basis of sex determination in the Lake Hora population remains unknown.

A custom ddPCR method for the detection of copy number variations in the nebulin triplicate region.

The human genome contains repetitive regions, such as segmental duplications, known to be prone to copy number variation. Segmental duplications are highly identical and homologous sequences, posing a specific challenge for most mutation detection methods. The giant nebulin gene is expressed in skeletal muscle. It harbors a large segmental duplication region composed of eight exons repeated three times, the so-called triplicate region. Mutations in nebulin are known to cause nemaline myopathy and other congenital myopathies. Using our custom targeted Comparative Genomic Hybridization arrays, we have previously shown that copy number variations in the nebulin triplicate region are pathogenic when the copy number of the segmental duplication block deviates two or more copies from the normal number, which is three per allele. To complement our Comparative Genomic Hybridization arrays, we have established a custom Droplet Digital PCR method for the detection of copy number variations within the nebulin triplicate region. The custom Droplet Digital PCR assays allow sensitive, rapid, high-throughput, and cost-effective detection of copy number variations within this region and is ready for implementation a screening method for disease-causing copy number variations of the nebulin triplicate region. We suggest that Droplet Digital PCR may also be used in the study and diagnostics of other segmental duplication regions of the genome.

Gene arrangement, phylogeny and divergence time estimation of mitogenomes in Thrips.

The metazoan mitogenomes usually display conserved gene arrangement while thrips are known for their extensive gene rearrangement, and duplication of the control region. We sequenced complete mitogenomes of eight species of thrips to determine the gene arrangement, phylogeny and divergence time estimation. All contain 37 genes and one control region, (CR) except four species with two CRs. Duplicated tRNAs were detected in Mycterothrips nilgiriensis and Thrips florum. nad4-nad4L were not found adjacent to each other in Phibalothrips peringueyi and Plicothrips apicalis. Both Bayesian and likelihood phylogenetic analyses of thrips mitogenomes supported the monophyly of two suborders (Terebrantia and Tubulifera) and the two largest families (Phlaeothripidae and Thripidae). Out of seven earlier proposed ancestral gene blocks, six are conserved in Panchaetothripinae, three in Thripinae and two in Phlaeothripidae. Additionally, eight Thrips Gene Blocks were identified, of which, three conserved in Tubulifera, four in Terebrantia, and one only in Aeolothripidae. Forty-two gene boundaries (15 from previous study + 27 new) were identified. The molecular divergence time is estimated for the order Thysanoptera and suggested that these insects may have been diversified from hemipterans in the late Permian period. The most recent ancestors belong to family Thripidae and Phlaeothripidae, which were diversified in upper Cretaceous period and showed higher rates of rearrangement from the ancestral gene order. The current study is the first largest effort to provide the new insights into the mitogenomic features, gene arrangement, phylogeny and divergence time estimation of thrips belonging to the order Thysanoptera.

CRISPR/Cas9-based genome editing for the modification of multiple duplications that cause Duchenne muscular dystrophy.

With the development of basic research, some genetic-based methods have been found to treat Duchenne muscular dystrophy (DMD) with large deletion mutations and nonsense mutations. Appropriate therapeutic approaches for repairing multiple duplications are limited. We used the CRISPR (clustered regularly interspaced short palindromic repeat)/Cas9 system with patient-derived primary myoblasts to correct multiple duplications of the dystrophin gene. Muscle tissues from a patient carrying duplications of dystrophin were obtained, and tissue-derived primary cells were cultured. Myoblasts were purified with an immunomagnetic sorting system using CD56 microbeads. After transduction by lentivirus with a designed single guide RNA (sgRNA) targeting a duplicated region, myoblasts were allowed to differentiate for 7 days. Copy number variations in the exons of the patient's myotubes were quantified by real-time PCR before and after genetic editing. Western blot analysis was performed to detect the full-length dystrophin protein before and after genetic editing. The ten sequences predicted to be the most likely off-targets were determined by Sanger sequencing. The patient carried duplications of exon 18-25, dystrophin protein expression was completely abrogated. Real-time PCR showed that the copy number of exon 25 in the patient's myotubes was 2.015 ± 0.079 compared with that of the healthy controls. After editing, the copy number of exon 25 in the patient's modified myotubes was 1.308 ± 0.083 compared with that of the healthy controls (P < 0.001). Western blot analysis revealed no expression of the dystrophin protein in the patient's myotubes before editing. After editing, the patient's myotubes expressed the full-length dystrophin protein at a level that was ~6.12% of that in the healthy control samples. Off-target analysis revealed no abnormal editing at the ten sites predicted to be the most likely off-target sites. The excision of multiple duplications by the CRISPR/Cas9 system restored the expression of full-length dystrophin. This study provides proof of evidence for future genome-editing therapy in patients with DMD caused by multiple duplication mutations.