Search

Background: Powdery mildew is a devastating fungal disease that poses a significant threat to wheat yield and quality worldwide. Identifying resistance genes is highly advantageous for the molecular breeding of resistant cultivars. GRAS proteins are important transcription factors that regulate plant development and stress responses. Nonetheless, their roles in wheat-pathogen interactions remain poorly understood. Results: In this study, we used bioinformatics tools to identify and analyze wheat GRAS family genes responsive to biotic stresses and elucidated the function of TaeSCL6 within this family. A total of 179 GRAS genes in wheat were unevenly distributed on 7 chromosomes, and classified into 12 subfamilies based on phylogenetic relationship analysis. Gene duplication analysis revealed 13 pairs of tandem repeats and 142 pairs of segmental duplications, which may account for the rapid expansion of the wheat GRAS family. Expression pattern analysis revealed that 75% of the expressed TaeGRAS genes are responsive to biotic stresses. Few studies have focused on the roles of HAM subfamily genes. Consequently, we concentrated our analysis on the members of the HAM subfamily. Fourteen motifs were identified in the HAM family proteins from both Triticeae species and Arabidopsis, indicating that these motifs were highly conserved during evolution. Promoter analysis indicated that the promoters of HAM genes contain several cis-regulatory elements associated with hormone response, stress response, light response, and growth and development. Both qRT-PCR and RNA-seq data analyses demonstrated that TaeSCL6 responds to Blumeria graminis infection. Therefore, we investigated the role of TaeSCL6 in regulating wheat resistance via RNA interference and barley stripe mosaic virus induced gene silencing. Wheat plants with silenced TaeSCL6 exhibited increased susceptibility to powdery mildew. Conclusions: In summary, this study not only validates the positive role of TaeSCL6 in wheat resistance to powdery mildew, but also provides candidate gene resources for future breeding of disease-resistance wheat cultivars. [ABSTRACT FROM AUTHOR]

Background: Microsatellites are increasingly realized to have biological significance in human genome and health in past decades, the assembled complete reference sequence of human genome T2T-CHM13 brought great help for a comprehensive study of short tandem repeats in the human genome. Results: Microsatellites density landscapes of all 24 chromosomes were built here for the first complete reference sequence of human genome T2T-CHM13. These landscapes showed that short tandem repeats (STRs) are prone to aggregate characteristically to form a large number of STRs density peaks. We classified 8,823 High Microsatellites Density Peaks (HMDPs), 35,257 Middle Microsatellites Density Peaks (MMDPs) and 199, 649 Low Microsatellites Density Peaks (LMDPs) on the 24 chromosomes; and also classified the motif types of every microsatellites density peak. These STRs density aggregation peaks are mainly composing of a single motif, and AT is the most dominant motif, followed by AATGG and CCATT motifs. And 514 genomic regions were characterized by microsatellite density feature in the full T2T-CHM13 genome. Conclusions: These landscape maps exhibited that microsatellites aggregate in many genomic positions to form a large number of microsatellite density peaks with composing of mainly single motif type in the complete reference genome, indicating that the local microsatellites density varies enormously along the every chromosome of T2T-CHM13. [ABSTRACT FROM AUTHOR]

Background: The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology. Results: The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress. Conclusions: Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species. [ABSTRACT FROM AUTHOR]

Background: Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns. Results: A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans. Conclusions: Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida. [ABSTRACT FROM AUTHOR]

Background: 'Taishuu' has a crisp texture, abundant juice, and sweet flavor with hints of cantaloupe. The availability of mitochondrial genome data of Diospyros species is far from the known number of species. Results: The sequencing data were assembled into a closed circular mitochondrial chromosome with a 421,308 bp length and a 45.79% GC content. The mitochondrial genome comprised 40 protein-coding, 24 tRNA, and three rRNA genes. The most common codons for arginine (Arg), proline (Pro), glycine (Gly), tryptophan (Trp), valine (Val), alanine (Ala), and leucine (Leu) were AGA, CCA, GGA, UGG, GUA, GCA, and CUA, respectively. The start codon for cox1 and nad4L protein-coding genes was ACG (ATG), whereas the remaining protein-coding genes started with ATG. There are four types of stop codons: CGA, TAA, TAG, and TGA, with TAA being the most frequently used stop codon (45.24%). In the D. kaki Thunb. 'Taishuu' mitochondrial genome, a total of 645 repeat sequences were identified, including 125 SSRs, 7 tandem repeats, and 513 dispersed repeats. Collinearity analysis revealed a close relationship between D. kaki Thunb. 'Taishuu' and Diospyros oleifera, with conserved homologous gene fragments shared among these species in large regions of the mitochondrial genome. The protein-coding genes ccmB and nad4L were observed to undergo positive selection. Analysis of homologous sequences between chloroplasts and mitochondria identified 28 homologous segments, with a total length of 24,075 bp, accounting for 5.71% of the mitochondrial genome. These homologous segments contain 8 annotated genes, including 6 tRNA genes and 2 protein-coding genes (rrn18 and ccmC). There are 23 homologous genes between chloroplasts and nuclei. Mitochondria, chloroplasts, and nuclei share two homologous genes, which are trnV-GAC and trnW-CCA. Conclusion: In conclusion, a high-quality chromosome-level draft genome for D. kaki was generated in this study, which will contribute to further studies of major economic traits in the genus Diospyros. [ABSTRACT FROM AUTHOR]

Background: Calmodulins (CaMs)/CaM-like proteins (CMLs) are crucial Ca2+-binding sensors that can decode and transduce Ca2+ signals during plant development and in response to various stimuli. The CaM/CML gene family has been characterized in many plant species, but this family has not yet been characterized and analyzed in peanut, especially for its functions in response to Ralstonia solanacearum. In this study, we performed a genome-wide analysis to analyze the CaM/CML genes and their functions in resistance to R. solanacearum. Results: Here, 67, 72, and 214 CaM/CML genes were identified from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. The genes were divided into nine subgroups (Groups I-IX) with relatively conserved exon‒intron structures and motif compositions. Gene duplication, which included whole-genome duplication, tandem repeats, scattered repeats, and unconnected repeats, produced approximately 81 pairs of homologous genes in the AhCaM/CML gene family. Allopolyploidization was the main reason for the greater number of AhCaM/CML members. The nonsynonymous (Ka) versus synonymous (Ks) substitution rates (less than 1.0) suggested that all homologous pairs underwent intensive purifying selection pressure during evolution. AhCML69 was constitutively expressed in different tissues of peanut plants and was involved in the response to R. solanacearum infection. The AhCML69 protein was localized in the cytoplasm and nucleus. Transient overexpression of AhCML69 in tobacco leaves increased resistance to R. solanacearum infection and induced the expression of defense-related genes, suggesting that AhCML69 is a positive regulator of disease resistance. Conclusions: This study provides the first comprehensive analysis of the AhCaM/CML gene family and potential genetic resources for the molecular design and breeding of peanut bacterial wilt resistance. [ABSTRACT FROM AUTHOR]

Background: Venoms have evolved independently over a hundred times in the animal kingdom to deter predators and/or subdue prey. Venoms are cocktails of various secreted toxins, whose origin and diversification provide an appealing system for evolutionary researchers. Previous studies of the ant venom of Tetramorium bicarinatum revealed several Myrmicitoxin (MYRTX) peptides that gathered into seven precursor families suggesting different evolutionary origins. Analysis of the T. bicarinatum genome enabling further genomic approaches was necessary to understand the processes underlying the evolution of these myrmicitoxins. Results: Here, we sequenced the genome of Tetramorium bicarinatum and reported the organisation of 44 venom peptide genes (vpg). Of the eleven chromosomes that make up the genome of T. bicarinatum, four carry the vpg which are organized in tandem repeats. This organisation together with the ML evolutionary analysis of vpg sequences, is consistent with evolution by local duplication of ancestral genes for each precursor family. The structure of the vpg into two or three exons is conserved after duplication events while the promoter regions are the least conserved parts of the vpg even for genes with highly identical sequences. This suggests that enhancer sequences were not involved in duplication events, but were recruited from surrounding regions. Expression level analysis revealed that most vpg are highly expressed in venom glands, although one gene or group of genes is much more highly expressed in each family. Finally, the examination of the genomic data revealed that several genes encoding transcription factors (TFs) are highly expressed in the venom glands. The search for binding sites (BS) of these TFs in the vpg promoters revealed hot spots of GATA sites in several vpg families. Conclusion: In this pioneering investigation on ant venom genes, we provide a high-quality assembly genome and the annotation of venom peptide genes that we think can fosters further genomic research to understand the evolutionary history of ant venom biochemistry. [ABSTRACT FROM AUTHOR]

