Beyond junk-variable tandem repeats as facilitators of rapid evolution of regulatory and coding sequences

Hytrosavirus genetic diversity and eco-regional spread in Glossina species

BACKGROUND

The management of the tsetse species Glossina pallidipes (Diptera; Glossinidae) in Africa by the sterile insect technique (SIT) has been hindered by infections of G. pallidipes production colonies with Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; Hytrosaviridae family). This virus can significantly decrease productivity of the G. pallidipes colonies. Here, we used three highly diverged genes and two variable number tandem repeat regions (VNTRs) of the GpSGHV genome to identify the viral haplotypes in seven Glossina species obtained from 29 African locations and determine their phylogenetic relatedness.

RESULTS

GpSGHV was detected in all analysed Glossina species using PCR. The highest GpSGHV prevalence was found in G. pallidipes colonized at FAO/IAEA Insect Pest Control Laboratory (IPCL) that originated from Uganda (100%) and Tanzania (88%), and a lower prevalence in G. morsitans morsitans from Tanzania (58%) and Zimbabwe (20%). Whereas GpSGHV was detected in 25-40% of G. fuscipes fuscipes in eastern Uganda, the virus was not detected in specimens of neighboring western Kenya. Most of the identified 15 haplotypes were restricted to specific Glossina species in distinct locations. Seven haplotypes were found exclusively in G. pallidipes. The reference haplotype H1 (GpSGHV-Uga; Ugandan strain) was the most widely distributed, but was not found in G. swynnertoni GpSGHV. The 15 haplotypes clustered into three distinct phylogenetic clades, the largest contained seven haplotypes, which were detected in six Glossina species. The G. pallidipes-infecting haplotypes H10, H11 and H12 (from Kenya) clustered with H7 (from Ethiopia), which presumably corresponds to the recently sequenced GpSGHV-Eth (Ethiopian) strain. These four haplotypes diverged the most from the reference H1 (GpSGHV-Uga). Haplotypes H1, H5 and H14 formed three main genealogy hubs, potentially representing the ancestors of the 15 haplotypes.

CONCLUSION

These data identify G. pallidipes as a significant driver for the generation and diversity of GpSGHV variants. This information may provide control guidance when new tsetse colonies are established and hence, for improved management of the virus in tsetse rearing facilities that maintain multiple Glossina species.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

GtTR: Bayesian estimation of absolute tandem repeat copy number using sequence capture and high throughput sequencing

BACKGROUND

Tandem repeats comprise significant proportion of the human genome including coding and regulatory regions. They are highly prone to repeat number variation and nucleotide mutation due to their repetitive and unstable nature, making them a major source of genomic variation between individuals. Despite recent advances in high throughput sequencing, analysis of tandem repeats in the context of complex diseases is still hindered by technical limitations. We report a novel targeted sequencing approach, which allows simultaneous analysis of hundreds of repeats. We developed a Bayesian algorithm, namely - GtTR - which combines information from a reference long-read dataset with a short read counting approach to genotype tandem repeats at population scale. PCR sizing analysis was used for validation.

RESULTS

We used a PacBio long-read sequenced sample to generate a reference tandem repeat genotype dataset with on average 13% absolute deviation from PCR sizing results. Using this reference dataset GtTR generated estimates of VNTR copy number with accuracy within 95% high posterior density (HPD) intervals of 68 and 83% for capture sequence data and 200X WGS data respectively, improving to 87 and 94% with use of a PCR reference. We show that the genotype resolution increases as a function of depth, such that the median 95% HPD interval lies within 25, 14, 12 and 8% of the its midpoint copy number value for 30X, 200X WGS, 395X and 800X capture sequence data respectively. We validated nine targets by PCR sizing analysis and genotype estimates from sequencing results correlated well with PCR results.

CONCLUSIONS

The novel genotyping approach described here presents a new cost-effective method to explore previously unrecognized class of repeat variation in GWAS studies of complex diseases at the population level. Further improvements in accuracy can be obtained by improving accuracy of the reference dataset.

Selection of the sex-linked inhibitor of apoptosis in mountain pine beetle (Dendroctonus ponderosae) driven by enhanced expression during early overwintering

The mountain pine beetle (Dendroctonus ponderosae) is an insect native to western North America; however, its geographical range has recently expanded north in BC and east into Alberta. To understand the population structure in the areas of expansion, 16 gene-linked microsatellites were screened and compared to neutral microsatellites using outlier analyses of F_st and F_ct values. One sex-linked gene, inhibitor of apoptosis (IAP), showed a strong signature of positive selection for neo-X alleles and was analyzed for evidence of adaptive variation. Alleles of IAP were sequenced, and differences between the neo-X and neo-Y alleles were consistent with neutral evolution suggesting that the neo-Y allele may not be under functional constraints. Neo-Y alleles were amplified from gDNA, but not effectively from cDNA, suggesting that there was little IAP expression from neo-Y alleles. There were no differences in overall IAP expression between males and females with the common northern neo-X allele suggesting that the neo-X allele in males compensates for the reduced expression of neo-Y alleles. However, males lacking the most common northern neo-X allele thought to be selected for in northern populations had reduced overall IAP expression in early October-at a time when beetles are preparing for overwintering. This suggests that the most common allele may have more rapid upregulation. The reduced function of neo-Y alleles of IAP suggested by both sequence differences and lower levels of expression may foster a highly selective environment for neo-X alleles such as the common northern allele with more efficient upregulation.

