Search

About half of the human genome is covered by repetitive sequences. These sequences have a large share in the other mammalian genomes, therefore studying this part of the genome can provide researchers valuable information on evolution. The aim of this study was to sequencing and assembly the whole genome of Iranian Bactrian camels to identify transposable elements and their distribution in the genome of this species. In addition, the results of Iranian Bactrian camels were compared with non-Iranian Bactrian camels and dromedary camels. Materials and methods In this study, the whole genome of six Iranian Bactrian camels was sequenced to transposable elements identification. Iranian Bactrian camel whole genome sequenced using Illumina HiSeq 2000 system in paired-end. FastQC and Trimmomatic software were used to quality control and quality filtering of raw sequencing reads, respectively. CLC Genomics Workbench (CLC Bio, Aarhus, Denmark) was used to de novo assembly of trimmed reads. Also, we used the RepeatMasker program to search for transposable elements using a homology-based method. Results Results of the assembling of sequenced genomes showed that the genome size in these samples ranged from 1.9 to 1.97 Gb. In the present study, the percentage of transposable elements for six Iranian Bactrian camels was 29.89% on average of the whole genome. The percentage of LINE sequences for the Iranian Bactrian camel was 17.58% on average. So, these sequences were considered as the largest group of transposable elements in the Bactrian camel in this study. SINE elements showed a lower number in comparison with LINEs. So that, only 3.45% of the total Bactria camel genome length was dedicated to the SINEs. In accordance with the results of Iranian dromedary camels, no Alu element was identified in the genome of Iranian Bactrian camels. Conclusion Shortage of genomic and biological information about camels is one of the inhibiting factors in advancing the breeding goals and programs. Although this study is not enough alone, it can be a step towards starting the production of genomic data for camels. Continuing this kind of study and integrating biological and genomic information will provide the ground for the start of modern breeding in Iranian camels.

Among the three genomes in plant cells, the mitochondrial genome (mitogenome) is the least studied due to complex recombination and intergenomic transfer. In gymnosperms only ∼20 mitogenomes have been released thus far, which hinders a systematic investigation into the tempo and mode of mitochondrial DNA evolution in seed plants. Here, we report the complete mitogenome sequence of Platycladus orientalis (Cupressaceae). This mitogenome is assembled as two circular-mapping chromosomes with a size of ∼2.6 Mb and which contains 32 protein-coding genes, three rRNA and seven tRNA genes, and 1,068 RNA editing sites. Repetitive sequences, including dispersed repeats, transposable elements (TEs), and tandem repeats, made up 23% of the genome. Comparative analyses with 17 other mitogenomes representing the five gymnosperm lineages revealed a 30-fold difference in genome size, 80-fold in repetitive content, and 230-fold in substitution rate. We found dispersed repeats are highly associated with mitogenome expansion (r = 0.99), and most of them were accumulated during recent duplication events. Syntenic blocks and shared sequences between mitogenomes decay rapidly with divergence time (r = 0.53), with the exceptions of Ginkgo and Cycads which retained conserved genome structure over long evolutionary time. Our phylogenetic analysis supports a sister group relationship of Cupressophytes and Gnetophytes; both groups are unique in that they lost 8–12 protein-coding genes, of which 4–7 intact genes are likely transferred to nucleus. These two clades also show accelerated and highly variable substitution rates relative to other gymnosperms. Our study highlights the dynamic and enigmatic evolution of gymnosperm mitogenomes.

MotivationRepeat proteins, which contain multiple repeats of short sequence motifs, form a large but seldom-studied group of proteins. Methods focusing on the analysis of 3D structures of such proteins identified many subtle effects in length distribution of individual motifs that are important for their functions. However, similar analysis was yet not applied to the vast majority of repeat proteins with unknown 3D structures, mostly because of the extreme diversity of the underlying motifs and the resulting difficulty to detect those.ResultsWe developed FAIT, a sequence-based algorithm for the precise assignment of individual repeats in repeat proteins and introduced a framework to classify and compare aperiodicity patterns for large protein families. FAIT extracts repeat positions by post-processing FFAS alignment matrices with image processing methods. On examples of proteins with Leucine Rich Repeat (LRR) domains and other solenoids like proteins, we show that the automated analysis with FAIT correctly identifies exact lengths of individual repeats based entirely on sequence information.Availability and implementationhttps://github.com/GodzikLab/FAIT CONTACT: adam@godziklab.orgSupplementary informationSupplementary data are available at Bioinformatics online.

MYC is a key driver of cellular transformation and is deregulated in most human cancers. Studies of MYC and its interactors have provided mechanistic insight into its role as a regulator of gene transcription. MYC has been previously linked to chromatin regulation through its interaction with INI1 (SMARCB1/hSNF5/BAF47), a core member of the SWI/SNF chromatin remodeling complex. INI1 is a potent tumor suppressor that is inactivated in several types of cancers, most prominently as the hallmark alteration in pediatric malignant rhabdoid tumors. However, the molecular and functional interaction of MYC and INI1 remains unclear. Here, we characterize the MYC-INI1 interaction in mammalian cells, mapping their minimal binding domains to functionally significant regions of MYC (leucine zipper) and INI1 (repeat motifs), and demonstrating that the interaction does not interfere with MYC-MAX interaction. Protein-protein interaction network analysis expands the MYC-INI1 interaction to the SWI/SNF complex and a larger network of chromatin regulatory complexes. Genome-wide analysis reveals that the DNA-binding regions and target genes of INI1 significantly overlap with those of MYC. In an INI1-deficient rhabdoid tumor system, we observe that with re-expression of INI1, MYC and INI1 bind to common target genes and have opposing effects on gene expression. Functionally, INI1 re-expression suppresses cell proliferation and MYC-potentiated transformation. Our findings thus establish the antagonistic roles of the INI1 and MYC transcriptional regulators in mediating cellular and oncogenic functions.

A selective sweep is the result of strong positive selection driving newly occurring or standing genetic variants to fixation, and can dramatically alter the pattern and distribution of allelic diversity in a population. Population-level sequencing data have enabled discoveries of selective sweeps associated with genes involved in recent adaptations in many species. In contrast, much debate but little evidence addresses whether selfish genes are capable of fixation-thereby leaving signatures identical to classical selective sweeps-despite being neutral or deleterious to organismal fitness. We previously described R2d2, a large copy-number variant that causes nonrandom segregation of mouse Chromosome 2 in females due to meiotic drive. Here we show population-genetic data consistent with a selfish sweep driven by alleles of R2d2 with high copy number (R2d2(HC)) in natural populations. We replicate this finding in multiple closed breeding populations from six outbred backgrounds segregating for R2d2 alleles. We find that R2d2(HC) rapidly increases in frequency, and in most cases becomes fixed in significantly fewer generations than can be explained by genetic drift. R2d2(HC) is also associated with significantly reduced litter sizes in heterozygous mothers, making it a true selfish allele. Our data provide direct evidence of populations actively undergoing selfish sweeps, and demonstrate that meiotic drive can rapidly alter the genomic landscape in favor of mutations with neutral or even negative effects on overall Darwinian fitness. Further study will reveal the incidence of selfish sweeps, and will elucidate the relative contributions of selfish genes, adaptation and genetic drift to evolution.

