eDNA-based detection reveals invasion risks of a biofouling bivalve in the world's largest water diversion project.

Environmental DNA (eDNA) has increasingly been used to detect rare species (e.g., newly introduced nonindigenous species) in both terrestrial and aquatic ecosystems, often with distinct advantages over traditional methods. However, whether water eDNA signals can be used to inform invasion risks remains debatable owing to inherent uncertainties associated with the methods used and the varying conditions among study systems. Here, we sampled eDNA from canals of the central route of the South-to-North Water Diversion Project (hereafter SNWDP) in China to investigate eDNA distribution and efficacy to inform invasion risks in a unique lotic system. We first conducted a total of 16 monthly surveys in this system (two sites in the source reservoir and four sites in the main canal) to test if eDNA could be applied to detect an invasive, biofouling bivalve, the golden mussel Limnoperna fortunei. Second, we initiated a one-time survey in a sub-canal of the SNWDP using refined sampling (12 sites in ~22 km canal) and considered a few environmental predictors. We found that detection of target eDNA in the main canal was achieved up to 1100 km from the putative source population but was restricted to the warmer months (May-November). Detection probability exhibited a significant positive relationship with average daily minimum air temperature and with water temperature, consistent with the expected spawning season. eDNA concentration in the main canal generally fluctuated across months and sites and was generally higher in warmer months. Golden mussel eDNA concentration in the sub-canal decreased significantly with distance from the source and with increasing water temperature and became almost undetectable at ~22 km distance. Given the enormity of the SNWDP, golden mussels may eventually expand their distribution in the main canal, with established "bridgehead" populations facilitating further spread. Our findings suggest an elevated invasion risk of golden mussels in the SNWDP in warm months, highlighting the critical period for spread and, possibly, management.

Development and clinical applications of an enclosed automated targeted NGS library preparation system.

The rapid development of next-generation sequencing (NGS) technology has promoted its wide clinical application in precision medicine for oncology. However, laborious and time-consuming manual operations, highly skilled personnel requirements, and cross-contamination are major challenges for the clinical implementation of NGS technology-based tests. The Automated NGS Diagnostic Solutions (ANDiS) 500 system is a fully enclosed cassette-dependent automated NGS library preparation system. This platform could produce qualified targeted amplicon library in three steps with only 15 min of hands-on time. Rigorous cross-contamination test using simulated contaminant plasmids confirmed that the design of disposable cassette guarantees zero sample cross-contamination. The BRCA1 and BRCA2 mutation detection panel and gastrointestinal cancer-related gene analysis panel for the ANDiS 500 platform showed 100% accuracy and precision in detecting germ-line mutations and somatic mutations respectively. Furthermore, those panels showed 100% concordance with verified methods in a prospective cohort study enrolling 363 patients and a cohort of 45 pan-cancer samples. In conclusion, the ANDiS 500 automated platform could overcome major challenges for implementing NGS assays clinically and is eligible for routine clinical tests.

Comprehensive Evaluation and Application of a Novel Method to Isolate Cell-Free DNA Derived From Bile of Biliary Tract Cancer Patients.

Cell-free DNA (cfDNA) exists in various types of bodily fluids, including plasma, urine, bile, and others. Bile cfDNA could serve as a promising liquid biopsy for biliary tract cancer (BTC) patients, as bile directly contacts tumors in the biliary tract system. However, there is no commercial kit or widely acknowledged method for bile cfDNA extraction. In this study, we established a silica-membrane-based method, namely 3D-BCF, for bile cfDNA isolation, exhibiting effective recovery of DNA fragments in the spike-in assay. We then compared the 3D-BCF method with four other commercial kits: the BIOG cfDNA Easy Kit (BIOG), QIAamp DNA Mini Kit (Qiagen), MagMAXTM Cell-Free DNA Isolation Kit (Thermo Fisher), and NORGEN Urine Cell-Free Circulating DNA Purification Mini Kit (Norgen Biotek). The proposed 3D-BCF method exhibited the highest cfDNA isolation efficiency (p < 0.0001) from patient bile samples, and bile cfDNA of short, medium or long fragments could all be extracted effectively. To test whether the extracted bile cfDNA from patients carries tumor-related genomic information, we performed next-generation sequencing on the cfDNA and verified the gene-mutation results by polymerase chain reaction (PCR)-Sanger chromatograms and copy-number-variation (CNV) detection by fluorescence in situ hybridization (FISH) of tumor tissues. The 3D-BCF method could efficiently extract cfDNA from bile samples, providing technical support for bile cfDNA as a promising liquid biopsy for BTC patient diagnosis and prognosis.