Background: Cytochrome P450 monooxygenases (CYP450s) play a crucial role in various biochemical reactions involved in the synthesis of antioxidants, pigments, structural polymers, and defense-related compounds in plants. As sweet potato (Ipomoea batatas L.) holds significant economic importance, a comprehensive analysis of CYP450 genes in this plant species can offer valuable insights into the evolutionary relationships and functional characteristics of these genes. Results: In this study, we successfully identified and categorized 95 CYP450 genes from the sweet potato genome into 5 families and 31 subfamilies. The predicted subcellular localization results indicate that CYP450s are distributed in the cell membrane system. The promoter region of the IbCYP450 genes contains various cis-acting elements related to plant hormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) have been identified in the IbCYP450 family proteins, with 5 genes lacking introns and only one exon. We observed extensive duplication events within the CYP450 gene family, which may account for its expansion. The gene duplication analysis results showed the presence of 15 pairs of genes with tandem repeats. Interaction network analysis reveals that IbCYP450 families can interact with multiple target genes and there are protein-protein interactions within the family. Transcription factor interaction analysis suggests that IbCYP450 families interact with multiple transcription factors. Furthermore, gene expression analysis revealed tissue-specific expression patterns of CYP450 genes in sweet potatoes, as well as their response to abiotic stress and plant hormones. Notably, quantitative real-time polymerase chain reaction (qRT‒PCR) analysis indicated the involvement of CYP450 genes in the defense response against nonbiological stresses in sweet potatoes. Conclusions: These findings provide a foundation for further investigations aiming to elucidate the biological functions of CYP450 genes in sweet potatoes. [ABSTRACT FROM AUTHOR]

Background: γ-glutamylcyclotransferase (GGCT), an enzyme to maintain glutathione homeostasis, plays a vital role in the response to plant growth and development as well as the adaptation to various stresses. Although the GGCT gene family analysis has been conducted in Arabidopsis and rice, the family genes have not yet been well identified and analyzed at the genome-wide level in wheat (Triticum aestivum L.). Results: In the present study, 20 TaGGCT genes were identified in the wheat genome and widely distributed on chromosomes 2A, 2B, 2D, 3A, 4A, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B, and 7D. Phylogenetic and structural analyses showed that these TaGGCT genes could be classified into three subfamilies: ChaC, GGGACT, and GGCT-PS. They exhibited similar motif compositions and distribution patterns in the same subgroup. Gene duplication analysis suggested that the expansion of TaGGCT family genes was facilitated by segmental duplications and tandem repeats in the wheat evolutionary events. Identification of diverse cis-acting response elements in TaGGCT promoters indicated their potential fundamental roles in response to plant development and abiotic stresses. The analysis of transcriptome data combined with RT-qPCR results revealed that the TaGGCTs genes exhibited ubiquitous expression across plant organs, with highly expressed in roots, stems, and developing grains. Most TaGGCT genes were up-regulated after 6 h under 20% PEG6000 and ABA treatments. Association analysis revealed that two haplotypes of TaGGCT20 gene displayed significantly different Thousand-kernel weight (TKW), Kernel length (KL), and Kernel width (KW) in wheat. The geographical and annual distribution of the two haplotypes of TaGGCT20 gene further revealed that the frequency of the favorable haplotype TaGGCT20-Hap-I was positively selected in the historical breeding process of wheat. Conclusion: This study investigated the genome-wide identification, structure, evolution, and expression analysis of TaGGCT genes in wheat. The motifs of TaGGCTs were highly conserved throughout the evolutionary history of wheat. Most TaGGCT genes were highly expressed in roots, stems, and developing grains, and involved in the response to drought stresses. Two haplotypes were developed in the TaGGCT20 gene, where TaGGCT20-Hap-I, as a favorable haplotype, was significantly associated with higher TKW, KL, and KW in wheat, suggesting that the haplotype is used as a function marker for the selection in grain yield in wheat breeding. [ABSTRACT FROM AUTHOR]

Background: As a population genetic tool, mitochondrial DNA is commonly divided into the ~ 1-kb control region (CR), in which single nucleotide variant (SNV) diversity is relatively high, and the coding region, in which selective constraint is greater and diversity lower, but which provides an informative phylogeny. In some species, the CR contains variable tandemly repeated sequences that are understudied due to heteroplasmy. Domestic cats (Felis catus) have a recent origin and therefore traditional CR-based analysis of populations yields only a small number of haplotypes. Results: To increase resolution we used Nanopore sequencing to analyse 119 cat mitogenomes via a long-amplicon approach. This greatly improves discrimination (from 15 to 87 distinct haplotypes in our dataset) and defines a phylogeny showing similar starlike topologies within all major clades (haplogroups), likely reflecting post-domestication expansion. We sequenced RS2, a CR tandem array of 80-bp repeat units, placing RS2 array structures within the phylogeny and increasing overall haplotype diversity. Repeat number varies between 3 and 12 (median: 4) with over 30 different repeat unit types differing largely by SNVs. Five SNVs show evidence of independent recurrence within the phylogeny, and seven are involved in at least 11 instances of rapid spread along repeat arrays within haplogroups. Conclusions: In defining mitogenome variation our study provides key information for the forensic genetic analysis of cat hair evidence, and for the first time a phylogenetically informed picture of tandem repeat variation that reveals remarkably dynamic mutation processes at work in the mitochondrion. [ABSTRACT FROM AUTHOR]

Background: The 1RS arm of wheat-rye 1BL.1RS translocations contains several subtelomeric tandem repeat families. To study the effect of the difference in the composition of these tandem repeats on the meiotic recombination of 1RS arms can help to enrich the genetic diversity of 1BL.1RS translocation chromosomes. Results: Five wheat-rye 1BL.1RS translocation cultivars/lines were used to build two cross combinations including group 1 (20T401 × Zhou 8425B, 20T401 × Lovrin 10 and 20T401 × Chuannong 17) and group 2 (20T360-2 × Zhou 8425B, 20T360-2 × Lovrin 10 and 20T360-2 × Chuannong 17). Oligonucleotide (oligo) probes Oligo-s120.3, Oligo-TR72, and Oligo-119.2-2 produced the same signal pattern on the 1RS arms in lines 20T401 and 20T360-2, and another signal pattern in the three cultivars Zhou 8425B, Lovrin 10 and Chuannong 17. The Oligo-pSc200 signal disappeared from the 1RS arms of the line 20T401, and the signal intensity of this probe on the 1RS arms of the line 20T360-2 was weaker than that of the three cultivars. The five cultivars/lines had the same signal pattern of the probe Oligo-pSc250. The recombination rate of 1RS arms in group 1 was significantly lower than that in group 2. In the progenies from group 1, unequal meiotic recombination in the subtelomeric pSc119.2 and pSc250 tandem repeat regions, and a 1BL.1RS with inversion of 1RS segment between the pSc200 and the nucleolar organizer region were found. Conclusions: This study provides a visual tool to detect the meiotic recombination of 1RS arms. The meiotic recombination rate of 1RS arms was affected by the variation of pSc200 tandem repeat, indicating the similar composition of subtelomeric tandem repeats on these arms could increase their recombination rate. These results indicate that the 1RS subtelomeric structure will affect its recombination, and thus the localization of genes on 1RS by means of meiotic recombination might also be affected. [ABSTRACT FROM AUTHOR]

Self-incompatibility (SI) is a reproductive protection mechanism that plants acquired during evolution to prevent self-recession. As the female determinant of SI specificity, SRK has been shown to be the only recognized gene on the stigma and plays important roles in SI response. Asteraceae is the largest family of dicotyledonous plants, many of which exhibit self-incompatibility. However, systematic studies on SRK gene family in Asteraceae are still limited due to lack of high-quality genomic data. In this study, we performed the first systematic genome-wide identification of S-locus receptor like kinases (SRLKs) in the self-incompatible Asteraceae species, Erigeron breviscapus, which is also a widely used perennial medicinal plant endemic to China.52 SRLK genes were identified in the E. breviscapus genome. Structural analysis revealed that the EbSRLK proteins in E. breviscapus are conserved. SRLK proteins from E. breviscapus and other SI plants are clustered into 7 clades, and the majority of the EbSRLK proteins are distributed in Clade I. Chromosomal and duplication analyses indicate that 65% of the EbSRLK genes belong to tandem repeats and could be divided into six tandem gene clusters. Gene expression patterns obtained in E. breviscapus multiple-tissue RNA-Seq data revealed differential temporal and spatial features of EbSRLK genes. Among these, two EbSRLK genes having high expression levels in tongue flowers were cloned. Subcellular localization assay demonstrated that both of their fused proteins are localized on the plasma membrane. All these results indicated that EbSRLK genes possibly involved in SI response in E. breviscapus. This comprehensive genome-wide study of the SRLK gene family in E. breviscapus provides valuable information for understanding the mechanism of SSI in Asteraceae. [ABSTRACT FROM AUTHOR]

