Search

MicroRNAs are conserved, endogenous small RNAs with critical post-transcriptional regulatory functions throughout eukaryota, including prominent roles in development and disease. Despite much effort, microRNA annotations still contain errors and are incomplete due especially to challenges related to identifying valid miRs that have small numbers of reads, to properly locating hairpin precursors and to balancing precision and recall. Here, we present miRWoods, which solves these challenges using a duplex-focused precursor detection method and stacked random forests with specialized layers to detect mature and precursor microRNAs, and has been tuned to optimize the harmonic mean of precision and recall. We trained and tuned our discovery pipeline on data sets from the well-annotated human genome, and evaluated its performance on data from mouse. Compared to existing approaches, miRWoods better identifies precursor spans, and can balance sensitivity and specificity for an overall greater prediction accuracy, recalling an average of 10% more annotated microRNAs, and correctly predicts substantially more microRNAs with only one read. We apply this method to the under-annotated genomes of Felis catus (domestic cat) and Bos taurus (cow). We identified hundreds of novel microRNAs in small RNA sequencing data sets from muscle and skin from cat, from 10 tissues from cow and also from human and mouse cells. Our novel predictions include a microRNA in an intron of tyrosine kinase 2 (TYK2) that is present in both cat and cow, as well as a family of mirtrons with two instances in the human genome. Our predictions support a more expanded miR-2284 family in the bovine genome, a larger mir-548 family in the human genome, and a larger let-7 family in the feline genome. [ABSTRACT FROM AUTHOR]

Cultivated watermelon (Citrullus lanatus) is one of the most important food crops in the Cucurbitaceae family. Diversification after domestication has led cultivated watermelons to exhibit diverse fruit flesh colors, including red, yellow, and orange. Recently, there has been increased interest in red-fleshed watermelons because they contain the antioxidant cis-isomeric lycopene. We performed whole genome resequencing (WGRS) of 24 watermelons with different flesh colors to identify single-nucleotide polymorphisms (SNPs) related to high lycopene content. The resequencing data revealed 203,894–279,412 SNPs from read mapping between inbred lines and the 97103 reference genome. In total, 295,065 filtered SNPs were identified, which had an average polymorphism information content of 0.297. Most of these SNPs were intergenic (90.1%) and possessed a transversion (Tv) rate of 31.64%. Overall, 2,369 SNPs were chosen at 0.5 Mb physical intervals to analyze genetic diversity across the 24 inbred lines. A neighbor-joining dendrogram and principal coordinate analysis (PCA) based on the 2,369 SNPs revealed that the 24 inbred lines could be grouped into high and low lycopene-type watermelons. In addition, we analyzed SNPs that could discriminate high lycopene content, red-fleshed watermelon from low lycopene, yellow or orange watermelon inbred lines. For validation, 19 SNPs (designated as WMHL1–19) were chosen randomly, and cleavage amplified polymorphic sequence (CAPS) markers were designed. Genotyping of the above 24 lines and 12 additional commercial cultivars using WMHL1–19 CAPS markers resulted in match rates of over 0.92 for most validated markers in correlation with the flesh color phenotypes. Our results provide valuable genomic information regarding the high lycopene content phenotype of red-fleshed cultivated watermelons, and the identified SNPs will be useful for the development of molecular markers in the marker-assisted breeding of watermelons with high lycopene content. [ABSTRACT FROM AUTHOR]

Two novel and related pathogenic variants of the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene were structurally and functionally characterized. These alterations have not been previously described in literature. Two patients with diagnosis of Cystic Fibrosis (CF) based on the presence of one mutated allele, p.Phe508del, pathological sweat test and clinical symptoms were studied. To complete the genotypes of both patients, an extensive genetic and functional analysis of the CFTR gene was performed. Extensive genetic characterization confirmed the presence of p.Phe508del pathogenic variant and revealed, in both patients, the presence of an insertion of part of intron 10 in intron 9 of the CFTR gene, within the (TG)m repeat, with a variable poly-T stretch. The molecular lesions resulted to be very similar in both patients, with only a difference in the number of T in the poly-T stretch. The functional characterization at RNA level revealed a complete anomalous splicing, without exon 10, from the allele with the insertion of both patients. Consequently, the alleles with the insertions are expected not to contribute to the formation of a functional CFTR protein. Molecular and functional features of these alterations are compatible with the definition of novel CF-causing variants of the CFTR gene. This also allowed the completion of the genetic characterization of both patients. [ABSTRACT FROM AUTHOR]

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5’ splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5’ splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing. [ABSTRACT FROM AUTHOR]

Plasmodium vivax invasion of human erythrocytes depends on the Duffy Binding Protein (PvDBP) which interacts with the Duffy antigen. PvDBP copy number has been recently shown to vary between P. vivax isolates in Sub-Saharan Africa. However, the extent of PvDBP copy number variation, the type of PvDBP multiplications, as well as its significance across broad samples are still unclear. We determined the prevalence and type of PvDBP duplications, as well as PvDBP copy number variation among 178 Ethiopian P. vivax isolates using a PCR-based diagnostic method, a novel quantitative real-time PCR assay and whole genome sequencing. For the 145 symptomatic samples, PvDBP duplications were detected in 95 isolates, of which 81 had the Cambodian and 14 Malagasy-type PvDBP duplications. PvDBP varied from 1 to >4 copies. Isolates with multiple PvDBP copies were found to be higher in symptomatic than asymptomatic infections. For the 33 asymptomatic samples, PvDBP was detected with two copies in two of the isolates, and both were the Cambodian-type PvDBP duplication. PvDBP copy number in Duffy-negative heterozygotes was not significantly different from that in Duffy-positives, providing no support for the hypothesis that increased copy number is a specific association with Duffy-negativity, although the number of Duffy-negatives was small and further sampling is required to test this association thoroughly. [ABSTRACT FROM AUTHOR]

Terpenes are responsible for most or all of the odor and flavor properties of Cannabis sativa, and may also impact effects users experience either directly or indirectly. We report the diversity of terpene profiles across samples bound for the Washington dispensary market. The remarkable degree of variation in terpene profiles ultimately results from action of a family of terpene synthase genes, only some of which have been described. Using a recently available genome assembly we describe 55 terpene synthases with genomic context, and tissue specific expression. The family is quite diverse from a protein similarity perspective, and subsets of the family are expressed in all tissues in the plant, including a set of root specific monoterpene synthases that could well have agronomic importance. Ultimately understanding and breeding for specific terpene profiles will require a good understanding of the gene family that underlies it. We intend for this work to serve as a foundation for that. [ABSTRACT FROM AUTHOR]

Background: Tyrosine kinase inhibitors targeted to anaplastic lymphoma kinase (ALK) have been demonstrated to be effective for lung cancer patients with an ALK fusion gene. Application of liquid biopsy, i.e., detection and quantitation of the fusion product in plasma cell-free DNA (cfDNA), could improve clinical practice. To detect ALK fusions, because fusion breakpoints occur somewhere in intron 19 of the ALK gene, sequencing of the entire intron is required to locate breakpoints. Results: We constructed a target sequencing system using an adapter and a set of primers that cover the entire ALK intron 19. This system can amplify fragments, including breakpoints, regardless of fusion partners. The data analysis pipeline firstly detected fusions by alignment to selected target sequences, and then quantitated the fusion alleles aligning to the identified breakpoint sequences. Performance was validated using 20 cfDNA samples from ALK-positive non-small cell lung cancer patients and samples from 10 healthy volunteers. Sensitivity and specificity were 50 and 100%, respectively. Conclusions: We demonstrated that PCR-based target sequencing using a tiling primer set and two-step mapping/alignment quantitatively detected ALK fusions in cfDNA from lung cancer patients. The system offers an alternative to existing approaches based on hybridization capture. [ABSTRACT FROM AUTHOR]

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms. [ABSTRACT FROM AUTHOR]

