Microsatellite markers for linkage and association studies

Rediscovering tandem repeat variation in schizophrenia: challenges and opportunities

Tandem repeats (TRs) are prevalent throughout the genome, constituting at least 3% of the genome, and often highly polymorphic. The high mutation rate of TRs, which can be orders of magnitude higher than single-nucleotide polymorphisms and indels, indicates that they are likely to make significant contributions to phenotypic variation, yet their contribution to schizophrenia has been largely ignored by recent genome-wide association studies (GWAS). Tandem repeat expansions are already known causative factors for over 50 disorders, while common tandem repeat variation is increasingly being identified as significantly associated with complex disease and gene regulation. The current review summarizes key background concepts of tandem repeat variation as pertains to disease risk, elucidating their potential for schizophrenia association. An overview of next-generation sequencing-based methods that may be applied for TR genome-wide identification is provided, and some key methodological challenges in TR analyses are delineated.

Genetic Diversity and Structure Revealed by Genomic Microsatellite Markers of Mytilus unguiculatus in the Coast of China Sea

The hard-shelled mussel Mytilus unguiculatus plays an important role in mussel aquaculture in China due to its characteristic and nutritive value. In this study, 10 microsatellite loci are used to study the genetic diversity and genetic structure of seven location populations of M. unguiculatus in coastal areas of China. The results of amplification and genotyping show that the observed heterozygosity (H_o) and the expected heterozygosity (H_e) are 0.61~0.71 and 0.72~0.83, respectively. M. unguiculatus has high genetic diversity. The inbreeding index (F_IS) of M. unguiculatus is significantly positive (F_IS: 0.14~0.19), indicating that inbreeding might exist within populations. The genetic structure of M. unguiculatus is weak within populations from the East China Sea All results showed that genetic differences existed between the Qingdao population from the Yellow Sea and other populations from the East China Sea. It does not detect a population bottleneck event or expansion event in the populations. The results from this study can be used to provide important insights in genetic management units and sustainable utilization of M. unguiculatus resources and provide a better understand of genetic structure of marine bivalve with similar planktonic larval stage in the China Sea.

LT-RPA: An Isothermal DNA Amplification Approach for Improved Microsatellite Genotyping and Microsatellite Instability Detection

Microsatellites are short tandem repeats of one to six nucleotides that are highly polymorphic and extensively used as genetic markers in numerous biomedical applications, including the detection of microsatellite instability (MSI) in cancer. The standard analytical method for microsatellite analysis relies on PCR amplification followed by capillary electrophoresis or, more recently, next-generation sequencing (NGS). However, their amplification during PCR generates undesirable frameshift products known as stutter peaks caused by polymerase slippage, complicating data analysis and interpretation, while very few alternative methods for microsatellite amplification have been developed to reduce the formation of these artifacts. In this context, the recently developed low-temperature recombinase polymerase amplification (LT-RPA) is an isothermal DNA amplification method at low temperature (32 °C) that drastically reduces and sometimes completely abolishes the formation of stutter peaks. LT-RPA greatly simplifies the genotyping of microsatellites and improves the detection of MSI in cancer. In this chapter, we describe in detail all the experimental steps necessary for the development of LT-RPA simplex and multiplex assays for microsatellite genotyping and MSI detection, including the design, optimization, and validation of the assays combined with capillary electrophoresis or NGS.

Genetic and Epigenetic Interplay Define Disease Onset and Severity in Repeat Diseases

Repeat diseases, such as fragile X syndrome, myotonic dystrophy, Friedreich ataxia, Huntington disease, spinocerebellar ataxias, and some forms of amyotrophic lateral sclerosis, are caused by repetitive DNA sequences that are expanded in affected individuals. The age at which an individual begins to experience symptoms, and the severity of disease, are partially determined by the size of the repeat. However, the epigenetic state of the area in and around the repeat also plays an important role in determining the age of disease onset and the rate of disease progression. Many repeat diseases share a common epigenetic pattern of increased methylation at CpG islands near the repeat region. CpG islands are CG-rich sequences that are tightly regulated by methylation and are often found at gene enhancer or insulator elements in the genome. Methylation of CpG islands can inhibit binding of the transcriptional regulator CTCF, resulting in a closed chromatin state and gene down regulation. The downregulation of these genes leads to some disease-specific symptoms. Additionally, a genetic and epigenetic interplay is suggested by an effect of methylation on repeat instability, a hallmark of large repeat expansions that leads to increasing disease severity in successive generations. In this review, we will discuss the common epigenetic patterns shared across repeat diseases, how the genetics and epigenetics interact, and how this could be involved in disease manifestation. We also discuss the currently available stem cell and mouse models, which frequently do not recapitulate epigenetic patterns observed in human disease, and propose alternative strategies to study the role of epigenetics in repeat diseases.

Genetic characterization of F0 outbred and F1 inbred Swiss albino mice using microsatellite markers and their performance evaluation

Swiss albino mice have been widely utilized in various biological researches worldwide. Phenotypic and fitness related traits of F0 and F1 inbred mice were estimated on 918 and 707 individual offsprings, respectively. The influence of fixed effects (litter size and sex) on birth weight (BW), weaning weight (WW) and adult body weight (ABW) in both the generations were found to be statistically significant. Genetic characterization of F0 outbred and the F1 inbred strain of Swiss albino mice were evaluated by using 10 microsatellites markers. The results indicated that total number of alleles per locus ranged from 3 (D2Mit61, D3Mit55, D8Mit14, D9Mit27, D10Mit180, D11Mit167) to 4 (D1Mit15, D2Mit51, D5Mit18, D7Mit323) in F0 and F1 inbred population, with a mean value of 3.4 indicating polymorphism in all 10 loci. The mean of effective number of alleles was 2.935 and 2.733 in F0 and F1 population, respectively. Estimates of the FIS ranged from 0.139 (D10Mit180) to 0.999 (D9Mit27); and from 0.109 (D3Mit55) to 0.679 (D2Mit51) in F0 and F1 inbred population, respectively. The estimated mean markerbased FIS was 0.294 and 0.372 in F0 and F1 populations, respectively. The mean values of observed heterozygosity (Ho) and expected heterozygosity (He) were 0.460 and 0.654, respectively for F0 and 0.390 and 0.627, respectively for F1 inbred mice population. Slight reduction in heterozygosity and 7.8% increase in inbreeding coefficient were observed in F1 inbred in comparison to F0 population. The results suggested that genome wide microsatellite genotyping might be more useful for accurate measuring and reliable estimation of population genetic parameters and inbreeding coefficient.