Abstract Background Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study. Results A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation. Conclusions The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation.

Background: Methyl-binding domain (MBD) is a class of methyl-CpG-binding domain proteins that affects the regulation of gene expression through epigenetic modifications. MBD genes are not only inseparable from DNA methylation but have also been identified and validated in various plants. Although MBD is involved in a group of physiological processes and stress regulation in these plants, MBD genes in Eleutherococcus senticosus remain largely unknown. Results: Twenty EsMBD genes were identified in E. senticosus. Among the 24 chromosomes of E. senticosus, EsMBD genes were unevenly distributed on 12 chromosomes, and only one tandem repeat gene existed. Collinearity analysis showed that the fragment duplication was the main motif for EsMBD gene expansion. As the species of Araliaceae evolved, MBD genes also evolved and gradually exhibited different functional differentiation. Furthermore, cis-acting element analysis showed that there were numerous cis-acting elements in the EsMBD promoter region, among which light response elements and anaerobic induction elements were dominant. The expression motif analysis revealed that 60% of the EsMBDs were up-regulated in the 30% water content group. Conclusions: By comparing the transcriptome data of different saponin contents of E. senticosus and integrating them with the outcomes of molecular docking analysis, we hypothesized that EsMBD2 and EsMBD5 jointly affect the secondary metabolic processes of E. senticosus saponins by binding to methylated CpG under conditions of drought stress. The results of this study laid the foundation for subsequent research on the E. senticosus and MBD genes. [ABSTRACT FROM AUTHOR]

Background: Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) is a frequently used typing method for identifying the Beijing genotype of Mycobacterium tuberculosis (Mtb), which is easily transformed into rifampicin (RIF) resistance. The RIF resistance of Mtb is considered to be highly related with the mutation of rpoB gene. Therefore, this study aimed to analyze the relationship between the repetitive number of MIRU loci and the mutation of rpoB gene. Methods: An open-source whole-genome sequencing data of Mtb was used to detect the mutation of rpoB gene and the repetitive number of MIRU loci by bioinformatics methods. Cochran-Armitage analysis was performed to analyze the trend of the rpoB gene mutation rate and the repetitive number of MIRU loci. Results: Among 357 rifampicin-resistant tuberculosis (RR-TB), 304 strains with mutated rpoB genes were detected, and 6 of 67 rifampicin susceptible strains were detected mutations. The rpoB gene mutational rate showed an upward trend with the increase of MIRU10, MIRU39, QUB4156 and MIRU16 repetitive number, but only the repetitive number of MIRU10, MRIU39 and QUB4156 were risk factors for rpoB gene mutation. The Hunter-Gaston discriminatory index (HGDI) of MIRU10 (0.65) and QUB4156 (0.62) was high in the overall sample, while MIRU39 (0.39) and MIRU16 (0.43) showed a moderate discriminatory Power. Conclusion: The mutation rate of rpoB gene increases with the addition of repetitive numbers of MIRU10, QUB4156 and MIRU39 loci. [ABSTRACT FROM AUTHOR]

Background: The subgenus Gynopodium belonging to genus Magnolia have high ornamental, economic, and ecological value. Subgenus Gynopodium contains eight species, but six of these species are threatened. No studies to date have characterized the characteristics of the chloroplast genomes (CPGs) within subgenus Gynopodium species. In this study, we compared the structure of CPGs, identified the mutational hotspots and resolved the phylogenetic relationship of subgenus Gynopodium. Results: The CPGs of six subgenus Gynopodium species ranged in size from 160,027 bp to 160,114 bp. A total of 131 genes were identified, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. We detected neither major expansions or contractions in the inverted repeat region, nor rearrangements or insertions in the CPGs of six subgenus Gynopodium species. A total of 300 large repeat sequences (forward, reverse, and palindrome repeats), 847 simple sequence repeats, and five highly variable regions were identified. One gene (ycf1) and four intergenic regions (psbA-trnH-GUG, petA-psbJ, rpl32-trnL-UAG, and ccsA-ndhD) were identified as mutational hotspots by their high nucleotide diversity (Pi) values (≥ 0.004), which were useful for species discrimination. Maximum likelihood and Bayesian inference trees were concordant and indicated that Magnoliaceae consisted of two genera Liriodendron and Magnolia. Six species of subgenus Gynopodium clustered as a monophyletic clade, forming a sister clade with subgenus Yulania (BS = 100%, PP = 1.00). Due to the non-monophyly of subgenus Magnolia, subgenus Gynopodium should be treated as a section of Magnolia. Within section Gynopodium, M. sinica diverged first (posterior probability = 1, bootstrap = 100), followed by M. nitida, M. kachirachirai and M. lotungensis. M. omeiensis was sister to M. yunnanensis (posterior probability = 0.97, bootstrap = 50). Conclusion: The CPGs and characteristics information provided by our study could be useful in species identification, conservation genetics and resolving phylogenetic relationships of Magnoliaceae species. [ABSTRACT FROM AUTHOR]

Background: Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study. Results: A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation. Conclusions: The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation. [ABSTRACT FROM AUTHOR]

Background: Ribosomally-synthesized cyclic peptides are widely found in plants and exhibit useful bioactivities for humans. The identification of cyclic peptide sequences and their precursor proteins is facilitated by the growing number of sequenced genomes. While previous research largely focused on the chemical diversity of these peptides across various species, there is little attention to a broader range of potential peptides that are not chemically identified. Results: A pioneering study was initiated to explore the genetic diversity of linusorbs, a group of cyclic peptides uniquely occurring in cultivated flax (Linum usitatissimum). Phylogenetic analysis clustered the 5 known linusorb precursor proteins into two clades and one singleton. Preliminary tBLASTn search of the published flax genome using the whole protein sequence as query could only retrieve its homologues within the same clade. This limitation was overcome using a profile-based mining strategy. After genome reannotation, a hidden Markov Model (HMM)-based approach identified 58 repeats homologous to the linusorb-embedded repeats in 8 novel proteins, implying that they share common ancestry with the linusorb-embedded repeats. Subsequently, we developed a customized profile composed of a random linusorb-like domain (LLD) flanked by 5 conserved sites and used it for string search of the proteome, which extracted 281 LLD-containing repeats (LLDRs) in 25 proteins. Comparative analysis of different repeat categories suggested that the 5 conserved flanking sites among the non-homologous repeats have undergone convergent evolution driven by functional selection. Conclusions: The profile-based mining approach is suitable for analyzing repetitive sequences. The 25 LLDR proteins identified herein represent the potential diversity of cyclic peptides within the flax genome and lay a foundation for further studies on the functions and evolution of these protein tandem repeats. [ABSTRACT FROM AUTHOR]

Background: Structural variants (SVs) play a crucial role in gene regulation, trait association, and disease in humans. SV genotyping has been extensively applied in genomics research and clinical diagnosis. Although a growing number of SV genotyping methods for long reads have been developed, a comprehensive performance assessment of these methods has yet to be done. Results: Based on one simulated and three real SV datasets, we performed an in-depth evaluation of five SV genotyping methods, including cuteSV, LRcaller, Sniffles, SVJedi, and VaPoR. The results show that for insertions and deletions, cuteSV and LRcaller have similar F1 scores (cuteSV, insertions: 0.69–0.90, deletions: 0.77–0.90 and LRcaller, insertions: 0.67–0.87, deletions: 0.74–0.91) and are superior to other methods. For duplications, inversions, and translocations, LRcaller yields the most accurate genotyping results (0.84, 0.68, and 0.47, respectively). When genotyping SVs located in tandem repeat region or with imprecise breakpoints, cuteSV (insertions and deletions) and LRcaller (duplications, inversions, and translocations) are better than other methods. In addition, we observed a decrease in F1 scores when the SV size increased. Finally, our analyses suggest that the F1 scores of these methods reach the point of diminishing returns at 20× depth of coverage. Conclusions: We present an in-depth benchmark study of long-read SV genotyping methods. Our results highlight the advantages and disadvantages of each genotyping method, which provide practical guidance for optimal application selection and prospective directions for tool improvement. [ABSTRACT FROM AUTHOR]