There is a need to investigate and better understand the inherited risk of cancer to ensure that clinical applications provide more accurate assessments and management strategies. Developing research-based next-generation sequencing gene panels that not only target (present-day) clinically actionable susceptibility genes but also genes that currently lack sufficient evidence for risk as well as candidate genes, such as those in DNA repair pathways, can help aid this effort. Therefore, gene panel B.O.P. (reast, varian, and rostate) was developed to evaluate the genetic risk of breast, ovarian and/or prostate cancer, and this manuscript serves as an introduction to B.O.P. and highlights its initial analytical validity assessment. B.O.P targets 87 genes that have been suggested, predicted, or clinically proven to be associated with breast, ovarian, and/or prostate cancer risk using Agilent Technologies Haloplex probes. The probes were designed for 100 base pair reads on an Illumina platform and target both coding and non-coding exons as well as 10 intronic base pairs flanking the intron-exon boundaries. The initial B.O.P screening involved 43 individuals from the Alabama Hereditary Cancer Cohort, and an average sequencing depth of 809X was obtained. Upon variant filtering and validation, true positives had an average sequencing depth of 659X and allele balance of 0.51. The average false positive sequencing depth was 34X and allele balance was 0.33. Although low sequencing depth was not always indicative of a false positive, high sequencing depths (>100X) signified a true positive. Furthermore, sensitivity and specificity of BRCA1/2 were calculated to be 100% and 92.3%, respectively. Overall, this screening enabled the establishment of criteria to alleviate future validation efforts and strongly supports the use of B.O.P. to further elucidate hereditary cancer susceptibility. Ultimately, continued B.O.P. screening will provide insights toward the genetic risk of and overlap between breast, ovarian, and/or prostate cancer. [ABSTRACT FROM AUTHOR]

Vesicoureteral reflux (VUR) is a complex, heritable disorder. Genome-wide linkage analyses of families affected by VUR have revealed multiple genomic loci linked to VUR. These loci normally harbor a number of genes whose biologically functional variant is yet to be identified. DNA copy number variations (CNVs) have not been extensively studied at high resolution in VUR patients. In this study, we performed array comparative genomic hybridization (aCGH) on a cohort of patients with a history of both VUR and urinary tract infection (UTI) with the objective of identifying genetic variations responsible for VUR and/or UTI susceptibility. UTI/VUR-associated CNVs were identified by aCGH results from the 192 Randomized Intervention for Children With Vesicoureteral Reflux (RIVUR) patients compared to 683 controls. Rare, large CNVs that are likely pathogenic and lead to VUR development were identified using stringent analysis criteria. Because UTI is a common affliction with multiple risk factors, we utilized standard analysis to identify potential disease-modifying CNVs that can contribute to UTI risk. Gene ontology analysis identified that CNVs in innate immunity and development genes were enriched in RIVUR patients. CNVs affecting innate immune genes may contribute to UTI susceptibility in VUR patients and may provide the first step in assisting clinical medicine in determining adverse outcome risk in children with VUR. [ABSTRACT FROM AUTHOR]

Adeno-associated viruses (AAV) are Dependoparvoviruses that have shown promise as recombinant vectors for gene therapy. While infectious pathways of AAV are well studied, gaps remain in our understanding of host factors affecting vector genome expression. Here, we map the role of ring finger protein 121 (RNF121), an E3 ubiquitin ligase, as a key regulator of AAV genome transcription. CRISPR-mediated knockout of RNF121 (RNF121 KO) in different cells markedly decreased AAV transduction regardless of capsid serotype or vector dose. Recombinant AAV transduction is partially rescued by overexpressing RNF121, but not by co-infection with helper Adenovirus. Major steps in the AAV infectious pathway including cell surface binding, cellular uptake, nuclear entry, capsid uncoating and second strand synthesis are unaffected. While gene expression from transfected plasmids or AAV genomes is unaffected, mRNA synthesis from AAV capsid-associated genomes is markedly decreased in RNF121 KO cells. These observations were attributed to transcriptional arrest as corroborated by RNAPol-ChIP and mRNA half-life measurements. Although AAV capsid proteins do not appear to be direct substrates of RNF121, the catalytic domain of the E3 ligase appears essential. Inhibition of ubiquitin-proteasome pathways revealed that blocking Valosin Containing Protein (VCP/p97), which targets substrates to the proteasome, can selectively and completely restore AAV-mediated transgene expression in RNF121 KO cells. Expanding on this finding, transcriptomic and proteomic analysis revealed that the catalytic subunit of DNA PK (DNAPK-Cs), a known activator of VCP, is upregulated in RNF121 KO cells and that the DNA damage machinery is enriched at sites of stalled AAV genome transcription. We postulate that a network of RNF121, VCP and DNA damage response elements function together to regulate transcriptional silencing and/or activation of AAV vector genomes. [ABSTRACT FROM AUTHOR]

Microvascular proliferation (MVP), an aberrant vascular structure containing multilayered mitotically active endothelial- and smooth-muscle cells/pericytes, is a histopathological hallmark of glioblastoma multiforme (GBM). Although MVP tends to be associated with high-grade glioma, it has also been detected in WHO grade I pilocytic astrocytoma (PA). However, little is known about the mechanism underlying its formation. Using TP53 point mutations as a marker for tumor-derived cells, we earlier reported that MVP was partially converted from tumor cells via mesenchymal transition. In the current study we used the KIAA1549-BRAF fusion gene as a marker to assess whether MVPs in PA contained tumor-derived cells and/or phenotypically distinct tumor cells expressing vascular markers. cDNA synthesized from frozen tissue of six PA patients operated at our institute was analyzed to detect the KIAA1549-BRAF fusion gene by reverse transcription polymerase chain reaction (RT-PCR) assay. The breakpoint in the fusion gene was identified by long and accurate PCR (LA-PCR) and Sanger sequencing of genomic DNA. Distinct tumor cells and cellular components of MVP were obtained by laser microdissection. For the qualitative and quantitative detection of the KIAA1549-BRAF fusion gene we performed genomic and digital PCR assays. Fluorescence in situ hybridization (FISH) was used to assess gene fusion in cellular components of MVP. Samples from three PA patients harbored the KIAA1549 exon 15, BRAF exon 9 fusion gene. In two patient samples with abundant MVP, RT-PCR assay detected strong bands arising from the KIAA1549-BRAF fusion gene in both tumor cells and cellular components of MVP. Digital PCR showed that vis-à-vis tumor tissue, its relative expression in cellular components of MVP was 42% in one- and 76% in another sample. FISH revealed amplified signals in both tumor cells and cellular components of MVP indicative of tandem duplication. Our findings suggest that in patients with PA, some cellular components of MVP contained tumor derived cell and/or phenotypically distinct tumor cells expressing vascular markers. [ABSTRACT FROM AUTHOR]

X-box binding protein 1 (XBP1) mRNA processing plays a crucial role in the unfolded protein response (UPR), which is activated in response to endoplasmic reticulum (ER) stress. Upon accumulation of the UPR-converted XBP1 mRNA splicing from an unspliced (u) XBP1 (inactive) isoform to the spliced (s) XBP1 (active) isoform, inositol-requiring enzyme 1 α (IRE1α) removes a 26-nucleotide intron from uXBP1 mRNA. Recent studies have reported the assessment of ER stress by examining the ratio of sXBP1 to uXBP1 mRNA (s/uXBP1 ratio) via densitometric analysis of PCR bands relative to increased levels of sXBP1 to uXBP1 using a housekeeping gene for normalization. However, this measurement is visualized by gel electrophoresis, making it very difficult to quantify differences between the two XBP1 bands and complicating data interpretation. Moreover, most commonly used housekeeping genes display an unacceptably high variable expression pattern of the s/uXBP1 ratio under different experimental conditions, such as various phases of development and different cell types, limiting their use as internal controls. For a more quantitative determination of XBP1 splicing activity, we measured the expression levels of total XBP1 (tXBP1: common region of s/uXBP1) and sXBP1 via real-time PCR using specific primer sets. We also designed universal real-time PCR primer sets capable of amplifying a portion of each u/s/tXBP1 mRNA that is highly conserved in eukaryotes, including humans, monkeys, cows, pigs, and mice. Therefore, we provide a more convenient and easily approachable quantitative real-time PCR method that can be used in various research fields to assess ER stress. [ABSTRACT FROM AUTHOR]