Biological Roles of Protein-Coding Tandem Repeats in the Yeast Candida Albicans

Tandem repeat (TR) DNA mutates faster than other DNA by insertion and deletion of repeats. Large parts of eukaryotic proteomes are encoded by ORFs containing protein-coding TRs (TR-ORFs, pcTRs) with largely unknown biological consequences. We explored these in the yeast Candida albicans, an opportunistic human pathogen. We found that almost half of C. albicans’ proteins are encoded by TR-ORFs. pcTR frequency differed only moderately between different gene (GO) categories. Bioinformatic predictions of genome-wide mutation rates and clade-specific differences in pcTR allele frequencies indicated that pcTRs (i) significantly increase the genome-wide mutation rate; (ii) significantly impact on fitness and (iii) allow the evolution of selectively advantageous clade-specific protein variants. Synonymous mutations reduced the repetitiveness of many amino acid repeat-encoding pcTRs. A survey, in 58 strains, revealed that in some pcTR regions in which repetitiveness was not significantly diminished by synonymous mutations the habitat predicted which alleles were present, suggesting roles of pcTR mutation in short-term adaptation and pathogenesis. In C. albicans pcTR mutation apparently is an important mechanism for mutational advance and possibly also rapid adaptation, with synonymous mutations providing a mechanism for adjusting mutation rates of individual pcTRs. Analyses of Arabidopsis and human pcTRs showed that the latter also occurs in other eukaryotes.

The Complete Chloroplast Genome of Heimia myrtifolia and Comparative Analysis within Myrtales

Heimia myrtifolia is an important medicinal plant with several pharmacologically active alkaloids and is also used as an ornamental landscape plant. The purpose of this study is to complete and characterize the chloroplast (cp) genome of H. myrtifolia and compare genomic features to other Myrtales species' cp genomes. The analysis showed that H. myrtifolia has a total length of 159,219 bp with a typical quadripartite structure containing two identical inverted repeats (IRs) of 25,643 bp isolated by one large single copy (LSC) of 88,571 bp and one small single copy (SSC) of 18,822 bp. The H. myrtifolia cp genome contains 129 genes with eight ribosomal RNAs, 30 transfer RNAs, and 78 protein coding genes, in which 17 genes are duplicated in two IR regions. The genome organization including gene type and number and guanine-cytosine (GC) content is analyzed among the 12 cp genomes in this study. Approximately 255 simple sequence repeats (SSRs) and 16 forward, two reverses, and two palindromic repeats were identified in the H. myrtifolia cp genome. By comparing the whole H. myrtifolia cp genome with 11 other Myrtales species, the results showed that the sequence similarity was high between coding regions while sequence divergence was high between intergenic regions. By employing the full cp genomes for phylogenetic analysis, structural and sequence differences were characterized between H. myrtifolia and 11 Myrtales species illustrating what patterns are common in the evolution of cp genomes within the Myrtales. The first entire cp genome in the genus Heimia provides a valuable resource for further studies in these medicinally and ornamentally important taxa.

The Complete Chloroplast Genome of Catha edulis: A Comparative Analysis of Genome Features with Related Species

Qat (Catha edulis, Celastraceae) is a woody evergreen species with great economic and cultural importance. It is cultivated for its stimulant alkaloids cathine and cathinone in East Africa and southwest Arabia. However, genome information, especially DNA sequence resources, for C. edulis are limited, hindering studies regarding interspecific and intraspecific relationships. Herein, the complete chloroplast (cp) genome of Catha edulis is reported. This genome is 157,960 bp in length with 37% GC content and is structurally arranged into two 26,577 bp inverted repeats and two single-copy areas. The size of the small single-copy and the large single-copy regions were 18,491 bp and 86,315 bp, respectively. The C. edulis cp genome consists of 129 coding genes including 37 transfer RNA (tRNA) genes, 8 ribosomal RNA (rRNA) genes, and 84 protein coding genes. For those genes, 112 are single copy genes and 17 genes are duplicated in two inverted regions with seven tRNAs, four rRNAs, and six protein coding genes. The phylogenetic relationships resolved from the cp genome of qat and 32 other species confirms the monophyly of Celastraceae. The cp genomes of C. edulis, Euonymus japonicus and seven Celastraceae species lack the rps16 intron, which indicates an intron loss took place among an ancestor of this family. The cp genome of C. edulis provides a highly valuable genetic resource for further phylogenomic research, barcoding and cp transformation in Celastraceae.