Background: Ziziphus hajarensis is an endemic plant species well-distributed in the Western Hajar mountains of Oman. Despite its potential medicinal uses, little is known regarding its genomic architecture, phylogenetic position, or evolution. Here we sequenced and analyzed the entire chloroplast (cp) genome of Z. hajarensis to understand its genetic organization, structure, and phylogenomic disposition among Rhamnaceae species. Results: The results revealed the genome of Z. hajarensis cp comprised 162,162 bp and exhibited a typical quadripartite structure, with a large single copy (LSC) region of 895,67 bp, a small single copy (SSC) region of 19,597 bp and an inverted repeat (IR) regions of 26,499 bp. In addition, the cp genome of Z. hajarensis comprises 126 genes, including 82 protein-coding genes, eight rRNA genes, and 36 tRNA genes. Furthermore, the analysis revealed 208 microsatellites, 96.6% of which were mononucleotides. Similarly, a total of 140 repeats were identified, including 11 palindromic, 24 forward, 14 reverse, and 104 tandem repeats. The whole cp genome comparison of Z. hajarensis and nine other species from family Rhamnaceae showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. Comparative phylogenetic analysis based on the complete cp genome, 66 shared genes and matK gene revealed that Z. hajarensis shares a clade with Z. jujuba and that the family Rhamnaceae is the closest family to Barbeyaceae and Elaeagnaceae. Conclusion: All the genome features such as genome size, GC content, genome organization and gene order were highly conserved compared to the other related genomes. The whole cp genome of Z. hajarensis gives fascinating insights and valuable data that may be used to identify related species and reconstruct the phylogeny of the species. [ABSTRACT FROM AUTHOR]

Background: The yield of many crop plants can be substantially reduced by plant-pathogenic Xanthomonas bacteria. The infection strategy of many Xanthomonas strains is based on transcription activator-like effectors (TALEs), which are secreted into the host cells and act as transcriptional activators of plant genes that are beneficial for the bacteria.The modular DNA binding domain of TALEs contains tandem repeats, each comprising two hyper-variable amino acids. These repeat-variable diresidues (RVDs) bind to their target box and determine the specificity of a TALE.All available tools for the prediction of TALE targets within the host plant suffer from many false positives. In this paper we propose a strategy to improve prediction accuracy by considering the epigenetic state of the host plant genome in the region of the target box. Results: To this end, we extend our previously published tool PrediTALE by considering two epigenetic features: (i) chromatin accessibility of potentially bound regions and (ii) DNA methylation of cytosines within target boxes. Here, we determine the epigenetic features from publicly available DNase-seq, ATAC-seq, and WGBS data in rice.We benchmark the utility of both epigenetic features separately and in combination, deriving ground-truth from RNA-seq data of infections studies in rice. We find an improvement for each individual epigenetic feature, but especially the combination of both.Having established an advantage in TALE target predicting considering epigenetic features, we use these data for promoterome and genome-wide scans by our new tool EpiTALE, leading to several novel putative virulence targets. Conclusions: Our results suggest that it would be worthwhile to collect condition-specific chromatin accessibility data and methylation information when studying putative virulence targets of Xanthomonas TALEs. [ABSTRACT FROM AUTHOR]

Background: F-box proteins represent a diverse class of adaptor proteins of the ubiquitin-proteasome system (UPS) that play critical roles in the cell cycle, signal transduction, and immune response by removing or modifying cellular regulators. Among closely related organisms of the Caenorhabditis genus, remarkable divergence in F-box gene copy numbers was caused by sizeable species-specific expansion and contraction. Although F-box gene number expansion plays a vital role in shaping genomic diversity, little is known about molecular evolutionary mechanisms responsible for substantial differences in gene number of F-box genes and their functional diversification in Caenorhabditis. Here, we performed a comprehensive evolution and underlying mechanism analysis of F-box genes in five species of Caenorhabditis genus, including C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei. Results: Herein, we identified and characterized 594, 192, 377, 39, 1426 F-box homologs encoding putative F-box proteins in the genome of C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. Our work suggested that extensive species-specific tandem duplication followed by a small amount of gene loss was the primary mechanism responsible for F-box gene number divergence in Caenorhabditis genus. After F-box gene duplication events occurred, multiple mechanisms have contributed to gene structure divergence, including exon/intron gain/loss, exonization/pseudoexonization, exon/intron boundaries alteration, exon splits, and intron elongation by tandem repeats. Based on high-throughput RNA sequencing data analysis, we proposed that F-box gene functions have diversified by sub-functionalization through highly divergent stage-specific expression patterns in Caenorhabditis species. Conclusions: Massive species-specific tandem duplications and occasional gene loss drove the rapid evolution of the F-box gene family in Caenorhabditis, leading to complex gene structural variation and diversified functions affecting growth and development within and among Caenorhabditis species. In summary, our findings outline the evolution of F-box genes in the Caenorhabditis genome and lay the foundation for future functional studies. [ABSTRACT FROM AUTHOR]

Background: Members of the basic helix-loop-helix (bHLH) transcription factor family perform indispensable functions in various biological processes, such as plant growth, seed maturation, and abiotic stress responses. However, the bHLH family in foxtail millet (Setaria italica), an important food and feed crop, has not been thoroughly studied. Results: In this study, 187 bHLH genes of foxtail millet (SibHLHs) were identified and renamed according to the chromosomal distribution of the SibHLH genes. Based on the number of conserved domains and gene structure, the SibHLH genes were divided into 21 subfamilies and two orphan genes via phylogenetic tree analysis. According to the phylogenetic tree, the subfamilies 15 and 18 may have experienced stronger expansion in the process of evolution. Then, the motif compositions, gene structures, chromosomal spread, and gene duplication events were discussed in detail. A total of sixteen tandem repeat events and thirty-eight pairs of segment duplications were identified in bHLH family of foxtail millet. To further investigate the evolutionary relationship in the SibHLH family, we constructed the comparative syntenic maps of foxtail millet associated with representative monocotyledons and dicotyledons species. Finally, the gene expression response characteristics of 15 typical SibHLH genes in different tissues and fruit development stages, and eight different abiotic stresses were analysed. The results showed that there were significant differences in the transcription levels of some SibHLH members in different tissues and fruit development stages, and different abiotic stresses, implying that SibHLH members might have different physiological functions. Conclusions: In this study, we identified 187 SibHLH genes in foxtail millet and further analysed the evolution and expression patterns of the encoded proteins. The findings provide a comprehensive understanding of the bHLH family in foxtail millet, which will inform further studies on the functional characteristics of SibHLH genes. [ABSTRACT FROM AUTHOR]

Background: Plant non-specific lipid transfer proteins (nsLTPs), a group of small, basic ubiquitous proteins to participate in lipid transfer, cuticle formation and stress response, are involved in the regulation of plant growth and development. To date, although the nsLTP gene family of barley (Hordeum vulgare L.) has been preliminarily identified, it is still unclear in the recently completed genome database of barley and Qingke, and its transcriptional profiling under abiotic stress has not been elucidated as well. Results: We identified 40 barley nsLTP (HvLTP) genes through a strict screening strategy based on the latest barley genome and 35 Qingke nsLTP (HtLTP) orthologues using blastp, and these LTP genes were divided into four types (1, 2, D and G). At the same time, a comprehensive analysis of the physical and chemical characteristics, homology alignment, conserved motifs, gene structure and evolution of HvLTPs and HtLTPs further supported their similar nsLTP characteristics and classification. The genomic location of HvLTPs and HtLTPs showed that these genes were unevenly distributed, and obvious HvLTP and HtLTP gene clusters were found on the 7 chromosomes including six pairs of tandem repeats and one pair of segment repeats in the barley genome, indicating that these genes may be co-evolutionary and co-regulated. A spatial expression analysis showed that most HvLTPs and HtLTPs had different tissue-specific expression patterns. Moreover, the upstream cis-element analysis of HvLTPs and HtLTPs showed that there were many different stress-related transcriptional regulatory elements, and the expression pattern of HvLTPs and HtLTPs under abiotic stress also indicated that numerous HvLTP and HtLTP genes were related to the abiotic stress response. Taken together, these results may be due to the differences in promoters rather than by genes themselves resulting in different expression patterns under abiotic stress. Conclusion: Due to a stringent screening and comprehensive analysis of the nsLTP gene family in barley and Qingke and its expression profile under abiotic stress, this study can be considered a useful source for the future studies of nsLTP genes in either barley or Qingke or for comparisons of different plant species. [ABSTRACT FROM AUTHOR]