Functional response and size-selective clearance of suspended matter by an invasive mussel.

Filter feeding activities link suspension feeders with their environment and underpin their impact on aquatic ecosystems. Despite their ecological and economic impacts, the functional response and size-selective capture of suspended particulates have not been well documented for the golden mussel Limnoperna fortunei. Here we demonstrated that golden mussels had a type I functional response, with an attack rate a = 0.085 and negligible handling time (h). Clearance rate ranged between 72.6 ± 27.0 and 305.5 ± 105.9 mL ind.-1h-1 (Mean ± S.E.), depending on food concentrations, which exhibited an inverse relationship with clearance rate. Presence of golden mussels suppressed chlorophyll a concentration in experimental mesocosms, the extent of which was dependent on mussel abundance. Concentration of suspended particles in experimental mesocosms experienced a sharp initial decline across all size categories (≤1->50 µm), though with increased final concentration of large particles (>25 µm), indicating packaging and egestion by golden mussels of fine particles (down to ≤1 µm). Capture efficiency of quantitatively-dominant suspended matter (≤1-50 µm) by golden mussels was inversely related to particle size. Animal abundance, particle size, and their interaction (abundance × particle size) determined the extent to which matter was removed from the water column. Presently L. fortunei occurs primarily in the southern end of the central route of South to North Water Diversion Project (China), but the species is spreading north; we anticipate that impacts associated with filtering of L. fortunei will correspond with local population abundance along this gradient.

Bile cell­free DNA as a novel and powerful liquid biopsy for detecting somatic variants in biliary tract cancer.

Tissue sampling of biliary tract carcinomas (BTCs) for molecular characterization is challenging. The aim of this study was to investigate the possibility of identifying individual actionable mutations derived from bile cell­free DNA (cfDNA) using targeted deep sequencing. Ten BTC patients, four with gallbladder carcinomas and six with cholangiocarcinomas, were enrolled in the present study. Using targeted deep sequencing with a panel of 150 tumor­related genes, paired bile cfDNA and tumor DNA were analyzed for mutational variants individually and then compared. The present study, to the best of our knowledge, is the first to reveal that bile cfDNA is predominantly comprised of long DNA fragments, which is not the case for plasma cfDNA. Herein, paired bile cfDNA and tumors from ten BTC patients were examined using targeted deep sequencing. When comparing bile cfDNA and tumor DNA for single nucleotide variation (SNV)/insertion and deletion (Indel), the results using targeted deep sequencing revealed high sensitivity (94.7%) and specificity (99.9%). Additionally, the sensitivity of detecting a copy number variation (CNV) was 75.0%, with a specificity of 98.9%. When comparing two bile extraction methods, including percutaneous transhepatic cholangial drainage and operation, no significant difference in SNV/Indel or CNV detection sensitivity was noted. Moreover, when examining the tumor stage and incidence site, AJCC stage II and the distal bile duct both had significantly decreased CNV detection sensitivities. The present study revealed that targeted deep sequencing can reliably detect mutational variants within bile cfDNA obtained from BTC patients. These preliminary results may shed light on bile cfDNA as a promising liquid biopsy for BTC patients.

Mutation spectrum of germline cancer susceptibility genes among unselected Chinese colorectal cancer patients.

Background: Genetic factors play an important role in colorectal cancer (CRC) risk, yet the prevalence and spectrum of germline cancer susceptibility gene mutations among unselected Chinese CRC patients is largely undetermined. Methods: We performed next-generation sequencing with a 73-genes panel and analyzed the prevalence and spectrum of germline mutations in 618 unselected Chinese CRC patients. We classified all identified germline alterations for pathogenicity and calculated the frequencies of pathogenic mutations. Clinical characteristics were assessed by age and mutation status. Protein expressions and interactions of MLH1 missense variants were evaluated by western blot and co- immunoprecipitation. Results: Overall, 112 (18.1%) of 618 unselected Chinese CRC patients were found to carry at least one pathogenic or likely pathogenic variant (totaling 97 variants), including 70 (11.3%) Lynch syndrome (LS) mutation carriers and 42 (6.8%) non-LS mutation carriers. LS mutation carriers were significantly younger at CRC diagnosis and were more likely to have right-sided, poorly differentiated, early stage, high-frequency microsatellite instability (MSI-H) or dMMR CRC and a family history of cancer compared with noncarriers. Non-LS mutation carriers were more likely to be proficient mismatch repair (pMMR) than noncarriers (p=0.039). We found four clinical variables (gender, tumor histological stage, cancer stage and mutation status) that showed significant differences between patients younger and older than 50 years old. Higher mutation rates were found in patients under 50 years old (p=0.017). Thirty-three novel variants were discovered and evaluated as pathogenic mutations by our study. Conclusion: Given the high frequency and wide spectrum of mutations, genetic testing with a multigene panel should be considered for all Chinese CRC patients under 50 years old and is also needed to determine whether a gene is associated with CRC susceptibility and to promote clinical translation.