Sugarcane (Saccharum spp.) has a complex genome with variable ploidy and frequent aneuploidy, which hampers the understanding of phenotype and genotype relations. Despite this complexity, genome-wide association studies (GWAS) may be used to identify favorable alleles for target traits in core collections and then assist breeders in better managing crosses and selecting superior genotypes in breeding populations. Therefore, in the present study, we used a diversity panel of sugarcane, called the Brazilian Panel of Sugarcane Genotypes (BPSG), with the following objectives: (i) estimate, through a mixed model, the adjusted means and genetic parameters of the five yield traits evaluated over two harvest years; (ii) detect population structure, linkage disequilibrium (LD) and genetic diversity using simple sequence repeat (SSR) markers; (iii) perform GWAS analysis to identify marker-trait associations (MTAs); and iv) annotate the sequences giving rise to SSR markers that had fragments associated with target traits to search for putative candidate genes. The phenotypic data analysis showed that the broad-sense heritability values were above 0.48 and 0.49 for the first and second harvests, respectively. The set of 100 SSR markers produced 1,483 fragments, of which 99.5% were polymorphic. These SSR fragments were useful to estimate the most likely number of subpopulations, found to be four, and the LD in BPSG, which was stronger in the first 15 cM and present to a large extension (65 cM). Genetic diversity analysis showed that, in general, the clustering of accessions within the subpopulations was in accordance with the pedigree information. GWAS performed through a multilocus mixed model revealed 23 MTAs, six, three, seven, four and three for soluble solid content, stalk height, stalk number, stalk weight and cane yield traits, respectively. These MTAs may be validated in other populations to support sugarcane breeding programs with introgression of favorable alleles and marker-assisted selection. [ABSTRACT FROM AUTHOR]

Testicular or ovotesticular disorders of sex development (DSD) in individuals with female karyotype (XX) lacking the SRY gene has been observed in several mammalian species, including dogs. A genetic background for this abnormality has been extensively sought, and the region harboring the SOX9 gene has often been considered key in canine DSD. Three types of polymorphism have been studied in this region to date: a) copy number variation (CNV) in a region about 400 kb upstream of SOX9, named CNVR1; b) duplication of SOX9; and c) insertion of a single G-nucleotide (rs852549625) approximately 2.2 Mb upstream of SOX9. The aim of this study was thus to comprehensively analyze these polymorphisms in a large multibreed case-control cohort containing 45 XX DSD dogs, representing 23 breeds. The control set contained 57 fertile females. Droplet digital PCR (ddPCR) was used to study CNVR1 and the duplication of SOX9. Fluorescent in situ hybridization (FISH) was used to visualize copy numbers on a cellular level. The Sanger sequencing approach was performed to analyze the region harboring the G-insertion. We confirmed that CNVR1 is highly polymorphic and that copy numbers varied between 0 and 7 in the case and control cohorts. Interestingly, the number of copies was significantly higher (P = 0.038) in XX DSD dogs (mean = 2.7) than in the control females (mean = 2.0) but not in all studied breeds. Duplication of the SOX9 gene was noted only in a single XX DSD dog (an American Bully), which had three copies of SOX9. Distribution of the G-nucleotide insertion was similar in the XX DSD (frequency 0.20) and control (frequency 0.14) cohorts. Concluding, our study showed that CNVR1, located upstream of SOX9, is associated with the XX DSD phenotype, though in a breed-specific manner. Duplication of the SOX9 gene is a rare cause of this disorder in dogs. Moreover, we did not observe any association of G-insertion with the DSD phenotype. We assume that the genetic background of XX DSD can be different in certain breeds. [ABSTRACT FROM AUTHOR]

In contrast to human cells, very few HSV-1 genes are known to be spliced, although the same pre-mRNA processing machinery is shared. Here, through global analysis of splice junctions in cells infected with HSV-1 and an HSV-1 mutant virus with deletion of infectious cell culture protein 27 (ICP27), one of two viral immediate early (IE) genes essential for viral replication, we identify hundreds of novel alternative splice junctions mapping to both previously known HSV-1 spliced genes and previously unknown spliced genes, the majority of which alter the coding potential of viral genes. Quantitative and qualitative splicing efficiency analysis of these novel alternatively spliced genes based on RNA-Seq and RT-PCR reveals that splicing at these novel splice sites is efficient only when ICP27 is absent; while in wildtype HSV-1 infected cells, the splicing of these novel splice junctions is largely silenced in a gene/sequence specific manner, suggesting that ICP27 not only promotes accumulation of ICP27 targeted transcripts but also ensures correctness of the functional coding sequences through inhibition of alternative splicing. Furthermore, ICP27 toggles expression of ICP34.5, the major viral neurovirulence factor, through inhibition of splicing and activation of a proximal polyadenylation signal (PAS) in the newly identified intron, revealing a novel regulatory mechanism for expression of a viral gene. Thus, through the viral IE protein ICP27, HSV-1 co-opts both splicing and polyadenylation machinery to achieve optimal viral gene expression during lytic infection. On the other hand, during latent infection when ICP27 is absent, HSV-1 likely takes advantages of host splicing machinery to restrict expression of randomly activated antigenic viral genes to achieve immune evasion. [ABSTRACT FROM AUTHOR]

The human dopamine transporter gene SLC6A3 is involved in substance use disorders (SUDs) among many other common neuropsychiatric illnesses but allelic association results including those with its classic genetic markers 3’VNTR or Int8VNTR remain mixed and unexplainable. To better understand the genetics for reproducible association signals, we report the presence of recombination hotspots based on sequencing of the entire 5’ promoter regions in two small SUDs cohorts, 30 African Americans (AAs) and 30 European Americans (EAs). Recombination rate was the highest near the transcription start site (TSS) in both cohorts. In addition, each cohort carried 57 different promoter haplotypes out of 60 and no haplotypes were shared between the two ethnicities. A quarter of the haplotypes evolved in an ethnicity-specific manner. Finally, analysis of five hundred subjects of European ancestry, from the 1000 Genome Project, confirmed the promoter recombination hotspots and also revealed several additional ones in non-coding regions only. These findings provide an explanation for the mixed results as well as guidance for selection of effective markers to be used in next generation association validation (NGAV), facilitating the delineation of pathogenic variation in this critical neuropsychiatric gene. [ABSTRACT FROM AUTHOR]

The computational reconstruction of genome sequences from shotgun sequencing data has been greatly simplified by the advent of sequencing technologies that generate long reads. In the case of relatively small genomes (e.g., bacterial or viral), complete genome sequences can frequently be reconstructed computationally without the need for further experiments. However, large and complex genomes, such as those of most animals and plants, continue to pose significant challenges. In such genomes, assembly software produces incomplete and fragmented reconstructions that require additional experimentally derived information and manual intervention in order to reconstruct individual chromosome arms. Recent technologies originally designed to capture chromatin structure have been shown to effectively complement sequencing data, leading to much more contiguous reconstructions of genomes than previously possible. Here, we survey these technologies and the algorithms used to assemble and analyze large eukaryotic genomes, placed within the historical context of genome scaffolding technologies that have been in existence since the dawn of the genomic era. [ABSTRACT FROM AUTHOR]