Nucleocytoplasmic transport is mediated by the interaction of transport factors (TFs) with disordered phenylalanine-glycine (FG) repeats that fill the central channel of the nuclear pore complex (NPC). However, the mechanism by which TFs rapidly diffuse through multiple FG repeats without compromising NPC selectivity is not yet fully understood. In this study, we build on our recent NMR investigations showing that FG repeats are highly dynamic, flexible, and rapidly exchanging among TF interaction sites. We use unbiased long timescale all-atom simulations on the Anton supercomputer, combined with extensive enhanced sampling simulations and NMR experiments, to characterize the thermodynamic and kinetic properties of FG repeats and their interaction with a model transport factor. Both the simulations and experimental data indicate that FG repeats are highly dynamic random coils, lack intrachain interactions, and exhibit significant entropically driven resistance to spatial confinement. We show that the FG motifs reversibly slide in and out of multiple TF interaction sites, transitioning rapidly between a strongly interacting state and a weakly interacting state, rather than undergoing a much slower transition between strongly interacting and completely noninteracting (unbound) states. In the weakly interacting state, FG motifs can be more easily displaced by other competing FG motifs, providing a simple mechanism for rapid exchange of TF/FG motif contacts during transport. This slide-and-exchange mechanism highlights the direct role of the disorder within FG repeats in nucleocytoplasmic transport, and resolves the apparent conflict between the selectivity and speed of transport.

In mammals, small RNAs are important players in post-transcriptional gene regulation. While their roles in mRNA destabilization and translational repression are well appreciated, their involvement in endonucleolytic cleavage of target RNAs is poorly understood. Very few microRNAs are known to guide RNA cleavage. Endogenous small interfering RNAs are expected to induce target cleavage, but their target genes remain largely unknown. We report a systematic study of small RNA-mediated endonucleolytic cleavage in mouse through integrative analysis of small RNA and degradome sequencing data without imposing any bias toward known small RNAs. Hundreds of small cleavage-inducing RNAs and their cognate target genes were identified, significantly expanding the repertoire of known small RNA-guided cleavage events. Strikingly, both small RNAs and their target sites demonstrated significant overlap with retrotransposons, providing evidence for the long-standing speculation that retrotransposable elements in mRNAs are leveraged as signals for gene targeting. Furthermore, our analysis showed that the RNA cleavage pathway is also present in human cells but affecting a different repertoire of retrotransposons. These results show that small RNA-guided cleavage is more widespread than previously appreciated. Their impact on retrotransposons in non-coding regions shed light on important aspects of mammalian gene regulation.

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

In somatic cells, BRCA2 is needed for RAD51-mediated homologous recombination. The meiosis-specific DNA strand exchange protein, DMC1, promotes the formation of DNA strand invasion products (joint molecules) between homologous molecules in a fashion similar to RAD51. BRCA2 interacts directly with both human RAD51 and DMC1; in the case of RAD51, this interaction results in stimulation of RAD51-promoted DNA strand exchange. However, for DMC1, little is known regarding the basis and functional consequences of its interaction with BRCA2. Here we report that human DMC1 interacts directly with each of the BRC repeats of BRCA2, albeit most tightly with repeats 1-3 and 6-8. However, BRC1-3 bind with higher affinity to RAD51 than to DMC1, whereas BRC6-8 bind with higher affinity to DMC1, providing potential spatial organization to nascent filament formation. With the exception of BRC4, each BRC repeat stimulates joint molecule formation by DMC1. The basis for this stimulation is an enhancement of DMC1-ssDNA complex formation by the stimulatory BRC repeats. Lastly, we demonstrate that full-length BRCA2 protein stimulates DMC1-mediated DNA strand exchange between RPA-ssDNA complexes and duplex DNA, thus identifying BRCA2 as a mediator of DMC1 recombination function. Collectively, our results suggest unique and specialized functions for the BRC motifs of BRCA2 in promoting homologous recombination in meiotic and mitotic cells.

Metazoans start embryogenesis with a relatively naïve genome. The transcriptionally inert, late-replicating heterochromatic regions, including the constitutive heterochromatin on repetitive sequences near centromeres and telomeres, need to be re-established during development. To explore the events initiating heterochromatin formation and examine their temporal control, sequence specificity, and immediate regulatory consequence, we established a live imaging approach that enabled visualization of steps in heterochromatin emergence on specific satellite sequences during the mid-blastula transition (MBT) in Drosophila. Unexpectedly, only a subset of satellite sequences, including the 359-base-pair (bp) repeat sequence, recruited HP1a at the MBT. The recruitment of HP1a to the 359-bp repeat was dependent on HP1a's chromoshadow domain but not its chromodomain and was guided by maternally provided signals. HP1a recruitment to the 359-bp repeat was required for its programmed shift to later replication, and ectopic recruitment of HP1a was sufficient to delay replication timing of a different repeat. Our results reveal that emergence of constitutive heterochromatin follows a stereotyped developmental program in which different repetitive sequences use distinct interactions and independent pathways to arrive at a heterochromatic state. This differential emergence of heterochromatin on various repetitive sequences changes their replication order and remodels the DNA replication schedule during embryonic development.

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F1 individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.

We here report statistical analyses of 76 families of integral outer membrane pore-forming proteins (OMPPs) found in bacteria and eukaryotic organelles. 47 of these families fall into one superfamily (SFI) which segregate into fifteen phylogenetic clusters. Families with members of the same protein size, topology and substrate specificities often cluster together. Virtually all OMPP families include only proteins that form transmembrane pores. Nine such families, all of which cluster together in the SFI phylogenetic tree, contain both α- and β-structures, are multi domain, multi subunit systems, and transport macromolecules. Most other SFI OMPPs transport small molecules. SFII and SFV homologues derive from Actinobacteria while SFIII and SFIV proteins derive from chloroplasts. Three families of actinobacterial OMPPs and two families of eukaryotic OMPPs apparently consist primarily of α-helices (α-TMSs). Of the 71 families of (putative) β-barrel OMPPs, only twenty could not be assigned to a superfamily, and these derived primarily from Actinobacteria (1), chloroplasts (1), spirochaetes (8), and proteobacteria (10). Proteins were identified in which two or three full length OMPPs are fused together. Family characteristic are described and evidence agrees with a previous proposal suggesting that many arose by adjacent β-hairpin structural unit duplications.

Repeat elements are important components of eukaryotic genomes. One limitation in our understanding of repeat elements is that most analyses rely on reference genomes that are incomplete and often contain missing data in highly repetitive regions that are difficult to assemble. To overcome this problem we develop a new method, REPdenovo, which assembles repeat sequences directly from raw shotgun sequencing data. REPdenovo can construct various types of repeats that are highly repetitive and have low sequence divergence within copies. We show that REPdenovo is substantially better than existing methods both in terms of the number and the completeness of the repeat sequences that it recovers. The key advantage of REPdenovo is that it can reconstruct long repeats from sequence reads. We apply the method to human data and discover a number of potentially new repeats sequences that have been missed by previous repeat annotations. Many of these sequences are incorporated into various parasite genomes, possibly because the filtering process for host DNA involved in the sequencing of the parasite genomes failed to exclude the host derived repeat sequences. REPdenovo is a new powerful computational tool for annotating genomes and for addressing questions regarding the evolution of repeat families. The software tool, REPdenovo, is available for download at https://github.com/Reedwarbler/REPdenovo.