In silico characterization of tandem repeats in Trichophyton rubrum and related dermatophytes provides new insights into their role in pathogenesis

Trichophyton rubrum is the most common etiological agent of dermatophytoses worldwide, which is able to degrade keratinized tissues. The sequencing of the genome of different dermatophyte species has provided a large amount of data, including tandem repeats that may play a role in genetic variability and in the pathogenesis of these fungi. Tandem repeats are adjacent DNA sequences of 2–200 nucleotides in length, which exert regulatory and adaptive functions. These repetitive DNA sequences are found in different classes of fungal proteins, especially those involved in cell adhesion, a determinant factor for the establishment of fungal infection. The objective of this study was to develop a Dermatophyte Tandem Repeat Database (DTRDB) for the storage and identification of tandem repeats in T. rubrum and six other dermatophyte species. The current version of the database contains 35 577 tandem repeats detected in 16 173 coding sequences. The repeats can be searched using entry parameters such as repeat unit length (nt—nucleotide), repeat number, variability score, and repeat sequence motif. These data were used to study the relative frequency and distribution of repeats in the sequences, as well as their possible functions in dermatophytes. A search of the database revealed that these repeats occur in 22–33% of genes transcribed in dermatophytes where they could be involved in the success of adaptation to the host tissue and establishment of infection. The repeats were detected in transcripts that are mainly related to three biological processes: regulation, adhesion, and metabolism. The database developed enables users to identify and analyse tandem repeat regions in target genes related to pathogenicity and fungal–host interactions in dermatophytes and may contribute to the discovery of new targets for the development of antifungal agents.

Database URL:http://comp.mch.ifsuldeminas.edu.br/dtrdb/

Comparative analysis of low complexity regions in Plasmodia

Low complexity regions (LCRs) are a common feature shared by many genomes, but their evolutionary and functional significance remains mostly unknown. At the core of the uncertainty is a poor understanding of the mechanisms that regulate their retention in genomes, whether driven by natural selection or neutral evolution. Applying a comparative approach of LCRs to multiple strains and species is a powerful approach to identify patterns of conservation in these regions. Using this method, we investigate the evolutionary history of LCRs in the genus Plasmodium based on orthologous protein coding genes shared by 11 species and strains from primate and rodent-infecting pathogens. We find multiple lines of evidence in support of natural selection as a major evolutionary force shaping the composition and conservation of LCRs through time and signatures that their evolutionary paths are species specific. Our findings add a comparative analysis perspective to the debate on the evolution of LCRs and harness the power of sequence comparisons to identify potential functionally important LCR candidates.

Optimization of PCR conditions for amplifying an AT-rich amino acid transporter promoter sequence with high number of tandem repeats from Arabidopsis thaliana

OBJECTIVE

The aim of the present study is to optimize the PCR conditions required to amplify the promoter sequence of an amino acid transporter having an AT-rich base composition with a high number of tandem repeats.

RESULT

Results show that successful amplification can be achieved by performing a 2-step PCR at a lower extension temperature of 65 °C for an increased extension period of 1.5 min/kb, with MgCl₂ concentration ranging from 2.5 to 3.0 mM. The results also suggest that the DNA concentration of about 25-30 ng/µl was essential to achieve this amplification.

Distinct patterns of simple sequence repeats and GC distribution in intragenic and intergenic regions of primate genomes

As the first systematic examination of simple sequence repeats (SSRs) and guanine-cytosine (GC) distribution in intragenic and intergenic regions of ten primates, our study showed that SSRs and GC displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation. Our results suggest that the majority of SSRs are distributed in non-coding regions, such as the introns, TEs, and intergenic regions. In these primates, trinucleotide perfect (P) SSRs were the most abundant repeats type in the 5'UTRs and CDSs, whereas, mononucleotide P-SSRs were the most in the intron, 3'UTRs, TEs, and intergenic regions. The GC-contents varied greatly among different intragenic and intergenic regions: 5'UTRs > CDSs > 3'UTRs > TEs > introns > intergenic regions, and high GC-content was frequently distributed in exon-rich regions. Our results also showed that in the same intragenic and intergenic regions, the distribution of GC-contents were great similarity in the different primates. Tri- and hexanucleotide P-SSRs had the most GC-contents in the 5'UTRs and CDSs, whereas mononucleotide P-SSRs had the least GC-contents in the six genomic regions of these primates. The most frequent motifs for different length varied obviously with the different genomic regions.

Signatures of selection in mammalian clock genes with coding trinucleotide repeats: Implications for studying the genomics of high-pace adaptation