Motivation: Recently, copy number variation (CNV) has gained considerable interest as a type of genomic variation that plays an important role in complex phenotypes and disease susceptibility. Since a number of CNV detection methods have recently been developed, it is necessary to help investigators choose suitable methods for CNV detection depending on their objectives. For this reason, this study compared ten commonly used CNV detection applications, including CNVnator, ReadDepth, RDXplorer, LUMPY and Control-FREEC, benchmarking the applications by sensitivity, specificity and computational demands. Taking the DGV gold standard variants as a standard dataset, we evaluated the ten applications with real sequencing data at sequencing depths from 5X to 50X. Among the ten methods benchmarked, LUMPY performs the best for both high sensitivity and specificity at each sequencing depth. For the purpose of high specificity, Canvas is also a good choice. If high sensitivity is preferred, CNVnator and RDXplorer are better choices. Additionally, CNVnator and GROM-RD perform well for low-depth sequencing data. Our results provide a comprehensive performance evaluation for these selected CNV detection methods and facilitate future development and improvement in CNV prediction methods. [ABSTRACT FROM AUTHOR]

The apolipoprotein E-C1-C4-C2 gene cluster at 19q13.32 encodes four amphipathic apolipoproteins. The influence of APOE common polymorphisms on plasma lipid/lipoprotein profile, especially on LDL-related traits, is well recognized; however, little is known about the role of other genes/variants in this gene cluster. In this study, we evaluated the role of common and uncommon/rare genetic variation in this gene region on inter-individual variation in plasma lipoprotein levels in non-Hispanic Whites (NHWs) and African blacks (ABs). In the variant discovery step, the APOE, APOC1, APOC4, APOC2 genes were sequenced along with their flanking and hepatic control regions (HCR1 and HCR2) in 190 subjects with extreme HDL-C/TG levels. The next step involved the genotyping of 623 NHWs and 788 ABs for the identified uncommon/rare variants and common tagSNPs along with additional relevant SNPs selected from public resources, followed by association analyses with lipid traits. A total of 230 sequence variants, including 15 indels were identified, of which 65 were novel. A total of 70 QC-passed variants in NHWs and 108 QC-passed variants in ABs were included in the final association analyses. Single-site association analysis of SNPs with MAF>1% revealed 20 variants in NHWs and 24 variants in ABs showing evidence of association with at least one lipid trait, including several variants exhibiting independent associations from the established APOE polymorphism even after multiple-testing correction. Overall, our study has confirmed known associations and also identified novel associations in this genomic region with various lipid traits. Our data also support the contribution of both common and uncommon/rare variation in this gene region in affecting plasma lipid profile in the general population. [ABSTRACT FROM AUTHOR]

Multiscale signal processing techniques such as wavelet filtering have proved to be particularly successful in predicting exon sequences. Traditional wavelet predictor is domain filtering, and enforces exon features by weighting nucleotide values with coefficients. Such a measure performs linear filtering and is not suitable for preserving the short coding exons and the exon-intron boundaries. This paper describes a prediction framework that is capable of non-linearly processing DNA sequences while achieving high prediction rates. There are two key contributions. The first is the introduction of a genomic-inspired multiscale bilateral filtering (MSBF) which exploits both weighting coefficients in the spatial domain and nucleotide similarity in the range. Similarly to wavelet transform, the MSBF is also defined as a weighted sum of nucleotides. The difference is that the MSBF takes into account the variation of nucleotides at a specific codon position. The second contribution is the exploitation of inter-scale correlation in MSBF domain to find the inter-scale dependency on the differences between the exon signal and the background noise. This favourite property is used to sharp the important structures while weakening noise. Three benchmark data sets have been used in the evaluation of considered methods. By comparison with four existing techniques, the prediction results demonstrate that: the proposed method reveals at least improvement of 4.1%, 50.5%, 25.6%, 2.5%, 10.8%, 15.5%, 11.1%, 12.3%, 9.2% and 2.4% on the exons length of 1–24, 25–49, 50–74, 75–99, 100–124, 125–149, 150–174, 175–199, 200–299 and 300–300+, respectively. The MSBF of its nonlinear nature is good at energy compaction, which makes it capable of locating the sharp variations around short exons. The direct scale multiplication of coefficients at several adjacent scales obviously enhanced exon features while the noise contents were suppressed. We show that the non-linear nature and correlation-based property achieved in proposed predictor is greater than that for traditional filtering, which leads to better exon prediction performance. There are some possible applications of this predictor. Its good localization and protection of sharp variations will make the predictor be suitable to perform fault diagnosis of aero-engine. [ABSTRACT FROM AUTHOR]

Mast cell tumours are the most common type of skin cancer in dogs, representing a significant concern in canine health. The molecular pathogenesis is largely unknown, but breed-predisposition for mast cell tumour development suggests the involvement of inherited genetic risk factors in some breeds. In this study, we aimed to identify germline risk factors associated with the development of mast cell tumours in Labrador Retrievers, a breed with an elevated risk of mast cell tumour development. Using a methodological approach that combined a genome-wide association study, targeted next generation sequencing, and TaqMan genotyping, we identified a synonymous variant in the DSCAM gene on canine chromosome 31 that is associated with mast cell tumours in Labrador Retrievers. DSCAM encodes a cell-adhesion molecule. We showed that the variant has no effect on the DSCAM mRNA level but is associated with a significant reduction in the level of the DSCAM protein, suggesting that the variant affects the dynamics of DSCAM mRNA translation. Furthermore, we showed that the variant is also associated with mast cell tumours in Golden Retrievers, a breed that is closely related to Labrador Retrievers and that also has a predilection for mast cell tumour development. The variant is common in both Labradors and Golden Retrievers and consequently is likely to be a significant genetic contributor to the increased susceptibility of both breeds to develop mast cell tumours. The results presented here not only represent an important contribution to the understanding of mast cell tumour development in dogs, as they highlight the role of cell adhesion in mast cell tumour tumourigenesis, but they also emphasise the potential importance of the effects of synonymous variants in complex diseases such as cancer. [ABSTRACT FROM AUTHOR]

Arabis alpina is a perennial arctic-alpine plant and an upcoming model organism for genetics and molecular biology for the Brassicaceae family. One essential method for most molecular approaches is the analysis of gene expression by reverse-transcription quantitative Real-Time PCR (RT-qPCR). For the normalisation of expression data in RT-qPCR experiments, it is essential to use reliable reference genes that are not affected under a wide range of conditions. In this study we establish a set of 15 A. alpina reference genes that were tested under different conditions including cold, drought, heat, salt and gibberellic acid treatments. Data analyses with geNORM, BestKeeper and NormFinder revealed the most stable reference genes for the tested conditions: RAN3, HCF and PSB33 are most suitable for cold treatments; UBQ10 and TUA5 for drought; RAN3, PSB33 and EIF4a for heat; CAC, TUA5, ACTIN 2 and PSB33 for salt and PSB33 and TUA5 for gibberellic acid treatments. CAC and ACTIN 2 showed the least variation over all tested samples. In addition, we show that two reference genes are sufficient to normalize RT-qPCR data under our treatment conditions. In future studies, these reference genes can be used for an adequate normalisation and thus help to generate high quality RT-qPCR data in A. alpina. [ABSTRACT FROM AUTHOR]

Transmission ratio distortion (TRD) by the mouse t-haplotype, a variant region on chromosome 17, is a well-studied model of non-Mendelian inheritance. It is characterized by the high transmission ratio (up to 99%) of the t-haplotype from t/+ males to their offspring. TRD is achieved by the exquisite ability of the responder (Tcr) to trigger non-Mendelian inheritance of homologous chromosomes. Several distorters (Tcd1-Tcd4), which act cumulatively, together promote the high transmission ratio of Tcr and the t-haplotype. Molecularly, TRD is brought about by deregulation of Rho signaling pathways via the distorter products, which impair sperm motility, and the t-sperm specific rescue of sperm motility by the responder. The t-sperm thus can reach the egg cells faster than +-sperm and fertilize them. Previously we have shown that the responder function is accomplished by a dominant negative form of sperm motility kinase (SMOKTCR), while the distorter functions are accomplished by the Rho G protein regulators TAGAP, FGD2 and NME3 proposed to function in two oppositely acting pathways. Here we identify the RAC1-specific guanine nucleotide exchange factor TIAM2 as modifier of t-haplotype TRD. Tiam2 is expressed in two isoforms, the full-length (Tiam2l) and a short transcript (Tiam2s). Tiam2s expression from the t-allele is strongly increased compared to the wild-type allele. By transgenic approaches we show that Tiam2s enhances t-haplotype transmission, while Tiam2l has the opposite effect. Our data show that a single modifier locus can encode different gene products exerting opposite effects on a trait. They also suggest that the expression ratio of the isoforms determines if the outcome is an enhancing or a suppressive effect on the trait. [ABSTRACT FROM AUTHOR]