BackgroundDowny mildews are the most speciose group of oomycetes and affect crops of great economic importance. So far, there is only a single deeply-sequenced downy mildew genome available, from Hyaloperonospora arabidopsidis. Further genomic resources for downy mildews are required to study their evolution, including pathogenicity effector proteins, such as RxLR effectors. Plasmopara halstedii is a devastating pathogen of sunflower and a potential pathosystem model to study downy mildews, as several Avr-genes and R-genes have been predicted and unlike Arabidopsis downy mildew, large quantities of almost contamination-free material can be obtained easily.ResultsHere a high-quality draft genome of Plasmopara halstedii is reported and analysed with respect to various aspects, including genome organisation, secondary metabolism, effector proteins and comparative genomics with other sequenced oomycetes. Interestingly, the present analyses revealed further variation of the RxLR motif, suggesting an important role of the conservation of the dEER-motif. Orthology analyses revealed the conservation of 28 RxLR-like core effectors among Phytophthora species. Only six putative RxLR-like effectors were shared by the two sequenced downy mildews, highlighting the fast and largely independent evolution of two of the three major downy mildew lineages. This is seemingly supported by phylogenomic results, in which downy mildews did not appear to be monophyletic.ConclusionsThe genome resource will be useful for developing markers for monitoring the pathogen population and might provide the basis for new approaches to fight Phytophthora and downy mildew pathogens by targeting core pathogenicity effectors.

BackgroundDNA cytosine methylation is an important epigenetic modification that has significant effects on a variety of biological processes in animals. Avian species hold a crucial position in evolutionary history. In this study, we used whole-genome bisulfite sequencing (MethylC-seq) to generate single base methylation profiles of lungs in two genetically distinct and highly inbred chicken lines (Fayoumi and Leghorn) that differ in genetic resistance to multiple pathogens, and we explored the potential regulatory role of DNA methylation associated with immune response differences between the two chicken lines.MethodsThe MethylC-seq was used to generate single base DNA methylation profiles of Fayoumi and Leghorn birds. In addition, transcriptome profiling using RNA-seq from the same chickens and tissues were obtained to interrogate how DNA methylation regulates gene transcription on a genome-wide scale.ResultsThe general DNA methylation pattern across different regions of genes was conserved compared to other species except for hyper-methylation of repeat elements, which was not observed in chicken. The methylation level of miRNA and pseudogene promoters was high, which indicates that silencing of these genes may be partially due to promoter hyper-methylation. Interestingly, the promoter regions of more recently evolved genes tended to be more highly methylated, whereas the gene body regions of evolutionarily conserved genes were more highly methylated than those of more recently evolved genes. Immune-related GO (Gene Ontology) terms were significantly enriched from genes within the differentially methylated regions (DMR) between Fayoumi and Leghorn, which implicates DNA methylation as one of the regulatory mechanisms modulating immune response differences between these lines.ConclusionsThis study establishes a single-base resolution DNA methylation profile of chicken lung and suggests a regulatory role of DNA methylation in controlling gene expression and maintaining genome transcription stability. Furthermore, profiling the DNA methylomes of two genetic lines that differ in disease resistance provides a unique opportunity to investigate the potential role of DNA methylation in host disease resistance. Our study provides a foundation for future studies on epigenetic modulation of host immune response to pathogens in chickens.

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

Echinococcosis poses a significant concern in the fields of public health and veterinary care as it can be transmitted between animals and humans. The primary endemic subtypes are cystic echinococcosis (CE) and alveolar echinococcosis (AE), which result from infestation by Echinococcus granulosus and Echinococcus multilocularis , respectively. A prominent epidemic of echinococcosis greatly affects the Tibet Autonomous Region (TAR) in China. A new technique called the loop-mediated isothermal amplification-lateral flow dipstick (LAMP-LFD) test is introduced in this research to differentiate between E. granulosus and E. multilocularis using their repetitive genetic sequences. The test is characterized by its portable nature, simple operation, quick result production, high sensitivity, and low susceptibility to aerosol contamination. The LAMP-LFD method demonstrated an exceptional minimal detection limit, reaching levels as low as approximately 1 fg/μL (femtogram per microliter) of genomic DNA. The assay's specificity was assessed, and no cross-reactivity was seen. A total of 982 dog fecal samples were collected from 54 counties in the TAR region between July 2021 and June 2022. The established method underwent validation using a commercially available ELISA kit. The agreement rate between the LAMP-LFD and ELISA methods was 97.25%, with a sensitivity of 96.05% and a specificity of 97.35%. The assay described in this study improves specificity by using a double-labeled probe, and it reduces the risk of false-positive results caused by aerosol contamination through the use of a sealed device. This makes it a suitable choice for quickly and accurately identifying the two main types of Echinococcus in field settings.

To study heterochromatin distribution differences among tamarins, we applied LINE‐1 probes using fluorescence in situ hybridization onto chromosomes of Saguinus mystax, Leontocebus fuscicollis, and Leontopithecus rosalia with the aim to investigate possible evolutionary implications. LINE‐1 repeats were shown to be involved in genome architecture and in the occurrence of chromosomal rearrangements in many vertebrates. We found bright LINE‐1 probe signals at centromeric or pericentromeric areas, GC rich, on almost all chromosomes in three tamarin species. We also found non‐centromeric signals along chromosome arms. In a phylogenetic perspective, we analyzed the pattern of LINE‐1 distribution considering human chromosomal homologies and C banding patterns. Our data indicate that LINE‐1 centromeric expansions and accumulation presumably arose in a common tamarin ancestor and that the presence of LINE‐1 at the junction of human chromosome associations is presumably linked to interchromosomal rearrangements. For example, we found bright centromeric signals as well as non‐centromeric signals on chromosomes 1 and 2, in all species analyzed, in correspondence to human chromosome associations 13/9/22 and 20/17/13, which are synapomorphic for all tamarins. Furthermore, we found other faint signals that could be apomorphisms linked both to intrachromosomal rearrangements as well as to retro‐transposition events. Our results confirm that the three species have similar karyotypes but small differences in LINE‐1 and heterochromatin amplification and distribution; in particular on chromosome pairs 19–22, where we show the occurrence of small inversions, in agreement with previous classic cytogenetic hypotheses. [ABSTRACT FROM AUTHOR]