Climate change is predicted to affect the reproductive ecology of wildlife; however, we have yet to understand if and how species can adapt to the rapid pace of change. Clock genes are functional genes likely critical for adaptation to shifting seasonal conditions through shifts in timing cues. Many of these genes contain coding trinucleotide repeats, which offer the potential for higher rates of change than single nucleotide polymorphisms (SNPs) at coding sites, and, thus, may translate to faster rates of adaptation in changing environments. We characterized repeats in 22 clock genes across all annotated mammal species and evaluated the potential for selection on repeat motifs in three clock genes (NR1D1,CLOCK, and PER1) in three congeneric species pairs with different latitudinal range limits: Canada lynx and bobcat (Lynx canadensis and L. rufus), northern and southern flying squirrels (Glaucomys sabrinus and G. volans), and white‐footed and deer mouse (Peromyscus leucopus and P. maniculatus). Signatures of positive selection were found in both the interspecific comparison of Canada lynx and bobcat, and intraspecific analyses in Canada lynx. Northern and southern flying squirrels showed differing frequencies at common CLOCK alleles and a signature of balancing selection. Regional excess homozygosity was found in the deer mouse at PER1 suggesting disruptive selection, and further analyses suggested balancing selection in the white‐footed mouse. These preliminary signatures of selection and the presence of trinucleotide repeats within many clock genes warrant further consideration of the importance of candidate gene motifs for adaptation to climate change.

Fish introductions in the former Soviet Union: The Sevan trout (Salmo ischchan) - 80 years later

The Soviet Union played the leading role in fish introductions in Eurasia. However, only 3% of all introductions prior to 1978 gave a commercial benefit. One of the noteworthy examples appears to be the Sevan trout (Salmo ischchan Kessler, 1877)-an endemic salmonid of Lake Sevan in Armenia. This species has been introduced to Kirghizstan, Kazakhstan, and Uzbekistan, however, only the Kirghiz population has persisted in relatively high numbers. In this paper we provide the first extensive molecular study of S. ischchan using samples from the native population from Lake Sevan and three hatcheries in Armenia, as well as from the population introduced to Lake Issyk Kul in Kirghizstan. The Kirghiz population has been isolated since the introductions took place in 1930 and 1936. Our results, based on 11 nuclear microsatellites and a 905 bp fragment of the mitochondrial control region suggest that hatcheries have maintained genetic variability by way of ongoing translocations of individuals from Lake Sevan. Simultaneously, significant Garza-Williamson M-values suggest that bottlenecks could have reduced the genetic variability of the wild populations in the past. This hypothesis is supported by historical data, indicating highly manipulated water-level regulations and poaching as two main factors that dramatically impact fish abundance in the lake. On the other hand, a similar situation has been observed in Kirghizstan, but this population likely rebounded from small population size faster than the other populations examined. The Kirghiz population is significantly genetically differentiated from the other groups and have morphological features and biological attributes not observed in the source population. Genetic data imply that the effective population size in the native population is lower than that found in the introduced population, suggesting that some active protection of the Lake Sevan population may be needed urgently.

Genetic variations at the human growth hormone receptor (GHR) gene locus are associated with idiopathic short stature

GH plays an essential role in the growing child by binding to the growth hormone receptor (GHR) on target cells and regulating multiple growth promoting and metabolic effects. Mutations in the GHR gene coding regions result in GH insensitivity (dwarfism) due to a dysfunctional receptor protein. However, children with idiopathic short stature (ISS) show growth impairment without GH or GHR defects. We hypothesized that decreased expression of the GHR gene may be involved. To test this, we investigated whether common genetic variants (microsatellites, SNPs) in regulatory regions of the GHR gene region were associated with the ISS phenotype. Genotyping of a GT‐repeat microsatellite in the GHR 5′UTR in a Montreal ISS cohort (n = 37 ISS, n = 105 controls) revealed that the incidence of the long/short (L/S) genotype was 3.3× higher in ISS children than controls (P = 0.04, OR = 3.85). In an Italian replication cohort (n = 143 ISS, n = 282 controls), the medium/short (M/S) genotype was 1.9× more frequent in the male ISS than controls (P = 0.017, OR = 2.26). In both ISS cohorts, logistic regression analysis of 27 SNPs showed an association of ISS with rs4292454, while haplotype analysis revealed specific risk haplotypes in the 3′ haploblocks. In contrast, there were no differences in GT genotype frequencies in a cohort of short stature (SS) adults versus controls (CARTaGENE: n = 168 SS, n = 207 controls) and the risk haplotype in the SS cohort was located in the most 5′ haploblock. These data suggest that the variants identified are potentially genetic markers specifically associated with the ISS phenotype.

Repetitive Elements in Mycoplasma hyopneumoniae Transcriptional Regulation