Strategic enhancement of immune checkpoint inhibition in refractory Colorectal Cancer: Trends and future prospective

Colorectal cancer (CRC), known as a frequently fatal disease, ranking as the third most common malignancy, is the second leading cause of cancer related mortality worldwide. Metastases are common in CRC patients which account for approximately 25% of the patients at diagnosis, 50% of patients during treatment which is associated closely with CRC mortality. Conventional therapies such as surgery, chemotherapy, and radiotherapy are standards of care for the treatment of CRC patients. However, primary tumor recurrence and secondary disease in patients receiving standard of care treatment modalities occur in 50% of patients so that new treatment modalities are needed. Immune checkpoint inhibition (ICI) has transformed the management of patients suffered from metastatic CRC (mCRC) with mismatch repair deficiency (dMMR) and microsatellite instability (MSI) -high (MSI-H) while manifests ineffectiveness in preserved mismatch repair (pMMR) or microsatellite stable (MSS) "cold" tumors which makes up the majority (95%) of mCRC. In this review, we mainly lay emphasis on the development of combinations in therapy strategies with ICIs with other immune based treatment approaches to increase the intra-tumoral immune response and render tumors 'immune-reactive', thereby increasing the efficacy of tumor immunotherapy.

Clear cell endometrial carcinomas with mismatch repair deficiency have a favorable prognosis: A systematic review and meta-analysis

INTRODUCTION

In the ESGO/ESTRO/ESP guidelines for endometrial carcinoma management, the risk category of clear cell carcinoma (CCC) is not well defined. In fact, while p53-abnormal (p53abn) CCC are known to be aggressive, the prognosis of mismatch repair-deficient (MMRd) and p53-wild-type (p53wt) CCCs is less clear.

OBJECTIVE

To assess the prognostic value of the MMRd and p53wt groups in CCC through a systematic review and meta-analysis.

METHODS

Electronic databases were searched from their inception to February 2021. All studies reporting p53 expression, MMR proteins expression and survival outcomes in endometrial CCC (either pure or mixed) were included. Kaplan-Meier and Cox regression survival analyses with hazard ratio (HR) for overall survival (OS) were performed by using the p53abn group as reference; a significant p-value<0.05 was adopted.

RESULTS

Six studies with 136 CCC (114 pure and 22 mixed) were included. Five-year OS was 95.7 ± 4.3% in the MMRd group, 48.4 ± 8.4% months in the p53wt group and 40.6 ± 10.4% in the p53abn group. The hazard of death was significantly lower in the MMRd group than in the p53abn group (HR = 0.062; p = 0.007), while it did not significantly differ between the p53wt and the p53abn group (HR = 0.673; p = 0.222). The POLEmut group could not be analyzed due to the absence of deaths. Similar results were observed in the pure CCC and mixed CCC subgroups.

CONCLUSION

MMRd CCCs seem to have a favorable prognosis and might be lumped together with MMRd endometrioid carcinoma for management purpose. On the other hand, p53wt CCCs appear prognostically more similar to p53abn CCCs.

Genetic evaluation of F2 inbred strain of Swiss albino mice by microsatellite markers

Genetic characterization of F2 inbred Swiss albino mice was done using 11 microsatellites markers (D1Mit15, D2Mit51, D2Mit61, D3Mit15, D3Mit55, D5Mit18, D7Mit323, D8Mit14, D9Mit27, D10Mit180, and D11Mit167). The results indicated that genotypic frequencies at all the loci were in Hardy-Weinberg disequilibrium (P<0.001). Total number of alleles per locus ranged from 3 (D2Mit61, D3Mit15, D8Mit14, D9Mit27, D10Mit180, D11Mit167) to 4 (D1Mit15, D2Mit51, D3Mit55, D5Mit18, D7Mit323), with a mean of 3.45. The numbers of effective alleles ranged between 1.847 and 3.814. The observed heterozygosity (Ho) was maximum for D1Mit15 (0.660) and minimum for D5Mit18 (0.080), with mean of 0.269. The FIS estimates was ranged from 0.0817 (D1Mit15) to 0.8799 (D5Mit18). The average inbreeding coefficient was 0.592, which indicates that parents were more related than expected under random mating. The range of PIC value (from 0.414 to 0.689) for various microsatellite loci was revealing that population under investigation was of high diversity maintaining a multiple allele.

A refined panel of 42 microsatellite loci to universally genotype catarrhine primates

Microsatellite genotyping is an important genetic method for a number of research questions in biology. Given that the traditional fragment length analysis using polyacrylamide gel or capillary electrophoresis has several drawbacks, microsatellite genotyping-by-sequencing (GBS) has arisen as a promising alternative. Although GBS mitigates many of the problems of fragment length analysis, issues with allelic dropout and null alleles often remain due to mismatches in primer binding sites and unnecessarily long PCR products. This is also true for GBS in catarrhine primates where cross-species amplification of loci (often human derived) is common.We therefore redesigned primers for 45 microsatellite loci based on 17 available catarrhine reference genomes. Next, we tested them in singleplex and different multiplex settings in a panel of species representing all major lineages of Catarrhini and further validated them in wild Guinea baboons (Papio papio) using fecal samples.The final panel of 42 microsatellite loci can efficiently be amplified with primers distributed into three amplification pools.With our microsatellite panel, we provide a tool to universally genotype catarrhine primates via GBS from different sample sources in a cost- and time-efficient way, with higher resolution, and comparability among laboratories and species.

Characterizing glycosyltransferases by a combination of sequencing platforms applied to the leaf tissues of Stevia rebaudiana

Background

Stevia rebaudiana (Bertoni) is considered one of the most valuable plants because of the steviol glycosides (SGs) that can be extracted from its leaves. Glycosyltransferases (GTs), which can transfer sugar moieties from activated sugar donors onto saccharide and nonsaccharide acceptors, are widely distributed in the genome of S. rebaudiana and play important roles in the synthesis of steviol glycosides.