Genome sizes of eukaryotic organisms vary substantially, with whole-genome duplications (WGD) and transposable element expansion acting as main drivers for rapid genome size increase. The two North American mudminnows, Umbra limi and Umbra pygmaea , feature genomes about twice the size of their sister lineage Esocidae (e.g. pikes and pickerels). However, it is unknown whether all Umbra species share this genome expansion and which causal mechanisms drive this expansion. Using flow cytometry, we find that the genome of the European mudminnow is expanded similarly to both North American species, ranging between 4.5 and 5.4 pg per diploid nucleus. Observed blocks of interstitially located telomeric repeats in U. limi suggest frequent Robertsonian rearrangements in its history. Comparative analyses of transcriptome and genome assemblies show that the genome expansion in Umbra is driven by the expansion of DNA transposon and unclassified repeat sequences without WGD. Furthermore, we find a substantial ongoing expansion of repeat sequences in the Alaska blackfish Dallia pectoralis , the closest relative to the family Umbridae, which might mark the beginning of a similar genome expansion. Our study suggests that the genome expansion in mudminnows, driven mainly by transposon expansion, but not WGD, occurred before the separation into the American and European lineage. [ABSTRACT FROM AUTHOR]

The loops of modular polyketide synthases (PKSs) serve diverse functions but are largely uncharacterized. They frequently contain amino acid repeats resulting from genetic events such as slipped‐strand mispairing. Determining the tolerance of loops to amino acid changes would aid in understanding and engineering these multidomain molecule factories. Here, tandem repeats in the DNA encoding 949 modules within 129 cis‐acyltransferase PKSs were cataloged, and the locations of the corresponding amino acids within the module were identified. The most frequently inserted interdomain loop corresponds with the updated module boundary immediately downstream of the ketosynthase (KS), while the loops bordering the dehydratase are nearly intolerant to such insertions. From the 949 modules, no repetitive sequence loop insertions are located within ACP, and only 2 reside within KS, indicating the sensitivity of these domains to alteration. [ABSTRACT FROM AUTHOR]

CRISPR-Cas adaptive immune systems protect bacteria and archaea against foreign genetic elements. In Escherichia coli, Cascade (CRISPR-associated complex for antiviral defense) is an RNA-guided surveillance complex that binds foreign DNA and recruits Cas3, a trans-acting nuclease helicase for target degradation. Here, we use single-molecule imaging to visualize Cascade and Cas3 binding to foreign DNA targets. Our analysis reveals two distinct pathways dictated by the presence or absence of a protospacer-adjacent motif (PAM). Binding to a protospacer flanked by a PAM recruits a nuclease-active Cas3 for degradation of short single-stranded regions of target DNA, whereas PAM mutations elicit an alternative pathway that recruits a nuclease-inactive Cas3 through a mechanism that is dependent on the Cas1 and Cas2 proteins. These findings explain how target recognition by Cascade can elicit distinct outcomes and support a model for acquisition of new spacer sequences through a mechanism involving processive, ATP-dependent Cas3 translocation along foreign DNA.

Vertebrates developed immunoglobulin heavy chain (IgH) class switch recombination (CSR) to express different IgH constant regions. Most double-strand breaks for Ig CSR occur within the repetitive portion of the switch regions located upstream of each set of constant domain exons for the Igγ, Igα or Igϵ heavy chain. Unlike mammalian switch regions, Xenopus switch regions do not have a high G-density on the non-template DNA strand. In previous studies, when Xenopus Sμ DNA was moved to the genome of mice, it is able to support substantial CSR when it is used to replace the murine Sγ1 region. Here, we tested both the 2kb repetitive portion and the 4.6 kb full-length portions of the Xenopus Sμ in both their natural (forward) orientation relative to the constant domain exons, as well as the opposite (reverse) orientation. Consistent with previous work, we find that the 4.6 kb full-length Sμ mediates similar levels of CSR in both the forward and reverse orientations. Whereas, the forward orientation of the 2kb portion can restore the majority of the CSR level of the 4.6 kb full-length Sμ, the reverse orientation poorly supports R-looping and no CSR. The forward orientation of the 2kb repetitive portion has more GG dinucleotides on the non-template strand than the reverse orientation. The correlation of R-loop formation with CSR efficiency, as demonstrated in the 2kb repetitive fragment of the Xenopus switch region, confirms a role played by R-looping in CSR that appears to be conserved through evolution.

Huntington disease (HD) is most widely known for its selective degeneration of striatal neurons but there is also growing evidence for white matter (WM) deterioration. The primary objective of this research was to conduct a large-scale analysis using multisite diffusion-weighted imaging (DWI) tractography data to quantify diffusivity properties along major prefrontal cortex WM tracts in prodromal HD. Fifteen international sites participating in the PREDICT-HD study collected imaging and neuropsychological data on gene-positive HD participants without a clinical diagnosis (i.e., prodromal) and gene-negative control participants. The anatomical prefrontal WM tracts of the corpus callosum (PFCC), anterior thalamic radiations (ATRs), inferior fronto-occipital fasciculi (IFO), and uncinate fasciculi (UNC) were identified using streamline tractography of DWI. Within each of these tracts, tensor scalars for fractional anisotropy, mean diffusivity, radial diffusivity, and axial diffusivity coefficients were calculated. We divided prodromal HD subjects into three CAG-age product (CAP) groups having Low, Medium, or High probabilities of onset indexed by genetic exposure. We observed significant differences in WM properties for each of the four anatomical tracts for the High CAP group in comparison to controls. Additionally, the Medium CAP group presented differences in the ATR and IFO in comparison to controls. Furthermore, WM alterations in the PFCC, ATR, and IFO showed robust associations with neuropsychological measures of executive functioning. These results suggest long-range tracts essential for cross-region information transfer show early vulnerability in HD and may explain cognitive problems often present in the prodromal stage. Hum Brain Mapp 36:3717-3732, 2015. © 2015 Wiley Periodicals, Inc.

Ferns are one of the few remaining major clades of land plants for which a complete genome sequence is lacking. Knowledge of genome space in ferns will enable broad-scale comparative analyses of land plant genes and genomes, provide insights into genome evolution across green plants, and shed light on genetic and genomic features that characterize ferns, such as their high chromosome numbers and large genome sizes. As part of an initial exploration into fern genome space, we used a whole genome shotgun sequencing approach to obtain low-density coverage (∼0.4X to 2X) for six fern species from the Polypodiales (Ceratopteris, Pteridium, Polypodium, Cystopteris), Cyatheales (Plagiogyria), and Gleicheniales (Dipteris). We explore these data to characterize the proportion of the nuclear genome represented by repetitive sequences (including DNA transposons, retrotransposons, ribosomal DNA, and simple repeats) and protein-coding genes, and to extract chloroplast and mitochondrial genome sequences. Such initial sweeps of fern genomes can provide information useful for selecting a promising candidate fern species for whole genome sequencing. We also describe variation of genomic traits across our sample and highlight some differences and similarities in repeat structure between ferns and seed plants.

Repetitive elements (REs) make up the vast majority of the mammalian genomes. We identified species-specific genomic libraries of RE arrays. The non-random configurations of RE arrays suggest their functions. We tested whether RE arrays undergo age- and tissue/cell-specific rearrangements. An RE array of C57BL/6J mice, containing tandem repeats of a mosaic of transposable REs, was selected to examine rearrangements in different ages and tissues. There were marked changes in the array configuration in the genomes of the skin and brain in all mice of six weeks and older, whereas the heart and liver had alterations at 29weeks. The temporal variations were confirmed by identifying putative rearrangement junctions. Temporal and spatial rearrangements of certain RE arrays may contribute to the acquired characteristics of the genome information system.