Transcriptional regulation, a multiple-step process, is still poorly understood in the important pig pathogen Mycoplasma hyopneumoniae. Basic motifs like promoters and terminators have already been described, but no other cis-regulatory elements have been found. DNA repeat sequences have been shown to be an interesting potential source of cis-regulatory elements. In this work, a genome-wide search for tandem and palindromic repetitive elements was performed in the intergenic regions of all coding sequences from M. hyopneumoniae strain 7448. Computational analysis demonstrated the presence of 144 tandem repeats and 1,171 palindromic elements. The DNA repeat sequences were distributed within the 5' upstream regions of 86% of transcriptional units of M. hyopneumoniae strain 7448. Comparative analysis between distinct repetitive sequences found in related mycoplasma genomes demonstrated different percentages of conservation among pathogenic and nonpathogenic strains. qPCR assays revealed differential expression among genes showing variable numbers of repetitive elements. In addition, repeats found in 206 genes already described to be differentially regulated under different culture conditions of M. hyopneumoniae strain 232 showed almost 80% conservation in relation to M. hyopneumoniae strain 7448 repeats. Altogether, these findings suggest a potential regulatory role of tandem and palindromic DNA repeats in the M. hyopneumoniae transcriptional profile.

Expansion of Lysine-rich Repeats in Plasmodium Proteins Generates Novel Localization Sequences That Target the Periphery of the Host Erythrocyte

Repetitive low complexity sequences, mostly assumed to have no function, are common in proteins that are exported by the malaria parasite into its host erythrocyte. We identify a group of exported proteins containing short lysine-rich tandemly repeated sequences that are sufficient to localize to the erythrocyte periphery, where key virulence-related modifications to the plasma membrane and the underlying cytoskeleton are known to occur. Efficiency of targeting is dependent on repeat number, indicating that novel targeting modules could evolve by expansion of short lysine-rich sequences. Indeed, analysis of fragments of GARP from different species shows that two novel targeting sequences have arisen via the process of repeat expansion in this protein. In the protein Hyp12, the targeting function of a lysine-rich sequence is masked by a neighboring repetitive acidic sequence, further highlighting the importance of repetitive low complexity sequences. We show that sequences capable of targeting the erythrocyte periphery are present in at least nine proteins from Plasmodium falciparum and one from Plasmodium knowlesi. We find these sequences in proteins known to be involved in erythrocyte rigidification and cytoadhesion as well as in previously uncharacterized exported proteins. Together, these data suggest that expansion and contraction of lysine-rich repeats could generate targeting sequences de novo as well as modulate protein targeting efficiency and function in response to selective pressure.

Recent evolution of the salivary mucin MUC7

Genomic structural variants constitute the majority of variable base pairs in primate genomes and affect gene function in multiple ways. While whole gene duplications and deletions are relatively well-studied, the biology of subexonic (i.e., within coding exon sequences), copy number variation remains elusive. The salivary MUC7 gene provides an opportunity for studying such variation, as it harbors copy number variable subexonic repeat sequences that encode for densely O-glycosylated domains (PTS-repeats) with microbe-binding properties. To understand the evolution of this gene, we analyzed mammalian and primate genomes within a comparative framework. Our analyses revealed that (i) MUC7 has emerged in the placental mammal ancestor and rapidly gained multiple sites for O-glycosylation; (ii) MUC7 has retained its extracellular activity in saliva in placental mammals; (iii) the anti-fungal domain of the protein was remodified under positive selection in the primate lineage; and (iv) MUC7 PTS-repeats have evolved recurrently and under adaptive constraints. Our results establish MUC7 as a major player in salivary adaptation, likely as a response to diverse pathogenic exposure in primates. On a broader scale, our study highlights variable subexonic repeats as a primary source for modular evolutionary innovation that lead to rapid functional adaptation.

Microsatellite markers: what they mean and why they are so useful

Microsatellites or Single Sequence Repeats (SSRs) are extensively employed in plant genetics studies, using both low and high throughput genotyping approaches. Motivated by the importance of these sequences over the last decades this review aims to address some theoretical aspects of SSRs, including definition, characterization and biological function. The methodologies for the development of SSR loci, genotyping and their applications as molecular markers are also reviewed. Finally, two data surveys are presented. The first was conducted using the main database of Web of Science, prospecting for articles published over the period from 2010 to 2015, resulting in approximately 930 records. The second survey was focused on papers that aimed at SSR marker development, published in the American Journal of Botany's Primer Notes and Protocols in Plant Sciences (over 2013 up to 2015), resulting in a total of 87 publications. This scenario confirms the current relevance of SSRs and indicates their continuous utilization in plant science.

Are all repeats created equal? Understanding DNA repeats at an individual level

Repetitive DNA sequences, comprising up to 50 % of the genome in all eukaryotes, play important roles in a wide range of cellular functions, such as transcriptional regulation, genome stability, and cellular differentiation. However, due to technical difficulties in differentiating their sequences, DNA repeats remain one of the most mysterious parts of eukaryotic genomes. Key questions, such as how repetitive entities behave at individual level and how the internal architecture of these repeats is organized, are still poorly understood. Recent advances from our group reveal unexpected position-dependent variation within tandem DNA repeats in fission yeast. Despite sharing identical DNA sequences, the peri-centromeric repeats are organized into diverse epigenetic states and chromatin structures. We demonstrate that this position-dependent variation requires key heterochromatin factors and condensin. Our works further suggest that the peri-centromeric repeats are organized into distinct higher order structures that ensure a proper positioning of CENP-A, the centromere-specific histone H3 variant, to centromeres. These most recent developments offer insights into the mechanisms underlying the position effect within tandem DNA arrays, and have broad implications in the field of epigenetics and chromatin biology.