Results

Six stevia genotypes with significantly different concentrations of SGs were obtained by induction through various mutagenic methods, and the contents of seven glycosides (stevioboside, Reb B, ST, Reb A, Reb F, Reb D and Reb M) in their leaves were considerably different. Then, NGS and single-molecule real-time (SMRT) sequencing were combined to analyse leaf tissue from these six different genotypes to generate a full-length transcriptome of S. rebaudiana. Two phylogenetic trees of glycosyltransferases (SrUGTs) were constructed by the neighbour-joining method and successfully predicted the functions of SrUGTs involved in SG biosynthesis. With further insight into glycosyltransferases (SrUGTs) involved in SG biosynthesis, the weighted gene co-expression network analysis (WGCNA) method was used to characterize the relationships between SrUGTs and SGs, and forty-four potential SrUGTs were finally obtained, including SrUGT85C2, SrUGT74G1, SrUGT76G1 and SrUGT91D2, which have already been reported to be involved in the glucosylation of steviol glycosides, illustrating the reliability of our results.

Conclusion

Combined with the results obtained by previous studies and those of this work, we systematically characterized glycosyltransferases in S. rebaudiana and forty-four candidate SrUGTs involved in the glycosylation of steviol glucosides were obtained. Moreover, the full-length transcriptome obtained in this study will provide valuable support for further research investigating S. rebaudiana.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07195-5.

Immune activation by a multigene family of lectins with variable tandem repeats in oriental river prawn (Macrobrachium nipponense)

Genomic regions with repeated sequences are unstable and prone to rapid DNA diversification. However, the role of tandem repeats within the coding region is not fully characterized. Here, we have identified a new hypervariable C-type lectin gene family with different numbers of tandem repeats (Rlecs; R means repeat) in oriental river prawn (Macrobrachium nipponense). Two types of repeat units (33 or 30 bp) are identified in the second exon, and the number of repeat units vary from 1 to 9. Rlecs can be classified into 15 types through phylogenetic analysis. The amino acid sequences in the same type of Rlec are highly conservative outside the repeat regions. The main differences among the Rlec types are evident in exon 5. A variable number of tandem repeats in Rlecs may be produced by slip mispairing during gene replication. Alternative splicing contributes to the multiplicity of forms in this lectin gene family, and different types of Rlecs vary in terms of tissue distribution, expression quantity and response to bacterial challenge. These variations suggest that Rlecs have functional diversity. The results of experiments on sugar binding, microbial inhibition and clearance, regulation of antimicrobial peptide gene expression and prophenoloxidase activation indicate that the function of Rlecs with the motif of YRSKDD in innate immunity is enhanced when the number of tandem repeats increases. Our results suggest that Rlecs undergo gene expansion through gene duplication and alternative splicing, which ultimately leads to functional diversity.

Genetic diversity of avocado from the southern highlands of Tanzania as revealed by microsatellite markers

BACKGROUND

Avocado is an important cash crop in Tanzania, however its genetic diversity is not thoroughly investigated. This study was undertaken to explore the genetic diversity of avocado in the southern highlands using microsatellite markers. A total of 226 local avocado trees originating from seeds were sampled in eight districts of the Mbeya, Njombe and Songwe regions. Each district was considered as a population. The diversity at 10 microsatellite loci was investigated.

RESULTS

A total of 167 alleles were detected across the 10 loci with an average of 16.7 ± 1.3 alleles per locus. The average expected and observed heterozygosity were 0.84 ± 0.02 and 0.65 ± 0.04, respectively. All but two loci showed a significant deviation from the Hardy-Weinberg principle. Analysis of molecular variance showed that about 6% of the variation was partitioned among the eight geographic populations. Population F_ST pairwise comparisons revealed lack of genetic differentiation for the seven of 28 population pairs tested. The principal components analysis (PCA) and hierarchical cluster analysis showed a mixing of avocado trees from different districts. The model-based STRUCTURE subdivided the trees samples into four major genetic clusters.

CONCLUSION

High diversity detected in the analysed avocado germplasm implies that this germplasm is a potentially valuable source of variable alleles that might be harnessed for genetic improvement of this crop in Tanzania. The mixing of avocado trees from different districts observed in the PCA and dendrogram points to strong gene flow among the avocado populations, which led to population admixture revealed in the STRUCTURE analysis. However, there is still significant differentiation among the tree populations from different districts that can be utilized in the avocado breeding program.

Characterization and analysis of full-length transcriptomes from two grasshoppers, Gomphocerus licenti and Mongolotettix japonicus

Acrididae are diverse in size, body shape, behavior, ecology and life history; widely distributed; easy to collect; and important to agriculture. They represent promising model candidates for functional genomics, but their extremely large genomes have hindered this research; establishing a reference transcriptome for a species is the primary means of obtaining genetic information. Here, two Acrididae species, Gomphocerus licenti and Mongolotettix japonicus, were selected for full-length (FL) PacBio transcriptome sequencing. For G. licenti and M. japonicus, respectively, 590,112 and 566,165 circular consensus sequences (CCS) were generated, which identified 458,131 and 428,979 full-length nonchimeric (FLNC) reads. After isoform-level clustering, next-generation sequencing (NGS) short sequences were used for error correction, and remove redundant sequences with CD-HIT, 17,970 and 16,766 unigenes were generated for G. licenti and M. japonicus. In addition, we obtained 17,495 and 16,373 coding sequences, 1,082 and 813 transcription factors, 11,840 and 10,814 simple sequence repeats, and 905 and 706 long noncoding RNAs by analyzing the transcriptomes of G. licenti and M. japonicus, respectively, and 15,803 and 14,846 unigenes were annotated in eight functional databases. This is the first study to sequence FL transcriptomes of G. licenti and M. japonicus, providing valuable genetic resources for further functional genomics research.

Alopecia areata susceptibility variant in MHC region impacts expressions of genes contributing to hair keratinization and is involved in hair loss

Background

Alopecia areata (AA) is considered a highly heritable, T-cell-mediated autoimmune disease of the hair follicle. However, no convincing susceptibility gene has yet been pinpointed in the major histocompatibility complex (MHC), a genome region known to be associated with AA as compared to other regions.