The Muller F element (4.2 Mb, ~80 protein-coding genes) is an unusual autosome of Drosophila melanogaster; it is mostly heterochromatic with a low recombination rate. To investigate how these properties impact the evolution of repeats and genes, we manually improved the sequence and annotated the genes on the D. erecta, D. mojavensis, and D. grimshawi F elements and euchromatic domains from the Muller D element. We find that F elements have greater transposon density (25-50%) than euchromatic reference regions (3-11%). Among the F elements, D. grimshawi has the lowest transposon density (particularly DINE-1: 2% vs. 11-27%). F element genes have larger coding spans, more coding exons, larger introns, and lower codon bias. Comparison of the Effective Number of Codons with the Codon Adaptation Index shows that, in contrast to the other species, codon bias in D. grimshawi F element genes can be attributed primarily to selection instead of mutational biases, suggesting that density and types of transposons affect the degree of local heterochromatin formation. F element genes have lower estimated DNA melting temperatures than D element genes, potentially facilitating transcription through heterochromatin. Most F element genes (~90%) have remained on that element, but the F element has smaller syntenic blocks than genome averages (3.4-3.6 vs. 8.4-8.8 genes per block), indicating greater rates of inversion despite lower rates of recombination. Overall, the F element has maintained characteristics that are distinct from other autosomes in the Drosophila lineage, illuminating the constraints imposed by a heterochromatic milieu.

MotivationStructural variation is common in human and cancer genomes. High-throughput DNA sequencing has enabled genome-scale surveys of structural variation. However, the short reads produced by these technologies limit the study of complex variants, particularly those involving repetitive regions. Recent 'third-generation' sequencing technologies provide single-molecule templates and longer sequencing reads, but at the cost of higher per-nucleotide error rates.ResultsWe present MultiBreak-SV, an algorithm to detect structural variants (SVs) from single molecule sequencing data, paired read sequencing data, or a combination of sequencing data from different platforms. We demonstrate that combining low-coverage third-generation data from Pacific Biosciences (PacBio) with high-coverage paired read data is advantageous on simulated chromosomes. We apply MultiBreak-SV to PacBio data from four human fosmids and show that it detects known SVs with high sensitivity and specificity. Finally, we perform a whole-genome analysis on PacBio data from a complete hydatidiform mole cell line and predict 1002 high-probability SVs, over half of which are confirmed by an Illumina-based assembly.

The boundaries of R-loops are well-documented at immunoglobulin heavy chain loci in mammalian B cells. Within primary B cells or B cell lines, the upstream boundaries of R-loops typically begin early in the repetitive portion of the switch regions. Most R-loops terminate within the switch repetitive zone, but the remainder can extend a few hundred base pairs further, where G-density on the non-template DNA strand gradually drops to the genome average. Whether the G-density determines how far the R-loops extend is an important question. We previously studied the role of G-clusters in initiating R-loop formation, but we did not examine the role of G-density in permitting the elongation of the R-loop, after it had initiated. Here, we vary the G-density of different portions of the switch region in a murine B cell line. We find that both class switch recombination (CSR) and R-loop formation decrease significantly when the overall G-density is reduced from 46% to 29%. Short 50 bp insertions with low G-density within switch regions do not appear to affect either CSR or R-loop elongation, whereas a longer (150 bp) insertion impairs both. These results demonstrate that G-density is an important determinant of the length over which mammalian genomic R-loops extend.

Recurrent deletions of chromosome 15q13.3 associate with intellectual disability, schizophrenia, autism and epilepsy. To gain insight into the instability of this region, we sequenced it in affected individuals, normal individuals and nonhuman primates. We discovered five structural configurations of the human chromosome 15q13.3 region ranging in size from 2 to 3 Mb. These configurations arose recently (∼0.5-0.9 million years ago) as a result of human-specific expansions of segmental duplications and two independent inversion events. All inversion breakpoints map near GOLGA8 core duplicons-a ∼14-kb primate-specific chromosome 15 repeat that became organized into larger palindromic structures. GOLGA8-flanked palindromes also demarcate the breakpoints of recurrent 15q13.3 microdeletions, the expansion of chromosome 15 segmental duplications in the human lineage and independent structural changes in apes. The significant clustering (P = 0.002) of breakpoints provides mechanistic evidence for the role of this core duplicon and its palindromic architecture in promoting the evolutionary and disease-related instability of chromosome 15.

BackgroundPeriodic proteins, characterized by the presence of multiple repeats of short motifs, form an interesting and seldom-studied group. Due to often extreme divergence in sequence, detection and analysis of such motifs is performed more reliably on the structural level. Yet, few algorithms have been developed for the detection and analysis of structures of periodic proteins.ResultsConSole recognizes modularity in protein contact maps, allowing for precise identification of repeats in solenoid protein structures, an important subgroup of periodic proteins. Tests on benchmarks show that ConSole has higher recognition accuracy as compared to Raphael, the only other publicly available solenoid structure detection tool. As a next step of ConSole analysis, we show how detection of solenoid repeats in structures can be used to improve sequence recognition of these motifs and to detect subtle irregularities of repeat lengths in three solenoid protein families.ConclusionsThe ConSole algorithm provides a fast and accurate tool to recognize solenoid protein structures as a whole and to identify individual solenoid repeat units from a structure. ConSole is available as a web-based, interactive server and is available for download at http://console.sanfordburnham.org.

Insects act as vectors for diseases of plants, animals, and humans. Replacement of wild insect populations with genetically modified individuals unable to transmit disease provides a potentially self-perpetuating method of disease prevention. Population replacement requires a gene drive mechanism in order to spread linked genes mediating disease refractoriness through wild populations. We previously reported the creation of synthetic Medea selfish genetic elements able to drive population replacement in Drosophila. These elements use microRNA-mediated silencing of myd88, a maternally expressed gene required for embryonic dorso-ventral pattern formation, coupled with early zygotic expression of a rescuing transgene, to bring about gene drive. Medea elements that work through additional mechanisms are needed in order to be able to carry out cycles of population replacement and/or remove existing transgenes from the population, using second-generation elements that spread while driving first-generation elements out of the population. Here we report the synthesis and population genetic behavior of two new synthetic Medea elements that drive population replacement through manipulation of signaling pathways involved in cellular blastoderm formation or Notch signaling, demonstrating that in Drosophila Medea elements can be generated through manipulation of diverse signaling pathways. We also describe the mRNA and small RNA changes in ovaries and early embryos associated from Medea-bearing females. Finally, we use modeling to illustrate how Medea elements carrying genes that result in diapause-dependent female lethality could be used to bring about population suppression.