The Complete Plastid Genome of Lagerstroemia fauriei and Loss of rpl2 Intron from Lagerstroemia (Lythraceae)

Lagerstroemia (crape myrtle) is an important plant genus used in ornamental horticulture in temperate regions worldwide. As such, numerous hybrids have been developed. However, DNA sequence resources and genome information for Lagerstroemia are limited, hindering evolutionary inferences regarding interspecific relationships. We report the complete plastid genome of Lagerstroemia fauriei. To our knowledge, this is the first reported whole plastid genome within Lythraceae. This genome is 152,440 bp in length with 38% GC content and consists of two single-copy regions separated by a pair of 25,793 bp inverted repeats. The large single copy and the small single copy regions span 83,921 bp and 16,933 bp, respectively. The genome contains 129 genes, including 17 located in each inverted repeat. Phylogenetic analysis of genera sampled from Geraniaceae, Myrtaceae, and Onagraceae corroborated the sister relationship between Lythraceae and Onagraceae. The plastid genomes of L. fauriei and several other Lythraceae species lack the rpl2 intron, which indicating an early loss of this intron within the Lythraceae lineage. The plastid genome of L. fauriei provides a much needed genetic resource for further phylogenetic research in Lagerstroemia and Lythraceae. Highly variable markers were identified for application in phylogenetic, barcoding and conservation genetic applications.

Profiles of low complexity regions in Apicomplexa

BACKGROUND

Low complexity regions (LCRs) are a ubiquitous feature in genomes and yet their evolutionary history and functional roles are unclear. Previous studies have shown contrasting evidence in favor of both neutral and selective mechanisms of evolution for different sets of LCRs suggesting that modes of identification of these regions may play a role in our ability to discern their evolutionary history. To further investigate this issue, we used a multiple threshold approach to identify species-specific profiles of proteome complexity and, by comparing properties of these sets, determine the influence that starting parameters have on evolutionary inferences.

RESULTS

We find that, although qualitatively similar, quantitatively each species has a unique LCR profile which represents the frequency of these regions within each genome. Inferences based on these profiles are more accurate in comparative analyses of genome complexity as they allow to determine the relative complexity of multiple genomes as well as the type of repetitiveness that is most common in each. Based on the multiple threshold LCR sets obtained, we identified predominant evolutionary mechanisms at different complexity levels, which show neutral mechanisms acting on highly repetitive LCRs (e.g., homopolymers) and selective forces becoming more important as heterogeneity of the LCRs increases.

CONCLUSIONS

Our results show how inferences based on LCRs are influenced by the parameters used to identify these regions. Sets of LCRs are heterogeneous aggregates of regions that include homo- and heteropolymers and, as such, evolve according to different mechanisms. LCR profiles provide a new way to investigate genome complexity across species and to determine the driving mechanism of their evolution.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10⁻⁸) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left–right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

A test of somatic mosaicism in the androgen receptor gene of Canada lynx (Lynx canadensis)

Background

The androgen receptor, an X-linked gene, has been widely studied in human populations because it contains highly polymorphic trinucleotide repeat motifs that have been associated with a number of adverse human health and behavioral effects. A previous study on the androgen receptor gene in carnivores reported somatic mosaicism in the tissues of a number of species including Eurasian lynx (Lynx lynx). We investigated this claim in a closely related species, Canada lynx (Lynx canadensis). The presence of somatic mosaicism in lynx tissues could have implications for the future study of exonic trinucleotide repeats in landscape genomic studies, in which the accurate reporting of genotypes would be highly problematic.

Methods

To determine whether mosaicism occurs in Canada lynx, two lynx individuals were sampled for a variety of tissue types (lynx 1) and tissue locations (lynx 1 and 2), and 1,672 individuals of known sex were genotyped to further rule out mosaicism.

Results

We found no evidence of mosaicism in tissues from the two necropsied individuals, or any of our genotyped samples.

Conclusions

Our results indicate that mosaicism does not manifest in Canada lynx. Therefore, the use of hide samples for further work involving trinucleotide repeat polymorphisms in Canada lynx is warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0284-y) contains supplementary material, which is available to authorized users.

A metalloprotease secreted by the type II secretion system links Vibrio cholerae with collagen