Methods

We engineered mice carrying AA risk allele identified by haplotype sequencing for the MHC region using allele-specific genome editing with the CRISPR/Cas9 system. Finally, we performed functional evaluations in the mice and AA patients with and without the risk allele.

Findings

We identified a variant (rs142986308, p.Arg587Trp) in the coiled-coil alpha-helical rod protein 1 (CCHCR1) gene as the only non-synonymous variant in the AA risk haplotype. Furthermore, mice engineered to carry the risk allele displayed a hair loss phenotype. Transcriptomics further identified CCHCR1 as a novel component interacting with hair cortex keratin in hair shafts. Both, these alopecic mice and AA patients with the risk allele displayed morphologically impaired hair and comparable differential expression of hair-related genes, including hair keratin and keratin-associated proteins (KRTAPs).

Interpretation

Our results implicate CCHCR1 with the risk allele in a previously unidentified subtype of AA based on aberrant keratinization in addition to autoimmune events.

Funding

This work was supported by JSPS KAKENHI (JP16K10177) and the NIHR UCLH Biomedical Research center (BRC84/CN/SB/5984).

Low temperature isothermal amplification of microsatellites drastically reduces stutter artifact formation and improves microsatellite instability detection in cancer

Microsatellites are polymorphic short tandem repeats of 1–6 nucleotides ubiquitously present in the genome that are extensively used in living organisms as genetic markers and in oncology to detect microsatellite instability (MSI). While the standard analysis method of microsatellites is based on PCR followed by capillary electrophoresis, it generates undesirable frameshift products known as ‘stutter peaks’ caused by the polymerase slippage that can greatly complicate the analysis and interpretation of the data. Here we present an easy multiplexable approach replacing PCR that is based on low temperature isothermal amplification using recombinase polymerase amplification (LT-RPA) that drastically reduces and sometimes completely abolishes the formation of stutter artifacts, thus greatly simplifying the calling of the alleles. Using HT17, a mononucleotide DNA repeat that was previously proposed as an optimal marker to detect MSI in tumor DNA, we showed that LT-RPA improves the limit of detection of MSI compared to PCR up to four times, notably for small deletions, and simplifies the identification of the mutant alleles. It was successfully applied to clinical colorectal cancer samples and enabled detection of MSI. This easy-to-handle, rapid and cost-effective approach may deeply improve the analysis of microsatellites in several biological and clinical applications.

A unified STR profiling system across multiple species with whole genome sequencing data

Background

Short tandem repeats (STRs) serve as genetic markers in forensic scenes due to their high polymorphism in eukaryotic genomes. A variety of STRs profiling systems have been developed for species including human, dog, cat, cattle, etc. Maintaining these systems simultaneously can be costly. These mammals share many high similar regions along their genomes. With the availability of the massive amount of the whole genomics data of these species, it is possible to develop a unified STR profiling system. In this study, our objective is to propose and develop a unified set of STR loci that could be simultaneously applied to multiple species.

Result

To find a unified STR set, we collected the whole genome sequence data of the concerned species and mapped them to the human genome reference. Then we extracted the STR loci across the species. From these loci, we proposed an algorithm which selected a subset of loci by incorporating the optimized combined power of discrimination. Our results show that the unified set of loci have high combined power of discrimination, >1−10⁻⁹, for both individual species and the mixed population, as well as the random-match probability, <10⁻⁷ for all the involved species, indicating that the identified set of STR loci could be applied to multiple species.

Conclusions

We identified a set of STR loci which shared by multiple species. It implies that a unified STR profiling system is possible for these species under the forensic scenes. The system can be applied to the individual identification or paternal test of each of the ten common species which are Sus scrofa (pig), Bos taurus (cattle), Capra hircus (goat), Equus caballus (horse), Canis lupus familiaris (dog), Felis catus (cat), Ovis aries (sheep), Oryctolagus cuniculus (rabbit), and Bos grunniens (yak), and Homo sapiens (human). Our loci selection algorithm employed a greedy approach. The algorithm can generate the loci under different forensic parameters and for a specific combination of species.

Complete chloroplast genomes of two Siraitia Merrill species: Comparative analysis, positive selection and novel molecular marker development

Siraitia grosvenorii fruit, known as Luo-Han-Guo, has been used as a traditional Chinese medicine for many years, and mogrosides are its primary active ingredients. Unfortunately, Siraitia siamensis, its wild relative, might be misused due to its indistinguishable appearance, not only threatening the reliability of the medication but also partly exacerbating wild resource scarcity. Therefore, high-resolution genetic markers must be developed to discriminate between these species. Here, the complete chloroplast genomes of S. grosvenorii and S. siamensis were assembled and analyzed for the first time; they were 158,757 and 159,190 bp in length, respectively, and possessed conserved quadripartite circular structures. Both contained 134 annotated genes, including 8 rRNA, 37 tRNA and 89 protein-coding genes. Twenty divergences (Pi > 0.03) were found in the intergenic regions. Nine protein-coding genes, accD, atpA, atpE, atpF, clpP, ndhF, psbH, rbcL, and rpoC2, underwent selection within Cucurbitaceae. Phylogenetic relationship analysis indicated that these two species originated from the same ancestor. Finally, four pairs of molecular markers were developed to distinguish the two species. The results of this study will be beneficial for taxonomic research, identification and conservation of Siraitia Merrill wild resources in the future.

Rare and common variant discovery in complex disease: the IBD case study

Complex diseases such as inflammatory bowel disease (IBD), which consists of ulcerative colitis and Crohn's disease, are a significant medical burden-70 000 new cases of IBD are diagnosed in the United States annually. In this review, we examine the history of genetic variant discovery in complex disease with a focus on IBD. We cover methods that have been applied to microsatellite, common variant, targeted resequencing and whole-exome and -genome data, specifically focusing on the progression of technologies towards rare-variant discovery. The inception of these methods combined with better availability of population level variation data has led to rapid discovery of IBD-causative and/or -associated variants at over 200 loci; over time, these methods have grown exponentially in both power and ascertainment to detect rare variation. We highlight rare-variant discoveries critical to the elucidation of the pathogenesis of IBD, including those in NOD2, IL23R, CARD9, RNF186 and ADCY7. We additionally identify the major areas of rare-variant discovery that will evolve in the coming years. A better understanding of the genetic basis of IBD and other complex diseases will lead to improved diagnosis, prognosis, treatment and surveillance.