Fibroblast growth factor receptor 2 (FGFR2) promotes osteoprogenitor proliferation and differentiation during bone development, yet how the receptor elicits these distinct cellular responses remains unclear. Analysis of the FGFR2-skeletal disorder bent bone dysplasia syndrome (BBDS) demonstrates that FGFR2, in addition to its canonical signaling activities at the plasma membrane, regulates bone formation from within the nucleolus. Previously, we showed that the unique FGFR2 mutations that cause BBDS reduce receptor levels at the plasma membrane and diminish responsiveness to extracellular FGF2. In this study, we find that these mutations, despite reducing canonical signaling, enhance nucleolar occupancy of FGFR2 at the ribosomal DNA (rDNA) promoter. Nucleolar FGFR2 activates rDNA transcription via interactions with FGF2 and UBF1 by de-repressing RUNX2. An increase in the nucleolar activity of FGFR2 in BBDS elevates levels of ribosomal RNA in the developing bone, consequently promoting osteoprogenitor cell proliferation and decreasing differentiation. Identifying FGFR2 as a transcriptional regulator of rDNA in bone unexpectedly reveals a nucleolar route for FGF signaling that allows for independent regulation of osteoprogenitor cell proliferation and differentiation.

Growth factor receptor levels are aberrantly high in diverse cancers, driving the proliferation and survival of tumor cells. Understanding the molecular basis for this aberrant elevation has profound clinical implications. Here we show that the pleckstrin homology domain leucine-rich repeat protein phosphatase (PHLPP) suppresses receptor tyrosine kinase (RTK) signaling output by a previously unidentified epigenetic mechanism unrelated to its previously described function as the hydrophobic motif phosphatase for the protein kinase AKT, protein kinase C, and S6 kinase. Specifically, we show that nuclear-localized PHLPP suppresses histone phosphorylation and acetylation, in turn suppressing the transcription of diverse growth factor receptors, including the EGF receptor. These data uncover a much broader role for PHLPP in regulation of growth factor signaling beyond its direct inactivation of AKT: By suppressing RTK levels, PHLPP dampens the downstream signaling output of two major oncogenic pathways, the PI3 kinase/AKT and the Rat sarcoma (RAS)/ERK pathways. Our data are consistent with a model in which PHLPP modifies the histone code to control the transcription of RTKs.

Fragile X-associated tremor/ataxia syndrome (FXTAS) is a late-onset neurodegenerative disorder affecting carriers of premutation forms of the FMR1 gene, resulting in a progressive development of tremor, ataxia and neuropsychological problems. The disease is caused by an expanded CGG repeat in the FMR1 gene, leading to an RNA gain-of-function toxicity mechanism. In order to study the pathogenesis of FXTAS, new inducible transgenic mouse models have been developed that expresses either 11CGGs or 90CGGs at the RNA level under control of a Tet-On promoter. When bred to an hnRNP-rtTA driver line, doxycycline (dox) induced expression of the transgene could be found in almost all tissues. Dox exposure resulted in loss of weight and death within 5 d for the 90CGG RNA expressing mice. Immunohistochemical examination of tissues of these mice revealed steatosis and apoptosis in the liver. Decreased expression of GPX1 and increased expression of cytochrome C is found. These effects were not seen in mice expressing a normal sized 11CGG repeat. In conclusion, we were able to show in vivo that expression of an expanded CGG-repeat rather than overexpression of a normal CGG-repeat causes pathology. In addition, we have shown that expanded CGG RNA expression can cause mitochondrial dysfunction by regulating expression levels of several markers. Although FTXAS patients do not display liver abnormalities, our findings contribute to understanding of the molecular mechanisms underlying toxicity of CGG repeat RNA expression in an animal model. In addition, the dox inducible mouse lines offer new opportunities to study therapeutic interventions for FXTAS.

Chromatin regulatory proteins affect diverse developmental and environmental response pathways via their influence on nuclear processes such as the regulation of gene expression. Through a genome-wide genetic screen, we implicate a novel protein called X-chromosome-associated protein 5 (Xap5) in chromatin regulation. We show that Xap5 is a chromatin-associated protein acting in a similar manner as the histone variant H2A.Z to suppress expression of antisense and repeat element transcripts throughout the fission yeast genome. Xap5 is highly conserved across eukaryotes, and a plant homolog rescues xap5 mutant yeast. We propose that Xap5 likely functions as a chromatin regulator in diverse organisms.

Many epidemiologic studies of environmental exposures and disease susceptibility measure DNA methylation in white blood cells (WBC). Some studies are also starting to use saliva DNA as it is usually more readily available in large epidemiologic studies. However, little is known about the correlation of methylation between WBC and saliva DNA. We examined DNA methylation in three repetitive elements, Sat2, Alu, and LINE-1, and in four CpG sites, including AHRR (cg23576855, cg05575921), cg05951221 at 2q37.1, and cg11924019 at CYP1A1, in 57 girls aged 6-15 years with blood and saliva collected on the same day. We measured all DNA methylation markers by bisulfite-pyrosequencing, except for Sat2 and Alu, which were measured by the MethyLight assay. Methylation levels measured in saliva DNA were lower than those in WBC DNA, with differences ranging from 2.8% for Alu to 14.1% for cg05575921. Methylation levels for the three repetitive elements measured in saliva DNA were all positively correlated with those in WBC DNA. However, there was a wide range in the Spearman correlations, with the smallest correlation found for Alu (0.24) and the strongest correlation found for LINE-1 (0.73). Spearman correlations for cg05575921, cg05951221, and cg11924019 were 0.33, 0.42, and 0.79, respectively. If these findings are replicated in larger studies, they suggest that, for selected methylation markers (e.g., LINE-1), methylation levels may be highly correlated between blood and saliva, while for others methylation markers, the levels may be more tissue specific. Thus, in studies that differ by DNA source, each interrogated site should be separately examined in order to evaluate the correlation in DNA methylation levels across DNA sources.

The transport of cargo across the nuclear membrane is highly selective and accomplished by a poorly understood mechanism involving hundreds of nucleoporins lining the inside of the nuclear pore complex (NPC). Currently, there is no clear picture of the overall structure formed by this collection of proteins within the pore, primarily due to their disordered nature. We perform coarse-grained simulations of both individual nucleoporins and grafted rings of nups mimicking the in vivo geometry of the NPC and supplement this with polymer brush modeling. Our results indicate that different regions or blocks of an individual NPC protein can have distinctly different forms of disorder and that this property appears to be a conserved functional feature. Furthermore, this block structure at the individual protein level is critical to the formation of a unique higher-order polymer brush architecture that can exist in distinct morphologies depending on the effective interaction energy between the phenylalanine glycine (FG) domains of different nups. Because the interactions between FG domains may be modulated by certain forms of transport factors, our results indicate that transitions between brush morphologies could play an important role in regulating transport across the NPC, suggesting novel forms of gated transport across membrane pores with wide biomimetic applicability.