Complex chromosomal rearrangements (CCRs) are rearrangements involving more than two chromosomes or more than two breakpoints. Whole genome sequencing (WGS) allows for outstanding high resolution characterization on the nucleotide level in unique sequences of such rearrangements, but problems remain for mapping breakpoints in repetitive regions of the genome, which are known to be prone to rearrangements. Hence, multiple complementary WGS experiments are sometimes needed to solve the structures of CCRs. We have studied three individuals with CCRs: Case 1 and Case 2 presented with de novo karyotypically balanced, complex interchromosomal rearrangements (46,XX,t(2;8;15)(q35;q24.1;q22) and 46,XY,t(1;10;5)(q32;p12;q31)), and Case 3 presented with a de novo, extremely complex intrachromosomal rearrangement on chromosome 1. Molecular cytogenetic investigation revealed cryptic deletions in the breakpoints of chromosome 2 and 8 in Case 1, and on chromosome 10 in Case 2, explaining their clinical symptoms. In Case 3, 26 breakpoints were identified using WGS, disrupting five known disease genes. All rearrangements were subsequently analyzed using optical maps, linked-read WGS, and short-read WGS. In conclusion, we present a case series of three unique de novo CCRs where we by combining the results from the different technologies fully solved the structure of each rearrangement. The power in combining short-read WGS with long-molecule sequencing or optical mapping in these unique de novo CCRs in a clinical setting is demonstrated. [ABSTRACT FROM AUTHOR]

Introns can be extraordinarily large and they account for the majority of the DNA sequence in human genes. However, little is known about their population patterns of structural variation and their functional implication. By combining the most extensive maps of CNVs in human populations, we have found that intronic losses are the most frequent copy number variants (CNVs) in protein-coding genes in human, with 12,986 intronic deletions, affecting 4,147 genes (including 1,154 essential genes and 1,638 disease-related genes). This intronic length variation results in dozens of genes showing extreme population variability in size, with 40 genes with 10 or more different sizes and up to 150 allelic sizes. Intronic losses are frequent in evolutionarily ancient genes that are highly conserved at the protein sequence level. This result contrasts with losses overlapping exons, which are observed less often than expected by chance and almost exclusively affect primate-specific genes. An integrated analysis of CNVs and RNA-seq data showed that intronic loss can be associated with significant differences in gene expression levels in the population (CNV-eQTLs). These intronic CNV-eQTLs regions are enriched for intronic enhancers and can be associated with expression differences of other genes showing long distance intron-promoter 3D interactions. Our data suggests that intronic structural variation of protein-coding genes makes an important contribution to the variability of gene expression and splicing in human populations. [ABSTRACT FROM AUTHOR]

Horizontal Gene Transfer was long thought to be marginal in Mycoplasma a large group of wall-less bacteria often portrayed as minimal cells because of their reduced genomes (ca. 0.5 to 2.0 Mb) and their limited metabolic pathways. This view was recently challenged by the discovery of conjugative exchanges of large chromosomal fragments that equally affected all parts of the chromosome via an unconventional mechanism, so that the whole mycoplasma genome is potentially mobile. By combining next generation sequencing to classical mating and evolutionary experiments, the current study further explored the contribution and impact of this phenomenon on mycoplasma evolution and adaptation using the fluoroquinolone enrofloxacin (Enro), for selective pressure and the ruminant pathogen Mycoplasma agalactiae, as a model organism. For this purpose, we generated isogenic lineages that displayed different combination of spontaneous mutations in Enro target genes (gyrA, gyrB, parC and parE) in association to gradual level of resistance to Enro. We then tested whether these mutations can be acquired by a susceptible population via conjugative chromosomal transfer knowing that, in our model organism, the 4 target genes are scattered in three distinct and distant loci. Our data shows that under antibiotic selective pressure, the time scale of the mutational pathway leading to high-level of Enro resistance can be readily compressed into a single conjugative step, in which several EnroR alleles were transferred from resistant to susceptible mycoplasma cells. In addition to acting as an accelerator for antimicrobial dissemination, mycoplasma chromosomal transfer reshuffled genomes beyond expectations and created a mosaic of resistant sub-populations with unpredicted and unrelated features. Our findings provide insights into the process that may drive evolution and adaptability of several pathogenic Mycoplasma spp. via an unconventional conjugative mechanism. [ABSTRACT FROM AUTHOR]

Background: The submission of DNA sequences to public sequence databases is an essential, but insufficiently automated step in the process of generating and disseminating novel DNA sequence data. Despite the centrality of database submissions to biological research, the range of available software tools that facilitate the preparation of sequence data for database submissions is low, especially for sequences generated via plant and fungal DNA barcoding. Current submission procedures can be complex and prohibitively time expensive for any but a small number of input sequences. A user-friendly software tool is needed that streamlines the file preparation for database submissions of DNA sequences that are commonly generated in plant and fungal DNA barcoding. Methods: A Python package was developed that converts DNA sequences from the common EMBL and GenBank flat file formats to submission-ready, tab-delimited spreadsheets (so-called ‘checklists’) for a subsequent upload to the annotated sequence section of the European Nucleotide Archive (ENA). The software tool, titled ‘EMBL2checklists’, automatically converts DNA sequences, their annotation features, and associated metadata into the idiosyncratic format of marker-specific ENA checklists and, thus, generates files that can be uploaded via the interactive Webin submission system of ENA. Results: EMBL2checklists provides a simple, platform-independent tool that automates the conversion of common DNA barcoding sequences into easily editable spreadsheets that require no further processing but their upload to ENA via the interactive Webin submission system. The software is equipped with an intuitive graphical as well as an efficient command-line interface for its operation. The utility of the software is illustrated by its application in four recent investigations, including plant phylogenetic and fungal metagenomic studies. Discussion: EMBL2checklists bridges the gap between common software suites for DNA sequence assembly and annotation and the interactive data submission process of ENA. It represents an easy-to-use solution for plant and fungal biologists without bioinformatics expertise to generate submission-ready checklists from common DNA sequence data. It allows the post-processing of checklists as well as work-sharing during the submission process and solves a critical bottleneck in the effort to increase participation in public data sharing. [ABSTRACT FROM AUTHOR]

Human Multipotent Stromal Cells (MSCs) are a valuable resource for regenerative medicine and are widely studied. They can be isolated from a variety of tissues and differentiate into multiple cell types (multi-potent). Many reports have been published using human MSCs and to be able to compare outcome, or be able to identify differences between MSCs, several cell surface markers have been proposed. Nevertheless, still many differences remain. Gene expression is known to be different between cell stage and origin. Furthermore, cells cultured on a culture dish (2D) show different gene expression profiles as compared to cells grown on scaffolds (3D). Even the RNA extraction method and the selection of genes used for normalisation have a role in gene expression profiling. To be able to compare gene expression data from samples cultured in different dimensions and RNA extracted using a variety of protocols we set out to define a set of reference genes suitable to normalise qPCR data from a very heterogeneous sample set. Hereto, Trizol was used to extract RNA from human MSCs cultured in 3D and 2D to validate newly designed and previously published primer sets. Subsequently, RNA from fresh human MSC samples and samples stored in RLT-buffer, Trizol or RNAlater was extracted using RNeasy and Trizol methods. All samples have been used to rank the candidate reference genes according to their stability after qPCR enabling identification of the most suitable reference gene(s) for normalisation of a heterogeneous sample set. The most stably expressed reference genes indicated superior normalisation of MSC marker gene expression over the least stable reference genes. [ABSTRACT FROM AUTHOR]