Comparison of chloroplast genomes of Gynura species: sequence variation, genome rearrangement and divergence studies

Background

Some Gynura species have been reported to be natural anti-diabetic plants. Improvement of their traits towards application relies on hybridization. Clearly, phylogenetic relationships could optimize compatible hybridizations. For flowerings plants, chloroplast genomes have been used to solve many phylogenetic relationships. To date, the chloroplast genome sequences of 4 genera of the tribe Senecioneae have been uploaded to GenBank. The internal relationships within the genus Gynura and the relationship of the genus Gynura with other genera in the tribe Senecioneae need further research.

Results

The chloroplast genomes of 4 Gynura species were sequenced, assembled and annotated. In comparison with those of 12 other Senecioneae species, the Gynura chloroplast genome features were analysed in detail. Subsequently, differences in the microsatellite and repeat types in the tribe were found. From the comparison, it was found that IR expansion and contraction are conserved in the genera Gynura, Dendrosenecio and Ligularia. Compared to other regions on the chloroplast genome, the region from 25,000 to 50,000 bp was not conserved. Seven ndh genes in this region are under purifying selection, with small changes in amino acids. The whole chloroplast genome sequences of 16 Senecioneae species were used to build a phylogenetic tree. Based on the oldest Artemisia pollen fossil, the divergence time was estimated.

Conclusions

Sequencing the chloroplast genomes of 4 Gynura species helps us to solve many problems. The phylogenetic relationships and divergence time among 4 Gynura and 16 Senecioneae species were evaluated by comparing their chloroplast genomes. The phylogenetic relationship of the genera Gynura and Ligularia was different from that observed previous work. In a previous phylogenetic tree, the genus Ligularia belonged to the Tussilagininae subtribe, which was in a lineage that diverged earlier than other genera. Further morphology and genome-wide analyses are needed to clarify the genus relationships.

Genetics without genes? The centrality of genetic markers in livestock genetics and genomics

In this paper, rather than focusing on genes as an organising concept around which historical considerations of theory and practice in genetics are elucidated, we place genetic markers at the heart of our analysis. This reflects their central role in the subject of our account, livestock genetics concerning the domesticated pig, Sus scrofa. We define a genetic marker as a (usually material) element existing in different forms in the genome, that can be identified and mapped using a variety (and often combination) of quantitative, classical and molecular genetic techniques. The conjugation of pig genome researchers around the common object of the marker from the early-1990s allowed the distinctive theories and approaches of quantitative and molecular genetics concerning the size and distribution of gene effects to align (but never fully integrate) in projects to populate genome maps. Critical to this was the nature of markers as ontologically inert, internally heterogeneous and relational. Though genes as an organising and categorising principle remained important, the particular concatenation of limitations, opportunities, and intended research goals of the pig genetics community, meant that a progressively stronger focus on the identification and mapping of markers rather than genes per se became a hallmark of the community. We therefore detail a different way of doing genetics to more gene-centred accounts. By doing so, we reveal the presence of practices, concepts and communities that would otherwise be hidden.

OSTRFPD: Multifunctional Tool for Genome-Wide Short Tandem Repeat Analysis for DNA, Transcripts, and Amino Acid Sequences with Integrated Primer Designer

Microsatellite mining is a common outcome of the in silico approach to genomic studies. The resulting short tandemly repeated DNA could be used as molecular markers for studying polymorphism, genotyping and forensics. The omni short tandem repeat finder and primer designer (OSTRFPD) is among the few versatile, platform-independent open-source tools written in Python that enables researchers to identify and analyse genome-wide short tandem repeats in both nucleic acids and protein sequences. OSTRFPD is designed to run either in a user-friendly fully featured graphical interface or in a command line interface mode for advanced users. OSTRFPD can detect both perfect and imperfect repeats of low complexity with customisable scores. Moreover, the software has built-in architecture to simultaneously filter selection of flanking regions in DNA and generate microsatellite-targeted primers implementing the Primer3 platform. The software has built-in motif-sequence generator engines and an additional option to use the dictionary mode for custom motif searches. The software generates search results including general statistics containing motif categorisation, repeat frequencies, densities, coverage, guanine–cytosine (GC) content, and simple text-based imperfect alignment visualisation. Thus, OSTRFPD presents users with a quick single-step solution package to assist development of microsatellite markers and categorise tandemly repeated amino acids in proteome databases. Practical implementation of OSTRFPD was demonstrated using publicly available whole-genome sequences of selected Plasmodium species. OSTRFPD is freely available and open-sourced for improvement and user-specific adaptation.

SMRT sequencing of the full-length transcriptome of the Sunda pangolin (Manis javanica)

It is widely known that transcriptional diversity contributes greatly to biological regulation in eukaryotes. With the development of next-generation sequencing (NGS) technologies, several studies on RNA sequencing have considerably improved our understanding of transcriptome complexity. However, obtaining full-length (FL) transcripts remains a considerable challenge because of difficulties in short read-based assembly. In the present study, single-molecule real-time (SMRT) sequencing and NGS were combined to generate the complete and FL transcriptome of Manis javanica. The results provide a comprehensive set of reference transcripts and hence contribute to the improved annotation of the M. javanica genome. We obtained 45,530 high-confidence transcripts from 19,109 genic loci, of which 8014 genes have not yet been annotated within the M. javanica genome. Furthermore, we revealed 8824 long-chain noncoding RNAs (lncRNAs). A total of 30,199 alternative splicing (AS) and 11,184 alternative polyadenylation (APA) events were identified in the sequencing data. The structure and expression level of 59 digestive enzyme genes, including 13 carbohydrase genes, 28 lipase genes and 18 protease genes, were analyzed, which might provide original data for further research on M. javanica.

GWAS analysis in spring barley (Hordeum vulgare L.) for morphological traits exposed to drought