Vibrio cholerae is autochthonous to various aquatic niches and is the etiological agent of the life-threatening diarrheal disease cholera. The persistence of V. cholerae in natural habitats is a crucial factor in the epidemiology of cholera. In contrast to the well-studied V. cholerae-chitin connection, scarce information is available about the factors employed by the bacteria for the interaction with collagens. Collagens might serve as biologically relevant substrates, because they are the most abundant protein constituents of metazoan tissues and V. cholerae has been identified in association with invertebrate and vertebrate marine animals, as well as in a benthic zone of the ocean where organic matter, including collagens, accumulates. Here, we describe the characterization of the V. cholerae putative collagenase, VchC, encoded by open reading frame VC1650 and belonging to the subfamily M9A peptidases. Our studies demonstrate that VchC is an extracellular collagenase degrading native type I collagen of fish and mammalian origin. Alteration of the predicted catalytic residues coordinating zinc ions completely abolished the protein enzymatic activity but did not affect the translocation of the protease by the type II secretion pathway into the extracellular milieu. We also show that the protease undergoes a maturation process with the aid of a secreted factor(s). Finally, we propose that V. cholerae is a collagenovorous bacterium, as it is able to utilize collagen as a sole nutrient source. This study initiates new lines of investigations aiming to uncover the structural and functional components of the V. cholerae collagen utilization program.

The peculiar landscape of repetitive sequences in the olive (Olea europaea L.) genome

Analyzing genome structure in different species allows to gain an insight into the evolution of plant genome size. Olive (Olea europaea L.) has a medium-sized haploid genome of 1.4 Gb, whose structure is largely uncharacterized, despite the growing importance of this tree as oil crop. Next-generation sequencing technologies and different computational procedures have been used to study the composition of the olive genome and its repetitive fraction. A total of 2.03 and 2.3 genome equivalents of Illumina and 454 reads from genomic DNA, respectively, were assembled following different procedures, which produced more than 200,000 differently redundant contigs, with mean length higher than 1,000 nt. Mapping Illumina reads onto the assembled sequences was used to estimate their redundancy. The genome data set was subdivided into highly and medium redundant and nonredundant contigs. By combining identification and mapping of repeated sequences, it was established that tandem repeats represent a very large portion of the olive genome (∼31% of the whole genome), consisting of six main families of different length, two of which were first discovered in these experiments. The other large redundant class in the olive genome is represented by transposable elements (especially long terminal repeat-retrotransposons). On the whole, the results of our analyses show the peculiar landscape of the olive genome, related to the massive amplification of tandem repeats, more than that reported for any other sequenced plant genome.

MIPSTR: a method for multiplex genotyping of germline and somatic STR variation across many individuals

Short tandem repeats (STRs) are highly mutable genetic elements that often reside in regulatory and coding DNA. The cumulative evidence of genetic studies on individual STRs suggests that STR variation profoundly affects phenotype and contributes to trait heritability. Despite recent advances in sequencing technology, STR variation has remained largely inaccessible across many individuals compared to single nucleotide variation or copy number variation. STR genotyping with short-read sequence data is confounded by (1) the difficulty of uniquely mapping short, low-complexity reads; and (2) the high rate of STR amplification stutter. Here, we present MIPSTR, a robust, scalable, and affordable method that addresses these challenges. MIPSTR uses targeted capture of STR loci by single-molecule Molecular Inversion Probes (smMIPs) and a unique mapping strategy. Targeted capture and our mapping strategy resolve the first challenge; the use of single molecule information resolves the second challenge. Unlike previous methods, MIPSTR is capable of distinguishing technical error due to amplification stutter from somatic STR mutations. In proof-of-principle experiments, we use MIPSTR to determine germline STR genotypes for 102 STR loci with high accuracy across diverse populations of the plant A. thaliana. We show that putatively functional STRs may be identified by deviation from predicted STR variation and by association with quantitative phenotypes. Using DNA mixing experiments and a mutant deficient in DNA repair, we demonstrate that MIPSTR can detect low-frequency somatic STR variants. MIPSTR is applicable to any organism with a high-quality reference genome and is scalable to genotyping many thousands of STR loci in thousands of individuals.

Dominant short repeated sequences in bacterial genomes

We use a novel multidimensional searching approach to present the first exhaustive search for all possible repeated sequences in 166 genomes selected to cover the bacterial domain. We found an overrepresentation of repeated sequences in all but one of the genomes. The most prevalent repeats by far were related to interspaced short palindromic repeats (CRISPRs)—conferring bacterial adaptive immunity. We identified a deep branching clade of thermophilic Firmicutes containing the highest number of CRISPR repeats. We also identified a high prevalence of tandem repeated heptamers. In addition, we identified GC-rich repeats that could potentially be involved in recombination events. Finally, we identified repeats in a 16322 amino acid mega protein (involved in biofilm formation) and inverted repeats flanking miniature transposable elements (MITEs). In conclusion, the exhaustive search for repeated sequences identified new elements and distribution of these, which has implications for understanding both the ecology and evolution of bacteria.

Evolution of invasion in a diverse set of Fusobacterium species