Transcriptomic profiling of complex tissues by single-nucleus RNA-sequencing (snRNA-seq) affords some advantages over single-cell RNA-sequencing (scRNA-seq). snRNA-seq provides less biased cellular coverage, does not appear to suffer cell isolation-based transcriptional artifacts, and can be applied to archived frozen specimens. We used well-matched snRNA-seq and scRNA-seq datasets from mouse visual cortex to compare cell type detection. Although more transcripts are detected in individual whole cells (~11,000 genes) than nuclei (~7,000 genes), we demonstrate that closely related neuronal cell types can be similarly discriminated with both methods if intronic sequences are included in snRNA-seq analysis. We estimate that the nuclear proportion of total cellular mRNA varies from 20% to over 50% for large and small pyramidal neurons, respectively. Together, these results illustrate the high information content of nuclear RNA for characterization of cellular diversity in brain tissues. [ABSTRACT FROM AUTHOR]

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1’s proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102–104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18–21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution. [ABSTRACT FROM AUTHOR]

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component and etiology characterized by chronic inflammation and autoantibody production. The purpose of this study was to ascertain copy number variation (CNV) in SLE using a case-control design in an admixed Brazilian population. The whole-genome detection of CNV was performed using Cytoscan HD array in SLE patients and healthy controls. The best CNV candidates were then evaluated by quantitative real-time PCR in a larger cohort or validated using droplet digital PCR. Logistic regression models adjusted for sex and ancestry covariates was applied to evaluate the association between CNV with SLE susceptibility. The data showed a synergistic effect between the FCGR3B and ADAM3A loci with the presence of deletions in both loci significantly increasing the risk to SLE (5.9-fold) compared to the deletion in the single FCGR3B locus (3.6-fold). In addition, duplications in these genes were indeed more frequent in healthy subjects, suggesting that high FCGR3B/ADAM3A gene copy numbers are protective factors against to disease development. Overall, 21 rare CNVs were identified in SLE patients using a four-step pipeline created for identification of rare variants. Furthermore, heterozygous deletions overlapping the CFHR4, CFHR5 and HLA-DPB2 genes were described for the first time in SLE patients. Here we present the first genome-wide CNV study of SLE patients in a tri-hybrid population. The results show that novel susceptibility loci to SLE can be found once the distribution of structural variants is analyzed throughout the whole genome. [ABSTRACT FROM AUTHOR]

Group II introns are ancient retroelements that significantly shaped the origin and evolution of contemporary eukaryotic genomes. These self-splicing ribozymes share a common ancestor with the telomerase enzyme, the spliceosome machinery as well as the highly abundant spliceosomal introns and non-LTR retroelements. More than half of the human genome thus consists of various elements that evolved from ancient group II introns, which altogether significantly contribute to key functions and genetic diversity in eukaryotes. Similarly, group II intron-related elements in bacteria such as abortive phage infection (Abi) retroelements, diversity generating retroelements (DGRs) and some CRISPR-Cas systems have evolved to confer important functions to their hosts. In sharp contrast, since bacterial group II introns are scarce, irregularly distributed and frequently spread by lateral transfer, they have mainly been considered as selfish retromobile elements with no beneficial function to their host. Here we unveil a new group II intron function that generates genetic diversity at the RNA level in bacterial cells. We demonstrate that Ll.LtrB, the model group II intron from Lactococcus lactis, recognizes specific sequence motifs within cellular mRNAs by base pairing, and invades them by reverse splicing. Subsequent splicing of ectopically inserted Ll.LtrB, through circularization, induces a novel trans-splicing pathway that generates exon 1-mRNA and mRNA-mRNA intergenic chimeras. Our data also show that recognition of upstream alternative circularization sites on intron-interrupted mRNAs release Ll.LtrB circles harboring mRNA fragments of various lengths at their splice junction. Intergenic trans-splicing and alternative circularization both produce novel group II intron splicing products with potential new functions. Overall, this work describes new splicing pathways in bacteria that generate, similarly to the spliceosome in eukaryotes, genetic diversity at the RNA level while providing additional functional and evolutionary links between group II introns, spliceosomal introns and the spliceosome. [ABSTRACT FROM AUTHOR]

CFTR modulators have revolutionized the treatment of individuals with cystic fibrosis (CF) by improving the function of existing protein. Unfortunately, almost half of the disease-causing variants in CFTR are predicted to introduce premature termination codons (PTC) thereby causing absence of full-length CFTR protein. We hypothesized that a subset of nonsense and frameshift variants in CFTR allow expression of truncated protein that might respond to FDA-approved CFTR modulators. To address this concept, we selected 26 PTC-generating variants from four regions of CFTR and determined their consequences on CFTR mRNA, protein and function using intron-containing minigenes expressed in 3 cell lines (HEK293, MDCK and CFBE41o-) and patient-derived conditionally reprogrammed primary nasal epithelial cells. The PTC-generating variants fell into five groups based on RNA and protein effects. Group A (reduced mRNA, immature (core glycosylated) protein, function <1% (n = 5)) and Group B (normal mRNA, immature protein, function <1% (n = 10)) variants were unresponsive to modulator treatment. However, Group C (normal mRNA, mature (fully glycosylated) protein, function >1% (n = 5)), Group D (reduced mRNA, mature protein, function >1% (n = 5)) and Group E (aberrant RNA splicing, mature protein, function > 1% (n = 1)) variants responded to modulators. Increasing mRNA level by inhibition of NMD led to a significant amplification of modulator effect upon a Group D variant while response of a Group A variant was unaltered. Our work shows that PTC-generating variants should not be generalized as genetic ‘nulls’ as some may allow generation of protein that can be targeted to achieve clinical benefit. [ABSTRACT FROM AUTHOR]

Root-knot nematodes (RKNs, Meloidogyne spp.) are destructive plant parasites with a wide host range. They severely reduce crop quality and yield worldwide. Tobacco is a versatile model plant organism for studying RKNs-host interactions and a key plant material for molecular research. Long noncoding RNAs (lncRNAs) play critical roles in post transcriptional and transcriptional regulation in a wide range of biological pathways, especially plant development and stress response. In the present study, we obtained 5,206 high-confidence lncRNAs based on RNA sequencing data. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that the target genes of these lncRNAs are mainly involved in plant biotic and abiotic stresses, plant hormone signal transduction, induced systemic resistance, plant-type hypersensitive response, plant-type cell wall organization or biogenesis. The 565 differentially expressed lncRNAs found to be involved in nematode stress response were validated by quantitative PCR using 15 randomly-selected lncRNA genes. Our study provides insights into the molecular mechanisms of RKNs-plant interactions that might help preventing nematode damages to crops. [ABSTRACT FROM AUTHOR]

Many eukaryotic genes play essential roles in multiple biological processes in several different tissues. Conditional mutants are needed to analyze genes with such pleiotropic functions. In vertebrates, conditional gene inactivation has only been feasible in the mouse, leaving other model systems to rely on surrogate experimental approaches such as overexpression of dominant negative proteins and antisense-based tools. Here, we have developed a simple and straightforward method to integrate loxP sequences at specific sites in the zebrafish genome using the CRISPR/Cas9 technology and oligonucleotide templates for homology directed repair. We engineered conditional (floxed) mutants of tbx20 and fleer, and demonstrate excision of exons flanked by loxP sites using tamoxifen-inducible CreERT2 recombinase. To demonstrate broad applicability of our method, we also integrated loxP sites into two additional genes, aldh1a2 and tcf21. The ease of this approach will further expand the use of zebrafish to study various aspects of vertebrate biology, especially post-embryonic processes such as regeneration. [ABSTRACT FROM AUTHOR]

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both “simple” and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements. [ABSTRACT FROM AUTHOR]