Assessment of Blood Tumor Mutational Burden as a Potential Biomarker for Immunotherapy in Patients With Non-Small Cell Lung Cancer With Use of a Next-Generation Sequencing Cancer Gene Panel.

IMPORTANCE: Tumor mutational burden (TMB), as measured by whole-exome sequencing (WES) or a cancer gene panel (CGP), is associated with immunotherapy responses. However, whether TMB estimated by circulating tumor DNA in blood (bTMB) is associated with clinical outcomes of immunotherapy remains to be explored. OBJECTIVES: To explore the optimal gene panel size and algorithm to design a CGP for TMB estimation, evaluate the panel reliability, and further validate the feasibility of bTMB as a clinical actionable biomarker for immunotherapy. DESIGN, SETTING, AND PARTICIPANTS: In this cohort study, a CGP named NCC-GP150 was designed and virtually validated using The Cancer Genome Atlas database. The correlation between bTMB estimated by NCC-GP150 and tissue TMB (tTMB) measured by WES was evaluated in matched blood and tissue samples from 48 patients with advanced NSCLC. An independent cohort of 50 patients with advanced NSCLC was used to identify the utility of bTMB estimated by NCC-GP150 in distinguishing patients who would benefit from anti-programmed cell death 1 (anti-PD-1) and anti-programmed cell death ligand 1 (anti-PD-L1) therapy. The study was performed from July 19, 2016, to April 20, 2018. MAIN OUTCOMES AND MEASURES: Assessment of the Spearman correlation coefficient between bTMB estimated by NCC-GP150 and tTMB calculated by WES. Evaluation of the association of bTMB level with progression-free survival and response to anti-PD-1 and anti-PD-L1 therapy. RESULTS: This study used 2 independent cohorts of patients with NSCLC (cohort 1: 48 patients; mean [SD] age, 60 [13] years; 15 [31.2%] female; cohort 2: 50 patients; mean [SD] age, 58 [8] years; 15 [30.0%] female). A CGP, including 150 genes, demonstrated stable correlations with WES for TMB estimation (median r2 = 0.91; interquartile range, 0.89-0.92), especially when synonymous mutations were included (median r2 = 0.92; interquartile range, 0.91-0.93), whereas TMB estimated by the NCC-GP150 panel found higher correlations with TMB estimated by WES than most of the randomly sampled 150-gene panels. Blood TMB estimated by NCC-GP150 correlated well with the matched tTMB calculated by WES (Spearman correlation = 0.62). In the anti-PD-1 and anti-PD-L1 treatment cohort, a bTMB of 6 or higher was associated with superior progression-free survival (hazard ratio, 0.39; 95% CI, 0.18-0.84; log-rank P = .01) and objective response rates (bTMB ≥6: 39.3%; 95% CI, 23.9%-56.5%; bTMB <6: 9.1%; 95% CI, 1.6%-25.9%; P = .02). CONCLUSIONS AND RELEVANCE: The findings suggest that established NCC-GP150 with an optimized gene panel size and algorithm is feasible for bTMB estimation, which may serve as a potential biomarker of clinical benefit in patients with NSCLC treated with anti-PD-1 and anti-PD-L1 agents.

Yield of trihalomethane, haloacetic acid and chloral upon chlorinating algae after coagulation-filtration: Is pre-oxidation necessarily negative for disinfection by-product control?