Association analysis based on linkage disequilibrium has become a common and powerful approach for detection of QTLs underlying complex agronomic traits including drought tolerance. To determine marker/trait association, 148 modern European spring barley cultivars were evaluated under drought stress. Associations of morphological traits with AFLP/SSR markers were investigated based on the mixed linear model using the TASSEL_3.0. Population structure was estimated using various methods including Bayesian clustering model by STRUCTURE software, PCoA analysis, NJ dendrogram and Hierarchical Clustering. Linkage disequilibrium patterns were explored among the whole genome and each chromosome separately. All the analysis for population structure divided the population into two sub-groups. Linkage disequilibrium analysis showed that by increasing genetic distance, LD decreases. Totally, 167 significant marker trait associations were found which delineated into 65 QTLs in both treatments. Two stable QTLs on 5H at 86.880 cM were detected for Internode Length and on 3H at 126.421 cM for flag leaf length in drought stress treatment. Fourteen QTLs were co-localized with previously reported QTLs and others were novel. The results indicate that these putative genomic regions contain genes that have pleiotropic effects on morphological traits in drought condition.

SMRT sequencing of full-length transcriptome of flea beetle Agasicles hygrophila (Selman and Vogt)

This study was aimed at generating the full-length transcriptome of flea beetle Agasicles hygrophila (Selman and Vogt) using single-molecule real-time (SMRT) sequencing. Four developmental stages of A. hygrophila, including eggs, larvae, pupae, and adults were harvested for isolating total RNA. The mixed samples were used for SMRT sequencing to generate the full-length transcriptome. Based on the obtained transcriptome data, alternative splicing event, simple sequence repeat (SSR) analysis, coding sequence prediction, transcript functional annotation, and lncRNA prediction were performed. Total 9.45 Gb of clean reads were generated, including 335,045 reads of insert (ROI) and 158,085 full-length non-chimeric (FLNC) reads. Transcript clustering analysis of FLNC reads identified 40,004 consensus isoforms, including 31,015 high-quality ones. After removing redundant reads, 28,982 transcripts were obtained. Total 145 alternative splicing events were predicted. Additionally, 12,753 SSRs and 16,205 coding sequences were identified based on SSR analysis. Furthermore, 24,031 transcripts were annotated in eight functional databases, and 4,198 lncRNAs were predicted. This is the first study to perform SMRT sequencing of the full-length transcriptome of A. hygrophila. The obtained transcriptome may facilitate further exploration of the genetic data of A. hygrophila and uncover the interactions between this insect and the ecosystem.

RXRB Is an MHC-Encoded Susceptibility Gene Associated with Anti-Topoisomerase I Antibody-Positive Systemic Sclerosis

Systemic sclerosis is a systemic autoimmune and connective tissue disorder associated with the human leukocyte antigen locus. However, the functional relationship between human leukocyte antigen gene(s) and disease development remains unknown. To elucidate major histocompatibility complex-linked systemic sclerosis genetics, we performed genotyping of major histocompatibility complex-borne microsatellites and HLA-DPB1 alleles using DNA obtained from 318 anti-topoisomerase I antibody-positive patients and 561 healthy controls, all of Japanese descent. Those results revealed two major histocompatibility complex haplotypes associated with systemic sclerosis. Exome sequencing and targeted analysis of these risk haplotypes identified rs17847931 in RXRB as a susceptibility variant (P = 1.3 × 10^-15; odds ratio [OR] = 9.4) with amino acid substitution p.V95A on the risk haplotype harboring HLA-DPB1∗13:01. No identical variant in the other haplotype including DPB1*09:01 was observed, though that haplotype also showed a significant association (P = 8.5 × 10^-22; OR = 4.3) with systemic sclerosis. Furthermore, the number of risk factors was shown to be a predominant factor, as individuals with two factors had elevated risk (P = 6.7 × 10^-13; OR = 30.2). We concluded that RXRB may be involved in antifibrotic activity in skin and chromatin remodeling.

New Repeat Polymorphism in the AKT1 Gene Predicts Striatal Dopamine D2/D3 Receptor Availability and Stimulant-Induced Dopamine Release in the Healthy Human Brain

The role of the protein kinase Akt1 in dopamine neurotransmission is well recognized and has been implicated in schizophrenia and psychosis. However, the extent to which variants in the AKT1 gene influence dopamine neurotransmission is not well understood. Here we investigated the effect of a newly characterized variant number tandem repeat (VNTR) polymorphism in AKT1 [major alleles: L- (eight repeats) and H- (nine repeats)] on striatal dopamine D2/D3 receptor (DRD2) availability and on dopamine release in healthy volunteers. We used PET and [¹¹C]raclopride to assess baseline DRD2 availability in 91 participants. In 54 of these participants, we also measured intravenous methylphenidate-induced dopamine release to measure dopamine release. Dopamine release was quantified as the difference in specific binding of [¹¹C]raclopride (nondisplaceable binding potential) between baseline values and values following methylphenidate injection. There was an effect of AKT1 genotype on DRD2 availability at baseline for the caudate (F_(2,90) = 8.2, p = 0.001) and putamen (F_(2,90) = 6.6, p = 0.002), but not the ventral striatum (p = 0.3). For the caudate and putamen, LL showed higher DRD2 availability than HH; HL were in between. There was also a significant effect of AKT1 genotype on dopamine increases in the ventral striatum (F_(2,53) = 5.3, p = 0.009), with increases being stronger in HH > HL > LL. However, no dopamine increases were observed in the caudate (p = 0.1) or putamen (p = 0.8) following methylphenidate injection. Our results provide evidence that the AKT1 gene modulates both striatal DRD2 availability and dopamine release in the human brain, which could account for its association with schizophrenia and psychosis. The clinical relevance of the newly characterized AKT1 VNTR merits investigation.SIGNIFICANCE STATEMENT The AKT1 gene has been implicated in schizophrenia and psychosis. This association is likely to reflect modulation of dopamine signaling by Akt1 kinase since striatal dopamine hyperstimulation is associated with psychosis and schizophrenia. Here, using PET with [¹¹C]raclopride, we identified in the AKT1 gene a new variable number tandem repeat (VNTR) marker associated with baseline striatal dopamine D2/D3 receptor availability and with methylphenidate-induced striatal dopamine increases in healthy volunteers. Our results confirm the involvement of the AKT1 gene in modulating striatal dopamine signaling in the human brain. Future studies are needed to assess the association of this new VNTR AKT1 variant in schizophrenia and drug-induced psychoses.