A nonsense variant (p.W358X) of human pancreatic lipase related protein 2 (PNLIPRP2) is present in different ethnic populations with a high allele frequency. In cell culture experiments, the truncated protein mainly accumulates inside the cells and causes endoplasmic reticulum stress. Here, we tested the hypothesis that variant p.W358X might increase risk for chronic pancreatitis through acinar cell stress. We sequenced exon 11 of PNLIPRP2 in a cohort of 256 subjects with chronic pancreatitis (152 alcoholic and 104 non-alcoholic) and 200 controls of Hungarian origin. We observed no significant difference in the distribution of the truncation variant between patients and controls. We analyzed mRNA expression in human pancreatic cDNA samples and found the variant allele markedly reduced. We conclude that the p.W358X truncation variant of PNLIPRP2 is expressed poorly and has no significant effect on the risk of chronic pancreatitis. [ABSTRACT FROM AUTHOR]

The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most serious insect pest species to evolve resistance against many insecticides from different chemical classes. This species has evolved resistance to the pyrethroid insecticides across its native range and is becoming a truly global pest after establishing in South America and having been recently recorded in North America. A chimeric cytochrome P450 gene, CYP337B3, has been identified as a resistance mechanism for resistance to fenvalerate and cypermethrin. Here we show that this resistance mechanism is common around the world with at least eight different alleles. It is present in South America and has probably introgressed into its closely related native sibling species, Helicoverpa zea. The different alleles of CYP337B3 are likely to have arisen independently in different geographic locations from selection on existing diversity. The alleles found in Brazil are those most commonly found in Asia, suggesting a potential origin for the incursion of H. armigera into the Americas. [ABSTRACT FROM AUTHOR]

During lytic replication of Kaposi’s sarcoma-associated herpesvirus (KSHV), a nuclear viral long noncoding RNA known as PAN RNA becomes the most abundant polyadenylated transcript in the cell. Knockout or knockdown of KSHV PAN RNA results in loss of late lytic viral gene expression and, consequently, reduction of progeny virion release from the cell. Here, we demonstrate that knockdown of PAN RNA from the related Rhesus macaque rhadinovirus (RRV) phenocopies that of KSHV PAN RNA. These two PAN RNA homologs, although lacking significant nucleotide sequence conservation, can functionally substitute for each other to rescue phenotypes associated with the absence of PAN RNA expression. Because PAN RNA is exclusively nuclear, previous studies suggested that it directly interacts with host and viral chromatin to modulate gene expression. We studied KSHV and RRV PAN RNA homologs using capture hybridization analysis of RNA targets (CHART) and observed their associations with host chromatin, but the loci differ between PAN RNA homologs. Accordingly, we find that KSHV PAN RNA is undetectable in chromatin following cell fractionation. Thus, modulation of gene expression at specific chromatin loci appears not to be the primary, nor the pertinent function of this viral long noncoding RNA. PAN RNA represents a cautionary tale for the investigation of RNA association with chromatin whereby cross-linking of DNA spatially adjacent to an abundant nuclear RNA gives the appearance of specific interactions. Similarly, PAN RNA expression does not affect viral transcription factor complex expression or activity, which is required for generation of the late lytic viral mRNAs. Rather, we provide evidence for an alternative model of PAN RNA function whereby knockdown of KSHV or RRV PAN RNA results in compromised nuclear mRNA export thereby reducing the cytoplasmic levels of viral mRNAs available for production of late lytic viral protein. [ABSTRACT FROM AUTHOR]

Mutations in interleukin-10 receptor (IL-10R) genes are one cause of very early-onset inflammatory bowel disease with perianal lesions, which can be cured by hematopoietic stem cell transplantation. Using a functional test, which assesses responsiveness of peripheral monocytes to IL-10, we identified three unrelated Portuguese patients carrying two novel IL-10RB mutations. In the three patients, sequencing of genomic DNA identified the same large deletion of exon 3 which precluded protein expression. This mutation was homozygous in two patients born from consanguineous families and heterozygous in the third patient born from unrelated parents. Microsatellite analysis of the IL10RB genomic region revealed a common haplotype in the three Portuguese families pointing to a founder deletion inherited from a common ancestor 400 years ago. In the third patient, surface expression of IL-10R was normal but signaling in response to IL-10 was impaired. Complementary DNA sequencing and next-generation sequencing of IL10RB locus with custom-made probes revealed a ≈ 6 Kb duplication encompassing the exon 6 which leads to a frameshift mutation and a loss of the TYK2-interacting Box 2 motif. Altogether, we describe two novel copy number variations in IL10RB, one with founder effect and one preserving cell surface expression but abolishing signaling. [ABSTRACT FROM AUTHOR]

Nuclear genes of euglenids and marine diplonemids harbor atypical, nonconventional introns which are not observed in the genomes of other eukaryotes. Nonconventional introns do not have the conserved borders characteristic for spliceosomal introns or the sequence complementary to U1 snRNA at the 5' end. They form a stable secondary structure bringing together both exon/intron junctions, nevertheless, this conformation does not resemble the form of self-splicing or tRNA introns. In the genes studied so far, frequent nonconventional introns insertions at new positions have been observed, whereas conventional introns have been either found at the conserved positions, or simply lost. In this work, we examined the order of intron removal from Euglena gracilis transcripts of the tubA and gapC genes, which contain two types of introns: nonconventional and spliceosomal. The relative order of intron excision was compared for pairs of introns belonging to different types. Furthermore, intermediate products of splicing were analyzed using the PacBio Next Generation Sequencing system. The analysis led to the main conclusion that nonconventional introns are removed in a rapid way but later than spliceosomal introns. Moreover, the observed accumulation of transcripts with conventional introns removed and nonconventional present may suggest the existence of a time gap between the two types of splicing. [ABSTRACT FROM AUTHOR]

The IGFN1 (Immunoglobulin-Like And Fibronectin Type III Domain Containing 1) gene has a role in skeletal muscle function and is also involved in metastatic breast cancer, and the isoforms with three N-terminal globular domains are sufficient for its function in skeletal muscle. Two novel splicing isoforms of IGFN1 have been identified in renal cell carcinoma (RCC), one with 5’exon extension and an isoform with a novel exon. The role of G-quadruplex, a non-B DNA, was explored for the splicing alteration of IGFN1 in RCC. G-quadruplexes are the secondary structures acquired by stacking of G-quartets by Hoogsteen hydrogen bonding in DNA and RNA. IGFN1 has intronic potential G-quadruplex forming sequence (PQS) folding into G-quadruplex and is studied for its involvement in aberrant splicing. A PQS in the intron 15 of IGFN1 gene was observed in our in silico analysis by QGRS mapper and non BdB web servers. We observed PQS folds into stable G-quadruplex structure in gel shift assay and circular dichroism (CD) spectroscopy in the presence of G-quadruplex stabilizing agents Pyridostatin (PDS) and KCl, respectively. G-quadruplex formation site with single base resolution was mapped by Sanger sequencing of the plasmid constructs harbouring the cloned PQS and its mutant. This stable G-quadruplex inhibits reverse transcriptase and taq polymerase in reverse transcriptase & PCR stop assays. PDS changes the different splicing isoforms of IGFN1 in UOK146 cell line, displaying involvement of intronic G-quadruplex in IGFN1 splicing. These results lead us to propose that a stable G-quadruplex structure is formed in IGFN1 intron and a reason behind IGFN1 aberrant splicing which could be targeted for therapeutic intervention. [ABSTRACT FROM AUTHOR]