Effect of pre-chlorination and pre-ozonation on Microcystis aeruginosa (MA) and Coccomyxa subellipsoidea (CS) as disinfection by-products (DBPs) precursors was investigated after coagulation-filtration. Pre-chlorination considerably decreased the autofluorescence of algae cells but barely influenced cell granularity. In comparison, after pre-ozonation more algae cells were associated with decreased cell size; yet less reduction in the autofluorescence was observed. In MA case, pre-chlorination increased the residual algae density after coagulation-filtration by 132%-146% while pre-ozonation enhanced the algae removal by 26%-28%. In CS case, algae removal was improved by pre-chlorination (32%-45%) and pre-ozonation (7%-45%). Pre-chlorination enhanced the removal of algogenic organic matters (AOM) by coagulation-filtration, especially for tryptophan-like and soluble microbial products. Effect of pre-ozonation on the fluorescence intensity of AOM after coauglation-filtration depended on AOM species and the ratio of [ozone dose]:[algae density]. In both MA and CS cases, chlorine increased the yields of trihalomethane (THM, 25%-78% and 51%-103%), haloacetic acid (HAA, 140%-360% and 167%-233%) and chloral (50%-161% and 68%-108%), respectively. Pre-ozonation decreased the total DBPs yields. For MA-added suspensions, ozone decreased the production of THM, HAA and chloral by 15%-37%, 28%-39% and 60%, respectively. In CS case, chloral yield was decreased by 12%-31% while THM formation was largely unchanged. HAA production varied by ± 1.5 µg/L.

The characteristics of ctDNA reveal the high complexity in matching the corresponding tumor tissues.

BACKGROUND: Next-generation sequencing (NGS) is an efficient and sensitive method to detect mutations from ctDNA. Many features and clinical conditions could significantly affect the concordance between ctDNA and corresponding tumor tissues. Our goal was to systematically investigate the critical factors contributing to different concordance between ctDNA and corresponding tumor tissues. METHODS: We recruited two groups of IIIB or IV lung cancer patients: The standard group to evaluate the accuracy of our method and the concordance between ctDNA and tumor tissues, and the study group with various clinical conditions. We applied our unique identification (UID) indexed capturing-based sequencing (UC-Seq) to ctDNA samples, and confirm the results by Droplet digital PCR (ddPCR). RESULTS: Considering mutations detected from NGS of tumor tissues as golden standard, UC-Seq achieved overall 93.6% sensitivity for SNVs and Indels, and 0.8 Pearson correlation between tumor TMB and bTMB. Efficacious treatments, long sampling date (more than 2 weeks) between tumor tissues and ctDNA and low concentrations of cfDNA (less than 9 ng/ml) could significantly decrease the concordance between ctDNA and tumor tissues. About 84% mutations showed shorter mutant fragment length than that of wild-type fragments, and the AFs of mutations could be significantly enriched in small-size ctDNA. CONCLUSIONS: In late-stage lung cancer patients, ctDNA generally has high concordance with tumor tissues. However it could be significantly affected by three clinical conditions which could dynamically change the content of ctDNA. Moreover, the detection limit could be further extended by enriching small-size ctDNA in the preparation of samples.

High performance of targeted next generation sequencing on variance detection in clinical tumor specimens in comparison with current conventional methods.

BACKGROUND: Next generation sequencing (NGS) is being increasingly applied for assisting cancer molecular diagnosis. However, it is still needed to validate NGS accuracy on detection of DNA alternations based on a large number of clinical samples, especially for DNA rearrangements and copy number variations (CNVs). This study is to set up basic parameters of targeted NGS for clinical diagnosis and to understand advantage of targeted NGS in comparison with the conventional methods of molecular diagnosis. METHODS: Genomic DNA from 1000 Genomes Project and DNA from cancer cell lines have been used to establish the basic parameters for targeted NGS. The following confirmation was conducted by clinical samples. The multiple variants tested by amplification-refractory mutation system (ARMS), fluorescence in situ hybridization (FISH) and immunohistochemistry (IHC) were evaluated by targeted NGS to determine the sensitivity. Furthermore, the multiple variants detected by targeted NGS were confirmed by current conventional methods to elucidate the specificity. RESULTS: At sequencing depth of 500×, the maximal sensitivities on detecting single nucletic variances (SNVs) and small insertions/deletions (Indels) can reach 99% and 98.7% respectively, and in 20% of cancer cells, CNV detection can reach to the maximal level. The following confirmation of the sensitivity and specificity was conducted by a large cohort of clinical samples. For SNV and indel detection in clinical samples, targeted NGS can identify all hotspot mutations with 100% sensitivity and specificity. On ALK fusion detection, about 86% IHC-identified cases could be identified by targeted NGS and all ALK fusion detected by targeted NGS were confirmed by IHC. For HER2-amplification, 14 HER2-amplification cases identified by target NGS were all confirmed by FISH and about 93.3% of Her-2 IHC (3+) cases were identified by targeted NGS. Finally, the targeted NGS platform developed here has accurately detected EGFR hotspot mutations in 215 NSCLC patients. CONCLUSIONS: DNA from cancer cell lines is better than standard DNA as a reference to establish basic parameters for targeted NGS. Comparison of the conventional methods using a large cohort of patient samples confirmed the high preformance of targeted NGS on detecting DNA alterations.