Background: GRAS, an important family of transcription factors, have played pivotal roles in regulating numerous intriguing biological processes in plant development and abiotic stress responses. Since the sequencing of the sorghum genome, a plethora of genetic studies were mainly focused on the genomic information. The indepth identification or genome-wide analysis of GRAS family genes, especially in Sorghum bicolor, have rarely been studied. Results: A total of 81 SbGRAS genes were identified based on the S. bicolor genome. They were named SbGRAS01 to SbGRAS81 and grouped into 13 subfamilies (LISCL, DLT, OS19, SCL4/7, PAT1, SHR, SCL3, HAM-1, SCR, DELLA, HAM-2, LAS and OS4). SbGRAS genes are not evenly distributed on the chromosomes. According to the results of the gene and motif composition, SbGRAS members located in the same group contained analogous intron/exon and motif organizations. We found that the contribution of tandem repeats to the increase in sorghum GRAS members was slightly greater than that of fragment repeats. By quantitative (q) RT-PCR, the expression of 13 SbGRAS members in different plant tissues and in plants exposed to six abiotic stresses at the seedling stage were quantified. We further investigated the relationship between DELLA genes, GAs and grain development in S. bicolor. The paclobutrazol treatment significantly increased grain weight, and affected the expression levels of all DELLA subfamily genes. SbGRAS03 is the most sensitive to paclobutrazol treatment, but also has a high response to abiotic stresses. Conclusions: Collectively, SbGRAs play an important role in plant development and response to abiotic stress. This systematic analysis lays the foundation for further study of the functional characteristics of GRAS genes of S. bicolor. [ABSTRACT FROM AUTHOR]

Abstract Background Microsatellites, or Simple Sequence Repeats (SSRs), are short tandem repeats of 1–6 nt motifs present in all genomes. Emerging evidence points to their role in cellular processes and gene regulation. Despite the huge resource of genomic information currently available, SSRs have been studied in a limited context and compared across relatively few species. Results We have identified ~ 685 million eukaryotic microsatellites and analyzed their genomic trends across 15 taxonomic subgroups from protists to mammals. The distribution of SSRs reveals taxon-specific variations in their exonic, intronic and intergenic densities. Our analysis reveals the differences among non-related species and novel patterns uniquely demarcating closely related species. We document several repeats common across subgroups as well as rare SSRs that are excluded almost throughout evolution. We further identify species-specific signatures in pathogens like Leishmania as well as in cereal crops, Drosophila, birds and primates. We also find that distinct SSRs preferentially exist as long repeating units in different subgroups; most unicellular organisms show no length preference for any SSR class, while many SSR motifs accumulate as long repeats in complex organisms, especially in mammals. Conclusions We present a comprehensive analysis of SSRs across taxa at an unprecedented scale. Our analysis indicates that the SSR composition of organisms with heterogeneous cell types is highly constrained, while simpler organisms such as protists, green algae and fungi show greater diversity in motif abundance, density and GC content. The microsatellite dataset generated in this work provides a large number of candidates for functional analysis and for studying their roles across the evolutionary landscape.

Background: Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively. Results: To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements. Conclusions: Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads and the generality of these causes and factors should be tested further in other species. [ABSTRACT FROM AUTHOR]

Background: The genus Ehrlichia consists of tick-borne obligatory intracellular bacteria that can cause deadly diseases of medical and agricultural importance. Ehrlichia sp. HF, isolated from Ixodes ovatus ticks in Japan [also referred to as I. ovatus Ehrlichia (IOE) agent], causes acute fatal infection in laboratory mice that resembles acute fatal human monocytic ehrlichiosis caused by Ehrlichia chaffeensis. As there is no small laboratory animal model to study fatal human ehrlichiosis, Ehrlichia sp. HF provides a needed disease model. However, the inability to culture Ehrlichia sp. HF and the lack of genomic information have been a barrier to advance this animal model. In addition, Ehrlichia sp. HF has several designations in the literature as it lacks a taxonomically recognized name. Results: We stably cultured Ehrlichia sp. HF in canine histiocytic leukemia DH82 cells from the HF strain-infected mice, and determined its complete genome sequence. Ehrlichia sp. HF has a single double-stranded circular chromosome of 1,148,904 bp, which encodes 866 proteins with a similar metabolic potential as E. chaffeensis. Ehrlichia sp. HF encodes homologs of all virulence factors identified in E. chaffeensis, including 23 paralogs of P28/OMP-1 family outer membrane proteins, type IV secretion system apparatus and effector proteins, two-component systems, ankyrin-repeat proteins, and tandem repeat proteins. Ehrlichia sp. HF is a novel species in the genus Ehrlichia, as demonstrated through whole genome comparisons with six representative Ehrlichia species, subspecies, and strains, using average nucleotide identity, digital DNA-DNA hybridization, and core genome alignment sequence identity. Conclusions: The genome of Ehrlichia sp. HF encodes all known virulence factors found in E. chaffeensis, substantiating it as a model Ehrlichia species to study fatal human ehrlichiosis. Comparisons between Ehrlichia sp. HF and E. chaffeensis will enable identification of in vivo virulence factors that are related to host specificity, disease severity, and host inflammatory responses. We propose to name Ehrlichia sp. HF as Ehrlichia japonica sp. nov. (type strain HF), to denote the geographic region where this bacterium was initially isolated. [ABSTRACT FROM AUTHOR]

Background: Telomeric DNA is typically comprised of G-rich tandem repeat motifs and maintained by telomerase (Greider CW, Blackburn EH; Cell 51:887–898; 1987). In eukaryotes lacking telomerase, a variety of DNA repair and DNA recombination based pathways for telomere maintenance have evolved in organisms normally dependent upon telomerase for telomere elongation (Webb CJ, Wu Y, Zakian VA; Cold Spring Harb Perspect Biol 5:a012666; 2013); collectively called Alternative Lengthening of Telomeres (ALT) pathways. By measuring (TTAGGG) n tract lengths from the same large DNA molecules that were optically mapped, we simultaneously analyzed telomere length dynamics and subtelomere-linked structural changes at a large number of specific subtelomeric loci in the ALT-positive cell lines U2OS, SK-MEL-2 and Saos-2. Results: Our results revealed loci-specific ALT telomere features. For example, while each subtelomere included examples of single molecules with terminal (TTAGGG) n tracts as well as examples of recombinant telomeric single molecules, the ratio of these molecules was subtelomere-specific, ranging from 33:1 (19p) to 1:25 (19q) in U2OS. The Saos-2 cell line shows a similar percentage of recombinant telomeres. The frequency of recombinant subtelomeres of SK-MEL-2 (11%) is about half that of U2OS and Saos-2 (24 and 19% respectively). Terminal (TTAGGG) n tract lengths and heterogeneity levels, the frequencies of telomere signal-free ends, and the frequency and size of retained internal telomere-like sequences (ITSs) at recombinant telomere fusion junctions all varied according to the specific subtelomere involved in a particular cell line. Very large linear extrachromosomal telomere repeat (ECTR) DNA molecules were found in all three cell lines; these are in principle capable of templating synthesis of new long telomere tracts via break-induced repair (BIR) long-tract DNA synthesis mechanisms and contributing to the very long telomere tract length and heterogeneity characteristic of ALT cells. Many of longest telomere tracts (both end-telomeres and linear ECTRs) displayed punctate CRISPR/Cas9-dependent (TTAGGG) n labeling patterns indicative of interspersion of stretches of non-canonical telomere repeats. Conclusion: Identifying individual subtelomeres and characterizing linked telomere (TTAGGG) n tract lengths and structural changes using our new single-molecule methodologies reveals the structural consequences of telomere damage, repair and recombination mechanisms in human ALT cells in unprecedented molecular detail and significant differences in different ALT-positive cell lines. [ABSTRACT FROM AUTHOR]

Background: Chloroplast (cp) genome information would facilitate the development and utilization of Taxodium resources. However, cp genome characteristics of Taxodium were poorly understood. Results: We determined the complete cp genome sequences of T. distichum, T. mucronatum, and T. ascendens. The cp genomes are 131,947 bp to 132,613 bp in length, encode 120 genes with the same order, and lack typical inverted repeat (IR) regions. The longest small IR, a 282 bp trnQ-containing IR, were involved in the formation of isomers. Comparative analysis of the 3 cp genomes showed that 91.57% of the indels resulted in the periodic variation of tandem repeat (TR) motifs and 72.46% single nucleotide polymorphisms (SNPs) located closely to TRs, suggesting a relationship between TRs and mutational dynamics. Eleven hypervariable regions were identified as candidates for DNA barcode development. Hypothetical cp open reading frame 1(Ycf1) was the only one gene that has an indel in coding DNA sequence, and the indel is composed of a long TR. When extended to cupressophytes, ycf1 genes have undergone a universal insertion of TRs accompanied by extreme length expansion. Meanwhile, ycf1 also located in rearrangement endpoints of cupressophyte cp genomes. All these characteristics highlight the important role of repeats in the evolution of cp genomes. Conclusions: This study added new evidence for the role of repeats in the dynamics mechanism of cp genome mutation and rearrangement. Moreover, the information of TRs and hypervariable regions would provide reliable molecular resources for future research focusing on the infrageneric taxa identification, phylogenetic resolution, population structure and biodiversity for the genus Taxodium and Cupressophytes. [ABSTRACT FROM AUTHOR]

Abstract Background Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. Methods We developed a hybrid assembly pipeline called “Alpaca” that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. Results Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. Conclusion Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