Genome-wide association study of IgA nephropathy using 23 465 microsatellite markers in a Japanese population

Immunoglobulin A nephropathy (IgAN) is the most common form of primary glomerulonephritis in many parts of the world. Although previous genome-wide association studies (GWAS) identified the major susceptibility loci for IgAN, the causal genes currently remain unknown. We performed a GWAS using 23 465 microsatellite (MS) markers to identify genes related to IgAN in a Japanese population. A pooled sample analysis was conducted in three-stage screenings of three independent case-control populations, and after the final step of individual typing, 11 markers survived. Of these, we focused on two regions on 6p21 and 12q21 because they (i) showed the strongest relationship with IgAN, and (ii) appeared to be highly relevant to IgAN in view of several previous studies. These regions contained the HLA, TSPAN8 and PTPRR genes. This study on GWAS, using >20 000 MS markers, provides a new approach regarding susceptible genes for IgAN for investigators seeking new tools for the prevention and treatment of IgAN.

Evidence for a genetic basis of urogenital carcinoma in the wild California sea lion

Although neoplasia is a major cause of mortality in humans and domestic animals, it has rarely been described in wildlife species. One of the few examples is a highly prevalent urogenital carcinoma in California sea lions (CSLs). Although the aetiology of this carcinoma is clearly multifactorial, inbreeding depression, as estimated using levels of microsatellite multilocus heterozygosity, is identified as predictive for this neoplasia. On further analysis, this relationship appears to be largely driven by one marker, suggesting that a single locus might be associated with the occurrence of this disease in CSLs. In a case–control study, carcinoma was significantly associated with homozygosity at the Pv11 microsatellite locus. Pv11 was mapped to intron 9 of the heparanase 2 gene (HPSE2) locus, a very large gene encoding heparanase 2, which in humans is associated with multiple carcinomas. Correspondingly, immunohistochemical labelling in tissues was present in carcinoma cases within a single homozygous Pv11 genotype. To our knowledge, this is the first report of an individual locus being associated with cancer in any wildlife species. This adds emphasis to the study of HPSE2 in other species, including humans and will guide future studies on this sentinel species that shares much of its diet and environment with humans

Chemotherapy induced microsatellite instability and loss of heterozygosity in chromosomes 2, 5, 10, and 17 in solid tumor patients

Background

The inevitable side effects of the currently used chemotherapy are associated with serious syndromes. Genotoxic effects and consequent genetic instability may play an important role in these syndromes. The aim of the study was to evaluate chemotherapy-related microsatellite instability (MSI), loss of heterozygosity (LOH), and loss of mismatch repair (MMR) expression in solid tumor patients.

Methods

Samples were collected from 117 de novo patients with solid tumors of different origins. Specimens, taken pre- and post-treatment, were screened for MSI and LOH in 10 microsatellite sequences in blood, and expression of five MMR proteins were analyzed in cancer tissues using immunohistochemistry. Statistical analysis included the use of; Fisher’s exact test, Chi Square, and an inter-rater reliability test using Cohen’s kappa coefficient.

Results

Microsatellite analysis showed that 66.7% of the patients had MSI, including 23.1% high-positive MSI and 43.6% low-positive MSI. A large portion (41%) of the patients exhibited LOH in addition to MSI. MSI and LOH were detected in seven loci in which incidence rates ranged from 3.8% positive for Bat-26 to 34.6% positive for Tp53-Alu. Immunohistochemistry revealed that human mutL homolog 1 (hMLH1) expression was deficient in 29.1% of the patients, whereas 18.8%, 23.9%, 13.4%, and 9.7% were deficient for human mutS homolog 2 (hMSH2), P53, human mutS homolog 6 (hMSH6) and human post-meiotic segregation increased 2 (hPMS2), respectively. There was a significant correlation between MSI and LOH incidence in Tp53-Alu, Mfd41, and APC with low or deficient expression of hMLH1, hMSH2, and P53. A significant association between MSI and LOH, and incidence of secondary tumors was also evident.

Conclusions

The negative correlation between MMR expression, MSI, and LOH and increased resistance to anti-cancer drugs and development of secondary cancers demonstrates a useful aid in early detection of potential chemotherapy-related side-effects. The diagnostic value demonstrated in our earlier study on breast cancer patients was confirmed for other solid tumors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12935-014-0118-4) contains supplementary material, which is available to authorized users.

Association of vWA and TPOX polymorphisms with venous thrombosis in Mexican mestizos

Objective. Venous thromboembolism (VTE) is a multifactorial disorder and, worldwide, the most important cause of morbidity and mortality. Genetic factors play a critical role in its aetiology. Microsatellites are the most important source of human genetic variation having more phenotypic effect than many single nucleotide polymorphisms. Hence, we evaluate a possible relationship between VTE and the genetic variants in von Willebrand factor, human alpha fibrinogen, and human thyroid peroxidase microsatellites to identify possible diagnostic markers. Methods. Genotypes were obtained from 177 patients with VTE and 531 nonrelated individuals using validated genotyping methods. The allelic frequencies were compared; Bayesian methods were used to correct population stratification to avoid spurious associations. Results. The vWA-18, TPOX-9, and TPOX-12 alleles were significantly associated with VTE. Moreover, subjects bearing the combination vWA-18/TPOX-12 loci exhibited doubled risk for VTE (95% CI = 1.02–3.64), whereas the combination vWA-18/TPOX-9 showed an OR = 10 (95% CI = 4.93–21.49). Conclusions. The vWA and TPOX microsatellites are good candidate biomarkers in venous thromboembolism diseases and could help to elucidate their origins. Additionally, these polymorphisms could become useful markers for genetic studies of VTE in the Mexican population; however, further studies should be done owing that this data only show preliminary evidence.

A systematic evaluation of short tandem repeats in lipid candidate genes: riding on the SNP-wave