With eukaryotic genome replication, incomplete telomere synthesis results in chromosome shortening and eventual compromise of genome stability. Telomerase counteracts this terminal sequence loss by synthesizing telomeric repeats through repeated cycles of reverse transcription of its internal RNA template. Using human telomerase domain-complementation assays for telomerase reverse transcriptase protein (TERT) and RNA in combination with the first direct footprinting assay for telomerase association with bound DNA, we resolve mechanisms by which TERT domains and RNA motifs direct repeat synthesis. Surprisingly, we find that product-template hybrid is sensed in a length- and sequence-dependent manner to set the template 5' boundary. We demonstrate that the TERT N-terminal (TEN) domain determines active-site use of the atypically short primer-template hybrid necessary for telomeric-repeat synthesis. Also against expectation, we show that the remainder of TERT (the TERT ring) supports functional recognition and physical protection of single-stranded DNA adjacent to the template hybrid. These findings establish unprecedented polymerase recognition specificities for DNA-RNA hybrid and single-stranded DNA and suggest a new perspective on the mechanisms of telomerase specialization for telomeric-repeat synthesis.

UnlabelledGenomic analysis of a large set of phages infecting the common host Mycobacterium smegmatis mc(2)155 shows that they span considerable genetic diversity. There are more than 20 distinct types that lack nucleotide similarity with each other, and there is considerable diversity within most of the groups. Three newly isolated temperate mycobacteriophages, Bongo, PegLeg, and Rey, constitute a new group (cluster M), with the closely related phages Bongo and PegLeg forming subcluster M1 and the more distantly related Rey forming subcluster M2. The cluster M mycobacteriophages have siphoviral morphologies with unusually long tails, are homoimmune, and have larger than average genomes (80.2 to 83.7 kbp). They exhibit a variety of features not previously described in other mycobacteriophages, including noncanonical genome architectures and several unusual sets of conserved repeated sequences suggesting novel regulatory systems for both transcription and translation. In addition to containing transfer-messenger RNA and RtcB-like RNA ligase genes, their genomes encode 21 to 24 tRNA genes encompassing complete or nearly complete sets of isotypes. We predict that these tRNAs are used in late lytic growth, likely compensating for the degradation or inadequacy of host tRNAs. They may represent a complete set of tRNAs necessary for late lytic growth, especially when taken together with the apparent lack of codons in the same late genes that correspond to tRNAs that the genomes of the phages do not obviously encode.ImportanceThe bacteriophage population is vast, dynamic, and old and plays a central role in bacterial pathogenicity. We know surprisingly little about the genetic diversity of the phage population, although metagenomic and phage genome sequencing indicates that it is great. Probing the depth of genetic diversity of phages of a common host, Mycobacterium smegmatis, provides a higher resolution of the phage population and how it has evolved. Three new phages constituting a new cluster M further expand the diversity of the mycobacteriophages and introduce novel features. As such, they provide insights into phage genome architecture, virion structure, and gene regulation at the transcriptional and translational levels.

BackgroundSpider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals.ResultsWe present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region.ConclusionsUsing the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

Evidence has been presented that 5+5 TMS and 7+7 TMS inverted repeat fold transporters are members of a single superfamily named the Amino acid-Polyamine-organoCation (APC) superfamily. However, the evolutionary relationship between the 5+5 and the 7+7 topological types has not been established. We have identified a common fold, consisting of a spiny membrane helix/sheet, followed by a U-like structure and a V-like structure that is recurrent between domain duplicated units of 5+5 and 7+7 inverted repeat folds. This fold is found in the following protein structures: AdiC, ApcT, LeuT, Mhp1, BetP, CaiT, and SglT (all 5+5 TMS repeats), as well as UraA and SulP (7+7 TMS repeats). AdiC, LeuT and Mhp1 have two extra TMSs after the second duplicated domain, SglT has four extra C-terminal TMSs, and BetP has two extra TMSs before the first duplicated domain. UraA and SulP on the other hand have two extra TMSs at the N-terminus of each duplicated domain unit. These observations imply that multiple hairpin and domain duplication events occurred during the evolution of the APC superfamily. We suggest that the five TMS architecture was primordial and that families gained two TMSs on either side of this basic structure via dissimilar hairpin duplications either before or after intragenic duplication. Evidence for homology between TMSs 1-2 of AdiC and TMSs 1-2 and 3-4 of UraA suggests that the 7+7 topology arose via an internal duplication of the N-terminal hairpin loop within the five TMS repeat unit followed by duplication of the 7 TMS domain.

BackgroundSpider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals.ResultsWe present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region.ConclusionsUsing the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

Microsatellites are repetitive DNA sequences usually found in non-coding regions of the genome. Their quantification and analysis have applications in fields from population genetics to evolutionary biology. As genome assemblies become commonplace, the need for software that can facilitate analyses has never been greater. In particular, R packages that can analyze genomic data are particularly important since this is one of the most popular software environments for biologists. We created an R package, micRocounter, to quantify microsatellites. We have optimized our package for speed, accessibility, and portability, making the automated analysis of large genomic data sets feasible. Computationally intensive algorithms were built in C++ to increase speed. Tests using benchmark datasets show a 200-fold improvement in speed over existing software. A moderately sized genome of 500 Mb can be processed in under 50 sec. Results are output as an object in R increasing accessibility and flexibility for practitioners.

Repetitive sequences represent about half of the human genome. They are actively transcribed and play a role during development and in epigenetic regulation. The altered activity of repetitive sequences can lead to genomic instability and they can contribute to the establishment or the progression of degenerative diseases and cancer transformation. In this work, we analyzed the expression profiles of DNA repetitive sequences in the breast cancer specimens of the HMUCC cohort. Satellite expression is generally upregulated in breast cancers, with specific families upregulated per histotype: in HER2-enriched cancers, they are the human satellite II (HSATII), in luminal A and B, they are part of the ALR family and in triple-negative, they are part of SAR and GSAT families, together with a perturbation in the transcription from endogenous retroviruses and their LTR sequences. We report that the background expression of repetitive sequences in healthy tissues of cancer patients differs from the tissues of non-cancerous controls. To conclude, peculiar patterns of expression of repetitive sequences are reported in each specimen, especially in the case of transcripts arising from satellite repeats.

Background and Purpose: Candida albicans is the most prevalent Candida species isolated from animals. Candidiasis can be systemic in animals or may affect a single organ, such as the mouth, urinary tract, and skin. The aim of the present study was to determine the genetic diversity of C. albicans isolated from different animals and investigate the presence of a relationship between host specificity and genetic typing of C. albicans. Materials and Methods: For the purpose of the study, DNA extraction was performed on 27 clinical isolates of C. albicans obtained from animals. Subsequently, they were subjected to 25S ribosomal DNA amplification and ALT repeats in repetitive sequences (RPSs). The minimum inhibitory concentrations of fluconazole, ketoconazole, clotrimazole, nystatin, amphotericin B, and caspofungin were determined using the microdilution method based on the Clinical and Laboratory Standards Institute M27-S4 standard. Results: Out of 27 C. albicans strains, 11, 6, 5, and 5 cases were recognized as genotypes A (40.8%), E (22.2%), B (18.5%), and C (18.5%), respectively, through amplification using AS-I, which revealed 17 different types of C. albicans. By combining the two typing methods, 27 C. albicans strains were finally divided into 22 genotypes. Conclusion: Different genotypes showed genetic diversity among the C. albicans strains isolated from animal sources. The results revealed no special genotype relationship according to the host, anatomical source of isolation, and antifungal susceptibility.