Background: The genus Hosta is a group of economically appreciated perennial herbs consisting of approximately 25 species that is endemic to eastern Asia. Due to considerable morphological variability, the genus has been well recognized as a group with taxonomic problems. Chloroplast is a cytoplasmic organelle with its own genome, which is the most commonly used for phylogenetic and genetic diversity analyses for land plants. To understand the genomic architecture of Hosta chloroplasts and examine the level of nucleotide and size variation, we newly sequenced four (H. clausa, H. jonesii, H. minor, and H. venusta) and analyzed six Hosta species (including the four, H. capitata and H. yingeri) distributed throughout South Korea. Results: The average size of complete chloroplast genomes for the Hosta taxa was 156,642 bp with a maximum size difference of ~ 300 bp. The overall gene content and organization across the six Hosta were nearly identical with a few exceptions. There was a single tRNA gene deletion in H. jonesii and four genes were pseudogenized in three taxa (H. capitata, H. minor, and H. jonesii). We did not find major structural variation, but there were a minor expansion and contractions in IR region for three species (H. capitata, H. minor, and H. venusta). Sequence variations were higher in non-coding regions than in coding regions. Four genic and intergenic regions including two coding genes (psbA and ndhD) exhibited the largest sequence divergence showing potential as phylogenetic markers. We found compositional codon usage bias toward A/T at the third position. The Hosta plastomes had a comparable number of dispersed and tandem repeats (simple sequence repeats) to the ones identified in other angiosperm taxa. The phylogeny of 20 Agavoideae (Asparagaceae) taxa including the six Hosta species inferred from complete plastome data showed well resolved monophyletic clades for closely related taxa with high node supports. Conclusions: Our study provides detailed information on the chloroplast genome of the Hosta taxa. We identified nucleotide diversity hotspots and characterized types of repeats, which can be used for developing molecular markers applicable in various research area. [ABSTRACT FROM AUTHOR]

Background: Increased availability of genome assemblies for non-model organisms has resulted in invaluable biological and genomic insight into numerous vertebrates, including teleosts. Sequencing of the Atlantic cod (Gadus morhua) genome and the genomes of many of its relatives (Gadiformes) demonstrated a shared loss of the major histocompatibility complex (MHC) II genes 100 million years ago. An improved version of the Atlantic cod genome assembly shows an extreme density of tandem repeats compared to other vertebrate genome assemblies. Highly contiguous assemblies are therefore needed to further investigate the unusual immune system of the Gadiformes, and whether the high density of tandem repeats found in Atlantic cod is a shared trait in this group. Results: Here, we have sequenced and assembled the genome of haddock (Melanogrammus aeglefinus) – a relative of Atlantic cod – using a combination of PacBio and Illumina reads. Comparative analyses reveal that the haddock genome contains an even higher density of tandem repeats outside and within protein coding sequences than Atlantic cod. Further, both species show an elevated number of tandem repeats in genes mainly involved in signal transduction compared to other teleosts. A characterization of the immune gene repertoire demonstrates a substantial expansion of MCHI in Atlantic cod compared to haddock. In contrast, the Toll-like receptors show a similar pattern of gene losses and expansions. For the NOD-like receptors (NLRs), another gene family associated with the innate immune system, we find a large expansion common to all teleosts, with possible lineage-specific expansions in zebrafish, stickleback and the codfishes. Conclusions: The generation of a highly contiguous genome assembly of haddock revealed that the high density of short tandem repeats as well as expanded immune gene families is not unique to Atlantic cod – but possibly a feature common to all, or most, codfishes. A shared expansion of NLR genes in teleosts suggests that the NLRs have a more substantial role in the innate immunity of teleosts than other vertebrates. Moreover, we find that high copy number genes combined with variable genome assembly qualities may impede complete characterization of these genes, i.e. the number of NLRs in different teleost species might be underestimates. [ABSTRACT FROM AUTHOR]

Background: Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. Methods: We developed a hybrid assembly pipeline called "Alpaca" that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. Results: Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. Conclusion: Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations. [ABSTRACT FROM AUTHOR]

Background: To date, mitochondrial genomes of more than one hundred flatworms (Platyhelminthes) have been sequenced. They show a high degree of similarity and a strong taxonomic bias towards parasitic lineages. The mitochondrial gene atp8 has not been confidently annotated in any flatworm sequenced to date. However, sampling of free-living flatworm lineages is incomplete. We addressed this by sequencing the mitochondrial genomes of the two small-bodied (about 1 mm in length) free-living flatworms Stenostomum sthenum and Macrostomum lignano as the first representatives of the earliest branching flatworm taxa Catenulida and Macrostomorpha respectively. Results: We have used high-throughput DNA and RNA sequence data and PCR to establish the mitochondrial genome sequences and gene orders of S. sthenum and M. lignano. The mitochondrial genome of S. sthenum is 16,944 bp long and includes a 1,884 bp long inverted repeat region containing the complete sequences of nad3, rrnS, and nine tRNA genes. The model flatworm M. lignano has the smallest known mitochondrial genome among freeliving flatworms, with a length of 14,193 bp. The mitochondrial genome of M. lignano lacks duplicated genes, however, tandem repeats were detected in a non-coding region. Mitochondrial gene order is poorly conserved in flatworms, only a single pair of adjacent ribosomal or protein-coding genes - nad4l-nad4 - was found in S. sthenum and M. lignano that also occurs in other published flatworm mitochondrial genomes. Unexpectedly, we unambiguously identified the full metazoan mitochondrial protein-coding gene complement including atp8 in S. sthenum and M. lignano. A subsequent search detected atp8 in all mitochondrial genomes of polyclad flatworms published to date, although the gene wasn't previously annotated in these species. Conclusions: Manual, but not automated genome annotation revealed the presence of atp8 in basally branching freeliving flatworms, signifying both the importance of manual data curation and of diverse taxon sampling. We conclude that the loss of atp8 within flatworms is restricted to the parasitic taxon Neodermata. [ABSTRACT FROM AUTHOR]

Background: The first Atlantic cod (Gadus morhua) genome assembly published in 2011 was one of the early genome assemblies exclusively based on high-throughput 454 pyrosequencing. Since then, rapid advances in sequencing technologies have led to a multitude of assemblies generated for complex genomes, although many of these are of a fragmented nature with a significant fraction of bases in gaps. The development of long-read sequencing and improved software now enable the generation of more contiguous genome assemblies. Results: By combining data from Illumina, 454 and the longer PacBio sequencing technologies, as well as integrating the results of multiple assembly programs, we have created a substantially improved version of the Atlantic cod genome assembly. The sequence contiguity of this assembly is increased fifty-fold and the proportion of gap-bases has been reduced fifteen-fold. Compared to other vertebrates, the assembly contains an unusual high density of tandem repeats (TRs). Indeed, retrospective analyses reveal that gaps in the first genome assembly were largely associated with these TRs. We show that 21% of the TRs across the assembly, 19% in the promoter regions and 12% in the coding sequences are heterozygous in the sequenced individual. Conclusions: The inclusion of PacBio reads combined with the use of multiple assembly programs drastically improved the Atlantic cod genome assembly by successfully resolving long TRs. The high frequency of heterozygous TRs within or in the vicinity of genes in the genome indicate a considerable standing genomic variation in Atlantic cod populations, which is likely of evolutionary importance. [ABSTRACT FROM AUTHOR]

Background: The program RepeatMasker and the database Repbase-ISB are part of the most widely used strategy for annotating repeats in animal genomes. They have been used to show that avian genomes have a lower repeat content (8-12 %) than the sequenced genomes of many vertebrate species (30-55 %). However, the efficiency of such a library-based strategies is dependent on the quality and completeness of the sequences in the database that is used. An alternative to these library based methods are methods that identify repeats de novo. These alternative methods have existed for a least a decade and may be more powerful than the library based methods. We have used an annotation strategy involving several complementary de novo tools to determine the repeat content of the model genome galGal4 (1.04 Gbp), including identifying simple sequence repeats (SSRs), tandem repeats and transposable elements (TEs). Results: We annotated over one Gbp. of the galGal4 genome and showed that it is composed of approximately 19 % SSRs and TEs repeats. Furthermore, we estimate that the actual genome of the red jungle fowl contains about 31-35 % repeats. We find that library-based methods tend to overestimate TE diversity. These results have a major impact on the current understanding of repeats distributions throughout chromosomes in the red jungle fowl. Conclusions: Our results are a proof of concept of the reliability of using de novo tools to annotate repeats in large animal genomes. They have also revealed issues that will need to be resolved in order to develop gold-standard methodologies for annotating repeats in eukaryote genomes. [ABSTRACT FROM AUTHOR]