Nucleosome eviction along with H3K9ac deposition enhances Sox2 binding during human neuroectodermal commitment.

Neuroectoderm is an important neural precursor. However, chromatin remodeling and its epigenetic regulatory roles during the differentiation of human neuroectodermal cells (hNECs) from human embryonic stem cells (hESCs) remain largely unexplored. Here, we obtained hNECs through directed differentiation from hESCs, and determined chromatin states in the two cell types. Upon differentiation, H2A.Z-mediated nucleosome depletion leads to an open chromatin structure in promoters and upregulates expression of neuroectodermal genes. Increase in H3K9ac signals and decrease in H3K27me3 signals in promoters result in an active chromatin state and activate neuroectodermal genes. Conversely, decrease in H3K9ac signals and increase in H3K27me3 signals in promoters repress pluripotency genes. Moreover, H3K9ac signals facilitate the pluripotency factor Sox2 binding to target sites unique to hNECs. Knockdown of the acetyltransferase Kat2b erases H3K9ac signals, disrupts Sox2 binding, and fails the differentiation. Our results demonstrate a hierarchy of epigenetic regulation of gene expression during the differentiation of hNECs from hESCs through chromatin remodeling.

Transcription factor ISL1 is essential for pacemaker development and function.

The sinoatrial node (SAN) maintains a rhythmic heartbeat; therefore, a better understanding of factors that drive SAN development and function is crucial to generation of potential therapies, such as biological pacemakers, for sinus arrhythmias. Here, we determined that the LIM homeodomain transcription factor ISL1 plays a key role in survival, proliferation, and function of pacemaker cells throughout development. Analysis of several Isl1 mutant mouse lines, including animals harboring an SAN-specific Isl1 deletion, revealed that ISL1 within SAN is a requirement for early embryonic viability. RNA-sequencing (RNA-seq) analyses of FACS-purified cells from ISL1-deficient SANs revealed that a number of genes critical for SAN function, including those encoding transcription factors and ion channels, were downstream of ISL1. Chromatin immunoprecipitation assays performed with anti-ISL1 antibodies and chromatin extracts from FACS-purified SAN cells demonstrated that ISL1 directly binds genomic regions within several genes required for normal pacemaker function, including subunits of the L-type calcium channel, Ank2, and Tbx3. Other genes implicated in abnormal heart rhythm in humans were also direct ISL1 targets. Together, our results demonstrate that ISL1 regulates approximately one-third of SAN-specific genes, indicate that a combination of ISL1 and other SAN transcription factors could be utilized to generate pacemaker cells, and suggest ISL1 mutations may underlie sick sinus syndrome.

Pwp1 is required for the differentiation potential of mouse embryonic stem cells through regulating Stat3 signaling.

Leukemia inhibitory factor/Stat3 signaling is critical for maintaining the self-renewal and differentiation potential of mouse embryonic stem cells (mESCs). However, the upstream effectors of this pathway have not been clearly defined. Here, we show that periodic tryptophan protein 1 (Pwp1), a WD-40 repeat-containing protein associated with histone H4 modification, is required for the exit of mESCs from the pluripotent state into all lineages. Knockdown (KD) of Pwp1 does not affect mESC proliferation, self-renewal, or apoptosis. However, KD of Pwp1 impairs the differentiation potential of mESCs both in vitro and in vivo. PWP1 chromatin immunoprecipitation-seq results revealed that the PWP1-occupied regions were marked with significant levels of H4K20me3. Moreover, Pwp1 binds to sites in the upstream region of Stat3. KD of Pwp1 decreases the level of H4K20me3 in the upstream region of Stat3 gene and upregulates the expression of Stat3. Furthermore, Pwp1 KD mESCs recover their differentiation potential through suppressing the expression of Stat3 or inhibiting the tyrosine phosphorylation of STAT3. Together, our results suggest that Pwp1 plays important roles in the differentiation potential of mESCs.

H2A.Z nucleosome positioning has no impact on genetic variation in Drosophila genome.