Structural genetic variants as short tandem repeats (STRs) are not targeted in SNP-based association studies and thus, their possible association signals are missed. We systematically searched for STRs in gene regions known to contribute to total cholesterol, HDL cholesterol, LDL cholesterol and triglyceride levels in two independent studies (KORA F4, n = 2553 and SAPHIR, n = 1648), resulting in 16 STRs that were finally evaluated. In a combined dataset of both studies, the sum of STR alleles was regressed on each phenotype, adjusted for age and sex. The association analyses were repeated for SNPs in a 200 kb region surrounding the respective STRs in the KORA F4 Study. Three STRs were significantly associated with total cholesterol (within LDLR, the APOA1/C3/A4/A5/BUD13 gene region and ABCG5/8), five with HDL cholesterol (3 within CETP, one in LPL and one inAPOA1/C3/A4/A5/BUD13), three with LDL cholesterol (LDLR, ABCG5/8 and CETP) and two with triglycerides (APOA1/C3/A4/A5/BUD13 and LPL). None of the investigated STRs, however, showed a significant association after adjusting for the lead or adjacent SNPs within that gene region. The evaluated STRs were found to be well tagged by the lead SNP within the respective gene regions. Therefore, the STRs reflect the association signals based on surrounding SNPs. In conclusion, none of the STRs contributed additionally to the SNP-based association signals identified in GWAS on lipid traits.

Investigation of phase shifts for different period lengths in the genomes of C. elegans, D. melanogaster and S. cerevisiae

We describe a new mathematical method for finding very diverged short tandem repeats containing a single indel. The method involves comparison of two frequency matrices: a first matrix for a subsequence before shift and a second one for a subsequence after it. A measure of comparison is based on matrix similarity. The approach developed was applied to analysis of the genomes of Caenorhabditis elegans, Drosophila melanogaster and Saccharomyces cerevisiae. They were investigated regarding the presence of tandem repeats having repeat length equal to 2 - 11 nucleotides except equal to 3, 6 and 9 nucleotides. A number of phase shift regions for these genomes was approximately 2.2 × 10(4), 1.5 × 10(4) and 1.7 × 10(2), respectively. Type I error was less than 5%. The mean length of fuzzy periodicity and phase shift regions was about 220 nucleotides. The regions of fuzzy periodicity having single insertion or deletion occupy substantial parts of the genomes: 5%, 3% and 0.3%, respectively. Only less than 10% of these regions have been detected previously. That is, the number of such regions in the genomes of C. elegans, D. melanogaster and S. cerevisiae is dramatically higher than it has been revealed by any known methods. We suppose that some found regions of fuzzy periodicity could be the regions for protein binding.

Large-scale analysis of tandem repeat variability in the human genome

Tandem repeats are short DNA sequences that are repeated head-to-tail with a propensity to be variable. They constitute a significant proportion of the human genome, also occurring within coding and regulatory regions. Variation in these repeats can alter the function and/or expression of genes allowing organisms to swiftly adapt to novel environments. Importantly, some repeat expansions have also been linked to certain neurodegenerative diseases. Therefore, accurate sequencing of tandem repeats could contribute to our understanding of common phenotypic variability and might uncover missing genetic factors in idiopathic clinical conditions. However, despite long-standing evidence for the functional role of repeats, they are largely ignored because of technical limitations in sequencing, mapping and typing. Here, we report on a novel capture technique and data filtering protocol that allowed simultaneous sequencing of thousands of tandem repeats in the human genomes of a three generation family using GS-FLX-plus Titanium technology. Our results demonstrated that up to 7.6% of tandem repeats in this family (4% in coding sequences) differ from the reference sequence, and identified a de novo variation in the family tree. The method opens new routes to look at this underappreciated type of genetic variability, including the identification of novel disease-related repeats.

STaRRRT: a table of short tandem repeats in regulatory regions of the human genome

Background

Tandem repeats (TRs) are unstable regions commonly found within genomes that have consequences for evolution and disease. In humans, polymorphic TRs are known to cause neurodegenerative and neuromuscular disorders as well as being associated with complex diseases such as diabetes and cancer. If present in upstream regulatory regions, TRs can modify chromatin structure and affect transcription; resulting in altered gene expression and protein abundance. The most common TRs are short tandem repeats (STRs), or microsatellites. Promoter located STRs are considerably more polymorphic than coding region STRs. As such, they may be a common driver of phenotypic variation. To study STRs located in regulatory regions, we have performed genome-wide analysis to identify all STRs present in a region that is 2 kilobases upstream and 1 kilobase downstream of the transcription start sites of genes.

Results

The Short Tandem Repeats in Regulatory Regions Table, STaRRRT, contains the results of the genome-wide analysis, outlining the characteristics of 5,264 STRs present in the upstream regulatory region of 4,441 human genes. Gene set enrichment analysis has revealed significant enrichment for STRs in cellular, transcriptional and neurological system gene promoters and genes important in ion and calcium homeostasis. The set of enriched terms has broad similarity to that seen in coding regions, suggesting that regulatory region STRs are subject to similar evolutionary pressures as STRs in coding regions and may, like coding region STRs, have an important role in controlling gene expression.

Conclusions

STaRRRT is a readily-searchable resource for investigating potentially polymorphic STRs that could influence the expression of any gene of interest. The processes and genes enriched for regulatory region STRs provide potential novel targets for diagnosing and treating disease, and support a role for these STRs in the evolution of the human genome.

Association of intercellular adhesion molecule 1 (ICAM1) with diabetes and diabetic nephropathy

Diabetes and diabetic nephropathy are complex diseases affected by genetic and environmental factors. Identification of the susceptibility genes and investigation of their roles may provide useful information for better understanding of the pathogenesis and for developing novel therapeutic approaches. Intercellular adhesion molecule 1 (ICAM1) is a cell surface glycoprotein expressed on endothelial cells and leukocytes in the immune system. The ICAM1 gene is located on chromosome 19p13 within the linkage region of diabetes. In the recent years, accumulating reports have implicated that genetic polymorphisms in the ICAM1 gene are associated with diabetes and diabetic nephropathy. Serum ICAM1 levels in diabetes patients and the icam1 gene expression in kidney tissues of diabetic animals are increased compared to the controls. Therefore, ICAM1 may play a role in the development of diabetes and diabetic nephropathy. In this review, we present genomic structure, variation, and regulation of the ICAM1 gene, summarized genetic and biological studies of this gene in diabetes and diabetic nephropathy and discussed about the potential application using ICAM1 as a biomarker and target for prediction and treatment of diabetes and diabetic nephropathy.