The diverse Fusobacterium genus contains species implicated in multiple clinical pathologies, including periodontal disease, preterm birth, and colorectal cancer. The lack of genetic tools for manipulating these organisms leaves us with little understanding of the genes responsible for adherence to and invasion of host cells. Actively invading Fusobacterium species can enter host cells independently, whereas passively invading species need additional factors, such as compromise of mucosal integrity or coinfection with other microbes. We applied whole-genome sequencing and comparative analysis to study the evolution of active and passive invasion strategies and to infer factors associated with active forms of host cell invasion. The evolution of active invasion appears to have followed an adaptive radiation in which two of the three fusobacterial lineages acquired new genes and underwent expansions of ancestral genes that enable active forms of host cell invasion. Compared to passive invaders, active invaders have much larger genomes, encode FadA-related adhesins, and possess twice as many genes encoding membrane-related proteins, including a large expansion of surface-associated proteins containing the MORN2 domain of unknown function. We predict a role for proteins containing MORN2 domains in adhesion and active invasion. In the largest and most comprehensive comparison of sequenced Fusobacterium species to date, we have generated a testable model for the molecular pathogenesis of Fusobacterium infection and illuminate new therapeutic or diagnostic strategies.

Fusobacterium species have recently been implicated in a broad spectrum of human pathologies, including Crohn’s disease, ulcerative colitis, preterm birth, and colorectal cancer. Largely due to the genetic intractability of member species, the mechanisms by which Fusobacterium causes these pathologies are not well understood, although adherence to and active invasion of host cells appear important. We examined whole-genome sequence data from a diverse set of Fusobacterium species to identify genetic determinants of active forms of host cell invasion. Our analyses revealed that actively invading Fusobacterium species have larger genomes than passively invading species and possess a specific complement of genes—including a class of genes of unknown function that we predict evolved to enable host cell adherence and invasion. This study provides an important framework for future studies on the role of Fusobacterium in pathologies such as colorectal cancer.

The overdue promise of short tandem repeat variation for heritability

Short tandem repeat (STR) variation has been proposed as a major explanatory factor in the heritability of complex traits in humans and model organisms. However, we still struggle to incorporate STR variation into genotype-phenotype maps. We review here the promise of STRs in contributing to complex trait heritability and highlight the challenges that STRs pose due to their repetitive nature. We argue that STR variants are more likely than single-nucleotide variants to have epistatic interactions, reiterate the need for targeted assays to genotype STRs accurately, and call for more appropriate statistical methods in detecting STR-phenotype associations. Lastly, we suggest that somatic STR variation within individuals may serve as a read-out of disease susceptibility, and is thus potentially a valuable covariate for future association studies.

Genome-wide analysis of tandem repeats in plants and green algae

Tandem repeats (TRs) extensively exist in the genomes of prokaryotes and eukaryotes. Based on the sequenced genomes and gene annotations of 31 plant and algal species in Phytozome version 8.0 (http://www.phytozome.net/), we examined TRs in a genome-wide scale, characterized their distributions and motif features, and explored their putative biological functions. Among the 31 species, no significant correlation was detected between the TR density and genome size. Interestingly, green alga Chlamydomonas reinhardtii (42,059 bp/Mbp) and castor bean Ricinus communis (55,454 bp/Mbp) showed much higher TR densities than all other species (13,209 bp/Mbp on average). In the 29 land plants, including 22 dicots, 5 monocots, and 2 bryophytes, 5′-UTR and upstream intergenic 200-nt (UI200) regions had the first and second highest TR densities, whereas in the two green algae (C. reinhardtii and Volvox carteri) the first and second highest densities were found in intron and coding sequence (CDS) regions, respectively. In CDS regions, trinucleotide and hexanucleotide motifs were those most frequently represented in all species. In intron regions, especially in the two green algae, significantly more TRs were detected near the intron–exon junctions. Within intergenic regions in dicots and monocots, more TRs were found near both the 5′ and 3′ ends of genes. GO annotation in two green algae revealed that the genes with TRs in introns are significantly involved in transcriptional and translational processing. As the first systematic examination of TRs in plant and green algal genomes, our study showed that TRs displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation in plants and green algae.

PCSK6 VNTR Polymorphism Is Associated with Degree of Handedness but Not Direction of Handedness

Although the left and right human cerebral hemispheres differ both functionally and anatomically, the mechanisms that underlie the establishment of these hemispheric specializations, as well as their physiological and behavioral implications, remain largely unknown. Since cerebral asymmetry is strongly correlated with handedness, and handedness is assumed to be influenced by a number of genetic and environmental factors, we performed an association study of LRRTM1 rs6733871 and a number of polymorphisms in PCSK6 and different aspects of handedness assessed with the Edinburgh handedness inventory in a sample of unrelated healthy adults (n = 1113). An intronic 33bp variable-number tandem repeat (VNTR) polymorphism in PCSK6 (rs10523972) shows a significant association (significance threshold: p<0.0025, adjusted for multiple comparisons) with a handedness category comparison (P = 0.0005) and degree of handedness (P = 0.001). These results provide further evidence for the role of PCSK6 as candidate for involvement in the biological mechanisms that underlie the establishment of normal brain lateralization and thus handedness and support the assumption that the degree of handedness, instead the direction, may be the more appropriate indicator of cerebral organization.