Nucleosome occupancy results in complex sequence variation rate heterogeneity by either increasing mutation rate or inhibiting DNA repair in yeast, fish, and human. H2A.Z nucleosome is extensively involved in gene transcription activation and regulation. To test whether H2A.Z nucleosome has the similar impact on sequence variability in the Drosophila genome, we profiled the H2A.Z nucleosome occupancy and sequence variation rate at gene ends and splicing sites. Consistent with previous studies, H2A.Z nucleosome positioning helps to demarcate the borders of exons. Nucleosome occupancy is anticorrelated with sequence divergence rate in the regions flanking transcription start sites and splicing sites. However, there is no rate heterogeneity between the linker DNA and H2A.Z nucleosomal DNA regardless of nucleosome occupancy, fuzziness, positioning in promoter, coding, and intergenic regions, young or old genes. But the rate at intergenic nucleosomes and the flanking linker regions is higher than that at the genic counterparts. Further analyses found that the high sequence divergence rate in the promoter regions that are usually nucleosome depleted regions may be likely resulted from the high mutation rate in the enriched tandem repeats. Interestingly, within nucleosomes spanning splicing sites, sequence variability of nucleosomal DNA significantly increases from the end within exons to the other end protruding into introns. The relaxed functional constraint in introns contributes to the high rate of nucleosomal DNA residing in introns while the strict functional constraint in exons maintains the low rate of nucleosomal DNA residing in exons. Taken together, H2A.Z nucleosome occupancy has no effect on sequence variability of Drosophila genome, which is likely determined by local sequence composition and the concomitant selection pressure.

[Degradation kinetics of naphthalene by anaerobic sludge and analysis of the bacterial biodiversity].

For efficient biological treatment of naphthalene in the industrial wastewater, activated anaerobic sludge was collected from a wastewater treatment plant of petroleum industry, and domesticated with naphthalene, naphthalene and lactate as electron donors, respectively. When the removal efficiency of naphthalene reached more than 90% in a domestication cycle, degradation kinetics were investigated in batch reactions with naphthalene, naphthalene and lactate as electron donors, respectively. Meanwhile, the microbial DNA was extracted from the sludge with high naphthalene removal efficiency, the 16S rDNA clone library was built up, and the bacterial community was analyzed. The results indicated that the degradation rate of naphthalene in reaction with naphthalene as the sole electron donor was much lower than that with naphthalene and lactate as electron donors. In both domestication modes, the naphthalene concentration and the time followed the first order reaction kinetics model and the kinetic constant K were 3.5 x 10(-3) h(-1) and 16 x 10(-3) h(-1), respectively. In addition, phylogenetic analysis indicated that the bacterial communities in naphthalene and lactate co-metabolism sludge were mainly composed of Deltaproteobacteria, Thermotogae, Bacteroidetes, Chloroflexi and Unclassified bacteria. Deltaproteobacteria was the main phylum in the sludge. In mature anaerobic activated sludge, Desulfobulbus sp. and Kosmotoga accounted for 24.2% and 21.0%, respectively. Smithella, Syntrophobacter and Levilinea were also found in the bioreactor. The study of the bacteria diversity in the anaerobic sludge is conducive to the optimization of reaction conditions for efficient removal of naphthalene.

MicroRNA-449 and microRNA-34b/c function redundantly in murine testes by targeting E2F transcription factor-retinoblastoma protein (E2F-pRb) pathway.

MicroRNAs (miRNAs) mainly function as post-transcriptional regulators and are involved in a wide range of physiological and pathophysiological processes such as cell proliferation, differentiation, apoptosis, and tumorigenesis. Mouse testes express a large number of miRNAs. However, the physiological roles of these testicular miRNAs remain largely unknown. Using microarray and quantitative real time PCR assays, we identified that miRNAs of the microRNA-449 (miR-449) cluster were preferentially expressed in the mouse testis, and their levels were drastically up-regulated upon meiotic initiation during testicular development and in adult spermatogenesis. The expression pattern of the miR-449 cluster resembled that of microRNA-34b/c (miR-34b/c) during spermatogenesis. Further analyses identified that cAMP-responsive element modulator τ and SOX5, two transcription factors essential for regulating male germ cell gene expression, acted as the upstream transactivators to stimulate the expression of the miR-449 cluster in mouse testes. Despite its abundant expression in testicular germ cells, miR-449-null male mice developed normally and exhibited normal spermatogenesis and fertility. Our data further demonstrated that miR-449 shared a cohort of target genes that belong to the E2F transcription factor-retinoblastoma protein pathway with the miR-34 family, and levels of miR-34b/c were significantly up-regulated in miR-449-null testes. Taken together, our data suggest that the miR-449 cluster and miR-34b/c function redundantly in the regulation of male germ cell development in murine testes.