Background: A prominent and distinctive feature of the rye (Secale cereale) chromosomes is the presence of massive blocks of subtelomeric heterochromatin, the size of which is correlated with the copy number of tandem arrays. The rapidity with which these regions have formed over the period of speciation remains unexplained. Results: Using a BAC library created from the short arm telosome of rye chromosome 1R we uncovered numerous arrays of the pSc200 and pSc250 tandem repeat families which are concentrated in subtelomeric heterochromatin and identified the adjacent DNA sequences. The arrays show significant heterogeneity in monomer organization. 454 reads were used to gain a representation of the expansion of these tandem repeats across the whole rye genome. The presence of multiple, relatively short monomer arrays, coupled with the mainly star-like topology of the monomer phylogenetic trees, was taken as indicative of a rapid expansion of the pSc200 and pSc250 arrays. The evolution of subtelomeric heterochromatin appears to have included a significant contribution of illegitimate recombination. The composition of transposable elements (TEs) within the regions flanking the pSc200 and pSc250 arrays differed markedly from that in the genome a whole. Solo-LTRs were strongly enriched, suggestive of a history of active ectopic exchange. Several DNA motifs were over-represented within the LTR sequences. Conclusion: The large blocks of subtelomeric heterochromatin have arisen from the combined activity of TEs and the expansion of the tandem repeats. The expansion was likely based on a highly complex network of recombination mechanisms. [ABSTRACT FROM AUTHOR]

Background: Colorectal cancer is a major contributor to cancer morbidity and mortality. Tandem repeat instability and its effect on cancer phenotypes remain so far poorly studied on a genome-wide scale. Results: Here we analyze the genomes of 35 colorectal tumors and their matched normal (healthy) tissues for two types of tandem repeat instability, de-novo repeat gain or loss and repeat copy number variation. Specifically, we study for the first time genome-wide repeat instability in the promoters and exons of 18,439 genes, and examine the association of repeat instability with genome-scale gene expression levels. We find that tumors with a microsatellite instable (MSI) phenotype are enriched in genes with repeat instability, and that tumor genomes have significantly more genes with repeat instability compared to healthy tissues. Genes in tumor genomes with repeat instability in their promoters are significantly less expressed and show slightly higher levels of methylation. Genes in well-studied cancer-associated signaling pathways also contain significantly more unstable repeats in tumor genomes. Genes with such unstable repeats in the tumor-suppressor p53 pathway have lower expression levels, whereas genes with repeat instability in the MAPK and Wnt signaling pathways are expressed at higher levels, consistent with the oncogenic role they play in cancer. Conclusions: Our results suggest that repeat instability in gene promoters and associated differential gene expression may play an important role in colorectal tumors, which is a first step towards the development of more effective molecular diagnostic approaches centered on repeat instability. [ABSTRACT FROM AUTHOR]

Background: VNTR (Variable Number of Tandem Repeats) composite retrotransposons - SVA (SINE-R-VNTR-Alu), LAVA (LINE-1-Alu-VNTR-Alu), PVA (PTGR2-VNTR-Alu) and FVA (FRAM-VNTR-Alu) - are specific to hominoid primates. Their assembly, the evolution of their 5' and 3' domains, and the functional significance of the shared 5' Alu-like region are well understood. The central VNTR domain, by contrast, has long been assumed to represent a more or less random collection of 30-50 bp GC-rich repeats. It is only recently that it attracted attention in the context of regulation of SVA expression. Results: Here we provide evidence that the organization of the VNTR is non-random, with conserved repeat unit (RU) arrays at both the 5' and 3' ends of the VNTRs of human, chimpanzee and orangutan SVA and gibbon LAVA. The younger SVA subfamilies harbour highly organized internal RU arrays. The composition of these arrays is specific to the human/chimpanzee and orangutan lineages, respectively. Tracing the development of the VNTR through evolution we show for the first time how tandem repeats evolve within the constraints set by a functional, non-autonomous non-LTR retrotransposon in two different families - LAVA and SVA - in different hominoid lineages. Our analysis revealed that a microhomology-driven mechanism mediates expansion/contraction of the VNTR domain at the DNA level. Elements of all four VNTR composite families have been shown to be mobilized by the autonomous LINE1 retrotransposon in trans. In case of SVA, key determinants of mobilization are found in the 5' hexameric repeat/ Alu-like region. We now demonstrate that in LAVA, by contrast, the VNTR domain determines mobilization efficiency in the context of domain swaps between active and inactive elements. Conclusions: The central domain of VNTR composites evolves in a lineage-specific manner which gives rise to distinct structures in gibbon LAVA, orangutan SVA, and human/chimpanzee SVA. The differences observed between the families and lineages are likely to have an influence on the expression and mobilization of the elements. [ABSTRACT FROM AUTHOR]

Background: Mycobacterium avium subsp. paratuberculosis (MAP), the causative bacterium of Johne's disease in dairy cattle, is widespread in the Canadian dairy industry and has significant economic and animal welfare implications. An understanding of the population dynamics of MAP can be used to identify introduction events, improve control efforts and target transmission pathways, although this requires an adequate understanding of MAP diversity and distribution between herds and across the country. Whole genome sequencing (WGS) offers a detailed assessment of the SNP-level diversity and genetic relationship of isolates, whereas several molecular typing techniques used to investigate the molecular epidemiology of MAP, such as variable number of tandem repeat (VNTR) typing, target relatively unstable repetitive elements in the genome that may be too unpredictable to draw accurate conclusions. The objective of this study was to evaluate the diversity of bovine MAP isolates in Canadian dairy herds using WGS and then determine if VNTR typing can distinguish truly related and unrelated isolates. Results: Phylogenetic analysis based on 3,039 SNPs identified through WGS of 124 MAP isolates identified eight genetically distinct subtypes in dairy herds from seven Canadian provinces, with the dominant type including over 80% of MAP isolates. VNTR typing of 527 MAP isolates identified 12 types, including "bison type" isolates, from seven different herds. At a national level, MAP isolates differed from each other by 1-2 to 239-240 SNPs, regardless of whether they belonged to the same or different VNTR types. A herd-level analysis of MAP isolates demonstrated that VNTR typing may both over-estimate and under-estimate the relatedness of MAP isolates found within a single herd. Conclusions: The presence of multiple MAP subtypes in Canada suggests multiple introductions into the country including what has now become one dominant type, an important finding for Johne's disease control. VNTR typing often failed to identify closely and distantly related isolates, limiting the applicability of using this typing scheme to study the molecular epidemiology of MAP at a national and herd-level. [ABSTRACT FROM AUTHOR]

Background By reshuffling genomes, structural genomic reorganizations provide genetic variation on which natural selection can work. Understanding the mechanisms underlying this process has been a long-standing question in evolutionary biology. In this context, our purpose in this study is to characterize the genomic regions involved in structural rearrangements between human and macaque genomes and determine their influence on meiotic recombination as a way to explore the adaptive role of genome shuffling in mammalian evolution. Results We first constructed a highly refined map of the structural rearrangements and evolutionary breakpoint regions in the human and rhesus macaque genomes based on orthologous genes and whole-genome sequence alignments. Using two different algorithms, we refined the genomic position of known rearrangements previously reported by cytogenetic approaches and described new putative micro-rearrangements (inversions and indels) in both genomes. A detailed analysis of the rhesus macaque genome showed that evolutionary breakpoints are in gene-rich regions, being enriched in GO terms related to immune system. We also identified defense-response genes within a chromosome inversion fixed in the macaque lineage, underlying the relevance of structural genomic changes in evolutionary and/or adaptation processes. Moreover, by combining in silico and experimental approaches, we studied the recombination pattern of specific chromosomes that have suffered rearrangements between human and macaque lineages. Conclusions Our data suggest that adaptive alleles - in this case, genes involved in the immune response - might have been favored by genome rearrangements in the macaque lineage. [ABSTRACT FROM AUTHOR]

Background Tandem repeats (TRs) are unstable regions commonly found within genomes that have consequences for evolution and disease. In humans, polymorphic TRs are known to cause neurodegenerative and neuromuscular disorders as well as being associated with complex diseases such as diabetes and cancer. If present in upstream regulatory regions, TRs can modify chromatin structure and affect transcription; resulting in altered gene expression and protein abundance. The most common TRs are short tandem repeats (STRs), or microsatellites. Promoter located STRs are considerably more polymorphic than coding region STRs. As such, they may be a common driver of phenotypic variation. To study STRs located in regulatory regions, we have performed genome-wide analysis to identify all STRs present in a region that is 2 kilobases upstream and 1 kilobase downstream of the transcription start sites of genes. Results The Short Tandem Repeats in Regulatory Regions Table, STaRRRT, contains the results of the genome-wide analysis, outlining the characteristics of 5,264 STRs present in the upstream regulatory region of 4,441 human genes. Gene set enrichment analysis has revealed significant enrichment for STRs in cellular, transcriptional and neurological system gene promoters and genes important in ion and calcium homeostasis. The set of enriched terms has broad similarity to that seen in coding regions, suggesting that regulatory region STRs are subject to similar evolutionary pressures as STRs in coding regions and may, like coding region STRs, have an important role in controlling gene expression. Conclusions STaRRRT is a readily-searchable resource for investigating potentially polymorphic STRs that could influence the expression of any gene of interest. The processes and genes enriched for regulatory region STRs provide potential novel targets for diagnosing and treating disease, and support a role for these STRs in the evolution of the human genome. [ABSTRACT FROM AUTHOR]