Visual rhodopsins are recognized members of the large and diverse family of G protein-coupled receptors (GPCRs), but their evolutionary origin and relationships to other proteins are not known. In a previous paper [Shlykov MA, Zheng WH, Chen JS & Saier MH Jr (2012) Biochim Biophys Acta 1818, 703-717], we characterized the 4-toluene sulfonate uptake permease (TSUP) family of transmembrane proteins, and showed that these 7-transmembrane segment (TMS) or 8-TMS proteins arose by intragenic duplication of a gene encoding a 4-TMS protein, sometimes followed by loss of a terminal TMS. In this study, we show that the TSUP, GPCR and microbial rhodopsin families are related to each other and to six other currently recognized transport protein families. We designate this superfamily the transporter/opsin/G protein-coupled receptor (TOG) superfamily. Despite their 8-TMS origins, the members of most constituent families exhibit 7-TMS topologies that are well conserved, and these arose by loss of either the N-terminal TMS (more frequent) or the C-terminal TMS (less frequent), depending on the family. Phylogenetic analyses revealed familial relationships within the superfamily and protein relationships within each of the nine families. The results of the statistical analyses leading to the conclusion of homology were confirmed using hidden Markov models, Pfam and 3D superimpositions. Proteins functioning by dissimilar mechanisms (channels, primary active transporters, secondary active transporters, group translocators and receptors) are interspersed on a phylogenetic tree of the TOG superfamily, suggesting that changes in the transport and energy-coupling mechanisms occurred multiple times during evolution of this superfamily.

We set out to determine whether expansions in the C9ORF72 repeat found in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) families are associated with Parkinson disease (PD). We determined the repeat size in a total of 889 clinically ascertained patients (including PD and essential tremor plus Parkinsonism (ETP)) and 1144 controls using a repeat-primed PCR assay. We found that large C9ORF72 repeat expansions (>30 repeats) were not contributing to PD risk. However, PD and ETP cases had a significant increase in intermediate (>20 to 30+) repeat copies compared to controls. Overall, 14 cases (13 PD, 1 ETP) and three controls had >20 repeat copies (Fisher's exact test p = 0.002). Further, seven cases and no controls had >23 repeat copies (p = 0.003). Our results suggest that intermediate copy numbers of the C9ORF72 repeat contribute to risk for PD and ETP. This also suggests that PD, ALS and FTD share some pathophysiological mechanisms of disease. Further studies are needed to elucidate the contribution of the C9ORF72 repeat in the overall PD population and to determine whether other common genetic risk factors exist between these neurodegenerative disorders.

BACKGROUND: ATP-dependent proteases translocate and unfold their substrates. RESULTS: A human virus sequence with only Gly and Ala residues causes similar dysfunctions of eukaryotic and prokaryotic protease motors: unfolding failure. CONCLUSION: Sequences with amino acids of simple shape and small size impair unfolding of contiguous stable domains. SIGNIFICANCE: Compartmented ATP-dependent proteases of diverse origin share conserved principles of interaction between translocase/effector and substrate/recipient. ATP-dependent proteases engage, translocate, and unfold substrate proteins. A sequence with only Gly and Ala residues (glycine-alanine repeat; GAr) encoded by the Epstein-Barr virus of humans inhibits eukaryotic proteasome activity. It causes the ATPase translocase to slip on its protein track, stalling unfolding and interrupting degradation. The bacterial protease ClpXP is structurally simpler than the proteasome but has related elements: a regulatory ATPase complex (ClpX) and associated proteolytic chamber (ClpP). In this study, GAr sequences were found to impair ClpXP function much as in proteasomes. Stalling depended on interaction between a GAr and a suitably spaced and positioned folded domain resistant to mechanical unfolding. Persistent unfolding failure results in the interruption of degradation and the production of partial degradation products that include the resistant domain. The capacity of various sequences to cause unfolding failure was investigated. Among those tested, a GAr was most effective, implying that viral selection had optimized processivity failure. More generally, amino acids of simple shape and small size promoted unfolding failure. The ClpX ATPase is a homohexamer. Partial degradation products could exit the complex through transient gaps between the ClpX monomers or, alternatively, by backing out. Production of intermediates by diverse topological forms of the hexamer was shown to be similar, excluding lateral escape. In principle, a GAr could interrupt degradation because 1) the translocase thrusts forward less effectively or because 2) the translocase retains substrate less well when resetting between forward strokes. Kinetic analysis showed that the predominant effect was through the second of these mechanisms.

Transcription activator-like effectors (TALEs) have revolutionized the field of genome engineering. We present here a systematic assessment of TALE DNA recognition, using quantitative electrophoretic mobility shift assays and reporter gene activation assays. Within TALE proteins, tandem 34-amino acid repeats recognize one base pair each and direct sequence-specific DNA binding through repeat variable di-residues (RVDs). We found that RVD choice can affect affinity by four orders of magnitude, with the relative RVD contribution in the order NG > HD ≈ NN >> NI > NK. The NN repeat preferred the base G over A, whereas the NK repeat bound G with 10(3)-fold lower affinity. We compared AvrBs3, a naturally occurring TALE that recognizes its target using some atypical RVD-base combinations, with a designed TALE that precisely matches 'standard' RVDs with the target bases. This comparison revealed unexpected differences in sensitivity to substitutions of the invariant 5'-T. Another surprising observation was that base mismatches at the 5' end of the target site had more disruptive effects on affinity than those at the 3' end, particularly in designed TALEs. These results provide evidence that TALE-DNA recognition exhibits a hitherto un-described polarity effect, in which the N-terminal repeats contribute more to affinity than C-terminal ones.

Spider silk fibers have impressive mechanical properties and are primarily composed of highly repetitive structural proteins (termed spidroins) encoded by a single gene family. Most characterized spidroin genes are incompletely known because of their extreme size (typically >9 kb) and repetitiveness, limiting understanding of the evolutionary processes that gave rise to their unusual gene architectures. The only complete spidroin genes characterized thus far form the dragline in the Western black widow, Latrodectus hesperus. Here, we describe the first complete gene sequence encoding the aciniform spidroin AcSp1, the primary component of spider prey-wrapping fibers. L. hesperus AcSp1 contains a single enormous (∼19 kb) exon. The AcSp1 repeat sequence is exceptionally conserved between two widow species (∼94% identity) and between widows and distantly related orb-weavers (∼30% identity), consistent with a history of strong purifying selection on its amino acid sequence. Furthermore, the 16 repeats (each 371-375 amino acids long) found in black widow AcSp1 are, on average, >99% identical at the nucleotide level. A combination of stabilizing selection on amino acid sequence, selection on silent sites, and intragenic recombination likely explains the extreme homogenization of AcSp1 repeats. In addition, phylogenetic analyses of spidroin paralogs support a gene duplication event occurring concomitantly with specialization of the aciniform glands and the tubuliform glands, which synthesize egg-case silk. With repeats that are dramatically different in length and amino acid composition from dragline spidroins, our L. hesperus AcSp1 expands the knowledge base for developing silk-based biomimetic technologies.