In plant breeding the use of molecular markers has resulted in tremendous improvement of the speed with which new crop varieties are introduced into the market. Single Nucleotide Polymorphism (SNP) genotyping is routinely used for association studies, Linkage Disequilibrium (LD) and Quantitative Trait Locus (QTL) mapping studies, marker-assisted backcrosses and validation of large numbers of novel SNPs. Here we present the KeyGene SNPSelect technology, a scalable and flexible multiplexed, targeted sequence-based, genotyping solution. The multiplex composition of SNPSelect assays can be easily changed between experiments by adding or removing loci, demonstrating their content flexibility. To demonstrate this versatility, we first designed a 1,056-plex maize assay and genotyped a total of 374 samples originating from an F2 and a Recombinant Inbred Line (RIL) population and a maize germplasm collection. Next, subsets of the most informative SNP loci were assembled in 384-plex and 768-plex assays for further genotyping. Indeed, selection of the most informative SNPs allows cost-efficient yet highly informative genotyping in a custom-made fashion, with average call rates between 88.1% (1,056-plex assay) and 99.4% (384-plex assay), and average reproducibility rates between duplicate samples ranging from 98.2% (1056-plex assay) to 99.9% (384-plex assay). The SNPSelect workflow can be completed from a DNA sample to a genotype dataset in less than three days. We propose SNPSelect as an attractive and competitive genotyping solution to meet the targeted genotyping needs in fields such as plant breeding. [ABSTRACT FROM AUTHOR]

The Prenatal Assessment of Genome and Exomes (PAGE) project is a UK-wide study aiming to gain a better understanding of genetic variants causing developmental problems during pregnancy. A further aim of the study is to provide an evidence-base for the introduction of prenatal whole genome and exome sequencing (PWGES) into prenatal diagnostics provided by the NHS, which is expected in 2018. This paper presents the findings of a qualitative interview study undertaken with 20 health professionals and researchers involved in the PAGE project, and explores their implications for understandings of ‘good practice’ in the uses of prenatal genomics clinically. A number of critical issues are identified that will need to be addressed in the development of a model of good ethical practice for prenatal genomics: consent, management of expectations, return of results, and professional duties in the context of PWGES. The analysis presented identifies and illustrates a great deal of complexity and qualitative richness in these issues as they arise in the day-to-day work of genomics professionals. Inclusive, critical discussion of these findings, together with the findings from other empirical studies, normative analysis and scientific discoveries resulting from PAGE, will be required to inform the development of appropriate guidelines of good ethical practice that address the needs and concerns to be encountered in daily clinical practice. [ABSTRACT FROM AUTHOR]

Sexually deceptive orchids of the genus Ophrys may rapidly evolve by adaptation to pollinators. However, understanding of the genetic basis of potential changes and patterns of relationships is hampered by a lack of genomic information. We report the complete plastid genome sequences of Ophrys iricolor and O. sphegodes, representing the two most species-rich lineages of the genus Ophrys. Both plastomes are circular DNA molecules (146754 bp for O. sphegodes and 150177 bp for O. iricolor) with the typical quadripartite structure of plastid genomes and within the average size of photosynthetic orchids. 213 Simple Sequence Repeats (SSRs) (31.5% polymorphic between O. iricolor and O. sphegodes) were identified, with homopolymers and dipolymers as the most common repeat types. SSRs were mainly located in intergenic regions but SSRs located in coding regions were also found, mainly in ycf1 and rpoC2 genes. The Ophrys plastome is predicted to encode 107 distinct genes, 17 of which are completely duplicated in the Inverted Repeat regions. 83 and 87 putative RNA editing sites were detected in 25 plastid genes of the two Ophrys species, all occurring in the first or second codon position. Comparing the rate of nonsynonymous (dN) and synonymous (dS) substitutions, 24 genes (including rbcL and ycf1) display signature consistent with positive selection. When compared with other members of the orchid family, the Ophrys plastome has a complete set of 11 functional ndh plastid genes, with the exception of O. sphegodes that has a truncated ndhF gene. Comparative analysis showed a large co-linearity with other related Orchidinae. However, in contrast to O. iricolor and other Orchidinae, O. sphegodes has a shift of the junction between the Inverted Repeat and Small Single Copy regions associated with the loss of the partial duplicated gene ycf1 and the truncation of the ndhF gene. Data on relative genomic coverage and validation by PCR indicate the presence, with a different ratio, of the two plastome types (i.e. with and without ndhF deletion) in both Ophrys species, with a predominance of the deleted type in O. sphegodes. A search for this deleted plastid region in O. sphegodes nuclear genome shows that the deleted region is inserted in a retrotransposon nuclear sequence. The present study provides useful genomic tools for studying conservation and patterns of relationships of this rapidly radiating orchid genus. [ABSTRACT FROM AUTHOR]

Genomic imprinting is the phenomena that leads to silencing of one copy of a gene inherited from a specific parent. Mutations in imprinted regions have been involved in diseases showing parent of origin effects. Identifying genes with evidence of parent of origin expression patterns in family studies allows the detection of more subtle imprinting. Here, we use allele specific expression in lymphoblastoid cell lines from 306 Hutterites related in a single pedigree to provide formal evidence for parent of origin effects. We take advantage of phased genotype data to assign parent of origin to RNA-seq reads in individuals with gene expression data. Our approach identified known imprinted genes, two putative novel imprinted genes, PXDC1 and PWAR6, and 14 genes with asymmetrical parent of origin gene expression. We used gene expression in peripheral blood leukocytes (PBL) to validate our findings, and then confirmed imprinting control regions (ICRs) using DNA methylation levels in the PBLs. [ABSTRACT FROM AUTHOR]

Recursive splicing (RS) is an evolutionarily conserved process of removing long introns via multiple steps of splicing. It was first discovered in Drosophila and recently proven to occur also in humans. The detailed mechanism of recursive splicing is not well understood, in particular, whether it is kinetically coupled with transcription. To investigate the dynamic process that underlies recursive splicing, we systematically characterized 342 RS sites in three human cell types using published time-series data that monitored synchronized Pol II elongation and nascent RNA production with 4-thiouridine labeling. We found that half of the RS events occurred post-transcriptionally with long delays. For at least 18–47% RS introns, we detected RS junction reads only after detecting canonical splicing junction reads, supporting the notion that these introns were removed by both recursive splicing and canonical splicing. Furthermore, the choice of which splicing mechanism was used showed cell type specificity. Our results suggest that recursive splicing supplements, rather than replaces, canonical splicing for removing long introns. [ABSTRACT FROM AUTHOR]

Recursive splicing, a process by which a single intron is removed from pre-mRNA transcripts in multiple distinct segments, has been observed in a small subset of Drosophila melanogaster introns. However, detection of recursive splicing requires observation of splicing intermediates that are inherently unstable, making it difficult to study. Here we developed new computational approaches to identify recursively spliced introns and applied them, in combination with existing methods, to nascent RNA sequencing data from Drosophila S2 cells. These approaches identified hundreds of novel sites of recursive splicing, expanding the catalog of recursively spliced fly introns by 4-fold. A subset of recursive sites were validated by RT-PCR and sequencing. Recursive sites occur in most very long (> 40 kb) fly introns, including many genes involved in morphogenesis and development, and tend to occur near the midpoints of introns. Suggesting a possible function for recursive splicing, we observe that fly introns with recursive sites are spliced more accurately than comparably sized non-recursive introns. [ABSTRACT FROM AUTHOR]

We present ggsashimi, a command-line tool for the visualization of splicing events across multiple samples. Given a specified genomic region, ggsashimi creates sashimi plots for individual RNA-seq experiments as well as aggregated plots for groups of experiments, a feature unique to this software. Compared to the existing versions of programs generating sashimi plots, it uses popular bioinformatics file formats, it is annotation-independent, and allows the visualization of splicing events even for large genomic regions by scaling down the genomic segments between splice sites. ggsashimi is freely available at . It is implemented in python, and internally generates R code for plotting. [ABSTRACT FROM AUTHOR]