Background: Microsatellite DNA is one of many powerful genetic markers used for the construction of genetic linkage maps and the study of population genetics. The biological databases in public domain hold vast numbers of microsatellite sequences for many organisms including fishes. The microsatellite data available in these data sources were extracted and managed into a database that facilitates sequences analysis and browsing relevant information. The system also helps to design primer sequences for flanking regions of repeat loci for PCR identification of polymorphism within populations. Description: FishMicrosat is a database of microsatellite sequences of fishes and shellfishes that includes important aquaculture species such as Lates calcarifer, Ctenopharyngodon idella, Hypophthalmichthys molitrix, Penaeus monodon, Labeo rohita, Oreochromis niloticus, Fenneropenaeus indicus and Macrobrachium rosenbergii. The database contains 4398 microsatellite sequences of 41 species belonging to 15 families from the Indian subcontinent. GenBank of NCBI was used as a prime data source for developing the database. The database presents information about simple and compound microsatellites, their clusters and locus orientation within sequences. The database has been integrated with different tools in a web interface such as primer designing, locus finding, mapping repeats, detecting similarities among sequences across species, and searching using motifs and keywords. In addition, the database has the ability to browse information on the top 10 families and the top 10 species, through record overview. Conclusions: FishMicrosat database is a useful resource for fish and shellfish microsatellite analyses and locus identification across species, which has important applications in population genetics, evolutionary studies and genetic relatedness among species. The database can be expanded further to include the microsatellite data of fishes and shellfishes from other regions and available information on genome sequencing project of species of aquaculture importance. [ABSTRACT FROM AUTHOR]

Background: Perennial ryegrass (Lolium perenne L.) is one of the most important forage and turf grass species of temperate regions worldwide. Its mitochondrial genome is inherited maternally and contains genes that can influence traits of agricultural importance. Moreover, the DNA sequence of mitochondrial genomes has been established and compared for a large number of species in order to characterize evolutionary relationships. Therefore, it is crucial to understand the organization of the mitochondrial genome and how it varies between and within species. Here, we report the first de novo assembly and annotation of the complete mitochondrial genome from perennial ryegrass. Results: Intact mitochondria from perennial ryegrass leaves were isolated and used for mtDNA extraction. The mitochondrial genome was sequenced to a 167-fold coverage using the Roche 454 GS-FLX Titanium platform, and assembled into a circular master molecule of 678,580 bp. A total of 34 proteins, 14 tRNAs and 3 rRNAs are encoded by the mitochondrial genome, giving a total gene space of 48,723 bp (7.2%). Moreover, we identified 149 open reading frames larger than 300 bp and covering 67,410 bp (9.93%), 250 SSRs, 29 tandem repeats, 5 pairs of large repeats, and 96 pairs of short inverted repeats. The genes encoding subunits of the respiratory complexes - nad1 to nad9, cob, cox1 to cox3 and atp1 to atp9 - all showed high expression levels both in absolute numbers and after normalization. Conclusions: The circular master molecule of the mitochondrial genome from perennial ryegrass presented here constitutes an important tool for future attempts to compare mitochondrial genomes within and between grass species. Our results also demonstrate that mitochondria of perennial ryegrass contain genes crucial for energy production that are well conserved in the mitochondrial genome of monocotyledonous species. The expression analysis gave us first insights into the transcriptome of these mitochondrial genes in perennial ryegrass. [ABSTRACT FROM AUTHOR]

Background: Tandem repeats are ubiquitous and abundant in higher eukaryotic genomes and constitute, along with transposable elements, much of DNA underlying centromeres and other heterochromatic domains. In maize, centromeric satellite repeat (CentC) and centromeric retrotransposons (CR), a class of Ty3/gypsy retrotransposons, are enriched at centromeres. Some satellite repeats have homology to retrotransposons and several mechanisms have been proposed to explain the expansion, contraction as well as homogenization of tandem repeats. However, the origin and evolution of tandem repeat loci remain largely unknown. Results: CRM1TR and CRM4TR are novel tandem repeats that we show to be entirely derived from CR elements belonging to two different subfamilies, CRM1 and CRM4. Although these tandem repeats clearly originated in at least two separate events, they are derived from similar regions of their respective parent element, namely the long terminal repeat (LTR) and untranslated region (UTR). The 5' ends of the monomer repeat units of CRM1TR and CRM4TR map to different locations within their respective LTRs, while their 3' ends map to the same relative position within a conserved region of their UTRs. Based on the insertion times of heterologous retrotransposons that have inserted into these tandem repeats, amplification of the repeats is estimated to have begun at least ~4 (CRM1TR) and ~1 (CRM4TR) million years ago. Distinct CRM1TR sequence variants occupy the two CRM1TR loci, indicating that there is little or no movement of repeats between loci, even though they are separated by only ~1.4 Mb. Conclusions: The discovery of two novel retrotransposon derived tandem repeats supports the conclusions from earlier studies that retrotransposons can give rise to tandem repeats in eukaryotic genomes. Analysis of monomers from two different CRM1TR loci shows that gene conversion is the major cause of sequence variation. We propose that successive intrastrand deletions generated the initial repeat structure, and gene conversions increased the size of each tandem repeat locus. [ABSTRACT FROM AUTHOR]

Background: Macrosatellite repeats (MSRs), usually spanning hundreds of kilobases of genomic DNA, comprise a significant proportion of the human genome. Because of their highly polymorphic nature, MSRs represent an extreme example of copy number variation, but their structure and function is largely understudied. Here, we describe a detailed study of six autosomal and two X chromosomal MSRs among 270 HapMap individuals from Central Europe, Asia and Africa. Copy number variation, stability and genetic heterogeneity of the autosomal macrosatellite repeats RS447 (chromosome 4p), MSR5p (5p), FLJ40296 (13q), RNU2 (17q) and D4Z4 (4q and 10q) and X chromosomal DXZ4 and CT47 were investigated. Results: Repeat array size distribution analysis shows that all of these MSRs are highly polymorphic with the most genetic variation among Africans and the least among Asians. A mitotic mutation rate of 0.4-2.2% was observed, exceeding meiotic mutation rates and possibly explaining the large size variability found for these MSRs. By means of a novel Bayesian approach, statistical support for a distinct multimodal rather than a uniform allele size distribution was detected in seven out of eight MSRs, with evidence for equidistant intervals between the modes. Conclusions: The multimodal distributions with evidence for equidistant intervals, in combination with the observation of MSR-specific constraints on minimum array size, suggest that MSRs are limited in their configurations and that deviations thereof may cause disease, as is the case for facioscapulohumeral muscular dystrophy. However, at present we cannot exclude that there are mechanistic constraints for MSRs that are not directly disease-related. This study represents the first comprehensive study of MSRs in different human populations by applying novel statistical methods and identifies commonalities and differences in their organization and function in the human genome. [ABSTRACT FROM AUTHOR]

Background: Tandem repeats (TRs) in the mitochondrial (mt) genome control region have been documented in a wide variety of vertebrate species. The mechanism by which repeated tracts originate and undergo duplication and deletion, however, remains unclear. Results: We analyzed DNA sequences of mt genome TRs (mtTRs) in the ridged-eye flounder (Pleuronichthys cornutus), and characterized DNA sequences of mtTRs from other vertebrates using the data available in GenBank. Tandem repeats are concentrated in the control regions; however, we found approximately 16.6% of the TRs elsewhere in the mt genome. The flounder mtTRs possess three motif types with hypervariable characteristics at the 30 end of the control region (CR). Conclusion: Based on our analysis of this larger dataset of mtTR sequences, we propose a novel model of Pause Melting Misalignment (PMM) to describe the birth and motif indel of tandem repeats. PMM is activated during a pause event in mitochondrial replication in which a dynamic competition between the nascent (N) heavy strand and the displaced (D) heavy strand may lead to the melting of the N-strand from the template (T) light strand. When mispairing occurs during rebinding of the N-strand, one or several motifs can be inserted or deleted in both strands during the next round of mt-replication or repair. This model can explain the characteristics of TRs in available vertebrate mt genomes. [ABSTRACT FROM AUTHOR]